Regulation of specialized metabolism by WRKY transcription factors

Co-Localization of Resistance and Metabolic Quantitative Trait Loci on Carrot Genome Reveals Fungitoxic Terpenes and Related Candidate Genes Associated with the Resistance to Alternaria dauci

Alternaria leaf blight, caused by the fungus Alternaria dauci, is the most damaging foliar disease of carrot. Some carrot genotypes exhibit partial resistance to this pathogen and resistance Quantitative Trait Loci (rQTL) have been identified. Co-localization of metabolic QTL and rQTL identified camphene, α-pinene, α-bisabolene, β-cubebene, caryophyllene, germacrene D and α-humulene as terpenes potentially involved in carrot resistance against ALB. By combining genomic and transcriptomic analyses, we identified, under the co-localization regions, terpene-related genes which are differentially expressed between a resistant and a susceptible carrot genotype. These genes include five terpene synthases and twenty transcription factors. In addition, significant mycelial growth inhibition was observed in the presence of α-humulene and caryophyllene.

Genomic profiling of WRKY transcription factors and functional analysis of CcWRKY7, CcWRKY29, and CcWRKY32 related to protoberberine alkaloids biosynthesis in Coptis chinensis Franch

Coptis chinensis Franch. (Huanglian in Chinese) is an important economic crop with medicinal value. Its rhizome has been used as a traditional herbal medicine for thousands of years in Asia. Protoberberine alkaloids, as the main bioactive component of Coptis chinensis, have a series of pharmacological activities. However, the protoberberine alkaloids content of C. chinensis is relatively low. Understanding the molecular mechanisms affecting the transcriptional regulation of protoberberine alkaloids would be crucial to increase their production via metabolic engineering. WRKY, one of the largest plant-specific gene families, regulates plant defense responses via the biosynthesis of specialized metabolites such as alkaloids. Totally, 41 WRKY transcription factors (TFs) related to protoberberine alkaloid biosynthesis were identified in the C. chinensis genome and classified into three groups based on phylogenetic and conserved motif analyses. Three WRKY genes (CcWRKY7, CcWRKY29, and CcWRKY32) may regulate protoberberine alkaloid biosynthesis, as suggested by gene-specific expression patterns, metabolic pathways, phylogenetic, and dual-luciferase analysis. Furthermore, the CcWRKY7, CcWRKY29, and CcWRKY32 proteins were specifically detected in the nucleus via subcellular localization. This study provides a basis for understanding the regulatory mechanisms of protoberberine alkaloid biosynthesis and valuable information for breeding C. chinensis varieties.

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

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

TaGSNE, a WRKY transcription factor, overcomes the trade-off between grain size and grain number in common wheat and is associated with root development

Wheat is one of the world's major staple food crops, and breeding for improvement of grain yield is a priority under the scenarios of climate change and population growth. WRKY transcription factors are multifaceted regulators in plant growth, development, and responses to environmental stimuli. In this study, we identify the WRKY gene TaGSNE (Grain Size and Number Enhancer) in common wheat, and find that it has relatively high expression in leaves and roots, and is induced by multiple abiotic stresses. Eleven single-nucleotide polymorphisms were identified in TaGSNE, forming two haplotypes in multiple germplasm collections, named as TaGSNE-Hap-1 and TaGSNE-Hap-2. In a range of different environments, TaGSNE-Hap-2 was significantly associated with increases in thousand-grain weight (TGW; 3.0%) and spikelet number per spike (4.1%), as well as with deeper roots (10.1%) and increased root dry weight (8.3%) at the mid-grain-filling stage, and these were confirmed in backcross introgression populations. Furthermore, transgenic rice lines overexpressing TaGSNE had larger panicles, more grains, increased grain size, and increased grain yield relative to the wild-type control. Analysis of geographic and temporal distributions revealed that TaGSNE-Hap-2 is positively selected in China and Pakistan, and TaGSNE-Hap-1 in Europe. Our findings demonstrate that TaGSNE overcomes the trade-off between TGW/grain size and grain number, leading us to conclude that these elite haplotypes and their functional markers could be utilized in marker-assisted selection for breeding high-yielding varieties.

A single-nucleotide polymorphism in WRKY33 promoter is associated with the cold sensitivity in cultivated tomato

Natural variations in cis-regulatory regions often affect crop phenotypes by altering gene expression. However, the mechanism of how promoter mutations affect gene expression and crop stress tolerance is still poorly understood. In this study, by analyzing RNA-sequencing (RNA-Seq) data and reverse transcription quantitative real-time PCR validation in the cultivated tomato and its wild relatives, we reveal that the transcripts of WRKY33 are almost unchanged in cold-sensitive cultivated tomato Solanum lycopersicum L. 'Ailsa Craig' but are significantly induced in cold-tolerant wild tomato relatives Solanum habrochaites LA1777 and Solanum pennellii LA0716 under cold stress. Overexpression of SlWRKY33 or ShWRKY33 positively regulates cold tolerance in tomato. Variant of the critical W-box in SlWRKY33 promoter results in the loss of self-transcription function of SlWRKY33 under cold stress. Analysis integrating RNA-Seq and chromatin immunoprecipitation sequencing data reveals that SlWRKY33 directly targets and induces multiple kinases, transcription factors, and molecular chaperone genes, such as CDPK11, MYBS3, and BAG6, thus enhancing cold tolerance. In addition, heat- and Botrytis-induced WRKY33s expression in both wild and cultivated tomatoes are independent of the critical W-box variation. Our findings suggest nucleotide polymorphism in cis-regulatory regions is crucial for different cold sensitivity between cultivated and wild tomato plants.

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

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

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.

Spatial Genomic Resource Reveals Molecular Insights into Key Bioactive-Metabolite Biosynthesis in Endangered Angelica glauca Edgew

Angelica glauca Edgew, which is an endangered medicinal and aromatic herb, is a rich source of numerous industrially important bioactive metabolites, including terpenoids, phenolics, and phthalides. Nevertheless, genomic interventions for the sustainable utilization and restoration of its genetic resources are greatly offset due to the scarcity of the genomic resources and key regulators of the underlying specialized metabolism. To unravel the global atlas of the specialized metabolism, the first spatial transcriptome sequencing of the leaf, stem, and root generated 109 million high-quality paired-end reads, assembled de novo into 81,162 unigenes, which exhibit a 61.53% significant homology with the six public protein databases. The organ-specific clustering grouped 1136 differentially expressed unigenes into four subclusters differentially enriched in the leaf, stem, and root tissues. The prediction of the transcriptional-interactome network by integrating enriched gene ontology (GO) and the KEGG metabolic pathways identified the key regulatory unigenes that correspond to terpenoid, flavonoid, and carotenoid biosynthesis in the leaf tissue, followed by the stem and root tissues. Furthermore, the stem and root-specific significant enrichments of phenylalanine ammonia lyase (PAL), cinnamate-4-hydroxylase (C4H), and caffeic acid 3-O-methyltransferase (COMT) indicate that phenylalanine mediated the ferulic acid biosynthesis in the stem and root. However, the root-specific expressions of NADPH-dependent alkenal/one oxidoreductase (NADPH-AOR), S-adenosyl-L-methionine-dependent methyltransferases (SDMs), polyketide cyclase (PKC), and CYP72A15 suggest the “root” as the primary site of phthalide biosynthesis. Additionally, the GC-MS and UPLC analyses corresponded to the organ-specific gene expressions, with higher contents of limonene and phthalide compounds in the roots, while there was a higher accumulation of ferulic acid in the stem, followed by in the root and leaf tissues. The first comprehensive genomic resource with an array of candidate genes of the key metabolic pathways can be potentially utilized for the targeted upscaling of aromatic and pharmaceutically important bioactive metabolites. This will also expedite genomic-assisted conservation and breeding strategies for the revival of the endangered A. glauca.

