Transcriptome-wide identification and screening of WRKY factors involved in the regulation of taxol biosynthesis in Taxus chinensis

Comprehensive characterization of the WRKY gene family and their potential roles in regulation phenylphenalenone biosynthesis in Musella lasiocarpa

Phenylphenalenone is an important phytoalexin for banana plant protection, yet the mechanisms governing its biosynthesis and regulation remain unclear in plant. WRKY transcription factors play essential roles in modulating plant growth, development, and the biosynthesis of secondary metabolites. In this study, we identified 158 WRKY genes (MlWRKYs) from a phenylphenalenone-rich plant species Musella lasiocarpa. Phylogenetic analysis classified the MlWRKY genes into three distinct subfamilies: type I, type II, and type III. Chromosomal distribution revealed that the MlWRKY genes are clustered on nine respective chromosomes. Additionally, synteny analysis between M. lasiocarpa and Musa balbisiana uncovered highly conserved collinear regions. MIWRKY15, MIWRKY111, MIWRKY122 were identified as candidate genes for regulating PhPNs biosynthesis by integration of multi-omics approaches. We further investigated the expression pattern of MIWRKY15, MIWRKY111, MIWRKY122 genes, as well as their putative target genes MlOMT22 and MlOMT27, the known phenylphenalenone biosynthesis genes in various tissues, including leaves, stems, roots, and seeds. MlWRKY15 and MlOMT22 showed similar expression patterns across tissues. MlWRKY122 and MlOMT27 also displayed consistent expression patterns, suggesting MlWRKY122 may regulate MlOMT27. Additionally, MlWRKY111's expression was inversely correlated with MlOMT27, indicating a potential negative regulation of MlOMT27 by MlWRKY111. This study provides valuable insights into the WRKY family in M. lasiocarpa and will serve as a useful genetic resource for elucidating the regulatory mechanisms of phenylphenalenone biosynthesis.

Regulation of important natural products biosynthesis by WRKY transcription factors in plants

BACKGROUND

Plants produce abundant natural products, among which are species-specific and diversified secondary metabolites that are essential for growth and development, as well as adaptation to adversity and ecology. Moreover, these secondary metabolites are extensively utilized in pharmaceuticals, fragrances, industrial materials, and more. WRKY transcription factors (TFs), as a family of TFs unique to plants, have significant functions in many plant life activities. Especially in recent years, their role in the field of secondary metabolite biosynthesis regulation has received much attention. However, very little comprehensive summarization has been done to review their research progress.

AIM OF REVIEW

The purpose of this work is not only to provide valuable insights into the regulation of WRKY TFs over metabolic pathways through compiling the WRKY TFs involved in these processes, but also to offer research directions for WRKY TFs by summarizing the regulatory modes of WRKY TFs in the biosynthesis of secondary metabolites, thereby increasing the yield of valuable natural products in the future.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Secondary metabolites can be categorized into three major classes-terpenoids, phenolic compounds, and nitrogen-containing compounds-based on their structural characteristics and biosynthetic pathways, and further subdivided into numerous subclasses. We review in detail the research progressregardingthe regulatory roles of WRKY TFs in plant secondary metabolitebiosynthesis and summarize more than 40 major related species. Additionally, we have presented the concepts of action modes of WRKY TFs involved in metabolic pathways, including direct regulation, indirect regulation, co-regulation, and self-regulation. It is helpful for others to investigate the molecular mechanisms of TF-mediated regulation. Furthermore, regarding future research prospects, we believe that research in this area lays the foundation for increasing the yield of important plant-derived natural products by molecular breeding, generating significant economic and social benefits.

Recent advances in paclitaxel biosynthesis and regulation

Paclitaxel (PTX) is a high value plant natural product derived from Taxus (yew) species. This plant specialized metabolite (PSM) and its derivatives constitute a cornerstone for the treatment of an increasing variety of cancers. New applications for PTX also continue to emerge, further promoting demand for this WHO-designated essential medicine. Here we review recent advances in our understanding of PTX biosynthesis and its cognate regulation, which have been enabled by the development of transcriptomic approaches and the recent sequencing and annotation of three Taxus genomes. Collectively, this has resulted in the elucidation of two functional gene sets for PTX biosynthesis, unlocking new potential for the use of heterologous hosts to produce PTX. Knowledge of the PTX pathway also provides a valuable resource for understanding the regulation of this key PSM. Epigenetic regulation of PSM in plant cell culture is a major concern for PTX production, given the loss of PSM production in long-term cell cultures. Recent developments aim to design tools for manipulating epigenetic regulation, potentially providing a means to reverse the silencing of PSM caused by DNA methylation. Exciting times clearly lie ahead for our understanding of this key PSM and improving its production potential.

Unlocking specialized metabolism in medicinal plant biotechnology through plant-microbiome interactions

Medicinal plants produce specialized metabolites (SM) that are used as drugs. However, due to low yields of field cultivation and the increasing market demand, this production method often failed to meet supply needs. Biotechnological alternatives, such as in vitro plant cultures, offer promising solutions. Nonetheless, SM production in these systems remains too low for industrial exploitation, necessitating an elicitation step to induce the plant defense metabolism. Traditional elicitation methods mimic environmental conditions that trigger plant-specialized metabolism, often with an artificial signal that mimics microbial interaction. Recent insights into the essential role of the plant microbiota, provides new opportunities for elicitation strategies by microbial coculture in a controlled environment. The successful co-culture of in vitro medicinal plants with synthetic microbial communities could enable sustainable production of pharmaceutically important SM.

Post-genomic illumination of paclitaxel biosynthesis

Paclitaxel rapidly became one of the most effective anticancer drugs. However, the production of paclitaxel is hindered by substantial challenges, particularly considering the significant quantities of drug required and the inherently low concentration of paclitaxel and its intermediates in plants. Paclitaxel is currently produced in a so-called semi-synthesis in which baccatin III is extracted from Taxus species and chemically converted to paclitaxel. Despite the fact that many of the intermediates of paclitaxel biosynthesis are yet to be experimentally determined, a set of recent papers-facilitated by the sequencing and assembly of three Taxus genomes-has uncovered the minimal gene sets for both baccatin III and paclitaxel biosynthesis. Here we summarize the key milestones towards our understanding of paclitaxel biosynthesis and highlight recent advancements made possible by genome-level analysis of potential key genes involved. We argue that these studies will ultimately pave the way towards the elucidation of the entire paclitaxel biosynthetic pathway and facilitate the industrial production of paclitaxel via synthetic biology approaches. However, several major challenges lie ahead before we can fully tap into the amazing curative potential that taxanes provide.

Transcriptomic analysis of codon usage patterns and gene expression characteristics in leafy spurge

Leafy spurge (Euphorbia esula) is an important herb and potential energy source with medicinal value. Codon usage bias (CUB) is a static feature of genes and genomes that results from adaptation and selection during long-term evolution and facilitates molecular breeding in transgenic plants. Here, we used TransDecoder to identify candidate coding regions from the downloaded leafy spurge transcriptome and generate coding region annotation files based on reference genomes. The whole genome showed A/T bias, especially at terminal positions, and seven high-frequency codons were identified. We compared codon usage frequencies to identify candidate exogenous expression receptor systems for leafy spurge. The identified factors affecting leafy spurge CUB included natural selection and other factors, mutation pressure and base composition, with natural selection and other factors being dominant. The observed CUB was significantly positively correlated with the gene expression levels. Systematic analysis of whole-genome leafy spurge revealed that highly expressed protein-coding genes presented greater CUB than did less expressed protein-coding genes. Furthermore, the highly expressed genes tended to have terminal G/C bases. In summary, we conducted a series of related studies based on the leafy spurge whole-genome sequence and laid a foundation for selecting suitable exogenous expression receptor systems and improving gene expression levels.

Genome-wide investigation of WRKY gene family in Lavandula angustifolia and potential role of LaWRKY57 and LaWRKY75 in the regulation of terpenoid biosynthesis

