Elucidating the evolutionary conserved DNA-binding specificities of WRKY transcription factors by molecular dynamics and in vitro binding assays

Transcription factor OsWRKY11 induces rice heading at low concentrations but inhibits rice heading at high concentrations

The heading date of rice is a crucial agronomic characteristic that influences its adaptability to different regions and its productivity potential. Despite the involvement of WRKY transcription factors in various biological processes related to development, the precise mechanisms through which these transcription factors regulate the heading date in rice have not been well elucidated. The present study identified OsWRKY11 as a WRKY transcription factor which exhibits a pivotal function in the regulation of the heading date in rice through a comprehensive screening of a clustered regularly interspaced palindromic repeats (CRISPR) ‒ CRISPR-associated nuclease 9 mutant library that specifically targets the WRKY genes in rice. The heading date of oswrky11 mutant plants and OsWRKY11-overexpressing plants was delayed compared with that of the wild-type plants under short-day and long-day conditions. Mechanistic investigation revealed that OsWRKY11 exerts dual effects on transcriptional promotion and suppression through direct and indirect DNA binding, respectively. Under normal conditions, OsWRKY11 facilitates flowering by directly inducing the expression of OsMADS14 and OsMADS15. The presence of elevated levels of OsWRKY11 protein promote formation of a ternary protein complex involving OsWRKY11, Heading date 1 (Hd1), and Days to heading date 8 (DTH8), and this complex then suppresses the expression of Ehd1, which leads to a delay in the heading date. Subsequent investigation revealed that a mild drought condition resulted in a modest increase in OsWRKY11 expression, promoting heading. Conversely, under severe drought conditions, a significant upregulation of OsWRKY11 led to the suppression of Ehd1 expression, ultimately causing a delay in heading date. Our findings uncover a previously unacknowledged mechanism through which the transcription factor OsWRKY11 exerts a dual impact on the heading date by directly and indirectly binding to the promoters of target genes.

Unlocking nature's secrets: The pivotal role of WRKY transcription factors in plant flowering and fruit development

The WRKY transcription factor family is a key player in the regulatory mechanisms of flowering plants, significantly influencing both their biotic and abiotic response systems as well as being vital to numerous physiological and biological functions. Over the past two decades, the functionality of WRKY proteins has been the subject of extensive research in over 50 plant species, with a strong focus on their roles in responding to various stresses. Despite this extensive research, there remains a notable gap in comprehensive studies aimed at understanding how specific WRKY genes directly influence the timing of flowering and fruit development. This review offers an up-to-date look at WRKY family genes and provides insights into the key genes of WRKY to control flowering, enhance fruit ripening and secondary metabolism synthesis, and maintain fruit quality of various plants, including annuals, perennials, medicinal, and crop plants. The WRKY transcription factors serve as critical regulators within the transcriptional regulatory network, playing a crucial role in the precise enhancement of flowering processes. It is also involved in the up-regulation of fruit ripening was strongly demonstrated by combined transcriptomics and metabolomic investigation. Therefore, we speculated that the WRKY family is known to be a key regulator of flowering and fruiting in plants. This detailed insight will enable the identification of the series of molecular occurrences featuring WRKY proteins throughout the stages of flowering and fruiting.

Genome-Wide Identification of the WRKY Gene Family and Functional Characterization of CpWRKY5 in Cucurbita pepo

The WRKY gene family is crucial for regulating plant growth and development. However, the WRKY gene is rarely studied in naked kernel formation in hull-less Cucurbita pepo L. (HLCP), a natural mutant that lacks the seed coat. In this research, 76 WRKY genes were identified through bioinformatics-based methods in C. pepo, and their phylogenetics, conserved motifs, synteny, collinearity, and temporal expression during seed coat development were analyzed. The results showed that 76 CpWRKYs were identified and categorized into three main groups (I-III), with Group II further divided into five subgroups (IIa-IIe). Moreover, 31 segmental duplication events were identified in 49 CpWRKY genes. A synteny analysis revealed that C. pepo shared more collinear regions with cucumber than with melon. Furthermore, quantitative RT-PCR (qRT-PCR) results indicated the differential expression of CpWRKYs across different varieties, with notable variations in seed coat development between HLCP and CP being attributed to differences in CpWRKY5 expression. To investigate this further, CpWRKY5-overexpression tobacco plants were generated, resulting in increased lignin content and an upregulation of related genes, as confirmed by qRT-PCR. This study offers valuable insights for future functional investigations of CpWRKY genes and presents novel information for understanding the regulation mechanism of lignin synthesis.

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

Differential heat-response characteristics of two plastid isoforms of triose phosphate isomerase in tomato

Heat stress causes dysfunction of the carbon-assimilation metabolism. As a member of Calvin-Benson-Bassham (CBB) cycle, the chloroplast triose phosphate isomerases (TPI) catalyse the interconversion of glyceraldehyde 3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP). The tomato (Solanum lycopersicum) genome contains two individual SlTPI genes, Solyc10g054870 and Solyc01g111120, which encode the chloroplast-located proteins SlTPI1 and SlTPI2, respectively. The tpi1 and tpi2 single mutants had no visible phenotypes, but the leaves of their double mutant lines tpi1tpi2 had obviously reduced TPI activity and displayed chlorotic variegation, dysplasic chloroplasts and lower carbon-assimilation efficiency. In addition to altering carbon metabolism, proteomic data showed that the loss of both SlTPI1 and SlTPI2 severely affected photosystem proteins, reducing photosynthetic capacity. None of these phenotypes was evident in the tpi1 or tpi2 single mutants, suggesting that SlTPI1 and SlTPI2 are functionally redundant. However, the two proteins differed in their responses to heat stress; the protein encoded by the heat-induced SlTPI2 showed a higher level of thermotolerance than that encoded by the heat-suppressed SlTPI1. Notably, heat-induced transcription factors, SlWRKY21 and SlHSFA2/7, which negatively regulated SlTPI1 expression and positively regulated SlTPI2 expression, respectively. Our findings thus reveal that SlTPI1 and SlTPI2 have different thermostabilities and expression patterns in response to heat stress, which have the potential to be applied in thermotolerance strategies in crops.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

Review: Structural-functional relationship of WRKY transcription factors: Unfolding the role of WRKY in plants

WRKY as the name suggests, are the transcription factors (TFs) that contain the signature WRKY domains, hence named after it. Since their discovery in 1994, they have been well studied in plants with exploration of approximately 74 WRKY genes in the model plant, Arabidopsis alone. However, the study of these transcription factors (TFs) is not just limited to model plant now. They have been studied widely in crop plants as well, because of their tremendous contribution in stress as well as in growth and development. Here, in this review, we describe the story of WRKY TFs from their identification to their origin, the binding mechanisms, structure and their contribution in regulating plant development and stress physiology. High throughput transcriptomics-based data also opened a doorway to understand the comprehensive and detailed functioning of WRKY TFs in plants. Indeed, the detailed functional role of each and every WRKY member in regulating the gene expression is required to pave the path to develop holistic understanding of their role in stress physiology and developmental processes in plants.

Genome-wide identification and expression profiling of the WRKY gene family reveals abiotic stress response mechanisms in Platycodon grandiflorus