Role of WRKY Transcription Factors in Regulation of Abiotic Stress Responses in Cotton

Environmental factors are the major constraints in sustainable agriculture. WRKY proteins are a large family of transcription factors (TFs) that regulate various developmental processes and stress responses in plants, including cotton. On the basis of Gossypium raimondii genome sequencing, WRKY TFs have been identified in cotton and characterized for their functions in abiotic stress responses. WRKY members of cotton play a significant role in the regulation of abiotic stresses, i.e., drought, salt, and extreme temperatures. These TFs either activate or repress various signaling pathways such as abscisic acid, jasmonic acid, salicylic acid, mitogen-activated protein kinases (MAPK), and the scavenging of reactive oxygen species. WRKY-associated genes in cotton have been genetically engineered in Arabidopsis, Nicotiana, and Gossypium successfully, which subsequently enhanced tolerance in corresponding plants against abiotic stresses. Although a few review reports are available for WRKY TFs, there is no critical report available on the WRKY TFs of cotton. Hereby, the role of cotton WRKY TFs in environmental stress responses is studied to enhance the understanding of abiotic stress response and further improve in cotton plants.

Functional analysis of Pogostemon cablin farnesyl pyrophosphate synthase gene and its binding transcription factor PcWRKY44 in regulating biosynthesis of patchouli alcohol

Farnesyl pyrophosphate synthase (FPPS) plays an important role in the synthesis of plant secondary metabolites, but its function and molecular regulation mechanism remain unclear in Pogostemon cablin. In this study, the full-length cDNA of the FPP synthase gene from P. cablin (PcFPPS) was cloned and characterized. The expressions of PcFPPS are different among different tissues (highly in P. cablin flowers). Subcellular localization analysis in protoplasts indicated that PcFPPS was located in the cytoplasm. PcFPPS functionally complemented the lethal FPPS deletion mutation in yeast CC25. Transient overexpression of PcFPPS in P. cablin leaves accelerated terpene biosynthesis, with an ~47% increase in patchouli alcohol. Heterologous overexpression of PcFPPS in tobacco plants was achieved, and it was found that the FPP enzyme activity was significantly up-regulated in transgenic tobacco by ELISA analysis. In addition, more terpenoid metabolites, including stigmasterol, phytol, and neophytadiene were detected compared with control by GC-MS analysis. Furthermore, with dual-LUC assay and yeast one-hybrid screening, we found 220 bp promoter of PcFPPS can be bound by the nuclear-localized transcription factor PcWRKY44. Overexpression of PcWRKY44 in P. cablin upregulated the expression levels of PcFPPS and patchoulol synthase gene (PcPTS), and then promote the biosynthesis of patchouli alcohol. Taken together, these results strongly suggest the PcFPPS and its binding transcription factor PcWRKY44 play an essential role in regulating the biosynthesis of patchouli alcohol.

Metabolic engineering to enhance the accumulation of bioactive flavonoids licochalcone A and echinatin in Glycyrrhiza inflata (Licorice) hairy roots

Echinatin and licochalcone A (LCA) are valuable chalcones preferentially accumulated in roots and rhizomes of licorice (Glycyrrhiza inflata). The licorice chalcones (licochalcones) are valued for their anti-inflammatory, antimicrobial, and antioxidant properties and have been widely used in cosmetic, pharmaceutical, and food industries. However, echinatin and LCA are accumulated in low quantities, and the biosynthesis and regulation of licochalcones have not been fully elucidated. In this study, we explored the potential of a R2R3-MYB transcription factor (TF) AtMYB12, a known regulator of flavonoid biosynthesis in Arabidopsis, for metabolic engineering of the bioactive flavonoids in G. inflata hairy roots. Overexpression of AtMYB12 in the hairy roots greatly enhanced the production of total flavonoids (threefold), echinatin (twofold), and LCA (fivefold). RNA-seq analysis of AtMYB12-overexpressing hairy roots revealed that expression of phenylpropanoid/flavonoid pathway genes, such as phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), and flavanone 3'-hydroxylase (F3'H), is significantly induced compared to the control. Transient promoter activity assay indicated that AtMYB12 activates the GiCHS1 promoter in plant cells, and mutation to the MYB-binding motif in the GiCHS1 promoter abolished activation. In addition, transcriptomic analysis revealed that AtMYB12 overexpression reprograms carbohydrate metabolism likely to increase carbon flux into flavonoid biosynthesis. Further, AtMYB12 activated the biotic defense pathways possibly by activating the salicylic acid and jasmonic acid signaling, as well as by upregulating WRKY TFs. The transcriptome of AtMYB12-overexpressing hairy roots serves as a valuable source in the identification of potential candidate genes involved in LCA biosynthesis. Taken together, our findings suggest that AtMYB12 is an effective gene for metabolic engineering of valuable bioactive flavonoids in plants.

Analysis of Tissue-Specific Defense Responses to Sclerotinia sclerotiorum in Brassica napus

Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum (S. sclerotiorum) is the main disease threat of oilseed rape (Brassica napus), resulting in huge economic losses every year. SSR resistance manifests as quantitative disease resistance (QDR), and no gene with complete SSR resistance has been cloned or reported so far. Transcriptome analysis has revealed a large number of defense-related genes and response processes. However, the similarities and differences in the defense responses of different tissues are rarely reported. In this study, we analyzed the similarities and differences of different tissues in response to S. sclerotiorum at 24 h post inoculation (hpi) by using the published transcriptome data for respective leaf and stem inoculation. At 24 hpi, large differences in gene expression exist in leaf and stem, and there are more differentially expressed genes and larger expression differences in leaf. The leaf is more sensitive to S. sclerotiorum and shows a stronger response than stem. Different defense responses appear in the leaf and stem, and the biosynthesis of lignin, callose, lectin, chitinase, PGIP, and PR protein is activated in leaf. In the stem, lipid metabolism-mediated defense responses are obviously enhanced. For the common defense responses in both leaf and stem, the chain reactions resulting from signal transduction and biological process take the primary responsibility. This research will be beneficial to exploit the potential of different tissues in plant defense and find higher resistance levels of genotypic variability in different environments. Our results are significant in the identification of resistance genes and analysis of defense mechanisms.

Integrated mRNA and Small RNA Sequencing Reveals microRNAs Associated With Xylem Development in Dalbergia odorifera

Dalbergia odorifera is a rare and precious rosewood specie, whose wood is a very high-quality material for valuable furniture and carving crafts. However, limited information is available about the process of wood formation in D. odorifera. To determine genes that might be closely associated with the xylem differentiation process, we analyzed the differentially expressed genes (DEGs) and microRNAs (miRNAs) from specific xylem tissues of D. odorifera by RNA sequencing (RNA-seq) and small RNA sequencing (small RNA-seq). In total, we obtained 134,221,955 clean reads from RNA-seq and 90,940,761 clean reads from small RNA-seq. By comparing the transition zone (Dotz) and sapwood (Dosw) samples, a total of 395 DEGs were identified. Further analysis revealed that DEGs encoded for WRKY transcription factors (eight genes), lignin synthesis (PER47, COMT, CCR2), cell wall composition (UXS2), gibberellin synthesis (KAO2, GA20OX1), jasmonic acid synthesis (OPR2, CYP74A), and synthesis of flavonoids (PAL2) and terpenoids (CYP71A1). Subsequently, a preliminary analysis by small RNA-seq showed that the expressions of 14 miRNAs (such as miR168a-5p, miR167f-5p, miR167h-5p, miR167e, miR390a, miR156g, novel_52, and novel_9) were significantly different between Dotz and Dosw. Further analysis revealed that the target genes of these differentially expressed miRNAs were enriched in the GO terms "amino acid binding," "cellulase activity," and "DNA beta-glucosyltransferase activity". Further, KEGG pathway annotation showed significant enrichment in "fatty acid elongation" and "biosynthesis of unsaturated fatty acids". These processes might be participating in the xylem differentiation of D. odorifera. Next, expression correlation analysis showed that nine differentially expressed miRNAs were significantly negatively associated with 21 target genes, which encoded for proteins such as pyrH, SPL6, SPL12, GCS1, and ARF8. Overall, this is the first study on miRNAs and their potential functions in the xylem development of D. odorifera, which provides a stepping stone for a detailed functional investigation of D. odorifera miRNAs.