The WRKY transcription factors are integral to plant biology, serving essential functions in growth, development, stress responses, and the control of secondary metabolism. Through the use of bioinformatics techniques, this research has effectively characterized 207 members of the WRKY family (LaWRKY) present in the entire genome of Lavandula angustifolia. Phylogenetic analysis classified LaWRKYs into three distinct categories based on conserved domains. Collinearity analysis revealed tandem repeats, segmental duplications, and whole genome duplications in LaWRKYs, especially for segmental duplication playing a significant role in gene family expansion. LaWRKYs displayed distinct tissue-specific expression profiles across six different tissues of L. angustifolia. Particularly noteworthy were 12 genes exhibiting high expression in flower buds and calyx, the primary sites of volatile terpenoid production, indicating their potential role in terpenoid biosynthesis in L. angustifolia. RT-qPCR analysis revealed a notable increase in the expression levels of most examined LaWRKY genes in flower buds in response to both intense light and low-temperature conditions, whereas the majority of these genes in leaves were primarily induced by drought stress. However, all genes exhibited downregulation following GA treatment in both flower buds and leaves. Overexpression of LaWRKY57 (La13G01665) and LaWRKY75 (La16G00030) in tobacco led to a reduction in the density of glandular trichomes on leaf surfaces, resulting in changes to the volatile terpenoid composition in the leaves. Specifically, the content of Neophytadiene was significantly elevated compared to wild-type tobacco, while compounds such as eucalyptol, cis-3-Hexenyl iso-butyrate, and D-Limonene were produced, which were absent in the wild type. These findings provide a valuable reference for future investigations into the biological functions of the WRKY gene family in L. angustifolia.

Genome-wide identification and expression analysis of the WRKY gene family in Rhododendron henanense subsp. lingbaoense

Background

This work explored the characteristics of the WRKY transcription factor family in Rhododendron henanense subsp. lingbaoense (Rhl) and the expression patterns of these genes under abiotic stress by conducting bioinformatics and expression analyses.

Methods

RhlWRKY genes were identified from a gene library of Rhl. Various aspects of these genes were analyzed, including genetic structures, conserved sequences, physicochemical properties, cis-acting elements, and chromosomal location. RNA-seq was employed to analyze gene expression in five different tissues of Rhl: roots, stems, leaves, flowers, and hypocotyls. Additionally, qRT-PCR was used to detect changes in the expression of five RhlWRKY genes under abiotic stress.

Result

A total of 65 RhlWRKY genes were identified and categorized into three subfamilies based on their structural characteristics: Groups I, II, and III. Group II was further divided into five subtribes, with shared similar genetic structures and conserved motifs among members of the same subtribe. The physicochemical properties of these proteins varied, but the proteins are generally predicted to be hydrophilic. Most proteins are predicted to be in the cell nucleus, and distributed across 12 chromosomes. A total of 84 cis-acting elements were discovered, with many related to responses to biotic stress. Among the identified RhlWRKY genes, there were eight tandem duplicates and 97 segmental duplicates. The majority of duplicate gene pairs exhibited Ka/Ks values <1, indicating purification under environmental pressure. GO annotation analysis indicated that WRKY genes regulate biological processes and participate in a variety of molecular functions. Transcriptome data revealed varying expression levels of 66.15% of WRKY family genes in all five tissue types (roots, stems, leaves, flowers, and hypocotyls). Five RhlWRKY genes were selected for further characterization and there were changes in expression levels for these genes in response to various stresses.

Conclusion

The analysis identified 65 RhlWRKY genes, among which the expression of WRKY_42 and WRKY_17 were mainly modulated by the drought and MeJA, and WRKY_19 was regulated by the low-temperature and high-salinity conditions. This insight into the potential functions of certain genes contributes to understanding the growth regulatory capabilities of Rhl.

Transcriptome and comparative chloroplast genome analysis of Taxus yunnanensis individuals with high and low paclitaxel yield

Paclitaxel is a potent anti-cancer drug that is mainly produced through semi-synthesis, which still requires plant materials as precursors. The content of paclitaxel and 10-deacetyl baccatin III (10-DAB) in Taxus yunnanensis has been found to differ from that of other Taxus species, but there is little research on the mechanism underlying the variation in paclitaxel content in T. yunnanensis of different provenances. In this experiment, the contents of taxoids and precursors in twigs between a high paclitaxel-yielding individual (TG) and a low paclitaxel-yielding individual (TD) of T. yunnanensis were compared, and comparative analyses of transcriptomes as well as chloroplast genomes were performed. High-performance liquid chromatography (HPLC) detection showed that 10-DAB and baccatin III contents in TG were 18 and 47 times those in TD, respectively. Transcriptomic analysis results indicated that genes encoding key enzymes in the paclitaxel biosynthesis pathway, such as taxane 10-β-hydroxylase (T10βH), 10-deacetylbaccatin III 10-O-acetyltransferase (DBAT), and debenzoyl paclitaxel N-benzoyl transferase (DBTNBT), exhibited higher expression levels in TG. Additionally, qRT-PCR showed that the relative expression level of T10βH and DBAT in TG were 29 and 13 times those in TD, respectively. In addition, six putative transcription factors were identified that may be involved in paclitaxel biosynthesis from transcriptome data. Comparative analysis of plastid genomes showed that the TD chloroplast contained a duplicate of rps12, leading to a longer plastid genome length in TD relative to TG. Fifteen mutation hotspot regions were identified between the two plastid genomes that can serve as candidate DNA barcodes for identifying high-paclitaxel-yield individuals. This experiment provides insight into the difference in paclitaxel accumulation among different provenances of T. yunnanensis individuals.

Genome and Transcriptome Analysis of the Torreya grandis WRKY Gene Family during Seed Development

Torreya grandis, an economically significant evergreen tree species exclusive to subtropical China, is highly valued for its seeds. However, the seed development process of T. grandis remains relatively unexplored. Given the pivotal role WRKY transcription factors (TFs) play in coordinating diverse cellular and biological activities, as well as crucial signaling pathways essential for plant growth and development, and the lack of comprehensive investigation into their specific functions in T. grandis, our study investigated its genome and successfully isolated 78 WRKY genes and categorized them into three distinct clades. A conserved motif analysis unveiled the presence of the characteristic WRKY domain in each identified TgWRKY protein. The examination of gene structures revealed variable numbers of introns (ranging from zero to eight) and exons (ranging from one to nine) among TgWRKY genes. A chromosomal distribution analysis demonstrated the presence of TgWRKY across eight chromosomes in T. grandis. Tissue-specific expression profiling unveiled distinctive patterns of these 78 TgWRKY genes across various tissues. Remarkably, a co-expression analysis integrating RNA-seq data and morphological assessments pinpointed the pronounced expression of TgWRKY25 during the developmental stages of T. grandis seeds. Moreover, a KEGG enrichment analysis, focusing on genes correlated with TgWRKY25 expression, suggested its potential involvement in processes such as protein processing in the endoplasmic reticulum, starch, and sucrose metabolism, thereby modulating seed development in T. grandis. These findings not only underscore the pivotal role of WRKY genes in T. grandis seed development but also pave the way for innovative breeding strategies.

Genomic characterization of WRKY transcription factors related to secoiridoid biosynthesis in Gentiana macrophylla