The WRKY family of transcription factors (TFs) is an important gene family involved in abiotic stress responses. Although the roles of WRKY TFs in plant abiotic stress responses are well studied, little is known about the stress-induced changes in WRKY family in Platycodon grandiflorus. 42 PgWRKY genes in seven subgroups were identified in the P. grandiflorus genome. The content of eight platycodins in P. grandiflorus was investigated under cold, heat, and drought stresses. Platycodin D levels significantly increased under three abiotic stresses, while the content changes of other platycodins varied. Transcriptome analysis showed that different WRKY family members exhibited varied expression patterns under different abiotic stresses. PgWRKY20, PgWRKY26, and PgWRKY39 were identified as three key candidates for temperature and drought stress responses, and were cloned and analysed for sequence characteristics, gene structure, subcellular localisation, and expression patterns. The RT-qPCR results showed that PgWRKY26 expression significantly increased after heat stress for 48 h, cold stress for 6 h, and drought stress for 2 d (DS_2 d). The PgWRKY39 expression level significantly increased at DS_2 d. This study provides a theoretical basis for clarifying the molecular mechanism of the abiotic stress responses of the WRKY gene family in P. grandiflorus.

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

Introduction

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

Methods

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

Results

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

Discussion

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

WRKY Transcription Factors in Jasminum sambac: An Insight into the Regulation of Aroma Synthesis

WRKY transcription factors are one of the largest families of transcription regulators that play essential roles in regulating the synthesis of secondary metabolites in plants. Jasmine (Jasminum sambac), renowned for its aromatic nature and fragrant blossoms, possesses a significant abundance of volatile terpene compounds. However, the role of the WRKY family in terpene synthesis in jasmine remains undetermined. In this study, 72 WRKY family genes of J. sambac were identified with their conserved WRKY domains and were categorized into three main groups based on their structural and phylogenetic characteristics. The extensive segmental duplications contributed to the expansion of the WRKY gene family. Expression profiles derived from the transcriptome data and qRT-PCR analysis showed that the majority of JsWRKY genes were significantly upregulated in fully bloomed flowers compared to buds. Furthermore, multiple correlation analyses revealed that the expression patterns of JsWRKYs (JsWRKY27/33/45/51/55/57) were correlated with both distinct terpene compounds (monoterpenes and sesquiterpenes). Notably, the majority of jasmine terpene synthase (JsTPS) genes related to terpene synthesis and containing W-box elements exhibited a significant correlation with JsWRKYs, particularly with JsWRKY51, displaying a strong positive correlation. A subcellular localization analysis showed that JsWRKY51 was localized in the nucleus. Moreover, transgenic tobacco leaves and jasmine calli experiments demonstrated that overexpression of JsWRKY51 was a key factor in enhancing the accumulation of β-ocimene, which is an important aromatic terpene component. Collectively, our findings suggest the roles of JsWRKY51 and other JsWRKYs in regulating the synthesis of aromatic compounds in J. sambac, providing a foundation for the potential utilization of JsWRKYs to facilitate the breeding of fragrant plant varieties with an improved aroma.

WRKY transcription factors and OBERON histone-binding proteins form complexes to balance plant growth and stress tolerance

WRKY transcription factors in plants are known to be able to mediate either transcriptional activation or repression, but the mechanism regulating their transcriptional activity is largely unclear. We found that group IId WRKY transcription factors interact with OBERON (OBE) proteins, forming redundant WRKY-OBE complexes in Arabidopsis thaliana. The coiled-coil domain of WRKY transcription factors binds to OBE proteins and is responsible for target gene selection and transcriptional repression. The PHD finger of OBE proteins binds to both histones and WRKY transcription factors. WRKY-OBE complexes repress the transcription of numerous stress-responsive genes and are required for maintaining normal plant growth. Several WRKY and OBE mutants show reduced plant size and increased drought tolerance, accompanied by increased expression of stress-responsive genes. Moreover, expression levels of most of these WRKY and OBE genes are reduced in response to drought stress, revealing a previously uncharacterized regulatory mechanism of the drought stress response. These results suggest that WRKY-OBE complexes repress transcription of stress-responsive genes, and thereby balance plant growth and stress tolerance.

KaScape: a sequencing-based method for global characterization of protein‒DNA binding affinity

It is difficult to exhaustively screen all possible DNA binding sequences for a given transcription factor (TF). Here, we developed the KaScape method, in which TFs bind to all possible DNA sequences in the same DNA pool where DNA sequences are prepared by randomized oligo synthesis and the random length can be adjusted to a length such as 4, 5, 6, or 7. After separating bound from unbound double-stranded DNAs (dsDNAs), their sequences are determined by next-generation sequencing. To demonstrate the relative binding affinities of all possible DNA sequences determined by KaScape, we developed three-dimensional KaScape viewing software based on a K-mer graph. We applied KaScape to 12 plant TF family AtWRKY proteins and found that all AtWRKY proteins bound to the core sequence GAC with similar profiles. KaScape can detect not only binding sequences consistent with the consensus W-box "TTGAC(C/T)" but also other sequences with weak affinity. KaScape provides a high-throughput, easy-to-operate, sensitive, and exhaustive method for quantitatively characterizing the relative binding strength of a TF with all possible binding sequences, allowing us to comprehensively characterize the specificity and affinity landscape of transcription factors, particularly for moderate- and low-affinity binding sites.

Senescence-related receptor kinase 1 functions downstream of WRKY53 in regulating leaf senescence in Arabidopsis

Receptor-like kinases (RLKs) are the most important class of cell surface receptors, and play crucial roles in plant development and stress responses. However, few studies have been reported about the biofunctions of RLKs in leaf senescence. Here, we characterized a novel Arabidopsis RLK-encoding gene, SENESCENCE-RELATED RECEPTOR KINASE 1 (SENRK1), which was significantly down-regulated during leaf senescence. Notably, the loss-of-function senrk1 mutants displayed an early leaf senescence phenotype, while overexpression of SENRK1 significantly delayed leaf senescence, indicating that SENRK1 negatively regulates age-dependent leaf senescence in Arabidopsis. Furthermore, the senescence-promoting transcription factor WRKY53 repressed the expression of SENRK1. While the wrky53 mutant showed a delayed senescence phenotype as previously reported, the wrky53 senrk1-1 double mutant exhibited precocious leaf senescence, suggesting that SENRK1 functions downstream of WRKY53 in regulating age-dependent leaf senescence in Arabidopsis.

Genome-Wide Identification of Bilberry WRKY Transcription Factors: Go Wild and Duplicate

WRKY transcription factor genes compose an important family of transcriptional regulators that are present in several plant species. According to previous studies, these genes can also perform important roles in bilberry (Vaccinium myrtillus L.) metabolism, making it essential to deepen our understanding of fruit ripening regulation and anthocyanin biosynthesis. In this context, the detailed characterization of these proteins will provide a comprehensive view of the functional features of VmWRKY genes in different plant organs and in response to different intensities of light. In this study, the investigation of the complete genome of the bilberry identified 76 VmWRKY genes that were evaluated and distributed in all twelve chromosomes. The proteins encoded by these genes were classified into four groups (I, II, III, and IV) based on their conserved domains and zinc finger domain types. Fifteen pairs of VmWRKY genes in segmental duplication and four pairs in tandem duplication were detected. A cis element analysis showed that all promoters of the VmWRKY genes contain at least one potential cis stress-response element. Differential expression analysis of RNA-seq data revealed that VmWRKY genes from bilberry show preferential or specific expression in samples. These findings provide an overview of the functional characterization of these proteins in bilberry.

Exploring the Binding Affinity of the ARR2 GARP DNA Binding Domain via Comparative Methods