SlWRKY35 positively regulates carotenoid biosynthesis by activating the MEP pathway in tomato fruit

Carotenoids are vital phytonutrients widely recognised for their health benefits. Therefore, it is vital to thoroughly investigate the metabolic regulatory network underlying carotenoid biosynthesis and accumulation to open new leads towards improving their contents in vegetables and crops. The outcome of our study defines SlWRKY35 as a positive regulator of carotenoid biosynthesis in tomato. SlWRKY35 can directly activate the expression of the 1-deoxy-d-xylulose 5-phosphate synthase (SlDXS1) gene to reprogramme metabolism towards the 2-C-methyl-d-erythritol 4-phosphate (MEP) pathway, leading to enhanced carotenoid accumulation. We also show that the master regulator SlRIN directly regulates the expression of SlWRKY35 during tomato fruit ripening. Compared with the SlLCYE overexpression lines, coexpression of SlWRKY35 and SlLCYE can further enhance lutein production in transgenic tomato fruit, indicating that SlWRKY35 represents a potential target towards designing innovative metabolic engineering strategies for carotenoid derivatives. In addition to providing new insights into the metabolic regulatory network associated with tomato fruit ripening, our data define a new tool for improving fruit content in specific carotenoid compounds.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Insight into the Phylogeny and Binding Ability of WRKY Transcription Factors

WRKY transcription factors (TFs), which make up one of the largest families of TFs in the plant kingdom, are key players in modulating gene expression relating to embryogenesis, senescence, pathogen resistance, and abiotic stress responses. However, the phylogeny and grouping of WRKY TFs and how their binding ability is affected by the flanking regions of W-box sequences remain unclear. In this study, we reconstructed the phylogeny of WRKY across the plant kingdom and characterized the DNA-binding profile of Arabidopsis thaliana WRKY (WRKY54) based on its W-box recognition sequence. We found that WRKY TFs could be separated into five clades, and that the functional zinc-finger motif at the C-terminal of WRKY appeared after several nucleotide substitutions had occurred at the 3′-end of the zinc-finger region in chlorophytes. In addition, we found that W-box flanking regions affect the binding ability of WRKY54 based on the results of a fluorescence-based electrophoretic mobility shift assay (fEMSA) and quartz crystal microbalance (QCM) analysis. The great abundance of WRKY TFs in plants implicates their involvement in diverse molecular regulatory networks, and the flanking regions of W-box sequences may contribute to their molecular recognition mechanism. This phylogeny and our findings on the molecular recognition mechanism of WRKY TFs should be helpful for further research in this area.

Comparative transcriptome mining for terpenoid biosynthetic pathway genes in wild and cultivated species of Plantago

Plantagos are important economical and medicinal plants that possess several bioactive secondary metabolites, such as phenolics, iridoids, triterpenes, and alkaloids. Triterpenoids are the ubiquitous and dynamic secondary metabolites that are deployed by plants for chemical interactions and protection under biotic/abiotic stress. Plantago ovata, a cultivated species, is the source of psyllium, while Plantago major, a wild species, has significant therapeutic potential. Wild species are considered more tolerant to stressful conditions in comparison to their cultivated allies. In view of this, the present study aimed to decipher the terpenoid biosynthetic pathway operative in P. ovata and P. major using a comparative transcriptomics approach. Majority of terpenoid biosynthetic genes were observed as upregulated in P. major including rate limiting genes of MVA (HMGR) and MEP (DXR) pathways and genes (α-AS, BAS, SM, and CYP716) involved in ursolic acid biosynthesis, an important triterpenoid prevalent in Plantago species. The HPLC output further confirmed the higher concentration of ursolic acid in P. major as compared to P. ovata leaf samples, respectively. In addition to terpenoid biosynthesis, KEGG annotation revealed the involvement of differentially expressed unigenes in several metabolic pathways, aminoacyl-tRNA biosynthesis, biosynthesis of antibiotics, and biosynthesis of secondary metabolites. MYB was found as the most abundant transcription factor family in Plantago transcriptome. We have been able to generate valuable information which can help in improving terpenoid production in Plantago. Additionally, the present study has laid a strong foundation for deciphering other important metabolic pathways in Plantago.

Grafting: a potential method to reveal the differential accumulation mechanism of secondary metabolites

Plant secondary metabolites make a great contribution to the agricultural and pharmaceutical industries. Their accumulation is determined by the integrated transport of target compounds and their biosynthesis-related RNA, protein, or DNA. However, it is hard to track the movement of these biomolecules in vivo. Grafting may be an ideal method to solve this problem. The differences in genetic and metabolic backgrounds between rootstock and scion, coupled with multiple omics approaches and other molecular tools, make it feasible to determine the movement of target compounds, RNAs, proteins, and DNAs. In this review, we will introduce methods of using the grafting technique, together with molecular biological tools, to reveal the differential accumulation mechanism of plant secondary metabolites at different levels. Details of the case of the transport of one diterpene alkaloid, fuziline, will be further illustrated to clarify how the specific accumulation model is shaped with the help of grafting and multiple molecular biological tools.

Transcription factor CsWRKY40 regulates L-theanine hydrolysis by activating the CsPDX2.1 promoter in tea leaves during withering

L-Theanine is a crucial secondary metabolite in tea and positively determines the potential quality and health benefits of tea products. Previous work found the content of L-theanine decreased during withering process, while the specific mechanism is still unknown. Here, weighted gene co-expression network analysis (WGCNA) was performed based on the transcriptome data obtained previously. The key hydrolysis gene CsPDX2.1 in L-theanine metabolism and seven candidate transcription factors were screened out. Among those transcription factors, CsWRKY40 presented the strongest activation on the CsPDX2.1 promoter (373.18-fold) by binding to W box element based on the dual luciferase assay and EMSA results. Meanwhile, CsWRKY40 protein was located in the nucleoplasm, while CsPDX2.1 was found in both the nucleoplasm and cytoplasm. Furthermore, it was confirmed that the water loss of tea leaves was the critical factor affecting the contents of ABA and L-theanine by activating the expression of CsPDX2.1 and CsPDX2.1 based on the analysis of the withering model, water-retention model and water-loss model. Our results provide a new insight into revealing the regulation mechanism of L-theanine hydrolysis metabolism.

Identification of a WRKY transcriptional activator from Camellia sinensis that regulates methylated EGCG biosynthesis

Naturally occurring methylated catechins, especially methylated EGCG in tea leaves are known to have many health benefits. Although the genes involved in methylated EGCG biosynthesis have been studied extensively, the transcriptional factors controlling methylated EGCG biosynthesis are still poorly understood. In the present study, a WRKY domain-containing protein termed CsWRKY57like was identified, which belongs to group IIc of the WRKY family, and contains one conserved WRKY motif. CsWRKY57like was found to localize in the nucleus, function as a transcriptional activator, and its expression positively correlated with methylated EGCG level. In addition, CsWRKY57like activated the transcription of three genes related to methylated EGCG biosynthesis, including CCoAOMT, CsLAR, and CsDFR by specifically interacting with their promoters via binding to the cis-acting element (C/T)TGAC(T/C). Further assays revealed that CsWRKY57like physically interacts with CsVQ4, and participates in the metabolic regulation of O-methylated catechin biosynthesis. Collectively, we conclude that CsWRKY57like may positively impact the biosynthesis of methylated EGCG in the tea plant, which comprehensively enriches the regulatory network of WRKY TFs associated with methylated EGCG and provide a potential strategy for the breeding of specific tea plant cultivars with high methylated EGCG .