Gentiana macrophylla is one of Chinese herbal medicines in which 4 kinds of iridoids or secoiridoids, such as loganic acid, sweroside, swertiamarin, and gentiopicroside, are identified as the dominant medicinal secondary metabolites. WRKY, as a large family of transcription factors (TFs), plays an important role in the synthesis of secondary metabolites in plants. Therefore, WRKY genes involved in the biosynthesis of secoiridoids in G. macrophylla were systematically studied. First, a comprehensive genome-wide analysis was performed, and 42 GmWRKY genes were identified, which were unevenly distributed in 12 chromosomes. Accordingly, gene structure, collinearity, sequence alignment, phylogenetic, conserved motif and promoter analyses were performed, and the GmWRKY proteins were divided into three subfamilies based on phylogenetic and multiple sequence alignment analyses. Moreover, the enzyme-encoding genes of the secoiridoid biosynthesis pathway and their promoters were then analysed, and the contents of the four secoiridoids were determined in different tissues. Accordingly, correlation analysis was performed using Pearson's correlation coefficient to construct WRKY gene-enzyme-encoding genes and WRKY gene-metabolite networks. Meanwhile, G. macrophylla seedlings were treated with methyl jasmonate (MeJA) to detect the dynamic change trend of GmWRKYs, biosynthetic genes, and medicinal ingredient accumulation. Thus, a total of 12 GmWRKYs were identified to be involved in the biosynthesis of secoiridoids, of which 8 (GmWRKY1, 6, 12, 17, 33, 34, 38 and 39) were found to regulate the synthesis of gentiopicroside, and 4 (GmWRKY7, 14, 26 and 41) were found to regulate the synthesis of loganic acid. Taken together, this study systematically identified WRKY transcription factors related to the biosynthesis of secoiridoids in G. macrophylla, which could be used as a cue for further investigation of WRKY gene functions in secondary metabolite accumulation.

WRKY transcription factors in passion fruit analysis reveals key PeWRKYs involved in abiotic stress and flavonoid biosynthesis

WRKY transcription factors (TFs) are a superfamily of regulators involved in plant responses to pathogens and abiotic stress. Passion fruit is famous for its unique flavor and nutrient-rich juice, but its growth is limited by environmental factors and pathogens. In this study, 55 WRKY genes were identified from the Passiflora edulis genome. The structure and evolutionary characteristics of PeWRKYs were analyzed using a bioinformatics approach. PeWRKYs were classified into seven subgroups (I, IIa, IIb, IIc, IId, IIe, III) according to their homologs in Arabidopsis thaliana. Group IIa PeWRKY48 gene was highly up-regulated under cold stress by RNA expression analysis, and transgenic PeWRKY48 in yeast and Arabidopsis showed resistance exposure to cold, salt, and drought stress. Metabolome and transcriptome co-expression analysis of two different disease resistance genotypes of P. edulis identified PeWRKY30 as a key TF co-expressed with flavonoid accumulation in yellow fruit P. edulis, which may contribute to biotic or abiotic resistance. The qRT-PCR verified the expression of key genes in different tissues of P. edulis and in different species of Passiflora. This study provides a set of WRKY candidate genes that will facilitate the genetic improvement of disease and abiotic tolerance in passion fruit.

Identification of WRKY gene family members in amaranth based on a transcriptome database and functional analysis of AtrWRKY42-2 in betalain metabolism

Introduction

WRKY TFs (WRKY transcription factors) contribute to the synthesis of secondary metabolites in plants. Betalains are natural pigments that do not coexist with anthocyanins within the same plant. Amaranthus tricolor ('Suxian No.1') is an important leaf vegetable rich in betalains. However, the WRKY family members in amaranth and their roles in betalain synthesis and metabolism are still unclear.

Methods

To elucidate the molecular characteristics of the amaranth WRKY gene family and its role in betalain synthesis, WRKY gene family members were screened and identified using amaranth transcriptome data, and their physicochemical properties, conserved domains, phylogenetic relationships, and conserved motifs were analyzed using bioinformatics methods.

Results

In total, 72 WRKY family members were identified from the amaranth transcriptome. Three WRKY genes involved in betalain synthesis were screened in the phylogenetic analysis of WRKY TFs. RT-qPCR showed that the expression levels of these three genes in red amaranth 'Suxian No.1' were higher than those in green amaranth 'Suxian No.2' and also showed that the expression level of AtrWRKY42 gene short-spliced transcript AtrWRKY42-2 in Amaranth 'Suxian No.1' was higher than that of the complete sequence AtrWRKY42-1, so the short-spliced transcript AtrWRKY42-2 was mainly expressed in 'Suxian No.2' amaranth. Moreover, the total expression levels of AtrWRKY42-1 and AtrWRKY42-2 were down-regulated after GA3 treatment, so AtrWRKY42-2 was identified as a candidate gene. Therefore, the short splice variant AtrWRKY42-2 cDNA sequence, gDNA sequence, and promoter sequence of AtrWRKY42 were cloned, and the PRI 101-AN-AtrWRKY42-2-EGFP vector was constructed to evaluate subcellular localization, revealing that AtrWRKY42-2 is located in the nucleus. The overexpression vector pRI 101-AN-AtrWRKY42-2-EGFP and VIGS (virus-induced gene silencing) vector pTRV2-AtrWRKY42-2 were transferred into leaves of 'Suxian No.1' by an Agrobacterium-mediated method. The results showed that AtrWRKY42-2 overexpression could promote the expression of AtrCYP76AD1 and increase betalain synthesis. A yeast one-hybrid assay demonstrated that AtrWRKY42-2 could bind to the AtrCYP76AD1 promoter to regulate betalain synthesis.

Discussion

This study lays a foundation for further exploring the function of AtrWRKY42-2 in betalain metabolism.

Functions of Representative Terpenoids and Their Biosynthesis Mechanisms in Medicinal Plants

Terpenoids are the broadest and richest group of chemicals obtained from plants. These plant-derived terpenoids have been extensively utilized in various industries, including food and pharmaceuticals. Several specific terpenoids have been identified and isolated from medicinal plants, emphasizing the diversity of biosynthesis and specific functionality of terpenoids. With advances in the technology of sequencing, the genomes of certain important medicinal plants have been assembled. This has improved our knowledge of the biosynthesis and regulatory molecular functions of terpenoids with medicinal functions. In this review, we introduce several notable medicinal plants that produce distinct terpenoids (e.g., Cannabis sativa, Artemisia annua, Salvia miltiorrhiza, Ginkgo biloba, and Taxus media). We summarize the specialized roles of these terpenoids in plant-environment interactions as well as their significance in the pharmaceutical and food industries. Additionally, we highlight recent findings in the fields of molecular regulation mechanisms involved in these distinct terpenoids biosynthesis, and propose future opportunities in terpenoid research, including biology seeding, and genetic engineering in medicinal plants.

miR5298b regulated taxol biosynthesis by acting on TcNPR3, resulting in an alleviation of the strong inhibition of the TcNPR3-TcTGA6 complex in Taxus chinensis

Taxol, a valuable but rare secondary metabolite of the genus Taxus, is an effective anticancer drug. Understanding the regulation of taxol biosynthesis may provide a means to increase taxol content. The microRNA miR5298b was found to promote the accumulation of taxol and upregulate several taxol biosynthesis genes, including DBAT, TASY, and T5H, as demonstrated by experiments using the overexpression and mimicry of transient leaves. Moreover, miR5298b cleaves the mRNA sequence of TcNPR3, a homolog of the salicylic acid receptor AtNPR3/4. Overexpression and knockdown by RNA interference of TcNPR3 confirmed that it repressed taxol biosynthesis. These results indicate that miR5298b enhances taxol biosynthesis via the cleavage of TcNPR3. Yeast two-hybrid bimolecular fluorescence complementation and pull-down assays revealed that TcTGA6, a TGA transcription factor, physically interacted with TcNPR3. Functional experiments showed that TcTGA6 negatively regulates taxol biosynthesis by directly combining with the TGACG motif in the promoters of TASY, T5H, and T10H. TcNPR3 enhances TcTGA6 inhibition Luciferase assays showed that miR5298b alleviated the repression of the TcNPR3-TcTGA6 complex. In summary, miR5298b can cleave TcNPR3, thereby alleviating the inhibition of the TcNPR3-TcTGA6 complex to upregulate taxol biosynthesis genes.

Physiological and metabolic traits of Taxol biosynthesis of endophytic fungi inhabiting plants: Plant-microbial crosstalk, and epigenetic regulators