Plants have evolved signaling mechanisms such as the multi-step phosphorelay (MSP) to respond to different internal and external stimuli. MSP responses often result in gene transcription regulation that is modulated through transcription factors such as B-type Arabidopsis response regulator (ARR) proteins. Among these proteins, ARR2 is a key component that is expressed ubiquitously and is involved in many aspects of plant development. Although it has been noted that B-type ARRs bind to their cognate genes through a DNA-binding domain termed the GARP domain, little is known about the structure and function of this type of DNA-binding domain; thus, how ARRs bind to DNA at a structural level is still poorly understood. In order to understand how the MSP functions in planta, it is crucial to unravel both the kinetics as well as the structural identity of the components involved in such interactions. For this reason, this work focusses on resolving how the GARP domain of ARR2 (GARP2) binds to the promoter region of ARR5, one of its native target genes in cytokinin signaling. We have established that GARP2 specifically binds to the ARR5 promoter with three different bi-molecular interaction systems-qDPI-ELISA, FCS, and MST-and we also determined the KD of this interaction. In addition, structural modeling of the GARP2 domain confirms that GARP2 entails a HTH motif, and that protein-DNA interaction most likely occurs via the α3-helix and the N-terminal arm of this domain since mutations in this region hinder ARR2's ability to activate transcription.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Rice WRKY13 TF protein binds to motifs in the promoter region to regulate downstream disease resistance-related genes

In plants, pathogen resistance is brought about by the binding of certain transcription factor (TF) proteins to the cis-elements of certain target genes. These cis-elements are present upstream in the motif of the promoters of each gene. This ensures the binding of a specific TF to a specific promoter, therefore regulating the expression of that gene. Therefore, the study of each promoter sequence of all the rice genes would help identify the target genes of a specific TF. Rice 1 kb upstream promoter sequences of 55,986 annotated genes were analyzed using the Perl program algorithm to detect WRKY13 binding motifs (bm). The resulting genes were grouped using Gene Ontology and gene set enrichment analysis. A gene with more than 4 TF bm in their promoter was selected. Ten genes reported to have a role in rice disease resistance were selected for further analysis. Cis-acting regulatory element analysis was carried out to find the cis-elements and confirm the presence of the corresponding motifs in the promoter sequences of these genes. The 3D structure of WRKY13 TF and the corresponding ten genes were built, and the interacting residues were determined. The binding capacity of WRKY13 to the promoter of these selected genes was analyzed using docking studies. WRKY13 was considered for docking analysis based on the prior reports of autoregulation. Molecular dynamic simulations provided more details regarding the interactions. Expression data revealed the expression of the genes that helped provide the mechanism of interaction. Further co-expression network helped to characterize the interaction of these selected disease resistance-related genes with the WRKY13 TF protein. This study suggests downstream target genes that are regulated by the WRKY13 TF. The molecular mechanism involving the gene network regulated by WRKY13 TF in disease resistance against rice fungal pathogens is explored.

Transcriptome-Wide Identification of the GRAS Transcription Factor Family in Pinus massoniana and Its Role in Regulating Development and Stress Response

Pinus massoniana is a species used in afforestation and has high economic, ecological, and therapeutic significance. P. massoniana experiences a variety of biotic and abiotic stresses, and thus presents a suitable model for studying how woody plants respond to such stress. Numerous families of transcription factors are involved in the research of stress resistance, with the GRAS family playing a significant role in plant development and stress response. Though GRASs have been well explored in various plant species, much research remains to be undertaken on the GRAS family in P. massoniana. In this study, 21 PmGRASs were identified in the P. massoniana transcriptome. P. massoniana and Arabidopsis thaliana phylogenetic analyses revealed that the PmGRAS family can be separated into nine subfamilies. The results of qRT-PCR and transcriptome analyses under various stress and hormone treatments reveal that PmGRASs, particularly PmGRAS9, PmGRAS10 and PmGRAS17, may be crucial for stress resistance. The majority of PmGRASs were significantly expressed in needles and may function at multiple locales and developmental stages, according to tissue-specific expression analyses. Furthermore, the DELLA subfamily members PmGRAS9 and PmGRAS17 were nuclear localization proteins, while PmGRAS9 demonstrated transcriptional activation activity in yeast. The results of this study will help explore the relevant factors regulating the development of P. massoniana, improve stress resistance and lay the foundation for further identification of the biological functions of PmGRASs.

Genome-Wide Identification and Characterization of WRKY Transcription Factors and Their Expression Profile in Loropetalum chinense var. rubrum

The WRKY gene family plays important roles in plant growth and development, as well as in the responses to biotic and abiotic stresses. Loropetalum chinense var. rubrum has high ornamental and medicinal value. However, few WRKY genes have been reported in this plant, and their functions remain unknown. To explore the roles that the WRKY genes play in L. chinense var. rubrum, we identified and characterized 79 LcWRKYs through BLAST homology analysis and renamed them (as LcWRKY1-79) based on their distribution on the chromosomes of L. chinense var. rubrum. In this way, according to their structural characteristics and phylogenetic analysis, they were divided into three groups containing 16 (Group I), 52 (Group II), and 11 (Group III) WRKYs, respectively. LcWRKYs in the same group have similar motifs and gene structures; for instance, Motifs 1, 2, 3, 4, and 10 constitute the WRKY domain and zinc-finger structure. The LcWRKY promoter region contains light response elements (ACE, G-box), stress response elements (TC-rich repeats), hormone response elements (TATC-box, TCA-element), and MYB binding sites (MBS, MBSI). Synteny analysis of LcWRKYs allowed us to establish orthologous relationships among the WRKY gene families of Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum L., Vitis vinifera L., Oryza sativa L., and Zea mays L.; furthermore, analysis of the transcriptomes of mature leaves and flowers from different cultivars demonstrated the cultivar-specific LcWRKY gene expression. The expression levels of certain LcWRKY genes also presented responsive changes from young to mature leaves, based on an analysis of the transcriptome in leaves at different developmental stages. White light treatment led to a significant decrease in the expression of LcWRKY6, 18, 24, 34, 36, 44, 48, 61, 62, and 77 and a significant increase in the expression of LcWRKY41, blue light treatment led to a significant decrease in the expression of LcWRKY18, 34, 50, and 77 and a significant increase in the expression of LcWRKY36 and 48. These results enable a better understanding of LcWRKYs, facilitating the further exploration of their genetic functions and the molecular breeding of L. chinense var. rubrum.

Solution structure of the DNA binding domain of Arabidopsis transcription factor WRKY11

The Arabidopsis WRKY11 (AtWRKY11) protein is an important transcription factor involved in plant response to biotic and abiotic stresses. Its DNA-binding domain specifically binds to gene promoter regions harboring the W-box consensus motif. Herein we report the high-resolution structure of the AtWRKY11 DNA-binding domain (DBD) determined by solution NMR spectroscopy. The results show that AtWRKY11-DBD adopts an all-β fold comprising five β-strands packed in an antiparallel topology, stabilized by a zinc-finger motif. Structural comparison reveals that the long β₁-β₂ loop shows the highest structural variation from other available WRKY domain structures. Moreover, this loop was further found to contribute to the binding between AtWRKY11-DBD and W-box DNA. Our current study provides atomic-level structural basis for further understanding the structure-function relationship of plant WRKY proteins.

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

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

LcWRKY17, a WRKY Transcription Factor from Litsea cubeba, Effectively Promotes Monoterpene Synthesis