PlantGF: an analysis and annotation platform for plant gene families

Gene families contain genes that come from the same ancestor and have similar sequences and structures. They perform certain specific functions within and among different species. Currently, there is no complete process or platform for the rapid analysis of plant gene families. In this study, a comprehensive query and analysis platform of plant gene families, the Plant Gene Family Platform (PlantGF), was constructed. The platform is composed of four main parts: Search, Tools, Statistics and Auxiliary. A total of 2 909 580 gene family members were identified from 138 plant species in PlantGF. The data can be queried in the Search section through a user-friendly interface. A general process for gene family analysis, having nine steps, is provided. The platform also includes four online tools (HMM-Search, BLAST, MAFFT and HMMER) in the Tools section for useful additional analyses. The statistical analysis of the relevant gene families is shown on the Statistics page. Auxiliary pages are provided for data downloading. The datasets for all 138 plant species' protein sequences and their gene families can be acquired on the Download page. A user's manual and some useful links are displayed on the Manual and Links pages, respectively. To the best of our knowledge, PlantGF is the first comprehensive platform for studying plant gene families, and it will make important contributions to plant gene family-related research. Database URL: http://biodb.sdau.edu.cn/PGF/index.html.

Comparative transcriptome analysis revealed the molecular mechanism of the effect of light intensity on the accumulation of rhynchophylline and isorhynchophylline in Uncaria rhynchophylla

Rhynchophylline (RIN) and isorhynchophylline (IRN), the main medicinal components in plant Uncaria rhynchophylla, have potential effects on Alzheimer's disease. Understanding the influence of environmental factors, especially light intensity, on the production of these active ingredients will help to improve cultivation techniques. Compared with the 100% light intensity (CK), the contents of RIN and IRN in U. rhynchophylla leaves significantly increased at 20% light intensity (HS) after 7 and 21 days. Short-term shading (21d) changed some morphological indicators of U. rhynchophylla, but did not affect its biomass. Transcriptome profile analysis was performed on data from two groups (7 and 21 days) of CK and HS samples and yielded 79,817 unigenes with an average length of 1023 bp. Concurrently, 2391 and 2136 differentially expressed genes were identified in the transcriptome data for, respectively, 7 and 21 days of shade treatment. Notably, unigenes known to be involved upstream in the biosynthesis of RIN and IRN, such as G8O, IO, 7-DLGT, LAMT, TDC, and STR, were mostly upregulated. In addition, 1065 putative transcription factors (TFs) were identified and grouped into 55 TF families; 26 TFs showed differential expression in the shade treatment after 7 and 21 days. HY5 and PIFs, two important TFs of the light signaling pathway, also showed differential expression. This study provides insight into how gene expression was affected by light intensity during RIN and IRN accumulation in U. rhynchophylla.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01142-2.

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.

Transcriptome Analysis Reveals Differentially Expressed Genes That Regulate Biosynthesis of the Active Compounds with Methyl Jasmonate in Rosemary Suspension Cells

To study the effects of Methyl jasmonates (MeJA) on rosemary suspension cells, the antioxidant enzymes’ change of activities under different concentrations of MeJA, including 0 (CK), 10 (M10), 50 (M50) and 100 μM MeJA (M100). The results demonstrated that MeJA treatments increased the activities of phenylalanine ammonla-lyase (PAL), superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) and polyphenol oxidase (PPO) and reduced the contents of hydrogen peroxide (H₂O₂) and malondialdehyde (MDA), thus accelerating the ROS scavenging. Comparative transcriptome analysis of different concentrations of MeJA showed that a total of 7836, 6797 and 8310 genes were differentially expressed in the comparisons of CKvsM10, CKvsM50, CKvsM100, respectively. The analysis of differentially expressed genes (DEGs) showed phenylpropanoid biosynthesis, vitamin B6, ascorbate and aldarate metabolism-related genes were significantly enriched. The transcripts of flavonoid and terpenoid metabolism pathways and plant hormone signal transduction, especially the jasmonic acid (JA) signal-related genes, were differentially expressed in CKvsM50 and CKvsM100 comparisons. In addition, the transcription factors (TFs), e.g., MYC2, DELLA, MYB111 played a key role in rosemary suspension cells under MeJA treatments. qRT-PCR of eleven DEGs showed a high correlation between the RNA-seq and the qRT-PCR result. Taken together, MeJA alleviated peroxidative damage of the rosemary suspension cells in a wide concentration range via concentration-dependent differential expression patterns. This study provided a transcriptome sequence resource responding to MeJA and a valuable resource for the genetic and genomic studies of the active compounds engineering in rosemary.

CaWRKY30 Positively Regulates Pepper Immunity by Targeting CaWRKY40 against Ralstonia solanacearum Inoculation through Modulating Defense-Related Genes

The WRKY transcription factors (TFs) network is composed of WRKY TFs’ subset, which performs a critical role in immunity regulation of plants. However, functions of WRKY TFs’ network remain unclear, particularly in non-model plants such as pepper (Capsicum annuum L.). This study functionally characterized CaWRKY30—a member of group III Pepper WRKY protein—for immunity of pepper against Ralstonia solanacearum infection. The CaWRKY30 was detected in nucleus, and its transcriptional expression levels were significantly upregulated by R. solanacearum inoculation (RSI), and foliar application ethylene (ET), abscisic acid (ABA), and salicylic acid (SA). Virus induced gene silencing (VIGS) of CaWRKY30 amplified pepper’s vulnerability to RSI. Additionally, the silencing of CaWRKY30 by VIGS compromised HR-like cell death triggered by RSI and downregulated defense-associated marker genes, like CaPR1, CaNPR1, CaDEF1, CaABR1, CaHIR1, and CaWRKY40. Conversely, transient over-expression of CaWRKY30 in pepper leaves instigated HR-like cell death and upregulated defense-related maker genes. Furthermore, transient over-expression of CaWRKY30 upregulated transcriptional levels of CaWRKY6, CaWRKY22, CaWRKY27, and CaWRKY40. On the other hand, transient over-expression of CaWRKY6, CaWRKY22, CaWRKY27, and CaWRKY40 upregulated transcriptional expression levels of CaWRKY30. The results recommend that newly characterized CaWRKY30 positively regulates pepper’s immunity against Ralstonia attack, which is governed by synergistically mediated signaling by phytohormones like ET, ABA, and SA, and transcriptionally assimilating into WRKY TFs networks, consisting of CaWRKY6, CaWRKY22, CaWRKY27, and CaWRKY40. Collectively, our data will facilitate to explicate the underlying mechanism of crosstalk between pepper’s immunity and response to RSI.

Nanoparticles as elicitors and harvesters of economically important secondary metabolites in higher plants: A review

Nanoparticles possess some unique properties which improve their biochemical reactivity. Plants, due to their stationary nature, are constantly exposed to nanoparticles present in the environment, which act as abiotic stress agents at sub-toxic concentrations and phytotoxic agents at higher concentrations. In general, nanoparticles exert their toxicological effect by the generation of reactive oxygen species to which plants respond by activating both enzymatic and non-enzymatic anti-oxidant defence mechanisms. One important manifestation of the defence response is the increased or de novo biosynthesis of secondary metabolites, many of which have commercial application. The present review extensively summarizes current knowledge about the application of different metallic, non-metallic and carbon-based nanoparticles as elicitors of economically important secondary metabolites in different plants, both in vivo and in vitro. Elicitation of secondary metabolites with nanoparticles in plant cultures, including hairy root cultures, is discussed. Another emergent technology is the ligand-harvesting of secondary metabolites using surface-functionalized nanoparticles, which is also mentioned. A brief explanation of the mechanism of action of nanoparticles on plant secondary metabolism is included. Optimum conditions and parameters to be evaluated and standardized for the successful commercial exploitation of this technology are also mentioned.