Attenuating the Taxol productivity of fungi with the subculturing and storage under axenic conditions is the challenge that halts the feasibility of fungi to be an industrial platform for Taxol production. This successive weakening of Taxol productivity by fungi could be attributed to the epigenetic down-regulation and molecular silencing of most of the gene clusters encoding Taxol biosynthetic enzymes. Thus, exploring the epigenetic regulating mechanisms controlling the molecular machinery of Taxol biosynthesis could be an alternative prospective technology to conquer the lower accessibility of Taxol by the potent fungi. The current review focuses on discussing the different molecular approaches, epigenetic regulators, transcriptional factors, metabolic manipulators, microbial communications and microbial cross-talking approaches on restoring and enhancing the Taxol biosynthetic potency of fungi to be industrial platform for Taxol production.

Genome-Wide Identification and Expression Analysis of WRKY Gene Family in Neolamarckia cadamba

The WRKY transcription factor family plays important regulatory roles in multiple biological processes in higher plants. They have been identified and functionally characterized in a number of plant species, but very little is known in Neolamarckia cadamba, a 'miracle tree' for its fast growth and potential medicinal resource in Southeast Asia. In this study, a total of 85 WRKY genes were identified in the genome of N. cadamba. They were divided into three groups according to their phylogenetic features, with the support of the characteristics of gene structures and conserved motifs of protein. The NcWRKY genes were unevenly distributed on 22 chromosomes, and there were two pairs of segmentally duplicated events. In addition, a number of putative cis-elements were identified in the promoter regions, of which hormone- and stress-related elements were shared in many NcWRKYs. The transcript levels of NcWRKY were analyzed using the RNA-seq data, revealing distinct expression patterns in various tissues and at different stages of vascular development. Furthermore, 16 and 12 NcWRKY genes were confirmed to respond to various hormone treatments and two different abiotic stress treatments, respectively. Moreover, the content of cadambine, the active metabolite used for the various pharmacological activities found in N. cadamba, significantly increased after Methyl jasmonate treatment. In addition, expression of NcWRKY64/74 was obviously upregulated, suggesting that they may have a potential function of regulating the biosynthesis of cadambine in response to MeJA. Taken together, this study provides clues into the regulatory roles of the WRKY gene family in N. cadamba.

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.

Evolution of the WRKY66 Gene Family and Its Mutations Generated by the CRISPR/Cas9 System Increase the Sensitivity to Salt Stress in Arabidopsis

Group Ⅲ WRKY transcription factors (TFs) play pivotal roles in responding to the diverse abiotic stress and secondary metabolism of plants. However, the evolution and function of WRKY66 remains unclear. Here, WRKY66 homologs were traced back to the origin of terrestrial plants and found to have been subjected to both motifs' gain and loss, and purifying selection. A phylogenetic analysis showed that 145 WRKY66 genes could be divided into three main clades (Clade A-C). The substitution rate tests indicated that the WRKY66 lineage was significantly different from others. A sequence analysis displayed that the WRKY66 homologs had conserved WRKY and C2HC motifs with higher proportions of crucial amino acid residues in the average abundance. The AtWRKY66 is a nuclear protein, salt- and ABA- inducible transcription activator. Simultaneously, under salt stress and ABA treatments, the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities, as well as the seed germination rates of Atwrky66-knockdown plants generated by the clustered, regularly interspaced, short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9) system, were all lower than those of wild type (WT) plants, but the relative electrolyte leakage (REL) was higher, indicating the increased sensitivities of the knockdown plants to the salt stress and ABA treatments. Moreover, RNA-seq and qRT-PCR analyses revealed that several regulatory genes in the ABA-mediated signaling pathway involved in stress response of the knockdown plants were significantly regulated, being evidenced by the more moderate expressions of the genes. Therefore, the AtWRKY66 likely acts as a positive regulator in the salt stress response, which may be involved in an ABA-mediated signaling pathway.

RNA-Seq Dissects Incomplete Activation of Phytoalexin Biosynthesis by the Soybean Transcription Factors GmMYB29A2 and GmNAC42-1

Glyceollins, isoflavonoid-derived antimicrobial metabolites, are the major phytoalexins in soybean (Glycine max). They play essential roles in providing resistance to the soil-borne pathogen Phytophthora sojae and have unconventional anticancer and neuroprotective activities that render them desirable for pharmaceutical development. Our previous studies revealed that the transcription factors GmMYB29A2 and GmNAC42-1 have essential roles in activating glyceollin biosynthesis, yet each cannot activate the transcription of all biosynthesis genes in the absence of a pathogen elicitor treatment. Here, we report that co-overexpressing both transcription factors is also insufficient to activate glyceollin biosynthesis. To understand this insufficiency, we compared the transcriptome profiles of hairy roots overexpressing each transcription factor with glyceollin-synthesizing roots treated with wall glucan elicitor (WGE) from P. sojae. GmMYB29A2 upregulated most of the WGE-regulated genes that encode enzymatic steps spanning from primary metabolism to the last step of glyceollin biosynthesis. By contrast, GmNAC42-1 upregulated glyceollin biosynthesis genes only when overexpressed in the presence of WGE treatment. This is consistent with our recent discovery that, in the absence of WGE, GmNAC42-1 is bound by GmJAZ1 proteins that inhibit its transactivation activity. WGE, and not GmMYB29A2 or GmNAC42-1, upregulated the heat shock family gene GmHSF6-1, the homolog of Arabidopsis HSFB2a that directly activated the transcription of several glyceollin biosynthesis genes. Our results provide important insights into what biosynthesis genes will need to be upregulated to activate the entire glyceollin biosynthetic pathway.

Transcription level differences in Taxus wallichiana var. mairei elicited by Ce3+, Ce4+ and methyl jasmonate

Taxol is a precious and effective anticancer drug. Cerium and methyl jasmonate (MJ) have been shown to increase the yield of taxol in taxus cells. However, the mechanisms of cerium-mediated and MJ-mediated taxol biosynthesis remain unknown. RNA-Seq was applied to study the overall regulation mechanism of cerium and MJ on taxol biosynthesis and analyze the differences among T. mairei cells elicited by Ce³⁺, Ce⁴⁺ and MJ on transcriptional level . Using sequence homology, 179 unigenes were identified as taxol synthesis genes. Under the condition of 100 μM MJ, taxol synthesis genes were up-regulated. Notably, taxol synthesis genes were down-regulated expression at 1 mM Ce³⁺ and 1 mM Ce⁴⁺. Differential expression genes involved in some related functions were analyzed, such as MAPK signaling pathway and plant-pathogen interaction. Sequence alignment and phylogenetic analysis of nine differentially expressed WRKYs in our data were carried out.

TcJAV3-TcWRKY26 Cascade Is a Missing Link in the Jasmonate-Activated Expression of Taxol Biosynthesis Gene DBAT in Taxus chinensis

Jasmonates (JAs) are the most effective inducers for the biosynthesis of various secondary metabolites. Currently, jasmonate ZIM domain (JAZ) and its interactors, such as MYC2, constitute the main JA signal transduction cascade, and such a cascade fails to directly regulate all the taxol biosynthesis genes, especially the rate-limit gene, DBAT. Another JA signaling branch, JAV and WRKY, would probably fill the gap. Here, TcJAV3 was the closest VQ-motif-containing protein in Taxus chinensis to AtJAV1. Although TcJAV3 was overexpressed in AtJAV1 knockdown mutant, JAVRi17, the enhanced disease resistance to Botrytis cinerea caused by silencing AtJAV1 was completely recovered. The results indicated that TcJAV3 indeed transduced JA signal as AtJAV1. Subsequently, TcWRKY26 was screened out to physically interact with TcJAV3 by using a yeast two-hybrid system. Furthermore, bimolecular fluorescence complementation and luciferase complementary imaging also confirmed that TcJAV3 and TcWRKY26 could form a protein complex in vivo. Our previous reports showed that transient TcWRKY26 overexpression could remarkably increase DBAT expression. Yeast one-hybrid and luciferase activity assays revealed that TcWRKY26 could directly bind with the wa-box of the DBAT promoter to activate downstream reporter genes. All of these results indicated that TcWRKY26 acts as a direct regulator of DBAT, and the TcJAV3−TcWRKY26 complex is actually another JA signal transduction mode that effectively regulates taxol biosynthesis in Taxus. Our results revealed that JAV−WRKY complexes directly regulated DBAT gene in response to JA stimuli, providing a novel model for JA-regulated secondary metabolism. Moreover, JAV could also transduce JA signal and function non-redundantly with JAZ during the regulation of secondary metabolisms.