The WRKY gene family is one of the most significant transcription factor (TF) families in higher plants and participates in many secondary metabolic processes in plants. Litsea cubeba (Lour.) Person is an important woody oil plant that is high in terpenoids. However, no studies have been conducted to investigate the WRKY TFs that regulate the synthesis of terpene in L. cubeba. This paper provides a comprehensive genomic analysis of the LcWRKYs. In the L. cubeba genome, 64 LcWRKY genes were discovered. According to a comparative phylogenetic study with Arabidopsis thaliana, these L. cubeba WRKYs were divided into three groups. Some LcWRKY genes may have arisen from gene duplication, but the majority of LcWRKY evolution has been driven by segmental duplication events. Based on transcriptome data, a consistent expression pattern of LcWRKY17 and terpene synthase LcTPS42 was found at different stages of L. cubeba fruit development. Furthermore, the function of LcWRKY17 was verified by subcellular localization and transient overexpression, and overexpression of LcWRKY17 promotes monoterpene synthesis. Meanwhile, dual-Luciferase and yeast one-hybrid (Y1H) experiments showed that the LcWRKY17 transcription factor binds to W-box motifs of LcTPS42 and enhances its transcription. In conclusion, this research provided a fundamental framework for future functional analysis of the WRKY gene families, as well as breeding improvement and the regulation of secondary metabolism in L. cubeba.

Regulation of the regulators: Transcription factors controlling biosynthesis of plant secondary metabolites during biotic stresses and their regulation by miRNAs

Biotic stresses threaten to destabilize global food security and cause major losses to crop yield worldwide. In response to pest and pathogen attacks, plants trigger many adaptive cellular, morphological, physiological, and metabolic changes. One of the crucial stress-induced adaptive responses is the synthesis and accumulation of plant secondary metabolites (PSMs). PSMs mitigate the adverse effects of stress by maintaining the normal physiological and metabolic functioning of the plants, thereby providing stress tolerance. This differential production of PSMs is tightly orchestrated by master regulatory elements, Transcription factors (TFs) express differentially or undergo transcriptional and translational modifications during stress conditions and influence the production of PSMs. Amongst others, microRNAs, a class of small, non-coding RNA molecules that regulate gene expression post-transcriptionally, also play a vital role in controlling the expression of many such TFs. The present review summarizes the role of stress-inducible TFs in synthesizing and accumulating secondary metabolites and also highlights how miRNAs fine-tune the differential expression of various stress-responsive transcription factors during biotic stress.

WRKY transcription factors: evolution, regulation, and functional diversity in plants

The recent advancements in sequencing technologies and informatic tools promoted a paradigm shift to decipher the hidden biological mysteries and transformed the biological issues into digital data to express both qualitative and quantitative forms. The transcriptomic approach, in particular, has added new dimensions to the versatile essence of plant genomics through the large and deep transcripts generated in the process. This has enabled the mining of super families from the sequenced plants, both model and non-model, understanding their ancestry, diversity, and evolution. The elucidation of the crystal structure of the WRKY proteins and recent advancement in computational prediction through homology modeling and molecular dynamic simulation has provided an insight into the DNA-protein complex formation, stability, and interaction, thereby giving a new dimension in understanding the WRKY regulation. The present review summarizes the functional aspects of the high volume of sequence data of WRKY transcription factors studied from different species, till date. The review focuses on the dynamics of structural classification and lineage in light of the recent information. Additionally, a comparative analysis approach was incorporated to understand the functions of the identified WRKY transcription factors subjected to abiotic (heat, cold, salinity, senescence, dark, wounding, UV, and carbon starvation) stresses as revealed through various sets of studies on different plant species. The review will be instrumental in understanding the events of evolution and the importance of WRKY TFs under the threat of climate change, considering the new scientific evidences to propose a fresh perspective.

Dioscorea composita WRKY12 is involved in the regulation of salt tolerance by directly activating the promoter of AtRCI2A

Dioscorea composita (D. composita) is an important medicinal plant worldwide with high economic value. However, its large-scale cultivation was limited by soil salinization. Identification of genes and their mechanisms of action in response to salt stress are critically important. In the present study, we isolated a classical WRKY transcription factor from D. composita, namely DcWRKY12, and analyzed its function in salt tolerance. Expression pattern analysis showed DcWRKY12 is mainly expressed in roots and significantly induced by NaCl, polyethylene glycol-6000 (PEG-6000), and abscisic acid (ABA). Phenotypic and physiological analyses revealed that heterologous expression of DcWRKY12 enhanced salt and osmotic stress tolerance by increasing antioxidant enzyme activity, osmoregulatory substance content, maintaining relative water content and ion homeostasis, decreasing reactive oxygen species and malondialdehyde content. Correspondingly, the overexpression of DcWRKY12 modulated the expression of salt stress-responsive and ion transport-related genes. Dual luciferase assay and Y1H were further confirmed that DcWRKY12 activates the promoter of AtRCI2A through directly binding to the specific W-box cis-acting elements. These results suggest that DcWRKY12 is a positive regulator of salt tolerance in D. composita and has potential applications in salt stress.

Genome-wide analysis of the Tritipyrum WRKY gene family and the response of TtWRKY256 in salt-tolerance

INTRODUCTION

The transcription factor WRKY is widespread in the plant kingdom and plays a crucial role in diverse abiotic stress responses in plant species. Tritipyrum, an octoploid derived from an intergeneric cross between Triticum aestivum (AABBDD) and Thinopyrum elongatum (EE), is a valuable germplasm resource for introducing superior traits of Th. elongatum into T. aestivum. The recent release of the complete genome sequences of T. aestivum and Th. elongatum enabled us to investigate the organization and expression profiling of Tritipyrum WRKY genes across the entire genome.

RESULTS

In this study, 346 WRKY genes, from TtWRKY1 to TtWRKY346, were identified in Tritipyrum. The phylogenetic analysis grouped these genes into three subfamilies (I-III), and members of the same subfamilies shared a conserved motif composition. The 346 TtWRKY genes were dispersed unevenly across 28 chromosomes, with 218 duplicates. Analysis of synteny suggests that the WRKY gene family may have a common ancestor. Expression profiles derived from transcriptome data and qPCR demonstrated that 54 TtWRKY genes exhibited relatively high levels of expression across various salt stresses and recovery treatments. Tel1E01T143800 (TtWRKY256) is extremely sensitive to salt stress and is on the same evolutionary branch as the salt-tolerant A. thaliana genes AtWRKY25 and AtWRKY33. From 'Y1805', the novel AtWRKY25 was cloned. The Pearson correlation analysis identified 181 genes that were positively correlated (R>0.9) with the expression of TtWRKY256, and these genes were mainly enriched in metabolic processes, cellular processes, response to stimulus, biological regulation, and regulation of biological. Subcellular localization and qRT-PCR analysis revealed that TtWRKY256 was located in the nucleus and was highly expressed in roots, stems, and leaves under salt stress.

DISCUSSION

The above results suggest that TtWRKY256 may be associated with salt stress tolerance in plants and may be a valuable alien gene for improving salt tolerance in wheat.

The colonization of land was a likely driving force for the evolution of mitochondrial retrograde signalling in plants