GhWRKY1-like, a WRKY transcription factor, mediates drought tolerance in Arabidopsis via modulating ABA biosynthesis

BACKGROUND

Drought stress has great negative effects on the plant growth and development. The tolerance of plants to such abiotic stress is triggered by complicated and multilayered signaling pathways to restore cellular homeostasis and to promote survival. The WRKY family is one of the largest transcription factor families in higher plants, and has been well recognized for the roles in regulating plants tolerance to abiotic and biotic stress. However, little is known about how the WRKY genes regulate drought resistance in cotton.

RESULTS

In this work, we identified the WRKY transcription factor GhWRKY1-like from upland cotton as a positive regulator of tolerance to drought that directly manipulates abscisic acid (ABA) biosynthesis. Overexpression of GhWRKY1-like in Arabidopsis constitutively activated ABA biosynthesis genes, signaling genes, responsive genes and drought related maker genes, and led to enhanced tolerance to drought. Further analysis has shown that GhWRKY1-like can interact with "W-box" cis-elements of the promoters of AtNCED2, AtNCED5, AtNCED6 and AtNCED9 which are essential enzymes for ABA biosynthesis, and promotes the expression of those target genes.

CONCLUSIONS

In summary, our findings suggest that GhWRKY1-like may act as a positive regulator in Arabidopsis tolerance to drought via directly interacting with the promoters of AtNCED2, AtNCED5, AtNCED6 and AtNCED9 to promote ABA biosynthesis.

Constitutive Defense Strategy of Coffee Under Field Conditions: A Comparative Assessment of Resistant and Susceptible Cultivars to Rust

Coffea arabica is the most economically important coffee species worldwide. However, its production is severely limited by diseases such as rust. The mechanisms underlying constitutive defense responses in coffee are still poorly understood, compared with induced defense mechanisms. We aimed to characterize constitutive defense responses of thirteen cultivars of C. arabica. Cultivars were classified under field conditions according to the level of resistance to rust: resistant (R), moderately resistant (MR), and susceptible (S). Based on this classification, the stability of eight reference genes (RGs) was evaluated. The most stable RGs were EF1α, APT1, and 24S. We also evaluated the expression of CaWRKY1, CaPAL1, CaCAD1, and CaPOX1, and activities of PAL, CAD, and POX, which are involved in lignin biosynthesis, and leaf content of total phenolic compounds and lignin. Gene expression and enzymatic activity were not correlated with defense metabolites in the R cultivar group but showed a negative correlation with phenolic compounds in MR cultivars. Cultivar S showed positive correlations of gene expression and enzyme activity with phenolic compounds. These results may assist coffee breeding programs regarding selection of genotypes and in optimization of rust resistance.

Integration of genome wide association studies and co-expression networks reveal roles of PtoWRKY 42-PtoUGT76C1-1 in trans-zeatin metabolism and cytokinin sensitivity in poplar

Cytokinins are important for in vitro shoot regeneration in plants. Cytokinin N-glucosides are produced via an irreversible glycosylation pathway, which regulates the endogenous cytokinin content. Although cytokinin N-glucoside pathways have been uncovered in higher plants, no regulator has been identified to date. We performed a metabolome genome-wide association study (mGWAS), weighted gene co-expression network analysis (WGCNA), and expression quantitative trait nucleotide (eQTN) mappings to build a core triple genetic network (mGWAS-gene expression-phenotype) for the trans-zeatin N-glucoside (ZNG) metabolite using data from 435 unrelated Populus tomentosa individuals. Variation of the ZNG level in poplar is attributed to the differential transcription of PtoWRKY42, a member of WRKY multigene family group IIb. Functional analysis revealed that PtoWRKY42 negatively regulated ZNG accumulation by binding directly to the W-box of the UDP-glycosyltransferase 76C 1-1 (PtoUGT761-1) promoter. Also, PtoWRKY42 was strongly induced by leaf senescence, 6-BA, wounding, and salt stress, resulting in a reduced ZNG level. We identified PtoWRKY42, a negative regulator of cytokinin N-glucosides, which contributes to the natural variation in ZNG level and mediates ZNG accumulation by directly modulating the key glycosyltransferase gene PtoUGT76C1-1.

An update on the function and regulation of methylerythritol phosphate and mevalonate pathways and their evolutionary dynamics

Isoprenoids are among the largest and most chemically diverse classes of organic compounds in nature and are involved in the processes of photosynthesis, respiration, growth, development, and plant responses to stress. The basic building block units for isoprenoid synthesis-isopentenyl diphosphate and its isomer dimethylallyl diphosphate-are generated by the mevalonate (MVA) and methylerythritol phosphate (MEP) pathways. Here, we summarize recent advances on the roles of the MEP and MVA pathways in plant growth, development and stress responses, and attempt to define the underlying gene networks that orchestrate the MEP and MVA pathways in response to developmental or environmental cues. Through phylogenomic analysis, we also provide a new perspective on the evolution of the plant isoprenoid pathway. We conclude that the presence of the MVA pathway in plants may be associated with the transition from aquatic to subaerial and terrestrial environments, as lineages for its core components are absent in green algae. The emergence of the MVA pathway has acted as a key evolutionary event in plants that facilitated land colonization and subsequent embryo development, as well as adaptation to new and varied environments.

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.

Transcription factors as key molecular target to strengthen the drought stress tolerance in plants

Amid apprehension of global climate change, crop plants are inevitably confronted with a myriad of abiotic stress factors during their growth that inflicts a serious threat to their development and overall productivity. These abiotic stresses comprise extreme temperature, pH, high saline soil, and drought stress. Among different abiotic stresses, drought is considered the most calamitous stressor with its serious impact on the crops' yield stability. The development of climate-resilient crops that withstands reduced water availability is a major focus of the scientific fraternity to ensure the food security of the sharply increasing population. Numerous studies aim to recognize the key regulators of molecular and biochemical processes associated with drought stress tolerance response. A few potential candidates are now considered as promising targets for crop improvement. Transcription factors act as a key regulatory switch controlling the gene expression of diverse biological processes and, eventually, the metabolic processes. Understanding the role and regulation of the transcription factors will facilitate the crop improvement strategies intending to develop and deliver agronomically-superior crops. Therefore, in this review, we have emphasized the molecular avenues of the transcription factors that can be exploited to engineer drought tolerance potential in crop plants. We have discussed the molecular role of several transcription factors, such as basic leucine zipper (bZIP), dehydration responsive element binding (DREB), DNA binding with one finger (DOF), heat shock factor (HSF), MYB, NAC, TEOSINTE BRANCHED1/CYCLOIDEA/PCF (TCP), and WRKY. We have also highlighted candidate transcription factors that can be used for the development of drought-tolerant crops.

A chloroplast heat shock protein modulates growth and abiotic stress response in creeping bentgrass

Small heat shock proteins (sHSPs), a family of the ubiquitous stress proteins in plants acting as molecular chaperones to protect other proteins from stress-induced damage, have been implicated in plant growth and development as well as plant response to environmental stress, especially heat stress. In this study, a chloroplast-localized sHSP, AsHSP26.8, was overexpressed in creeping bentgrass (Agrostis stolonifera L.) to study its role in regulating plant growth and stress response. Transgenic (TG) creeping bentgrass plants displayed arrested root development, slow growth rate, twisted leaf blades and are more susceptible to heat and salt but less sensitive to drought stress compared to wild-type (WT) controls. RNA-seq analysis revealed that AsHSP26.8 modulated the expression of genes in auxin signalling and stress-related genes such as those encoding HSPs, heat shock factors and other transcription factors. Our results provide new evidence demonstrating that AsHSP26.8 negatively regulates plant growth and development and plays differential roles in plant response to a plethora of diverse abiotic stresses.