Identification of WRKY transcription factor family genes in Pinus massoniana Lamb. and their expression patterns and functions in response to drought stress

BACKGROUND

Pinus massoniana Lamb. is the timber species with the widest distribution and the largest afforestation area in China, providing a large amount of timber, turpentine and ecological products. Seasonal drought caused by climate warming severely constrains the quality and growth of P. massoniana forests. WRKY transcription factors play an important role in plant responses to abiotic stress. In this study, the molecular mechanisms by which P. massoniana responds to drought stress were analysed based on the P. massoniana WRKY (PmWRKY) family of genes.

RESULTS

Forty-three PmWRKYs are divided into three major families, 7 sub-families, and the conserved motifs are essentially the same. Among these 43 PmWRKYs express under drought stress but with different expression patterns in response to stress. PmWRKYs respond to drought stress induced by exogenous hormones of SA, ABA, and MeJA. The expression of PmWRKY6, PmWRKY10, and PmWRKY30 up-regulate in different families and tissues under drought stress, while PmWRKY22 down-regulate. Transgenetic tobaccos of PmWRKY31 are with lower malondialdehyde (MDA) content and higher proline (Pro) content than wild type (WT) tobaccos. In transgenic tobaccos of PmWRKY31, expression levels of related genes significantly improve, and drought tolerance enhance.

CONCLUSIONS

This study analysed the molecular biological characteristics of PmWRKYs and investigated the expression patterns and functions of PmWRKYs in response to drought stress in P. massoniana. The results of this study provide a basis for in-depth research of the molecular functions of PmWRKYs in response to drought stress.

Genome-wide identification of Cymbidium sinense WRKY gene family and the importance of its Group III members in response to abiotic stress

Transcription factors (TFs) of the WRKY family play pivotal roles in defense responses and secondary metabolism of plants. Although WRKY TFs are well documented in numerous plant species, no study has performed a genome-wide investigation of the WRKY gene family in Cymbidium sinense. In the present work, we found 64 C. sinense WRKY (CsWRKY) TFs, and they were further divided into eight subgroups. Chromosomal distribution of CsWRKYs revealed that the majority of these genes were localized on 16 chromosomes, especially on Chromosome 2. Syntenic analysis implied that 13 (20.31%) genes were derived from segmental duplication events, and 17 orthologous gene pairs were identified between Arabidopsis thaliana WRKY (AtWRKY) and CsWRKY genes. Moreover, 55 of the 64 CsWRKYs were detectable in different plant tissues in response to exposure to plant hormones. Among them, Group III members were strongly induced in response to various hormone treatments, indicating their potential essential roles in hormone signaling. We subsequently analyzed the function of CsWRKY18 in Group III. The CsWRKY18 was localized in the nucleus. The constitutive expression of CsWRKY18 in Arabidopsis led to enhanced sensitivity to ABA-mediated seed germination and root growth and elevated plant tolerance to abiotic stress within the ABA-dependent pathway. Overall, our study represented the first genome-wide characterization and functional analysis of WRKY TFs in C. sinense, which could provide useful clues about the evolution and functional description of CsWRKY genes.

Genome-Wide Identification and Characterization of the WRKY Gene Family in Scutellaria baicalensis Georgi under Diverse Abiotic Stress

The WRKY gene family is an important inducible regulatory factor in plants, which has been extensively studied in many model plants. It has progressively become the focus of investigation for the secondary metabolites of medicinal plants. Currently, there is no systematic analysis of the WRKY gene family in Scutellaria baicalensis Georgi. For this study, a systematic and comprehensive bioinformatics analysis of the WRKY gene family was conducted based on the genomic data of S. baicalensis. A total of 77 WRKY members were identified and 75 were mapped onto nine chromosomes, respectively. Their encoded WRKY proteins could be classified into three subfamilies: Group I, Group II (II-a, II-b, II-c, II-d, II-e), and Group III, based on the characteristics of the amino acid sequences of the WRKY domain and genetic structure. Syntenic analysis revealed that there were 35 pairs of repetitive fragments. Furthermore, the transcriptome data of roots, stems, leaves, and flowers showed that the spatial expression profiles of WRKYs were different. qRT-PCR analysis revealed that 11 stress-related WRKYs exhibited specific expression patterns under diverse treatments. In addition, sub cellular localization analysis indicated that SbWRKY26 and SbWRKY41 were localized in nucleus. This study is the first to report the identification and characterization of the WRKY gene family in S. baicalensis, which is valuable for the further exploration of the biological function of SbWRKYs. It also provides valuable bioinformatics data for S. baicalensis and provides a reference for assessing the medicinal properties of the genus.

TcMYB29a, an ABA-Responsive R2R3-MYB Transcriptional Factor, Upregulates Taxol Biosynthesis in Taxus chinensis

Paclitaxel (Taxol), a highly modified diterpene agent mainly obtained from Taxus species, is the most widely used anticancer drug. Abscisic acid (ABA) is a well-known stress hormone that plays important roles in the secondary metabolism of plants, and it can also induce the accumulation of taxol in Taxus cell suspension cultures. However, the mechanism behind the regulation of taxol biosynthesis by ABA remains largely unknown. In previous research, a R2R3 MYB transcription factor (TF) TcMYB29a was observed to show a significant correlation with taxol biosynthesis, indicative of its potential role in the taxol biosynthesis. In this study, the TcMYB29a encoded by its gene was further characterized. An expression pattern analysis revealed that TcMYB29a was highly expressed in the needles and roots. Overexpression of TcMYB29a in Taxus chinensis cell suspension cultures led to an increased accumulation of taxol, and upregulated expression of taxol-biosynthesis-related genes, including the taxadiene synthase (TS) gene, the taxane 5α-hydroxylase (T5OH) gene, and the 3'-N-debenzoyl-2'-deoxytaxol-N-benzoyltransferase (DBTNBT) gene as compared to the controls. Chromatin immunoprecipitation (ChIP) assays, yeast one-hybrid (Y1H) assays, electrophoretic mobility shift assays (EMSAs), and dual-luciferase reporter assays verified that TcMYB29a could bind and activate the promoter of TcT5OH. Promoter sequence analysis of TcMYB29a revealed that its promoter containing an AERB site from -313 to -319 was a crucial ABA-responsive element. Subsequently, the ABA treatment assay showed that TcMYB29a was strongly upregulated at 6 h after ABA pretreatment. Furthermore, TcMYB29a was strongly suppressed at 3 h after the methyl jasmonate (MeJA) treatment and was depressed to the platform at 12 h. Taken together, these results reveal that TcMYB29a is an activator that improves the accumulation of taxol in Taxus chinensis cells through an ABA-medicated signaling pathway which is different from JA-medicated signaling pathways for the accumulation of taxol. These findings provide new insights into the potential regulatory roles of MYBs on the expression of taxol biosynthetic genes in Taxus.

Role of female-predominant MYB39-bHLH13 complex in sexually dimorphic accumulation of taxol in Taxus media

Taxus trees are major natural sources for the extraction of taxol, an anti-cancer agent used worldwide. Taxus media is a dioecious woody tree with high taxol yield. However, the sexually dimorphic accumulation of taxoids in T. media is largely unknown. Our study revealed high accumulation of taxoids in female T. media trees using a UPLC-MS/MS method. Thereafter, many differential metabolites and genes between female and male T. media trees were identified using metabolomic and transcriptomic analyses, respectively. Most of the taxol-related genes were predominantly expressed in female trees. A female-specific R2R3-MYB transcription factor gene, TmMYB39, was identified. Furthermore, bimolecular fluorescence complementation and yeast two-hybrid assays suggested the potential interaction between TmMYB39 and TmbHLH13. Several taxol biosynthesis-related promoter sequences were isolated and used for the screening of MYB recognition elements. The electrophoretic mobility shift assay indicated that TmMYB39 could bind to the promoters of the GGPPS, T10OH, T13OH, and TBT genes. Interaction between TmMYB39 and TmbHLH13 transactivated the expression of the GGPPS and T10OH genes. TmMYB39 might function in the transcriptional regulation of taxol biosynthesis through an MYB-bHLH module. Our results give a potential explanation for the sexually dimorphic biosynthesis of taxol in T. media.

