A sunflower WRKY transcription factor stimulates the mobilization of seed-stored reserves during germination and post-germination growth

Differential metabolic networks in three energy substances of flaxseed (Linum usitatissimum L.) during germination

Germinated flaxseed (Linum usitatissimum L.) is an essential potential food ingredient, but the major energy substances (proteins, lipids, and carbohydrates) metabolites and metabolic pathways are unknown. Comprehensive metabolomic analyses were performed using Fourier transform infrared spectroscopy and high-performance liquid chromatography mass spectrometry on flaxseed from 0 to 7 d. Additionally, the critical metabolites pathways networks of three energy substances metabolites during flaxseed germination were exhibited. The results showed that arginine was the most active metabolite during germination, strongly associated with the arginine biosynthesis and arginine and proline metabolism pathways. Carbohydrates predominantly comprised sucrose on 0-3 d, which participated in galactose metabolism and starch and sucrose metabolism. The main flaxseed phospholipid molecules were phosphatidic acid, phosphatidylethanolamine, lysophosphatidic acid, and lysophosphatidylcholine during germination. This study underscores the paramount metabolic pathways in proteins, lipids and carbohydrates were arginine and proline metabolism, linoleic acid metabolism, arachidonic acid metabolism, and ascorbate and aldarate metabolism during germination.

Integrating multiple regulations on enzyme activity: the case of phosphoenolpyruvate carboxykinases

Data on protein post-translational modifications (PTMs) increased exponentially in the last years due to the refinement of mass spectrometry techniques and the development of databases to store and share datasets. Nevertheless, these data per se do not create comprehensive biochemical knowledge. Complementary studies on protein biochemistry are necessary to fully understand the function of these PTMs at the molecular level and beyond, for example, designing rational metabolic engineering strategies to improve crops. Phosphoenolpyruvate carboxykinases (PEPCKs) are critical enzymes for plant metabolism with diverse roles in plant development and growth. Multiple lines of evidence showed the complex regulation of PEPCKs, including PTMs. Herein, we present PEPCKs as an example of the integration of combined mechanisms modulating enzyme activity and metabolic pathways. PEPCK studies strongly advanced after the production of the recombinant enzyme and the establishment of standardized biochemical assays. Finally, we discuss emerging open questions for future research and the challenges in integrating all available data into functional biochemical models.

Bioinformatics Analysis of WRKY Family Genes in Flax (Linum usitatissimum)

WRKY gene family is one of the largest transcription factor families involved in various physiological processes of plants. Flax (Linum usitatissimum) is an important stem fiber crop, and it is also an economically important crop in natural fiber and textile industries around the world. In this study, 105 WRKY genes were obtained by screening the whole genome of flax. There were 26 in group I, 68 in group II, 8 in group III and 3 in group UN. The characteristics of the WRKY motif and gene structure in each group are similar. The promoter sequence of WRKY genes includes photoresponsive elements, core regulatory elements and 12 cis-acting elements under abiotic stress. Similar to A. thaliana and Compositae plants, WRKY genes are evenly distributed on each chromosome, with segmental and tandem repeated events, which play a major role in the evolution of WRKY genes. The flax WRKY gene family is mainly concentrated in group I and group II. This study is mainly based on genome-wide information to classify and analyze the flax WRKY gene family, laying a foundation for further understanding the role of WRKY transcription factors in species evolution and functional analysis.

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.

The Arabidopsis transcription factors AtPHL1 and AtHB23 act together promoting carbohydrate transport from pedicel-silique nodes to seeds