The bZIP transcription factor GmbZIP15 facilitates resistance against Sclerotinia sclerotiorum and Phytophthora sojae infection in soybean

Soybean, one of the most valuable oilseed crops, is under constant pressure from pathogens. bZIP transcription factors (TFs) composing one of the largest TF families in plants have diverse functions. Biochemical and physiological analyses were performed to characterize the regulatory roles of soybean bZIP TF GmbZIP15 in response to pathogens. We found that transgenic soybean plants overexpressing GmbZIP15 has increased resistance against Sclerotinia sclerotiorum and Phytophthora sojae. Besides, GmbZIP15 regulates pathogen response by modulating the antioxidant defense system and phytohormone signaling. In addition, we performed chromatin immunoprecipitation sequencing to identify the downstream genes of GmbZIP15 in response to S. sclerotiorum and found that GmbZIP15 can activate or repress the expression of defense-related genes through direct promoter binding. Taken together, these results indicate that GmbZIP15 plays a positive role in pathogen resistance in soybean, and this activity may be dependent on phytohormone signaling.

•

GmbZIP15 improves resistance against pathogen

•

GmbZIP15 modulates the antioxidant defense system

•

GmbZIP15 regulates phytohormone signaling

•

GmbZIP15 can direct bind to G-box

Transcriptome analysis reveals gene responses to herbicide, tribenuron methyl, in Brassica napus L. during seed germination

Background

Tribenuron methyl (TBM) is an herbicide that inhibits sulfonylurea acetolactate synthase (ALS) and is one of the most widely used broad-leaved herbicides for crop production. However, soil residues or drifting of the herbicide spray might affect the germination and growth of rapeseed, Brassica napus, so it is imperative to understand the response mechanism of rape to TBM during germination. The aim of this study was to use transcriptome analysis to reveal the gene responses in herbicide-tolerant rapeseed to TBM stress during seed germination.

Results

2414, 2286, and 1068 differentially expressed genes (DEGs) were identified in TBM-treated resistant vs sensitive lines, treated vs. control sensitive lines, treated vs. control resistant lines, respectively. GO analysis showed that most DEGs were annotated to the oxidation-reduction pathways and catalytic activity. KEGG enrichment was mainly involved in plant-pathogen interactions, α-linolenic acid metabolism, glucosinolate biosynthesis, and phenylpropanoid biosynthesis. Based on GO and KEGG enrichment, a total of 137 target genes were identified, including genes involved in biotransferase activity, response to antioxidant stress and lipid metabolism. Biotransferase genes, CYP450, ABC and GST, detoxify herbicide molecules through physical or biochemical processes. Antioxidant genes, RBOH, WRKY, CDPK, MAPK, CAT, and POD regulate plant tolerance by transmitting ROS signals and triggering antioxidant enzyme expression. Lipid-related genes and hormone-related genes were also found, such as LOX3, ADH1, JAZ6, BIN2 and ERF, and they also played an important role in herbicide resistance.

Conclusions

This study provides insights for selecting TBM-tolerant rapeseed germplasm and exploring the molecular mechanism of TBM tolerance during germination.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07614-1.

Transcriptome-wide analysis of introgression-resistant regions reveals genetic divergence genes under positive selection in Populus trichocarpa

Comparing gene expression patterns and genetic polymorphisms between populations is of central importance for understanding the origin and maintenance of biodiversity. Based on population-specific gene expression levels and allele frequency differences, we sought to identify population divergence (PD) genes across the introgression-resistant genomic regions of Populus trichocarpa. Genes containing highly diverged loci [i.e., genetic divergence (GD)] or showing expression divergence (ED) between populations were widely distributed in the genome and substantially enriched in functional categories related to stress responses, disease resistance, timing of flowering, cell cycle regulation, plant growth, and development. Nine genomic regions showing evidence of strong positive selection were overlapped with GD genes, which had significant differences between Oregon (a southernmost peripheral deme) and the other demes. However, we did not find evidence that genes under positive selection show an enrichment for ED. PD genes and genes under selection pertained to the same gene classes, such as SERINE/CYSTEINE PROTEASE, ABC TRANSPORTER, GLYCOSYLTRANSFERASE and other transferases. Our analysis also revealed that GD genes were polymorphic within the species (41.9 ± 3.66 biallelic variants per gene), as previously reported in herbaceous plants. By contrast, ED genes contained less genetic variants (10.73 ± 1.14) and were likely highly expressed. In addition, we found that trans- rather than cis-acting variants considerably contribute to the evolution of >90% PD genes. Overall, this study elucidates that cohorts of PD genes agree with the general attributes of known speciation genes and GD genes will provide substrates for positive selection to operate on.

A Novel WRKY Transcription Factor HmoWRKY40 Associated with Betalain Biosynthesis in Pitaya (Hylocereus monacanthus) through Regulating HmoCYP76AD1

Betalains are water-soluble nitrogen-containing pigments with multiple bioactivities. Pitaya is the only large-scale commercially grown fruit containing abundant betalains for consumers. However, the upstream regulators in betalain biosynthesis are still not clear. In this study, HmoWRKY40, a novel WRKY transcription factor, was obtained from the transcriptome data of pitaya (Hylocereus monacanthus). HmoWRKY40 is a member of the Group IIa WRKY family, containing a conserved WRKY motif, and it is located in the nucleus. The betalain contents and expression levels of HmoWRKY40 increased rapidly during the coloration of pitaya and reached their maximums on the 23rd day after artificial pollination (DAAP). Yeast one-hybrid and transient expression assays showed that HmoWRKY40 could bind and activate the promoter of HmoCYP76AD1. Silencing the HmoWRKY40 gene resulted in a significant reduction of betacyanin contents. These results indicate that HmoWRKY40 transcriptionally activates HmoCYP76AD, which is involved in the regulation of pitaya betalain biosynthesis. The results of the present study provide new regulatory networks related to betalain biosynthesis in pitaya.

Transcriptomics Integrated with Free and Bound Terpenoid Aroma Profiling during "Shine Muscat" (Vitis labrusca × V. vinifera) Grape Berry Development Reveals Coordinate Regulation of MEP Pathway and Terpene Synthase Gene Expression

Terpenes and their derivatives are important biomarkers of grape quality as they contribute to the flavor and aroma of grapes. However, the molecular basis of terpene biosynthesis throughout the grapevine phenological developmental cycle remains elusive. Our current study investigates the free and bound terpene biosynthesis of berries at different phenological stages from preveraison to harvest. Detailed gene expression (transcriptomics) analysis, terpenoid volatile production by gas chromatography-mass spectrometry (GC-MS), and in planta transient expression were employed. Our results show that concentrations of most individual terpenes at different stages are distinctive and increase from preveraison to the veraison stage followed by a decrease from veraison to maturity. The combined transcriptomic analysis and terpene profiling revealed that 22 genes belonging to the MEP pathway and multiple classes of transcription factor family members including bHLH and several hormone biosynthesis- or signaling-related genes likely participate in the regulation of terpenoid biosynthesis according to their specific expression patterns in berries. Quantitative real-time polymerase chain expression analysis of 8 key differentially expressed genes in MEP pathways and further 12 randomly selected genes was performed during 8 sampling stages and validated the RNA-seq-derived expression profiles. To further confirm the function of a subset of the differentially expressed genes, we investigated the effects of combined overexpression of 1-deoxy-d-xylulose-5-phosphate synthase (VvDXS1-LOC100249323), 1-deoxy-d-xylulose-5-phosphate reductoisomerase (VvDXR-LOC100248516), and terpene synthase (VvTPS56-LOC100266449) on the production of terpenes by transient overexpression in Nicotiana benthamiana leaves. The overall developmental patterns of total terpenes and gene expression profiles will help guide the functional analyses of further candidate genes important for terpene biosynthesis of grape as well as identifying the master transcriptional and hormonal regulators of this pathway in the future.