Molecular cloning and expression analysis of a WRKY transcription factor gene, GbWRKY20, from Ginkgo biloba

WRKY transcription factors are important regulators of diverse plant life processes. Our aim was to clone and characterize GbWRKY20, a WRKY gene of group IIc, derived from Ginkgo biloba. The cDNA sequence of GbWRKY20 was 818 bp long, encoding a 271-amino acid proteins and containing two introns and three exons. The proteinic molecular weight was 30.99 kDa, with a relevant theoretical isoelectric point of 8.15. Subcellular localization analysis confirmed that the GbWRKY20 protein localized to the nucleus. In total, 75 cis-regulatory elements of 19 different types were identified in the GbWRKY20 promoter sequence, including some elements involved in light responsiveness, anaerobic induction and circadian control, low-temperature responsiveness, as well as salicylic acid (SA) and auxin responsiveness. Expression pattern analysis of plant samples from different developmental stages and tissue types, revealed differential GbWRKY20 expression. The GbWRKY20 transcript was downregulated 12 h after heat treatment and at 4-12 h after drought treatment, but was upregulated 12 h after NaCl, cold and methyl jasmonate treatments. For abscisic acid and SA treatments, the GbWRKY20 transcript was upregulated at 24 h. In summary, GbWRKY20 encoded a newly cloned WRKY transcription factor of G. biloba that might be involved in plant growth and plant responses to abiotic stresses and hormones treatments.

Salicylic Acid-Responsive Factor TcWRKY33 Positively Regulates Taxol Biosynthesis in Taxus chinensis in Direct and Indirect Ways

Taxol is a rare secondary metabolite that accumulates considerably in Taxus species under salicylic acid (SA) and methyl jasmonate treatment. However, the molecular mechanism of its accumulation remains unclear. We investigated TcWRKY33, a nuclear-localized group I WRKY transcription factor, as an SA-responsive regulator of taxol biosynthesis. Overexpression and RNA interference of TcWRKY33 confirmed that TcWRKY33 regulates the expression of most taxol biosynthesis genes, especially 10-deacetylbaccatin III-10-O-acetyltransferase (DBAT) and taxadiene synthase (TASY), which were considered as key enzymes in taxol biosynthesis. Transient overexpression of TcWRKY33 in Taxus chinensis leaves resulted in increased taxol and 10-deacetylbaccatin accumulation by 1.20 and 2.16 times compared with the control, respectively. Furthermore, TcWRKY33, DBAT, and TASY were confirmed to respond positively to SA signals. These results suggested that TcWRKY33 was the missing component of taxol biosynthesis that responds to SA. The sequence analysis identified two W-box motifs in the promoter of DBAT but not in the TASY. Yeast one-hybrid and dual-luciferase activity assays confirmed that TcWRKY33 can bind to the two W-boxes in the promoter of DBAT, upregulating its expression level. Hence, DBAT is a direct target of TcWRKY33. Furthermore, TcERF15, encoding a TASY activator, also contains two W-boxes in its promoter. Yeast one-hybrid and dual-luciferase activity assays further confirmed that TcWRKY33 can upregulate TASY expression through the activation of TcERF15. In summary, TcWRKY33 transmits SA signals and positively regulates taxol biosynthesis genes in two ways: directly and through the activation of other activators. Therefore, TcWRKY33 is an excellent candidate for genetically engineering regulation of taxol biosynthesis in Taxus plants.

The Taxus genome provides insights into paclitaxel biosynthesis

The ancient gymnosperm genus Taxus is the exclusive source of the anticancer drug paclitaxel, yet no reference genome sequences are available for comprehensively elucidating the paclitaxel biosynthesis pathway. We have completed a chromosome-level genome of Taxus chinensis var. mairei with a total length of 10.23 gigabases. Taxus shared an ancestral whole-genome duplication with the coniferophyte lineage and underwent distinct transposon evolution. We discovered a unique physical and functional grouping of CYP725As (cytochrome P450) in the Taxus genome for paclitaxel biosynthesis. We also identified a gene cluster for taxadiene biosynthesis, which was formed mainly by gene duplications. This study will facilitate the elucidation of paclitaxel biosynthesis and unleash the biotechnological potential of Taxus.

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.

The Research Progress of Taxol in Taxus

BACKGROUND

Taxus is a valuable woody species with important medicinal value. The bark of Taxus can produce taxol, a natural antineoplastic drug that is widely used in the treatment of breast, ovarian and lung cancers. However, the low content of taxol in the bark of Taxus can not meet the growing clinical demands, so the current research aims at finding ways to increase taxol production.

OBJECTIVE

In this review, the research progress of taxol including the factors affecting the taxol content, biosynthesis pathway of taxol, production of taxol in vitro and the application of multi-omics approaches in Taxus as well as future research prospects will be discussed.

RESULTS

The taxol content is not only dependent on the species, age and tissues but is also affected by light, moisture levels, temperature, soil fertility and microbes. Most of the enzymes in the taxol biosynthesis pathway have been identified and characterized. Total chemical synthesis, semi-synthesis, plant cell culture and biosynthesis in endophytic fungi have been explored to product taxol. Multi-omics have been used to study Taxus and taxol.

CONCLUSION

Further efforts in the identification of unknown enzymes in the taxol biosynthesis pathway, establishment of the genetic transformation system in Taxus and the regulatory mechanism of taxol biosynthesis and Taxus cell growth will play a significant role in improving the yield of taxol in Taxus cells and plants.

Proteomic insights into protostane triterpene biosynthesis regulatory mechanism after MeJA treatment in Alisma orientale (Sam.) Juz

Protostane triterpenes in Alisma orientale (Sam.) Juz. have unique structural features with distinct pharmacological activities. Previously we have demonstrated that protostane triterpene biosynthesis could be regulated by methyl jasmonate (MeJA) induction in A. orientale. Here, proteomic investigation reveals the MeJA mediated regulation of protostane triterpene biosynthesis. In our study, 281 differentially abundant proteins were identified from MeJA-treated compared to control groups, while they were mainly associated with triterpene biosynthesis, α-linolenic acid metabolism, carbohydrate metabolism and response to stress/defense. Key enzymes 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), squalene epoxidase (SE), oxidosqualene cyclase (OSC) and cytochrome P450s which potentially involved in protostane triterpene biosynthesis were significantly enriched in MeJA-treated group. Basic Helix-loop-helix (bHLH), MYB, and GRAS transcription factors were enhanced after MeJA treatment, and they also improved the expressions of key enzymes in Mevalonate pathway and protostane triterpene. Then, MeJA also could increase the expression of α-galactosidase (α-GAL), thereby promoting carbohydrate decomposition, and providing energy and carbon skeletons for protostane triterpene precursor biosynthesis. As well, exogenous MeJA treatment upregulated 13-lipoxygenase (13-LOX), allene oxide synthase (AOS) and allene oxide cyclase (AOC) involved in α-linolenic acid metabolism, leading to the accumulation of endogenous MeJA and activation of the protostane triterpene biosynthesis transduction. Finally, MeJA upregulated stress/defence-related proteins, as to enhance the defence responses activity of plants. These results were further verified by quantitative real-time PCR analysis of 19 selected genes and content analysis of protostane triterpene. The results provide some new insights into the role of MeJA in protostane triterpene biosynthesis.

Omic analysis of the endangered Taxaceae species Pseudotaxus chienii revealed the differences in taxol biosynthesis pathway between Pseudotaxus and Taxus yunnanensis trees

BACKGROUND