Most retrograde signalling research in plants was performed using Arabidopsis, so an evolutionary perspective on mitochondrial retrograde regulation (MRR) is largely missing. Here, we used phylogenetics to track the evolutionary origins of factors involved in plant MRR. In all cases, the gene families can be traced to ancestral green algae or earlier. However, the specific subfamilies containing factors involved in plant MRR in many cases arose during the transition to land. NAC transcription factors with C-terminal transmembrane domains, as observed in the key regulator ANAC017, can first be observed in non-vascular mosses, and close homologs to ANAC017 can be found in seed plants. Cyclin-dependent kinases (CDKs) are common to eukaryotes, but E-type CDKs that control MRR also diverged in conjunction with plant colonization of land. AtWRKY15 can be traced to the earliest land plants, while AtWRKY40 only arose in angiosperms and AtWRKY63 even more recently in Brassicaceae. Apetala 2 (AP2) transcription factors are traceable to algae, but the ABI4 type again only appeared in seed plants. This strongly suggests that the transition to land was a major driver for developing plant MRR pathways, while additional fine-tuning events have appeared in seed plants or later. Finally, we discuss how MRR may have contributed to meeting the specific challenges that early land plants faced during terrestrialization.

Comprehensive Genomic Survey, Structural Classification, and Expression Analysis of WRKY Transcription Factor Family in Rhododendron simsii

(1) Rhododendron is one of the top ten traditional flowers in China, with both high ornamental and economic values. However, with the change of the environment, Rhododendron suffers from various biological stresses. The WRKY transcription factor is a member of the most crucial transcription factor families, which plays an essential regulatory role in a variety of physiological processes and developmental stresses. (2) In this study, 57 RsWRKYs were identified using genome data and found to be randomly distributed on 13 chromosomes. Based on gene structure and phylogenetic relationships, 57 proteins were divided into three groups: I, II, and III. Multiple alignments of RsWRKYs with Arabidopsis thaliana homologous genes revealed that WRKY domains in different groups had different conserved sites. RsWRKYs have a highly conserved domain, WRKYGQK, with three variants, WRKYGKK, WRKYGEK, and WRKYGRK. Furthermore, cis-acting elements analysis revealed that all of the RsWRKYs had stress and plant hormone cis-elements, with figures varying by group. Finally, the expression patterns of nine WRKY genes treated with gibberellin acid (GA), methyl jasmonate (MeJA), heat, and drought in Rhododendron were also measured using quantitative real-time PCR (qRT-PCR). The results showed that the expression levels of the majority of RsWRKY genes changed in response to multiple phytohormones and abiotic stressors. (3) This current study establishes a theoretical basis for future studies on the response of RsWRKY transcription factors to various hormone and abiotic stresses as well as a significant foundation for the breeding of new stress-tolerant Rhododendron varieties.

Genome-Wide Identification and Characterization of the Oat (Avena sativa L.) WRKY Transcription Factor Family

The WRKY family is widely involved in the regulation of plant growth and stress response and is one of the largest gene families related to plant environmental adaptation. However, no systematic studies on the WRKY family in oat (Avena sativa L.) have been conducted to date. The recently published complete genome sequence of oat enables the systematic analysis of the AsWRKYs. Based on a genome-wide study of oat, we identified 162 AsWRKYs that were unevenly distributed across 21 chromosomes; a phylogenetic tree of WRKY domains divided these genes into three groups (I, II, and III). We also analyzed the gene duplication events and identified a total of 111 gene pairs that showed strong purifying selection during the evolutionary process. Surprisingly, almost all genes evolved after the completion of subgenomic differentiation of hexaploid oat. Further studies on the functional analysis indicated that AsWRKYs were widely involved in various biological processes. Notably, expression patterns of 16 AsWRKY genes revealed that the response of AsWRKYs were affected by stress level and time. In conclusion, this study provides a reference for further analysis of the role of WRKY transcription factors in species evolution and functional differentiation.

Genomic distribution and context dependent functionality of novel WRKY transcription factor binding sites

Background

The WT-boxes NGACTTTN are novel microbe-associated molecular pattern (MAMP)-responsive cis-regulatory sequences. Many of them are uncommon WRKY transcription factor (TF) binding sites.

Results

To understand their functional relevance, a genomic distribution analysis of the 16 possible WT-boxes and a functional analysis of a WT-box rich promoter was done. The genomic distribution analysis shows an enrichment of specific WT-boxes within 500 bp upstream of all Arabidopsis thaliana genes. Those that harbour a T 5′ to the core sequence GACTTT can also be part of the classic WRKY binding site the W-box TTGACT/C. The MAMP-responsive gene ATEP3, a class IV chitinase, harbours seven WT-boxes within its 1000 bp upstream region. In the context of synthetic promoters, the four proximal WT-boxes confer MAMP responsivity while the three WT-boxes further upstream have no effect. Rendering the nucleotides adjacent and in the vicinity of the WT-box core sequence reveals their functional importance for gene expression. A 158 bp long ATEP3 minimal promoter harbouring the two WT-boxes CGACTTTT, confers WT-box-dependent basal and MAMP-responsive reporter gene expression. The ATEP3 gene is a proposed target of WRKY50 and WRKY70. WRKY50 negatively regulates MAMP responsivity of the two WT-boxes CGACTTTT, while WRKY70 activates gene expression in a WT-box dependent manner. Both WRKY factors bind directly to the WT-box CGACTTTT.

Conclusion

In summary, WT-boxes are enriched in promoter regions and comprise novel and uncommon WRKY binding sites required for basal and MAMP-induced gene expression. WT-boxes not being part of a W-box may be a missing link for WRKY target gene prediction when these genes do not harbour a W-box.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08877-y.

The TaWRKY22–TaCOPT3D Pathway Governs Cadmium Uptake in Wheat

Cadmium (Cd) is a heavy metal nonessential for plants; this toxic metal accumulation in crops has significant adverse effects on human health. The crosstalk between copper (Cu) and Cd has been reported; however, the molecular mechanisms remain unknown. The present study investigated the function of wheat Cu transporter 3D (TaCOPT3D) in Cd tolerance. The TaCOPT3D transcripts significantly accumulated in wheat roots under Cd stress. Furthermore, TaCOPT3D-overexpressing lines were compared to wildtype (WT) plants to test the role of TaCOPT3D in Cd stress response. Under 20 mM Cd treatment, TaCOPT3D-overexpressing lines exhibited more biomass and lower root, shoot, and grain Cd accumulation than the WT plants. In addition, overexpression of TaCOPT3D decreased the reactive oxygen species (ROS) levels and increased the active antioxidant enzymes under Cd conditions. Moreover, the transcription factor (TF) TaWRKY22, which targeted the TaCOPT3D promoter, was identified in the regulatory pathway of TaCOPT3D under Cd stress. Taken together, these results show that TaCOPT3D plays an important role in regulating plant adaptation to cadmium stress through bound by TaWRKY22. These findings suggest that TaCOPT3D is a potential candidate for decreasing Cd accumulation in wheat through genetic engineering.

Identification of the WRKY Gene Family and Characterization of Stress-Responsive Genes in Taraxacum kok-saghyz Rodin

WRKY transcription factors present unusual research value because of their critical roles in plant physiological processes and stress responses. Taraxacum kok-saghyz Rodin (TKS) is a perennial herb of dandelion in the Asteraceae family. However, the research on TKS WRKY TFs is limited. In this study, 72 TKS WRKY TFs were identified and named. Further comparison of the core motifs and the structure of the WRKY motif was analyzed. These TFs were divided into three groups through phylogenetic analysis. Genes in the same group of TkWRKY usually exhibit a similar exon-intron structure and motif composition. In addition, virtually all the TKS WRKY genes contained several cis-elements related to stress response. Expression profiling of the TkWRKY genes was assessed using transcriptome data sets and Real-Time RT-PCR data in tissues during physiological development, under abiotic stress and hormonal treatments. For instance, the TkWRKY18, TkWRKY23, and TkWRKY38 genes were significantly upregulated during cold stress, whereas the TkWRKY21 gene was upregulated under heat-stress conditions. These results could provide a basis for further studies on the function of the TKS WRKY gene family and genetic amelioration of TKS germplasm.