Transcriptome-wide identification of WRKY family genes and their expression profiling toward salicylic acid in Camellia japonica

The ornamental plant Camellia japonica is widely distributed worldwide and is susceptible to various environmental stresses. The WRKY transcription factor (TF) is an important node of plant tolerance. However, WRKY TFs from C. japonica have not been reported yet. In this study, 48 CjWRKYs, namely, CjWRKY1 to CjWRKY48, were identified. Protein structure analysis revealed that CjWRKY proteins contain a highly conserved motif (WRKYGQK) and two variant motifs (WRKYGKK and WRKYGRK). Phylogenetic analysis indicated that the 48 CjWRKYs can be divided into three groups, which are further classified into six subgroups, namely, I-C, II-a, II-b, II-c, II-e, and III, which contain 10, 6, 8, 13, 7, and 4 members, respectively. The expression patterns of 15 CjWRKYs under salicylic acid (SA) treatment were investigated by real-time quantitative PCR (qRT-PCR). Results showed that the 15 CjWRKYs could be induced by SA treatment. This study is the first to screen CjWRKYs and identify the expression profile of CjWRKYs under SA treatment and provides a theoretical basis for analyzing the function of CjWRKY genes to SA stress tolerance in C. japonica.

HbWRKY82, a novel IIc WRKY transcription factor from Hevea brasiliensis associated with abiotic stress tolerance and leaf senescence in Arabidopsis

WRKY group transcription factors of model plants and major crops are confirmed to play essential roles in stress responses, senescence, secondary metabolism processes and hormone signal transduction. Previous studies have identified 81 HbWRKY genes from Hevea brasiliensis (the Pará rubber tree), but the functions of HbWRKYs in response to abiotic stresses and leaf senescence are unclear. In this study, one novel group IIc WRKY transcription factor named HbWRKY82 was identified and characterized as a stress-associated WRKY in rubber tree. Transient expression and transcriptional activation analyses indicated that HbWRKY82 encoded a nuclear protein and functioned as a transcription activator. The transcription levels of HbWRKY82 were induced by exogenous Ethrel (ET) (ethylene releaser) and abscisic acid (ABA) stimulations, down-regulated in tapping panel dryness rubber trees, and also exhibits significant decrease during the progression of leaf senescence. Overexpression of HbWRKY82 in Arabidopsis improved the tolerance to dehydration and salinity, and decreased the sensitivity to exogenous ABA. Moreover, real-time quantitative PCR analysis demonstrated that HbWRKY82 regulated the transcriptional expression of several stress-responsive genes (DREB1A, ERD10, HKT1, P5CS, RD22, RD29B, SKOR), leaf senescence marker genes (EIN3, WRKY53, NAP), ROS-related genes (RbohD, CSD1, CSD2, FSD3) and hormone signaling genes (EIN3, ABF3, ABF4). Collectively, our findings suggested that HbWRKY82 might function as an important transcriptional regulator in ET- and ABA-mediated leaf senescence and abiotic stress responses, and also be involved in tapping panel dryness, latex flow and regeneration processes of rubber trees via participating in the ET and reactive oxygen species signaling pathways.

Strawberry FaWRKY25 Transcription Factor Negatively Regulated the Resistance of Strawberry Fruits to Botrytis cinerea

WRKY genes and jasmonic acid (JA) play a crucial role in plants’ responses against biotic and abiotic stress. However, the regulating mechanism of WRKY genes on strawberry fruits’ resistance against Botrytis cinerea is largely unknown, and few studies have been performed on their effect on the JA-mediated defense mechanism against B. cinerea. This study explored the effect of FaWRKY25 on the JA-mediated strawberry resistance against B. cinerea. Results showed that the JA content decreased significantly as the fruits matured, whereas the FaWRKY25 expression rose substantially, which led to heightened susceptibility to B. cinerea and in strawberries. External JA treatment significantly increased the JA content in strawberries and reduced the FaWRKY25 expression, thereby enhancing the fruits’ resistance against B. cinerea. FaWRKY25 overexpression significantly lowered the fruits’ resistance against B. cinerea, whereas FaWRKY25 silencing significantly increased resistance. Moreover, FaWRKY25 overexpression significantly lowered the JA content, whereas FaWRKY25 silencing significantly increased it. FaWRKY25 expression level substantially affects the expression levels of genes related to JA biosynthesis and metabolism, other members of the WRKY family, and defense genes. Accordingly, FaWRKY25 plays a crucial role in regulating strawberries’ resistance against B. cinerea and may negatively regulate their JA-mediated resistance mechanism against B. cinerea.

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.

A WRKY transcription factor, TaWRKY40-D, promotes leaf senescence associated with jasmonic acid and abscisic acid pathways in wheat

Leaf senescence is a complex and precise regulatory process that is correlated with numerous internal and environmental factors. Leaf senescence is tightly related to the redistribution of nutrients, which significantly affects productivity and quality, especially in crops. Evidence shows that the mediation of transcriptional regulation by WRKY transcription factors is vital for the fine-tuning of leaf senescence. However, the underlying mechanisms of the involvement of WRKY in leaf senescence are still unclear in wheat. Using RNA sequencing data, we isolated a novel WRKY transcription factor, TaWRKY40-D, which localizes in the nucleus and is basically induced by the progression of leaf senescence. TaWRKY40-D is a promoter of natural and dark-induced leaf senescence in transgenic Arabidopsis thaliana and wheat. We also demonstrated a positive response of TaWRKY40-D in wheat upon jasmonic acid (JA) and abscisic acid (ABA) treatment. Consistent with this, the detached leaves of TaWRKY40-D VIGS (virus-induced gene silencing) wheat plants showed a stay-green phenotype, while TaWRKY40-D overexpressing Arabidopsis plants showed premature leaf senescence after JA and ABA treatment. Moreover, our results revealed that TaWRKY40-D positively regulates leaf senescence, possibly by altering the biosynthesis and signalling of JA and ABA pathway genes. Together, our results suggest a new regulator of JA- and ABA-related leaf senescence, as well as a new candidate gene that can be used for molecular breeding in wheat.

Molecular characterization and functional analysis of a novel WRKY transcription factor HbWRKY83 possibly involved in rubber production of Hevea brasiliensis

WRKY transcription factors play important roles in plant growth and developmental processes and various stress responses, and are also associated with jasmonic acid (JA) signaling in the regulation of secondary metabolite biosynthesis in plants. The regulatory networks mediated by WRKY proteins in the latex production of Hevea brasiliensis (the Pará rubber tree) are poorly understood. In this study, one novel WRKY gene (designated HbWRKY83) was identified from the latex of H. brasiliensis, and its functions were characterized via gene expression analysis in both the latex and HbWRKY83-overexpressing transgenic Arabidopsis. HbWRKY83 gene contains an open reading frame (ORF) of 921 bp encoding a 306-amino-acid protein which is clustered with group IIc WRKY TF. HbWRKY83 is a nuclear-localized protein with transcriptional activity. Real-time quantitative PCR analysis demonstrated that the transcription level of HbWRKY83 was up-regulated by exogenous methyl jasmonate, Ethrel (ethylene releaser) stimulation, and bark tapping (mechanical wounding). Compared with the wild-type plants, overexpression of HbWRKY83 improved the tolerance of transgenic Arabidopsis lines to drought and salt stresses by enhancing the expression levels of ethylene-insensitive3 transcription factors (EIN3s) and several stress-responsive genes, including Cu/Zn superoxide dismutases CSD1 (Cu/Zn-SOD1) and CSD2 (Cu/Zn-SOD2), related to reactive oxygen species scavenging. Additionally, these genes were also significantly up-regulated by bark tapping. In combination, these results suggest that HbWRKY83 might act as a positive regulator of rubber production by activating the expression of JA-, ethylene-, and wound-responsive genes in the laticiferous cells of rubber trees.