Taxol is an efficient anticancer drug accumulated in Taxus species. Pseudotaxus chienii is an important member of Taxaceae, however, the level of six taxoids in P. chienii is largely unknown.

RESULTS

High accumulation of 10-DAB, taxol, and 7-E-PTX suggested that P. chienii is a good taxol-yielding species for large-scale cultivation. By the omics approaches, a total of 3,387 metabolites and 61,146 unigenes were detected and annotated. Compared with a representative Taxus tree (Taxus yunnanensis), most of the differentially accumulated metabolites and differential expressed genes were assigned into 10 primary and secondary metabolism pathways. Comparative analyses revealed the variations in the precursors and intermediate products of taxol biosynthesis between P. chienii and T. yunnanensis. Taxusin-like metabolites highly accumulated in P. chienii, suggesting a wider value of P. chienii in pharmaceutical industry.

CONCLUSIONS

In our study, the occurrence of taxoids in P. chienii was determined. The differential expression of key genes involved in the taxol biosynthesis pathway is the major cause of the differential accumulation of taxoids. Moreover, identification of a number of differentially expressed transcription factors provided more candidate regulators of taxol biosynthesis. Our study may help to reveal the differences between Pseudotaxus and Taxus trees, and promote resource utilization of the endangered and rarely studied P. chienii.

Recent Advances in Metabolic Engineering, Protein Engineering, and Transcriptome-Guided Insights Toward Synthetic Production of Taxol

The diterpenoid paclitaxel (Taxol®) is a blockbuster anticancer agent that was originally isolated from the Pacific yew (Taxus brevifolia) five decades ago. Despite the wealth of information gained over the years on Taxol research, there still remains supply issues to meet increasing clinical demand. Although alternative Taxol production methods have been developed, they still face several drawbacks that cause supply shortages and high production costs. It is highly desired to develop biotechnological production platforms for Taxol, however, there are still gaps in our understanding of the biosynthetic pathway, catalytic enzymes, regulatory and control mechanisms that hamper production of this critical drug by synthetic biology approaches. Over the past 5 years, significant advances were made in metabolic engineering and optimization of the Taxol pathway in different hosts, leading to accumulation of taxane intermediates. Computational and experimental approaches were leveraged to gain mechanistic insights into the catalytic cycle of pathway enzymes and guide rational protein engineering efforts to improve catalytic fitness and substrate/product specificity, especially of the cytochrome P450s (CYP450s). Notable breakthroughs were also realized in engineering the pathway in plant hosts that are more promising in addressing the challenging CYP450 chemistry. Here, we review these recent advances and in addition, we summarize recent transcriptomic data sets of Taxus species and elicited culture cells, and give a bird's-eye view of the information that can be gleaned from these publicly available resources. Recent mining of transcriptome data sets led to discovery of two putative pathway enzymes, provided many lead candidates for the missing steps and provided new insights on the regulatory mechanisms governing Taxol biosynthesis. All these inferences are relevant to future biotechnological production of Taxol.

Tissue-specific study across the stem of Taxus media identifies a phloem-specific TmMYB3 involved in the transcriptional regulation of paclitaxel biosynthesis

Taxus stem barks can be used for extraction of paclitaxel. However, the composition of taxoids across the whole stem and the stem tissue-specificity of paclitaxel biosynthesis-related enzymes remain largely unknown. We used cultivated Taxus media trees for analyses of the chemical composition and protein of major stem tissues by an integrated metabolomic and proteomic approach, and the role of TmMYB3 in paclitaxel biosynthesis was investigated. The metabolomic landscape analysis showed differences in stem tissue-specific accumulation of metabolites. Phytochemical analysis revealed that there is high accumulation of paclitaxel in the phloem. Ten key enzymes involved in paclitaxel biosynthesis were identified, most of which are predominantly produced in the phloem. The full-length sequence of TmMYB3 and partial promoter sequences of five paclitaxel biosynthesis-related genes were isolated. Several MYB recognition elements were found in the promoters of TBT, DBTNBT and TS. Further in vitro and in vivo investigations indicated that TmMYB3 is involved in paclitaxel biosynthesis by activating the expression of TBT and TS. Differences in the taxoid composition of different stem tissues suggest that the whole stem of T. media has potential for biotechnological applications. Phloem-specific TmMYB3 plays a role in the transcriptional regulation of paclitaxel biosynthesis, and may explain the phloem-specific accumulation of paclitaxel.

Transcriptional reprogramming strategies and miRNA-mediated regulation networks of Taxus media induced into callus cells from tissues

Background

Taxus cells are a potential sustainable and environment-friendly source of taxol, but they have low survival ratios and slow grow rates. Despite these limitations, Taxus callus cells induced through 6 months of culture contain more taxol than their parent tissues. In this work, we utilized 6-month-old Taxus media calli to investigate their regulatory mechanisms of taxol biosynthesis by applying multiomics technologies. Our results provide insights into the adaptation strategies of T. media by transcriptional reprogramming when induced into calli from parent tissues.

Results

Seven out of 12 known taxol, most of flavonoid and phenylpropanoid biosynthesis genes were significantly upregulated in callus cells relative to that in the parent tissue, thus indicating that secondary metabolism is significantly strengthened. The expression of genes involved in pathways metabolizing biological materials, such as amino acids and sugars, also dramatically increased because all nutrients are supplied from the medium. The expression level of 94.1% genes involved in photosynthesis significantly decreased. These results reveal that callus cells undergo transcriptional reprogramming and transition into heterotrophs. Interestingly, common defense and immune activities, such as “plant–pathogen interaction” and salicylic acid- and jasmonic acid-signaling transduction, were repressed in calli. Thus, it’s an intelligent adaption strategy to use secondary metabolites as a cost-effective defense system. MiRNA- and degradome-sequencing results showed the involvement of a precise regulatory network in the miRNA-mediated transcriptional reprogramming of calli. MiRNAs act as direct regulators to enhance the metabolism of biological substances and repress defense activities. Given that only 17 genes of secondary metabolite biosynthesis were effectively regulated, miRNAs are likely to play intermediate roles in the biosynthesis of secondary metabolites by regulating transcriptional factors (TFs), such as ERF, WRKY, and SPL.

Conclusion

Our results suggest that increasing the biosynthesis of taxol and other secondary metabolites is an active regulatory measure of calli to adapt to heterotrophic culture, and this alteration mainly involved direct and indirect miRNA-induced transcriptional reprogramming. These results expand our understanding of the relationships among the metabolism of biological substances, the biosynthesis of secondary metabolites, and defense systems. They also provide a series of candidate miRNAs and transcription factors for taxol biosynthesis.

New different origins and evolutionary processes of AP2/EREBP transcription factors in Taxus chinensis

BACKGROUND

Taxus spp. produces the anticancer drug, taxol, and hence is planted as an industrial crop in China. APETALA2/ethylene response element binding proteins (AP2/EREBPs) are the key regulators of plant development, growth, and stress responses. Several homologues control taxol biosynthesis. Identifying the AP2/EREBP proteins from Taxus is important to increase breeding and production and clarify their evolutionary processes.

RESULTS

Among the 90 genes from multi Taxus chinensis transcriptome datasets, 81 encoded full-length AP2-containing proteins. A domain structure highly similar to that of angiosperm AP2/EREBPs was found in 2 AP2, 2 ANT, 1 RAV, 28 dehydration-responsive element-binding proteins, and 47 ethylene-responsive factors contained, indicating that they have extremely conservative evolution processes. A new subgroup protein, TcA3Bz1, contains three conserved AP2 domains and, a new domain structure of AP2/EREBPs that is different from that of known proteins. The new subtype AP2 proteins were also present in several gymnosperms (Gingko biloba) and bryophytes (Marchantia polymorpha). However, no homologue was found in Selaginella moellendorffii, indicating unknown evolutionary processes accompanying this plant's evolution. Moreover, the structures of the new subgroup AP2/EREBPs have different conserved domains, such as B3, zf-C3Hc3H, and agent domains, indicating their divergent evolution in bryophytes and gymnosperms. Interestingly, three repeats of AP2 domains have separately evolved from mosses to gymnosperms for most of the new proteins, but the AP2 domain of Gb_11937 has been replicated.