Carbohydrates are produced in green tissues through photosynthesis and then transported to sink tissues. Carbon partitioning is a strategic process, fine regulated, involving specific sucrose transporters in each connecting tissue. Here we report that a screening of an Arabidopsis transcription factor (TF) library using the homeodomain-leucine zipper I member AtHB23 as bait, allowed identifying the TF AtPHL1 interacting with the former. An independent Y2H assay, and in planta by BiFC, confirmed such interaction. AtHB23 and AtPHL1 coexpressed in the pedicel-silique nodes and the funiculus. Mutant plants (phl1, and amiR23) showed a marked reduction of lipid content in seeds, although lipid composition did not change compared to the wild type. While protein and carbohydrate contents were not significantly different between mutants and control mature seeds, we observed a reduced carbohydrate content in mutant plants young siliques (7 days after pollination). Moreover, using a CFDA probe, we revealed an impaired transport to the seeds, and the gene encoding the carbohydrate transporters SWEET10 and SWEET11, usually expressed in connecting tissues, was repressed in the amiR23 and phl1 mutant plants. Altogether, the results indicated that AtHB23 and AtPHL1 act together, promoting sucrose transport, and the lack of any of them provoked a reduction in seeds lipid content.

A salicylic acid inducible mulberry WRKY transcription factor, MiWRKY53 is involved in plant defence response

MiWRKY53 is expressed in response to various stresses and hormones. Although it is localized in the nucleus, it shows no transcriptional activation. Role of SA-mediated plant defence response is demonstrated. WRKY transcription factors are one the largest gene families in plants involved in almost every process in plants including development, physiological processes, and stress response. Salicylic acid (SA) is key regulator of biotic stress against various pathogens in plants acting via its multiple mechanisms to induce defence response. Herein, we have identified and functionally validated WRKY53 from mulberry (Morus indica var. K2). MiWRKY53 expressed differentially in response to different stress and hormonal treatments. MiWRKY53 belongs to group III of WKRY gene family, localized in nucleus, and lacks transcriptional activation activity in yeast. Hormone responsive behaviour of MiWRKY53 Arabidopsis overexpression (OE) transgenics preferentially was noted in root growth assay in response to Salicylic acid (SA). Arabidopsis overexpression plants also displayed alteration in leaf phenotype having wider leaves than the wild-type plants. PR-1 transcripts were higher in MiWRKY53 Arabidopsis OE plants and they displayed resistance towards biotrophic pathogen Pseudomonas syringae PstDC3000. MiWRKY53 Mulberry OE transgenics also depicted SA-responsive behaviour. Several hormones and stress-related cis-acting elements were also identified in the 1.2-Kb upstream regulatory region (URR) of MiWRKY53. Functional characterization of full-length promoter region revealed that it is induced by SA and further analysis of deletion constructs helped in the identification of minimal promoter responsible for its inducibility by SA. Altogether, the findings from this study point towards the SA preferential behaviour of MiWRKY53 and its function as regulator of plant defence response through SA-mediated mechanisms.

Tyr-Asp inhibition of glyceraldehyde 3-phosphate dehydrogenase affects plant redox metabolism

How organisms integrate metabolism with the external environment is a central question in biology. Here, we describe a novel regulatory small molecule, a proteogenic dipeptide Tyr-Asp, which improves plant tolerance to oxidative stress by directly interfering with glucose metabolism. Specifically, Tyr-Asp inhibits the activity of a key glycolytic enzyme, glyceraldehyde 3-phosphate dehydrogenase (GAPC), and redirects glucose toward pentose phosphate pathway (PPP) and NADPH production. In line with the metabolic data, Tyr-Asp supplementation improved the growth performance of both Arabidopsis and tobacco seedlings subjected to oxidative stress conditions. Moreover, inhibition of Arabidopsis phosphoenolpyruvate carboxykinase (PEPCK) activity by a group of branched-chain amino acid-containing dipeptides, but not by Tyr-Asp, points to a multisite regulation of glycolytic/gluconeogenic pathway by dipeptides. In summary, our results open the intriguing possibility that proteogenic dipeptides act as evolutionarily conserved small-molecule regulators at the nexus of stress, protein degradation, and metabolism.

Comprehensive RNA-seq analysis revealed molecular pathways and genes associated with drought tolerance in wild soybean (Glycine soja Sieb. and Zucc.)