The roles of WRKY transcription factors in Malus spp. and Pyrus spp

The WRKY transcription factor gene family is known to be involved in plant defense against pathogens and in tolerance to different environmental stresses at different stages of development. The response mechanisms through which these genes act can be influenced by different phytohormones as well as by many trans- and cis-acting elements, making this network an important topic for analysis, but still something complex to fully understand. According to available reports, these genes can also perform important roles in pome species (Malus spp. and Pyrus spp.) metabolism, especially in adaptation of these plants to stressful conditions. Here, we present a quick review of what is known about WRKY genes in Malus and Pyrus genomes offering a simple way to understand what is already known about this topic. We also add information connecting the evolution of these transcription factors with others that can also be found in pomes.

Panax notoginseng WRKY Transcription Factor 9 Is a Positive Regulator in Responding to Root Rot Pathogen Fusarium solani

Panax notoginseng (Burk) F.H. Chen is a rare and valuable Chinese herb, but root rot mainly caused by Fusarium solani severely affects the yield and quality of P. notoginseng herbal materials. In this study, we isolated 30 P. notoginseng WRKY transcription factors (TFs), which were divided into three groups (I, II, and III) on the basis of a phylogenetic analysis. The expression levels of 10 WRKY genes, including PnWRKY9, in P. notoginseng roots increased in response to a methyl jasmonate (MeJA) treatment and the following F. solani infection. Additionally, PnWRKY9 was functionally characterized. The PnWRKY9 protein was localized to the nucleus. The overexpression of PnWRKY9 in tobacco (Nicotiana tabacum) considerably increased the resistance to F. solani, whereas an RNAi-mediated decrease in the PnWRKY9 expression level in P. notoginseng leaves increased the susceptibility to F. solani. The RNA sequencing and hormone content analyses of PnWRKY9-overexpression tobacco revealed that PnWRKY9 and the jasmonic acid (JA) signaling pathway synergistically enhance disease resistance. The PnWRKY9 recombinant protein was observed to bind specifically to the W-box sequence in the promoter of a JA-responsive and F. solani resistance-related defensin gene (PnDEFL1). A yeast one-hybrid assay indicated that PnWRKY9 can activate the transcription of PnDEFL1. Furthermore, a co-expression assay in tobacco using β-glucuronidase (GUS) as a reporter further verified that PnWRKY9 positively regulates PnDEFL1 expression. Overall, in this study, we identified P. notoginseng WRKY TFs and demonstrated that PnWRKY9 positively affects plant defenses against the root rot pathogen. The data presented herein provide researchers with fundamental information regarding the regulatory mechanism mediating the coordinated activities of WRKY TFs and the JA signaling pathway in P. notoginseng responses to the root rot pathogen.

The WRKY10-VQ8 module safely and effectively regulates rice thermotolerance

WRKY transcription factors (TFs) play crucial roles in biotic and abiotic stress responses. However, their roles in thermal response are still largely elusive, especially in rice. In this study, we revealed the functions of WRKY10 TF and VQ8 protein containing VQ motif in rice thermotolerance. Overexpression of WRKY10 or loss of VQ8 function increases thermosensitivity, whereas conversely, overexpression of VQ8 or loss of WRKY10 function enhances thermotolerance. Overexpression of WRKY10 accelerates reactive oxygen species (ROS) accumulation in chloroplasts and apoplasts, and it also induces the expression of heat shock TF and protein genes. We also found that WRKY10 regulates nuclear DNA fragmentation and hypersensitive response by modulating NAC4 TF expression. The balance between destructive and protective responses in WRKY10-overexpression plant is more fragile and more easily broken by heat stress compared with wild type. In vitro and in vivo assays revealed that VQ8 interacts with WRKY10 and inhibits the transcription activity via repressing its DNA-binding activity. Our study demonstrates that WRKY10 negatively regulates thermotolerance by modulating the ROS balance and the hypersensitive response and that VQ8 functions antagonistically to positively regulate thermotolerance. The functional module of WRKY10-VQ8 provides safe and effective regulatory mechanisms in the heat stress response.

New aspects of DNA recognition by group II WRKY transcription factor revealed by structural and functional study of AtWRKY18 DNA binding domain

WRKY transcription factors (TFs) constitute one of the largest families of plant TFs. Based on the organization of domains and motifs, WRKY TFs are divided into three Groups (I-III). The WRKY subgroup IIa includes three representatives in A. thaliana, AtWRKY18, AtWRKY40, and AtWRKY60, that participate in biotic and abiotic stress responses. Here we present crystal structures of the DNA binding domain (DBD) of AtWRKY18 alone and in the complex with a DNA duplex containing the WRKY-recognition sequence, W-box. Subgroup IIa WRKY TFs are known to form homo and heterodimers. Our data suggest that the dimerization interface of the full-length AtWRKY18 involves contacts between the DBD subunits. DNA binding experiments and structural analysis point out novel aspects of DNA recognition by WRKY TFs. In particular, AtWRKY18-DBD preferentially binds an overlapping tandem of W-boxes accompanied by a quasi-W-box motif. The binding of DNA deforms the B-type double helix, which suggests that the DNA fragment must be prone to form a specific structure. This can explain why despite the short W-box consensus, WRKY TFs can precisely control gene expression. Finally, this first experimental structure of a Group II WRKY TF allowed us to compare Group I-III representatives.

The WRKY transcription factor family in cowpea: Genomic characterization and transcriptomic profiling under root dehydration

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most tolerant legume crops to drought and salt stresses. WRKY transcription factor (TF) family members stand out among plant transcriptional regulators related to abiotic stress tolerance. However, little information is currently available on the expression of the cowpea WRKY gene family (VuWRKY) in response to water deficit. Thus, we analyzed genomic and transcriptomic data from cowpea to identify VuWRKY members and characterize their structure and transcriptional response under root dehydration stress. Ninety-two complete VuWRKY genes were found in the cowpea genome based on their domain characteristics. They were clustered into three groups: I (15 members), II (58), and III (16), while three genes were unclassified. Domain analysis of the encoded proteins identified four major variants of the conserved heptapeptide motif WRKYGQK. In silico analysis of VuWRKY gene promoters identified eight candidate binding motifs of cis-regulatory elements, regulated mainly by six TF families associated with abiotic stress responses. Ninety-seven VuWRKY modulated splicing variants associated with 55 VuWRKY genes were identified via RNA-Seq analysis available at the Cowpea Genomics Consortium (CpGC) database. qPCR analyses showed that 22 genes are induced under root dehydration, with VuWRKY18, 21, and 75 exhibiting the most significant induction levels. Given their central role in activating signal transduction cascades in abiotic stress response, the data provide a foundation for the targeted modification of specific VuWRKY family members to improve drought tolerance in this important climate-resilient legume in the developing world and beyond.

Genome-wide identification and expression analysis response to GA3 stresses of WRKY gene family in seed hemp (Cannabis sativa L)