CONCLUSION

The new subtype AP2/EREBPs have different origins and would enrich our knowledge of the molecular structure, origin, and evolutionary processes of AP2/EREBP transcription factors in plants.

Transcriptome Sequencing Reveals Regulatory Mechanisms of Taxol Synthesis in Taxus wallichiana var. Mairei

Taxol is one of the most potent and effective anticancer drugs and is originally isolated from Taxus species. To investigate the specific regulatory mechanisms of taxol synthesis in Taxus wallichiana var. mairei, RNA-seq was conducted to reveal the differences in transcriptional levels between wild type (WT) and “Jinxishan” (JXS), a cultivar selected from a population of Taxus mairei that shows about 3-fold higher taxol content in the needles than WT. Our results indicated that high expressions of the genes taxadienol acetyltransferase (TAT), taxadiene 5-alpha hydroxylase (T5H), 5-alpha-taxadienol-10-beta-hydroxylase (T10OH), and 2-debenzoyl-7,13-diacetylbaccatin III-2-O-benzoyl-transferase (DBBT), which catalyze a series of key acetylation and hydroxylation steps, are the main cause of high taxol content in JXS. Moreover, in the present study, the activation of jasmonic acid (JA) signal transduction and its crosstalk with gibberellin (GA), auxin, and ethylene (ET) explained the elevation of differentially expressed genes (DEGs) from the taxol biosynthesis pathway. This also indicates that taxol biosynthesis in T. mairei is associated with the balance of cell development and defense. TF-encoding (transcriptional factor) genes, represented by the ethylene-responsive transcription factor (ERF), basic/helix-loop-helix (bHLH), MYB, and WRKY families, were detected as differentially expressed between JXS and WT, further indicating that the regulation of hormone signaling on taxol biosynthesis genes was mediated by transcription factors (TFs). To our knowledge, this is the first study to illustrate the regulatory mechanisms of taxol synthesis in a new cultivar of T. mairei with a high taxol content in its needles. These transcriptome data provide reasonable explanations for the variation of taxol content between WT and JXS.

Transcriptome analyses provide insights into the expression pattern and sequence similarity of several taxol biosynthesis-related genes in three Taxus species

BACKGROUND

Taxol is an efficient anticancer drug; however, the accumulation of taxoids can vary hugely among Taxus species. The mechanism underlying differential accumulation of taxoids is largely unknown. Thus, comparative analysis of the transcriptomes in three Taxus species, including T. media, T. mairei and T. cuspidata, was performed.

RESULTS

KEGG enrichment analysis revealed that the diterpenoid biosynthesis and cytochrome P450 pathways were significantly enriched in different comparisons. Differential expressions of these taxol biosynthesis related genes might be a potential explanation for the interspecific differential accumulation of taxol and its derivatives. Besides, the sequences of several MEP pathway-associated genes, such as DXS, DXR, MCT, CMK, MDS, HDS, HDR, IPPI, and GGPPS, were re-assembled based on independent transcriptomes from the three Taxus species. Phylogenetic analysis of these MEP pathway-associated enzymes also showed a high sequence similarity between T. media and T. cuspidata. Moreover, 48 JA-related transcription factor (TF) genes, including 10 MYBs, 5 ERFs, 4 RAPs, 3 VTCs, and 26 other TFs, were analyzed. Differential expression of these JA-related TF genes suggested distinct responses to exogenous JA applications in the three Taxus species.

CONCLUSIONS

Our results provide insights into the expression pattern and sequence similarity of several taxol biosynthesis-related genes in three Taxus species. The data give us an opportunity to reveal the mechanism underlying the variations in the taxoid contents and to select the highest-yielding Taxus species.

Genome-wide identification and expression analyses of WRKY transcription factor family members from chickpea (Cicer arietinum L.) reveal their role in abiotic stress-responses

BACKGROUND

WRKY proteins play a vital role in the regulation of several imperative plant metabolic processes and pathways, especially under biotic and abiotic stresses. Although WRKY genes have been characterized in various major crop plants, their identification and characterization in pulse legumes is still in its infancy. Chickpea (Cicer arietinum L.) is the most important pulse legume grown in arid and semi-arid tropics.

OBJECTIVE

In silico identification and characterization of WRKY transcription factor-encoding genes in chickpea genome.

METHODS

For this purpose, a systematic genome-wide analysis was carried out to identify the non-redundant WRKY transcription factors in the chickpea genome.

RESULTS

We have computationally identified 70 WRKY-encoding non-redundant genes which were randomly distributed on all the chickpea chromosomes except chromosome 8. The evolutionary phylogenetic analysis classified the WRKY proteins into three major groups (I, II and III) and seven sub-groups (IN, IC, IIa, IIb, IIc, IId and IIe). The gene structure analysis revealed the presence of 2-7 introns among the family members. Along with the presence of absolutely conserved signatory WRKY domain, 19 different domains were also found to be conserved in a group-specific manner. Insights of gene duplication analysis revealed the predominant role of segmental duplications for the expansion of WRKY genes in chickpea. Purifying selection seems to be operated during the evolution and expansion of paralogous WRKY genes. The transcriptome data-based in silico expression analysis revealed the differential expression of CarWRKY genes in root and shoot tissues under salt, drought, and cold stress conditions. Moreover, some of these genes showed identical expression pattern under these stresses, revealing the possibility of involvement of these genes in conserved abiotic stress-response pathways.

CONCLUSION

This genome-wide computational analysis will serve as a base to accelerate the functional characterization of WRKY TFs especially under biotic and abiotic stresses.

Transcriptome Profiles Reveal the Crucial Roles of Auxin and Cytokinin in the "Shoot Branching" of Cremastra appendiculata

Cremastra appendiculata has become endangered due to reproductive difficulties. Specifically, vegetative reproduction is almost its only way to reproduce, and, under natural conditions, it cannot grow branches, resulting in an extremely low reproductive coefficient (reproductive percentage). Here, we performed RNA-Seq and a differentially expressed gene (DEG) analysis of the three stages of lateral bud development in C. appendiculata after decapitation-dormancy (D2), transition (TD2), and emergence (TG2)-and the annual axillary bud natural break (G1) to gain insight into the molecular regulatory network of shoot branching in this plant. Additionally, we applied the auxin transport inhibitors N-1-naphthylphthalamic acid (NPA) and 2,3,5-triiodibenzoic acid (TIBA) to a treated pseudobulb string of C. appendiculata to verify the conclusions obtained by the transcriptome. RNA-Seq provided a wealth of valuable information. Successive pairwise comparative transcriptome analyses revealed 5988 genes as DEGs. GO (Gene Ontology) and KEGG (Kyoto encyclopedia of genes and genomes) analyses of DEGs showed significant enrichments in phytohormone biosynthesis and metabolism, regulation of hormone levels, and a hormone-mediated signaling pathway. qRT-PCR validation showed a highly significant correlation (p < 0.01) with the RNA-Seq generated data. High-performance liquid chromatography (HPLC) and qRT-PCR results showed that, after decapitation, the NPA- and TIBA-induced lateral buds germinated due to rapidly decreasing auxin levels, caused by upregulation of the dioxygenase for auxin oxidation gene (DAO). Decreased auxin levels promoted the expression of isopentenyl transferase (IPT) and cytochrome P450 monooxygenase, family 735, subfamily A (CYP735A) genes and inhibited two carotenoid cleavage dioxygenases (CCD7 and CCD8). Zeatin levels significantly increased after the treatments. The increased cytokinin levels promoted the expression of WUSCHEL (WUS) and inhibited expression of BRANCHED1 (BRC1) in the cytokinin signal transduction pathway and initiated lateral bud outgrowth. Our data suggest that our theories concerning the regulation of shoot branching and apical dominance is really similar to those observed in annual plants. Auxin inhibits bud outgrowth and tends to inhibit cytokinin levels. The pseudobulb in the plant behaves in a similar manner to that of a shoot above the ground.