Drought stress at the germination stage is an important environmental stress limiting crop yield. Hence, our study investigated comparative root transcriptome profiles of four contrasting soybean genotypes viz., drought-tolerant (PI342618B/DTP and A214/DTL) and drought-sensitive (NN86-4/DSP and A195/DSL) under drought stress using RNA-Seq approach. A total of 4850 and 6272 differentially expressed genes (DEGs) were identified in tolerant (DTP and DTL) and sensitive (DSP and DSL) genotypes, respectively. Principle component analysis (PCA) and correlation analysis revealed higher correlation between DTP and DTL. Both gene ontology (GO) and MapMan analyses showed that the drought response was enriched in DEGs associated with water and auxin transport, cell wall/membrane, antioxidant activity, catalytic activity, secondary metabolism, signaling and transcription factor (TF) activities. Out of 981 DEGs screened from above terms, only 547 showed consistent opposite expression between contrasting genotypes. Twenty-eight DEGs of 547 were located on Chr.08 rich in QTLs and "Hotspot regions" associated with drought stress, and eight of them showed non-synonymous single nucleotide polymorphism. Hence, 10 genes (including above eight genes plus two hub genes) were predicated as possible candidates regulating drought tolerance, which needs further functional validation. Overall, the transcriptome profiling provided in-depth understanding about the genetic mechanism and candidate genes underlying drought tolerance in soybean.

Genome-Wide Characterization of WRKY Transcription Factors Revealed Gene Duplication and Diversification in Populations of Wild to Domesticated Barley

The WRKY transcription factors (WRKYs) are known for their crucial roles in biotic and abiotic stress responses, and developmental and physiological processes. In barley, early studies revealed their importance, whereas their diversity at the population scale remains hardly estimated. In this study, 98 HsWRKYs and 103 HvWRKYs have been identified from the reference genome of wild and cultivated barley, respectively. The tandem duplication and segmental duplication events from the cultivated barley were observed. By taking advantage of early released exome-captured sequencing datasets in 90 wild barley accessions and 137 landraces, the diversity analysis uncovered synonymous and non-synonymous variants instead of loss-of-function mutations that had occurred at all WRKYs. For majority of WRKYs, the haplotype and nucleotide diversity both decreased in cultivated barley relative to the wild population. Five WRKYs were detected to have undergone selection, among which haplotypes of WRKY9 were enriched, correlating with the geographic collection sites. Collectively, profiting from the state-of-the-art barley genomic resources, this work represented the characterization and diversity of barley WRKY transcription factors, shedding light on future deciphering of their roles in barley domestication and adaptation.

Aethionema arabicum genome annotation using PacBio full-length transcripts provides a valuable resource for seed dormancy and Brassicaceae evolution research