WRKY transcription factor is one of the largest transcription factor families in higher plants. However, the investigations of the WRKY gene family have not yet been reported in seed hemp. In the present study, we identified 39 CasWRKYs at the genome-wide level and analyzed phylogenetic relationship, chromosome location, cis-acting elements, gene structure, conserved motif, and expression pattern. Based on the gene structure and phylogenetic analyses, CasWRKY proteins were divided into 3 groups and 7 subgroups. The gene duplication investigation revealed that 6 and 5 pairs of CasWRKY genes underwent tandem and segmental duplication events, respectively. These events may contribute to the diversity and expansion of the CasWRKY gene family. The regulatory elements in the promoter regions of CasWRKYs contained diverse cis-regulatory elements, among which P-box cis-regulatory elements showed high frequency, indicating that CasWRKYs can respond to the regulation of gibberellin. The expression profiles derived from RNA-seq and qRT-PCR showed that 13 CasWRKY genes could respond to GA₃ stress and affect fiber development, as well as play significant roles in stem growth and development. This study will serve as molecular basis and practical reference for further exploring the genetic evolution and biological function of CasWRKY genes in seed hemp.

Characterization of WRKY Gene Family in Whole-Genome and Exploration of Flowering Improvement Genes in Chrysanthemum lavandulifolium

Chrysanthemum is a well-known ornamental plant with numerous uses. WRKY is a large family of transcription factors known for a variety of functions ranging from stress resistance to plant growth and development. Due to the limited research on the WRKY family in chrysanthemums, we examined them for the first time in Chrysanthemum lavandulifolium. A total of 138 ClWRKY genes were identified, which were classified into three groups. Group III in C. lavandulifolium contains 53 members, which is larger than group III of Arabidopsis. The number of introns varied from one to nine in the ClWRKY gene family. The "WRKYGQK" motif is conserved in 118 members, while other members showed slight variations. AuR and GRE responsive cis-acting elements were located in the promoter region of WRKY members, which are important for plant development and flowering induction. In addition, the W box was present in most genes; the recognition site for the WRKY gene may play a role in autoregulation and cross-regulation. The expression of the most variable 19 genes in terms of different parameters was observed at different stages. Among them, 10 genes were selected due to the presence of CpG islands, while nine genes were selected based on their close association with important Arabidopsis genes related to floral traits. ClWRKY36 and ClWRKY45 exhibit differential expression at flowering stages in the capitulum, while methylation is detected in three genes, including ClWRKY31, ClWRKY100, and ClWRKY129. Our results provide a basis for further exploration of WRKY members to find their functions in plant growth and development, especially in flowering traits.

A tau class GST, OsGSTU5, interacts with VirE2 and modulates the Agrobacterium-mediated transformation in rice

OsGSTU5 interacts and glutathionylates the VirE2 protein of Agrobacterium and its (OsGSTU5) overexpression and downregulation showed a low and high AMT efficiency in rice, respectively. During Agrobacterium-mediated transformation (AMT), T-DNA along with several virulence proteins such as VirD2, VirE2, VirE3, VirD5, and VirF enter the plant cytoplasm. VirE2 serves as a single-stranded DNA binding (SSB) protein that assists the cytoplasmic trafficking of T-DNA inside the host cell. Though the regulatory roles of VirE2 have been established, the cellular reaction of their host, especially in monocots, has not been characterized in detail. This study identified a cellular interactor of VirE2 from the cDNA library of rice. The identified plant protein encoded by the gene cloned from rice was designated OsGSTU5, it interacted specifically with VirE2 in the host cytoplasm. OsGSTU5 was upregulated during Agrobacterium infection and involved in the post-translational glutathionylation of VirE2 (gVirE2). Interestingly, the in silico analysis showed that the 'gVirE2 + ssDNA' complex was structurally less stable than the 'VirE2 + ssDNA' complex. The gel shift assay also confirmed the attenuated SSB property of gVirE2 over VirE2. Moreover, knock-down and overexpression of OsGSTU5 in rice showed increased and decreased T-DNA expression, respectively after Agrobacterium infection. The present finding establishes the role of OsGSTU5 as an important target for modulation of AMT efficiency in rice.

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

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

Identification of WRKY transcription factors responding to abiotic stresses in Brassica napus L

A total of 278 BnWRKYs were identified and analyzed. Ectopic expression of BnWRKY149 and BnWRKY217 suggests that they function in the ABA signaling pathway. WRKY transcription factors play an important role in plant development, however, their function in Brassica napus L. abiotic stress response is still unclear. In this study, a total of 278 BnWRKY transcription factors were identified from the B. napus genome data, and they were subsequently distributed in three main groups. The protein motifs and classification of BnWRKY transcription factors were analyzed, and the locations of their corresponding encoding genes were mapped on the chromosomes of B. napus. Transcriptome analysis of rapeseed seedlings exposed to drought, salt, heat, cold and abscisic acid treatment revealed that 99 BnWRKYs responded to at least one of these stresses. The expression profiles of 12 BnWRKYs were examined with qPCR and the result coincided with RNA-seq analysis. Two genes of interest, BnWRKY149 and BnWRKY217 (homologs of AtWRKY40), were overexpressed in Arabidopsis, and the corresponding proteins were located to the nucleus. Transgene plants of BnWRKY149 and BnWRKY217 were less sensitive to ABA than Arabidopsis Col-0 plants, suggesting they might play important roles in the responses of rapeseed to abiotic stress.

Genome-wide investigation of the GRAS transcription factor family in foxtail millet (Setaria italica L.)

BACKGROUND

GRAS transcription factors perform indispensable functions in various biological processes, such as plant growth, fruit development, and biotic and abiotic stress responses. The development of whole-genome sequencing has allowed the GRAS gene family to be identified and characterized in many species. However, thorough in-depth identification or systematic analysis of GRAS family genes in foxtail millet has not been conducted.

RESULTS

In this study, 57 GRAS genes of foxtail millet (SiGRASs) were identified and renamed according to the chromosomal distribution of the SiGRAS genes. Based on the number of conserved domains and gene structure, the SiGRAS genes were divided into 13 subfamilies via phylogenetic tree analysis. The GRAS genes were unevenly distributed on nine chromosomes, and members of the same subfamily had similar gene structures and motif compositions. Genetic structure analysis showed that most SiGRAS genes lacked introns. Some SiGRAS genes were derived from gene duplication events, and segmental duplications may have contributed more to GRAS gene family expansion than tandem duplications. Quantitative polymerase chain reaction showed significant differences in the expression of SiGRAS genes in different tissues and stages of fruits development, which indicated the complexity of the physiological functions of SiGRAS. In addition, exogenous paclobutrazol treatment significantly altered the transcription levels of DELLA subfamily members, downregulated the gibberellin content, and decreased the plant height of foxtail millet, while it increased the fruit weight. In addition, SiGRAS13 and SiGRAS25 may have the potential for genetic improvement and functional gene research in foxtail millet.

CONCLUSIONS

Collectively, this study will be helpful for further analysing the biological function of SiGRAS. Our results may contribute to improving the genetic breeding of foxtail millet.

Dynamic differential evolution schemes of WRKY transcription factors in domesticated and wild rice