Maleic hydrazide elicits global transcriptomic changes in chemically topped tobacco to influence shoot bud development

MAIN CONCLUSION

Transcriptomic analysis revealed maleic hydrazide suppresses apical and axillary bud development by altering the expression of genes related to meristem development, cell division, DNA replication, DNA damage and recombination, and phytohormone signaling. Topping (removal of apical buds) is a common agricultural practice for some crop plants including cotton, cannabis, and tobacco. Maleic hydrazide (MH) is a systemic suckercide, a chemical that inhibits shoot bud growth, used to control the growth of apical (ApB) and axillary buds (AxB) following topping. However, the influence of MH on gene expression and the underlying molecular mechanism of controlling meristem development are not well studied. Our RNA sequencing analysis showed that MH significantly influences the transcriptomic landscape in ApB and AxB of chemically topped tobacco. Gene ontology (GO) enrichment analysis revealed that upregulated genes in ApB were enriched for phosphorelay signal transduction, and the regulation of transition timing from vegetative to reproductive phase, whereas downregulated genes were largely associated with meristem maintenance, cytokinin metabolism, cell wall synthesis, photosynthesis, and DNA metabolism. In MH-treated AxB, GO terms related to defense response and oxylipin metabolism were overrepresented in upregulated genes. GO terms associated with cell cycle, DNA metabolism, and cytokinin metabolism were enriched in downregulated genes. Expression of KNOX and MADS transcription factor (TF) family genes, known to be involved in meristem development, were affected in ApB and AxB by MH treatment. The promoters of MH-responsive genes are enriched for several known cis-acting elements, suggesting the involvement of a subset of TF families. Our findings suggest that MH affects shoot bud development in chemically topped tobacco by altering the expression of genes related to meristem development, DNA repair and recombination, cell division, and phytohormone signaling.

Overexpression of transcription factor SlWRKY28 improved the tolerance of Populus davidiana × P. bolleana to alkaline salt stress

BACKGROUND

WRKY transcription factors (TFs) have been suggested to play crucial roles in the response to biotic and abiotic stresses. This study is the first to report the alkaline salt regulation of the WRKY gene.

RESULTS

In this study, we cloned a WRKY gene (SlWRKY28) from the Salix linearistipularis and then transferred to the Populus davidiana × P. bolleana for expression. Sequence analysis on the transcriptome of Salix linearistipular showed the significant up-regulation of WRKY gene expression in response to salt-alkali stress in seedlings. Our data showed that SlWRKY28 localized to the nucleus. Furthermore, the expression of the SlWRKY28 from female plants increased with saline-alkali stress according to the northern blot analysis results. The results of 3,3'-Diaminobenzidine (DAB) staining showed that hydrogen peroxide (H2O2) concentration was lower under stress, but ascorbate peroxidase (APX) enzyme activity was significantly higher in the overexpressed plants than that in non-transgenic (NT) plants.

CONCLUSIONS

We found out the SlWRKY28 induced regulation of the enzyme gene in the reactive oxygen species (ROS) scavenging pathway is a potential mechanism for transgenic lines to improve their resistance to alkaline salt. This study shows theoretical and practical significance in determining SlWRKY28 transcription factors involved in the regulation of alkaline salt tolerance.

PkGPPS.SSU interacts with two PkGGPPS to form heteromeric GPPS in Picrorhiza kurrooa: Molecular insights into the picroside biosynthetic pathway

Geranyl geranyl pyrophosphate synthase (GGPPS) is known to form an integral subunit of the heteromeric GPPS (geranyl pyrophosphate synthase) complex and catalyzes the biosynthesis of monoterpene in plants. Picrorhiza kurrooa Royle ex Benth., a medicinally important high altitude plant is known for picroside biomolecules, the monoterpenoids. However, the significance of heteromeric GPPS in P. kurrooa still remains obscure. Here, transient silencing of PkGGPPS was observed to reduce picroside-I (P-I) content by more than 60% as well as picroside-II (P-II) by more than 75%. Thus, PkGGPPS was found to be involved in the biosynthesis of P-I and P-II besides other terpenoids. To unravel the mechanism, small subunit of GPPS (PkGPPS.SSU) was identified from P. kurrooa. Protein-protein interaction studies in yeast as well as bimolecular fluorescence complementation (BiFC) in planta have indicated that large subunit of GPPS PkGPPS.LSUs (PkGGPPS1 and PkGGPPS2) and PkGPPS.SSU form a heteromeric GPPS. Presence of similar conserved domains such as light responsive motifs, low temperature responsive elements (LTRE), dehydration responsive elements (DREs), W Box and MeJA responsive elements in the promoters of PkGPPS.LSU and PkGPPS.SSU documented their involvement in picroside biosynthesis. Further, the tissue specific transcript expression analysis vis-à-vis epigenetic regulation (DNA methylation) of promoters as well as coding regions of PkGPPS.LSU and PkGPPS.SSU has strongly suggested their role in picroside biosynthesis. Taken together, the newly identified PkGPPS.SSU formed the heteromeric GPPS by interacting with PkGPPS.LSUs to synthesize P-I and P-II in P. kurrooa.

Origin and evolution of jasmonate signaling

Jasmonate (JA) signaling is a key mediator of plant development and defense which arose during plants transition from an aqueous to terrestrial environment. Elucidating the evolution of JA signaling is important for understanding plant development, defense, and production of specialized metabolites. The lineage of key protein domains characterizing JA signaling factors was traced to identify the origins of CORONITINE INSENSITIVE 1 (COI1), JASMONATE ZIM-DOMAIN (JAZ), NOVEL INTERACTOR OF JAZ, MYC2, TOPLESS, and MEDIATOR SUBUNIT 25. Charophytes do not possess genes encoding key JA signaling components, including COI1, JAZ, MYC2, and the JAZ-interacting bHLH factors, yet their orthologs are present in bryophytes. TIFY family genes were found in charophyta and chlorophya algae. JAZs evolved from ZIM genes of the TIFY family through changes to several key amino acids. Dating placed the origin of JA signaling 515 to 473 million years ago during the middle Cambrian to early Ordovician periods. This time is known for rapid biodiversification and mass extinction events. An increased predation from the diversifying and changing fauna may have driven evolution of JA signaling and plant defense.

Two light responsive WRKY genes exhibit positive and negative correlation with picroside content in Picrorhiza kurrooa Royle ex Benth, an endangered medicinal herb

Picrorhiza kurrooa is an endangered herb known to produce the medicinally important picrosides through isoprenoid pathway. The present work showed the functionality of WRKY motifs (TGAC cis-acting elements) present in the promoters of regulatory genes 3-hydroxy-3-methylglutaryl coenzyme A reductase (Pkhmgr) and 1-deoxy-d-xylulose-5-phosphate synthase (Pkdxs) of the picrosides biosynthetic pathway by electrophoretic mobility shift assay. Also, the two WRKY genes, PkdWRKY and PksWRKY, were characterized and found to contain double and single characteristic WRKY domains, respectively along with a zinc-finger motif in each domain. Expression analysis revealed that PkdWRKY and PksWRKY exhibited a positive and negative correlation, respectively, with picrosides content under the environment of light and in different tissues. Functional evaluation in yeast showed DNA binding ability of both PksWRKY and PkdWRKY; however, only PkdWRKY exhibited transcriptional activation ability. Transient overexpression of PkdWRKY and PksWRKY in tobacco modulated the expression of selected native genes of tobacco involved in MVA and MEP pathway suggesting functionality of PkdWRKY and PksWRKY in planta. Collectively, data suggested that PkdWRKY and PksWRKY might be positive and negative regulators, respectively in the picrosides biosynthetic pathway.