Aethionema arabicum is an important model plant for Brassicaceae trait evolution, particularly of seed (development, regulation, germination, dormancy) and fruit (development, dehiscence mechanisms) characters. Its genome assembly was recently improved but the gene annotation was not updated. Here, we improved the Ae. arabicum gene annotation using 294 RNA-seq libraries and 136 307 full-length PacBio Iso-seq transcripts, increasing BUSCO completeness by 11.6% and featuring 5606 additional genes. Analysis of orthologs showed a lower number of genes in Ae. arabicum than in other Brassicaceae, which could be partially explained by loss of homeologs derived from the At-α polyploidization event and by a lower occurrence of tandem duplications after divergence of Aethionema from the other Brassicaceae. Benchmarking of MADS-box genes identified orthologs of FUL and AGL79 not found in previous versions. Analysis of full-length transcripts related to ABA-mediated seed dormancy discovered a conserved isoform of PIF6-β and antisense transcripts in ABI3, ABI4 and DOG1, among other cases found of different alternative splicing between Turkey and Cyprus ecotypes. The presented data allow alternative splicing mining and proposition of numerous hypotheses to research evolution and functional genomics. Annotation data and sequences are available at the Ae. arabicum DB (https://plantcode.online.uni-marburg.de/aetar_db).

Proteolytic cleavage of Arabidopsis thaliana phosphoenolpyruvate carboxykinase-1 modifies its allosteric regulation

Phosphoenolpyruvate carboxykinase (PEPCK) plays a crucial role in gluconeogenesis. In this work, we analyze the proteolysis of Arabidopsis thaliana PEPCK1 (AthPEPCK1) in germinating seedlings. We found that the amount of AthPEPCK1 protein peaks at 24-48 h post-imbibition. Concomitantly, we observed shorter versions of AthPEPCK1, putatively generated by metacaspase-9 (AthMC9). To study the impact of AthMC9 cleavage on the kinetic and regulatory properties of AthPEPCK1, we produced truncated mutants based on the reported AthMC9 cleavage sites. The Δ19 and Δ101 truncated mutants of AthPEPCK1 showed similar kinetic parameters and the same quaternary structure as the wild type. However, activation by malate and inhibition by glucose 6-phosphate were abolished in the Δ101 mutant. We propose that proteolysis of AthPEPCK1 in germinating seedlings operates as a mechanism to adapt the sensitivity to allosteric regulation during the sink-to-source transition.

Transcription Factor TaWRKY51 Is a Positive Regulator in Root Architecture and Grain Yield Contributing Traits

Wheat is one of the staple food crops. The utilization of elite genetic resources to develop resource-efficient wheat varieties is an effective approach to deal with the challenges of climate change and population growth. WRKY transcription factors (TFs) are multifaceted regulators of plant growth and development and response to environmental stress. The previous studies have shown that TaWRKY51 positively regulates the development of lateral roots, while its roles in agronomic trait development are not clear, and there is no functional marker for molecular breeding. To bridge the gap, we cloned the three members of TaWRKY51 and found they were highly expressed in the roots and flag leaves at the flowering stage and were induced by the multiple abiotic stresses and phytohormones. The highest expression level was observed in TaWRKY51-2D, followed by TaWRKY51-2A and -2B. The two haplotypes/alleles for each member were identified in the natural populations, and functional markers were developed accordingly. The association assays revealed that Hap-2A-I was an elite haplotype for the large spike, Hap-2B-II and allele-G were favorable haplotypes/alleles for long root. However, only Hap-2A-I was selected for wheat breeding in China. The results of transgenic experiments showed that the rice lines overexpressing TaWRKY51 had large panicle, high thousand-grain-weight, and more crown and lateral roots, which further confirmed the results of association analysis. In short, TaWRKY51 is a positive regulator of the root architecture and grain yield (GY) contributing traits. The elite gene resources and functional markers may be utilized in the marker-assisted selection for high-yield breeding in wheat.

Genome wide association mapping of epi-cuticular wax genes in Sorghum bicolor

Sorghum accumulates epi-cuticular wax (EW) in leaves, sheaths, and culms. EW reduces the transpirational and nontranspirational (nonstomatal) water loss and protects the plant from severe drought stress in addition to imparting resistance against insect pests. Results presented here are from the analysis of EW content of 387 diverse sorghum accessions and its genome-wide association study (GWAS). EW content in sorghum leaves ranged from 0.1 to 29.7 mg cm^-2 with a mean value of 5.1 mg cm^-2. GWAS using 265,487 single nucleotide polymorphisms identified thirty-seven putative genes associated (P < 9.89E-06) with EW biosynthesis and transport in sorghum. Major EW biosynthetic genes identified included 3-Oxoacyl-[acyl-carrier-protein (ACP)] synthase III, an Ankyrin repeat protein, a bHLH-MYC, and an R2R3-MYB transcription factor. Genes involved in EW regulation or transport included an ABC transporter, a Lipid exporter ABCA1, a Multidrug resistance protein, Inositol 1, 3, 4-trisphosphate 5/6-kinase, and a Cytochrome P450. This GWA study thus demonstrates the potential for genetic manipulation of EW content in sorghum for better adaptation to biotic and abiotic stress.

Seed Germination in Oil Palm (Elaeis guineensis Jacq.): A Review of Metabolic Pathways and Control Mechanisms

Oil palm is an oil-producing crop of major importance at the global scale. Oil palm mesocarp lipids are used for myriads industrial applications, and market demand has been growing for decades. In addition, oil palm seeds are oleaginous, and the oil extracted therefrom can be used for several purposes, from food to cosmetics. As such, there is a huge need in oil palm seeds to maintain the global cohort of more than 2 billion trees. However, oil palm seed germination is a rather difficult process, not only to break dormancy, but also because it is long and often reaches lower-than-expected germination rates. Surprisingly, despite the crucial importance of germination for oil palm plantation management, our knowledge is still rather limited, in particular about germinating oil palm seed metabolism. The present review incorporates different pieces of information that have been obtained in the past few years, in oil palm and in other palm species, in order to provide an overview of germination metabolism and its control. Further insights can also be gained from other oleaginous model plants, such as Arabidopsis or canola, however, palm seeds have peculiarities that must be accounted for, to gain a better understanding of germinating seed metabolism.

Cellular alteration and differential protein profile explain effects of GA3 and ABA and their inhibitor on Trichocline catharinensis (Asteraceae) seed germination

Seed physiology of wild species has not been studied as deeply as that of domesticated crop species. Trichocline catharinensis (Asteraceae) is an endemic wildflower species from the high-altitude fields of southern Brazil. This species is of interest as a source of genes to improve cultivated Asteraceae because of its ornamental features, disease resistance and ability to tolerate drought and poor soil conditions. We studied the effects of abscisic acid (ABA) and gibberellic acid (GA₃ ) and their inhibitors, fluridone (FLU) and paclobutrazol (PAC), on seed germination. We individually assessed ultrastructural changes and differential protein accumulation. The principal component analysis explained 69.66% of differential accumulation for 32 proteins at phase II of seed germination in response to hormone and inhibitor treatment. GA₃ -imbibed seed germination (98.75%) resulted in increased protein accumulation to meet energy demand, redox regulation, and reserve metabolism activation. FLU-imbibed seeds showed a higher germination speed index as a consequence of metabolism activation. ABA-imbibed seeds (58.75%) showed osmotolerance and flattened cells in the hypocotyl-radicular axis, suggesting that ABA inhibits cell expansion. PAC-imbibed seeds remained at phase II for 300 h, and germination was suppressed (7.5%) because of the increased signaling proteins and halted reserve mobilization. Therefore, our findings provide insight into the behavior of Asteraceae non-dormant seed germination, which broadens our knowledge of seed germination in a wild and endemic plant species from a threatened ecosystem.

Comprehensive dynamic transcriptome analysis at two seed germination stages in maize (Zea mays L.)

Seed germination, as an integral stage of crop production, directly affects Zea mays (maize) yield and grain quality. However, the molecular mechanisms of seed germination remain unclear in maize. We performed comparative transcriptome analysis of two maize inbred lines, Yu82 and Yu537A, at two stages of seed germination. Expression profile analysis during seed germination revealed that a total of 3381 and 4560 differentially expressed genes (DEGs) were identified in Yu82 and Yu537A at the two stages. Transcription factors were detected from several families, such as the bZIP, ERF, WRKY, MYB and bHLH families, which indicated that these transcription factor families might be involved in driving seed germination in maize. Prominent DEGs were submitted for KEGG enrichment analysis, which included plant hormones, amino acid mechanism, nutrient reservoir, metabolic pathways and ribosome. Of these pathways, genes associated with plant hormones, especially gibberellins, abscisic acid and auxin may be important for early germination in Yu82. In addition, DEGs involved in amino acid mechanism showed significantly higher expression levels in Yu82 than in Yu537A, which indicated that energy supply from soluble sugars and amino acid metabolism may contribute to early germination in Yu82. This results provide novel insights into transcriptional changes and gene interactions in maize during seed germination.

Genome-Wide Identification of WRKY Transcription Factors in the Asteranae

The WRKY transcription factors family, which participates in many physiological processes in plants, constitutes one of the largest transcription factor families. The Asterales and the Apiales are two orders of flowering plants in the superorder Asteranae. Among the members of the Asterales, globe artichoke (Cynara cardunculus var. scolymus L.), sunflower (Helianthus annuus L.), and lettuce (Lactuca sativa L.) are important economic crops worldwide. Within the Apiales, ginseng (Panax ginseng C. A. Meyer) and Panax notoginseng (Burk.) F.H. Chen are important medicinal plants, while carrot (Daucus carota subsp. carota L.) has significant economic value. Research involving genome-wide identification of WRKY transcription factors in the Asterales and the Apiales has been limited. In this study, 490 WRKY genes, 244 from three species of the Apiales and 246 from three species of the Asterales, were identified and categorized into three groups. Within each group, WRKY motif characteristics and gene structures were similar. WRKY gene promoter sequences contained light responsive elements, core regulatory elements, and 12 abiotic stress cis-acting elements. WRKY genes were evenly distributed on each chromosome. Evidence of segmental and tandem duplication events was found in all six species in the Asterales and the Apiales, with segmental duplication inferred to play a major role in WRKY gene evolution. Among the six species, we uncovered 54 syntenic gene pairs between globe artichoke and lettuce. The six species are thus relatively closely related, consistent with their traditional taxonomic placement in the Asterales. This study, based on traditional species classifications, was the first to identify WRKY transcription factors in six species from the Asteranae. Our results lay a foundation for further understanding of the role of WRKY transcription factors in species evolution and functional differentiation.

Reactive Oxygen Species and Gibberellin Acid Mutual Induction to Regulate Tobacco Seed Germination

Seed germination is a complex process controlled by various mechanisms. To examine the potential contribution of reactive oxygen species (ROS) and gibberellin acid (GA) in regulating seed germination, diphenylene iodonium chloride (DPI) and uniconazole (Uni), as hydrogen peroxide (H2O2) and GA synthesis inhibitor, respectively, were exogenously applied on tobacco seeds using the seed priming method. Seed priming with DPI or Uni decreased germination percentage as compared with priming with H2O, especially the DPI + Uni combination. H2O2 and GA completely reversed the inhibition caused by DPI or Uni. The germination percentages with H2O2 + Uni and GA + DPI combinations kept the same level as with H2O. Meanwhile, GA or H2O2 increased GA content and deceased ABA content through corresponding gene expressions involving homeostasis and signal transduction. In addition, the activation of storage reserve mobilization and the enhancement of soluble sugar content and isocitrate lyase (ICL) activity were also induced by GA or H2O2. These results strongly suggested that H2O2 and GA were essential for tobacco seed germination and by downregulating the ABA/GA ratio and inducing reserve composition mobilization mutually promoted seed germination. Meanwhile, ICL activity was jointly enhanced by a lower ABA/GA ratio and a higher ROS concentration.

Genome-wide analysis of WRKY transcription factors in wheat (Triticum aestivum L.) and differential expression under water deficit condition

BACKGROUND

WRKY proteins, which comprise one of the largest transcription factor (TF) families in the plant kingdom, play crucial roles in plant development and stress responses. Despite several studies on WRKYs in wheat (Triticum aestivum L.), functional annotation information about wheat WRKYs is limited.

RESULTS

Here, 171 TaWRKY TFs were identified from the whole wheat genome and compared with proteins from 19 other species representing nine major plant lineages. A phylogenetic analysis, coupled with gene structure analysis and motif determination, divided these TaWRKYs into seven subgroups (Group I, IIa-e, and III). Chromosomal location showed that most TaWRKY genes were enriched on four chromosomes, especially on chromosome 3B. In addition, 85 (49.7%) genes were either tandem (5) or segmental duplication (80), which suggested that though tandem duplication has contributed to the expansion of TaWRKY family, segmental duplication probably played a more pivotal role. Analysis of cis-acting elements revealed putative functions of WRKYs in wheat during development as well as under numerous biotic and abiotic stresses. Finally, the expression of TaWRKY genes in flag leaves, glumes, and lemmas under water-deficit condition were analyzed. Results showed that different TaWRKY genes preferentially express in specific tissue during the grain-filling stage.

CONCLUSION

Our results provide a more extensive insight on WRKY gene family in wheat, and also contribute to the screening of more candidate genes for further investigation on function characterization of WRKYs under various stresses.