WRKY transcription factors play key roles in stress responses, growth, and development. We previously reported on the evolution of WRKYs from unicellular green algae to land plants. To address recent evolution events, we studied three domesticated and eight wild species in the genus Oryza, an ideal model due to its long history of domestication, economic importance, and central role as a model system. We have identified prevalence of Group III WRKYs despite differences in breeding of cultivated and wild species. Same groups of WRKY genes tend to cluster together, suggesting recent, multiple duplication events. Duplications followed by divergence may result in neofunctionalizations of co-expressed WRKY genes that finely tune the regulation of target genes in a same metabolic or response pathway. WRKY genes have undergone recent rearrangements to form novel genes. Group Ib WRKYs, unique to AA genome type Oryza species, are derived from Group III genes dated back to 6.76 million years ago. Gene tree reconciliation analysis with the species tree revealed details of duplication and loss events in the 11 genomes. Selection analysis on single copy orthologs reveals the highly conserved nature of the WRKY domain and clusters of fast evolving sites under strong positive selection pressure. Also, the numbers of single copy orthologs under positive or negative selection almost evenly split. Our results provide valuable insights into the preservation and diversification of an important gene family under strong selective pressure for biotechnological improvements of the world's most valued food crop.

The DUB family in Populus: identification, characterization, evolution and expression patterns

Background

The deubiquitinase (DUB) family constitutes a group of proteases that regulate the stability or reverse the ubiquitination of many proteins in the cell. These enzymes participate in cell-cycle regulation, cell division and differentiation, diverse physiological activities such as DNA damage repair, growth and development, and response to stress. However, limited information is available on this family of genes in woody plants.

Results

In the present study, 88 DUB family genes were identified in the woody model plant Populus trichocarpa, comprising 44 PtrUBP, 3 PtrUCH, 23 PtrOTU, 4 PtrMJD, and 14 PtrJAMM genes with similar domains. According to phylogenetic analysis, the PtrUBP genes were classified into 16 groups, the PtrUCH genes into two, the PtrOTU genes into eight, the PtrMJD genes into two, and the PtrJAMM genes into seven. Members of same subfamily had similar gene structure and motif distribution characteristics. Synteny analysis of the DUB family genes from P. thrchocarpa and four other plant species provided insight into the evolutionary traits of DUB genes. Expression profiles derived from previously published transcriptome data revealed distinct expression patterns of DUB genes in various tissues. On the basis of the results of analysis of promoter cis-regulatory elements, we selected 16 representative PtrUBP genes to treatment with abscisic acid, methyl jasmonate, or salicylic acid applied as a foliar spray. The majority of PtrUBP genes were upregulated in response to the phytohormone treatments, which implied that the genes play potential roles in abiotic stress response in Populus.

Conclusions

The results of this study broaden our understanding of the DUB family in plants. Analysis of the gene structure, conserved elements, and expression patterns of the DUB family provides a solid foundation for exploration of their specific functions in Populus and to elucidate the potential role of PtrUBP gene in abiotic stress response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07844-3.

Identification, evolution and expression analysis of WRKY gene family in Eucommia ulmoides

The WRKY transcription factors is one of the largest families of transcription factors (TFs) in plants and involved in multiple biological processes. However, the role of the WRKY family had not been reported in Eucommia ulmoides. In this study, 45 WRKY genes (EuWRKY1-45) with conserved WRKY domain were identified in E. ulmoides and classified into three groups. The group II was further divided into five subgroups based on phylogenetic analysis, and each clade was well supported by the conserved motifs. All the genes were located on 34 different scaffolds respectively. A number of development-, light-, hormone-, and stress-related elements were randomly distributed in the promoter sequences of EuWRKYs. Expression profiles indicated that EuWRKY genes were involved in leaf development, and majority of EuWRKYs genes were highly expressed in leaf buds. Co-expression analysis of WRKYs suggested an intricate interplay of growth-related responses. EuWRKY4 was involved in a complex proteins interaction network. Collectively, our results provide extensive insights into the WRKY gene family, thereby contributing to the screening of additional candidate genes in E. ulmoides.

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.

Genome-Wide Analysis of Tubulin Gene Family in Cassava and Expression of Family Member FtsZ2-1 during Various Stress

Filamentous temperature-sensitive protein Z (Tubulin/FtsZ) family is a group of conserved GTP-binding (guanine nucleotide-binding) proteins, which are closely related to plant tissue development and organ formation as the major component of the cytoskeleton. According to the published genome sequence information of cassava (Manihot esculenta Crantz), 23 tubulin genes (MeTubulins) were identified, which were divided into four main groups based on their type and phylogenetic characteristics. The same grouping generally has the same or similar motif composition and exon–intron structure. Collinear analysis showed that fragment repetition event is the main factor in amplification of cassava tubulin superfamily gene. The expression profiles of MeTubulin genes in various tissue were analyzed, and it was found that MeTubulins were mainly expressed in leaf, petiole, and stem, while FtsZ2-1 was highly expressed in storage root. The qRT-PCR results of the FtsZ2-1 gene under hormone and abiotic stresses showed that indole-3-acetic acid (IAA) and gibberellin A3 (GA3) stresses could significantly increase the expression of the FtsZ2-1 gene, thereby revealing the potential role of FtsZ2-1 in IAA and GA3 stress-induced responses.

Identification and expression analysis of the WRKY gene family in Isatis indigotica

The WRKY proteins, which represent one of the largest families of transcriptional regulators in plants, play pivotal roles in regulating multiple processes of growth and development, particularly in diverse stress responses. Isatis indigotica is widely used in Traditional Chinese Medicine and is famous for its use as a dye for the color indigo. However, reports of the WRKY gene family in I. indigotica are limited. In this study, 64 IiWRKY genes encoding proteins with the complete WRKY domain were identified from genome of I. indigotica. Based on their structure and phylogenetic relationships of this gene family in I. indigotica, the IiWRKY genes were classified into three groups: Group I (n = 13), Group II (n = 35) and Group III (n = 16). Sequence alignment revealed that IiWRKY proteins harbored two variants, WRKYRQK and WRKYGKK, of the highly conserved WRKYGQK motif. The number of exons in IiWRKY genes varied from two to 14, with most of IiWRKY genes containing three exons. Investigation of gene duplication demonstrated that 10 and 14 IiWRKY genes were incorporated in tandem and segmental duplication events, respectively. Finally, the expression profiles derived from transcriptome data and quantitative real-time PCR analysis showed distinct expression patterns of these IiWRKY gene in five different organs or in response to four abiotic stresses. Taken together, our results will contribute to functional analysis of IiWRKY genes, and also provide a basis for further clarification of the molecular mechanism of stress responses in this important herb.

Highlighting the potential utility of MBP crystallization chaperone for Arabidopsis BIL1/BZR1 transcription factor-DNA complex

The maltose-binding protein (MBP) fusion tag is one of the most commonly utilized crystallization chaperones for proteins of interest. Recently, this MBP-mediated crystallization technique was adapted to Arabidopsis thaliana (At) BRZ-INSENSITIVE-LONG (BIL1)/BRASSINAZOLE-RESISTANT (BZR1), a member of the plant-specific BZR TFs, and revealed the first structure of AtBIL1/BZR1 in complex with target DNA. However, it is unclear how the fused MBP affects the structural features of the AtBIL1/BZR1-DNA complex. In the present study, we highlight the potential utility of the MBP crystallization chaperone by comparing it with the crystallization of unfused AtBIL1/BZR1 in complex with DNA. Furthermore, we assessed the validity of the MBP-fused AtBIL1/BZR1-DNA structure by performing detailed dissection of crystal packings and molecular dynamics (MD) simulations with the removal of the MBP chaperone. Our MD simulations define the structural basis underlying the AtBIL1/BZR1-DNA assembly and DNA binding specificity by AtBIL1/BZR1. The methodology employed in this study, the combination of MBP-mediated crystallization and MD simulation, demonstrates promising capabilities in deciphering the protein-DNA recognition code.