Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors under Multiple Stresses in Brassica napus

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.

The WRKY gene family in the halophyte Limonium bicolor: identification, expression analysis, and regulation of salt stress tolerance

KEY MESSAGE

63 L. bicolor WRKY genes were identified and their informatics was analyzed. The results suggested that the LbWRKY genes involved in the development and salt secretion of salt glands in L. bicolor. Salt stress, as a universal abiotic stress, severely inhibits the growth and development of plants. WRKY transcription factors play a vital role in plant growth and development, as well as in response to various stresses. Nevertheless, little is known of systematic genome-wide analysis of the WRKY genes in Limonium bicolor, a model recretohalophyte. In this study, 63 L. bicolor WRKY genes were identified (LbWRKY1-63), which were unevenly distributed across seven chromosomes and one scaffold. Based on the structural and phylogenetic characteristics, 63 LbWRKYs are divided into three main groups. Cis-elements in the LbWRKY promoters were related to growth and development, phytohormone responses, and stress responses. Colinearity analysis showed strong colinearity between LbWRKYs and GmWRKYs from soybean (Glycine max). Therefore, LbWRKY genes maybe have similar functions to GmWRKY genes. Expression analysis showed that 28 LbWRKY genes are highly expressed in roots, 9 in stems, 26 in leaves, and 12 in flowers and most LbWRKY genes responded to NaCl, ABA, and PEG6000. Silencing LbWRKY10 reduced salt gland density and salt secretion ability of leaves, and the salt tolerance of the species. Consistent with this, genes associated with salt gland development were markedly down-regulated in the LbWRKY10-silenced lines. Our findings suggested that the LbWRKY genes involved in the development and salt secretion of salt glands in L. bicolor. Our research provides new insights into the functions of the WRKY family in halophytes.

Genome-Wide Identification, Characterisation, and Evolution of the Transcription Factor WRKY in Grapevine (Vitis vinifera): New View and Update

WRKYs are a multigenic family of transcription factors that are plant-specific and involved in the regulation of plant development and various stress response processes. However, the evolution of WRKY genes is not fully understood. This family has also been incompletely studied in grapevine, and WRKY genes have been named with different numbers in different studies, leading to great confusion. In this work, 62 Vitis vinifera WRKY genes were identified based on six genomes of different cultivars. All WRKY genes were numbered according to their chromosomal location, and a complete revision of the numbering was performed. Amino acid variability between different cultivars was assessed for the first time and was greater than 5% for some WRKYs. According to the gene structure, all WRKYs could be divided into two groups: more exons/long length and fewer exons/short length. For the first time, some chimeric WRKY genes were found in grapevine, which may play a specific role in the regulation of different processes: VvWRKY17 (an N-terminal signal peptide region followed by a non-cytoplasmic domain) and VvWRKY61 (Frigida-like domain). Five phylogenetic clades A-E were revealed and correlated with the WRKY groups (I, II, III). The evolution of WRKY was studied, and we proposed a WRKY evolution model where there were two dynamic phases of complexity and simplification in the evolution of WRKY.

The translocation of a chloride channel from the Golgi to the plasma membrane helps plants adapt to salt stress

A key mechanism employed by plants to adapt to salinity stress involves maintaining ion homeostasis via the actions of ion transporters. While the function of cation transporters in maintaining ion homeostasis in plants has been extensively studied, little is known about the roles of their anion counterparts in this process. Here, we describe a mechanism of salt adaptation in plants. We characterized the chloride channel (CLC) gene AtCLCf, whose expression is regulated by WRKY transcription factor under salt stress in Arabidopsis thaliana. Loss-of-function atclcf seedlings show increased sensitivity to salt, whereas AtCLCf overexpression confers enhanced resistance to salt stress. Salt stress induces the translocation of GFP-AtCLCf fusion protein to the plasma membrane (PM). Blocking AtCLCf translocation using the exocytosis inhibitor brefeldin-A or mutating the small GTPase gene AtRABA1b/BEX5 (RAS GENES FROM RAT BRAINA1b homolog) increases salt sensitivity in plants. Electrophysiology and liposome-based assays confirm the Cl-/H+ antiport function of AtCLCf. Therefore, we have uncovered a mechanism of plant adaptation to salt stress involving the NaCl-induced translocation of AtCLCf to the PM, thus facilitating Cl- removal at the roots, and increasing the plant's salinity tolerance.

Analysis of sugar components and identification of SPS genes in citrus fruit development

Sugar is a primary determinant of citrus fruit flavour, but undergoes varied accumulation processes across different citrus varieties owing to high genetic variability. Sucrose phosphate synthase (SPS), a key enzyme in glucose metabolism, plays a crucial role in this context. Despite its significance, there is limited research on sugar component quality and the expression and regulatory prediction of SPS genes during citrus fruit development. Therefore, we analysed the sugar quality formation process in 'Kiyomi' and 'Succosa', two citrus varieties, and performed a comprehensive genome-wide analysis of citrus CsSPSs. We observed that the accumulation of sugar components significantly differs between the two varieties, with the identification of four CsSPSs in citrus. CsSPS sequences were highly conserved, featuring typical SPS protein domains. Expression analysis revealed a positive correlation between CsSPS expression and sugar accumulation in citrus fruits. However, CsSPS expression displays specificity to different citrus tissues and varieties. Transcriptome co-expression network analysis suggests the involvement of multiple transcription factors in shaping citrus fruit sugar quality through the regulation of CsSPSs. Notably, the expression levels of four CsWRKYs (CsWRKY2, CsWRKY20, CsWRKY28, CsWRKY32), were significantly positively correlated with CsSPSs and CsWRKY20 might can activate sugar accumulation in citrus fruit through CsSPS2. Collectively, we further emphasize the potential importance of CsWRKYs in citrus sugar metabolism, our findings serve as a reference for understanding sugar component formation and predicting CsSPS expression and regulation during citrus fruit development.

Transcription Factors and Their Regulatory Roles in the Male Gametophyte Development of Flowering Plants

Male gametophyte development in plants relies on the functions of numerous genes, whose expression is regulated by transcription factors (TFs), non-coding RNAs, hormones, and diverse environmental stresses. Several excellent reviews are available that address the genes and enzymes associated with male gametophyte development, especially pollen wall formation. Growing evidence from genetic studies, transcriptome analysis, and gene-by-gene studies suggests that TFs coordinate with epigenetic machinery to regulate the expression of these genes and enzymes for the sequential male gametophyte development. However, very little summarization has been performed to comprehensively review their intricate regulatory roles and discuss their downstream targets and upstream regulators in this unique process. In the present review, we highlight the research progress on the regulatory roles of TF families in the male gametophyte development of flowering plants. The transcriptional regulation, epigenetic control, and other regulators of TFs involved in male gametophyte development are also addressed.

Genome-wide identification and comprehensive analysis of WRKY transcription factor family in safflower during drought stress

The WRKY family is an important family of transcription factors in plant development and stress response. Currently, there are few reports on the WRKY gene family in safflower (Carthamus tinctorius L.). In this study, a total of 82 CtWRKY genes were identified from the safflower genome and could be classified into 3 major groups and 5 subgroups based on their structural and phylogenetic characteristics. The results of gene structure, conserved domain and motif analyses indicated that CtWRKYs within the same subfamily maintained a consistent exon/intron organization and composition. Chromosomal localization and gene duplication analysis results showed that CtWRKYs were randomly localized on 12 chromosomes and that fragment duplication and purification selection may have played an important role in the evolution of the WRKY gene family in safflower. Promoter cis-acting element analysis revealed that the CtWRKYs contain many abiotic stress response elements and hormone response elements. Transcriptome data and qRT-PCR analyses revealed that the expression of CtWRKYs showed tissue specificity and a strong response to drought stress. Notably, the expression level of the CtWRKY55 gene rapidly increased more than eightfold under drought treatment and rehydration, indicating that it may be a key gene in response to drought stress. These results provide useful insights for investigating the regulatory function of the CtWRKY gene in safflower growth and development, as well as identifying key genes for future molecular breeding programmes.

Plant salinity stress, sensing, and its mitigation through WRKY

Salinity or salt stress has deleterious effects on plant growth and development. It imposes osmotic, ionic, and secondary stresses, including oxidative stress on the plants and is responsible for the reduction of overall crop productivity and therefore challenges global food security. Plants respond to salinity, by triggering homoeostatic mechanisms that counter salt-triggered disturbances in the physiology and biochemistry of plants. This involves the activation of many signaling components such as SOS pathway, ABA pathway, and ROS and osmotic stress signaling. These biochemical responses are accompanied by transcriptional modulation of stress-responsive genes, which is mostly mediated by salt-induced transcription factor (TF) activity. Among the TFs, the multifaceted significance of WRKY proteins has been realized in many diverse avenues of plants' life including regulation of plant stress response. Therefore, in this review, we aimed to highlight the significance of salinity in a global perspective, the mechanism of salt sensing in plants, and the contribution of WRKYs in the modulation of plants' response to salinity stress. This review will be a substantial tool to investigate this problem in different perspectives, targeting WRKY and offering directions to better manage salinity stress in the field to ensure food security.

Comprehensive analyses of the citrus WRKY gene family involved in the metabolism of fruit sugars and organic acids

Sugars and organic acids are the main factors determining the flavor of citrus fruit. The WRKY transcription factor family plays a vital role in plant growth and development. However, there are still few studies about the regulation of citrus WRKY transcription factors (CsWRKYs) on sugars and organic acids in citrus fruit. In this work, a genome-wide analysis of CsWRKYs was carried out in the citrus genome, and a total of 81 CsWRKYs were identified, which contained conserved WRKY motifs. Cis-regulatory element analysis revealed that most of the CsWRKY promoters contained several kinds of hormone-responsive and abiotic-responsive cis-elements. Furthermore, gene expression analysis and fruit quality determination showed that multiple CsWRKYs were closely linked to fruit sugars and organic acids with the development of citrus fruit. Notably, transcriptome co-expression network analysis further indicated that three CsWRKYs, namely, CsWRKY3, CsWRKY47, and CsWRKY46, co-expressed with multiple genes involved in various pathways, such as Pyruvate metabolism and Citrate cycle. These CsWRKYs may participate in the metabolism of fruit sugars and organic acids by regulating carbohydrate metabolism genes in citrus fruit. These findings provide comprehensive knowledge of the CsWRKY family on the regulation of fruit quality.

The G2-Like gene family in Populus trichocarpa: identification, evolution and expression profiles

The Golden2-like (GLK) transcription factors are plant-specific transcription factors (TFs) that perform extensive and significant roles in regulating chloroplast development. Here, genome-wide identification, classification, conserved motifs, cis-elements, chromosomal locations, evolution and expression patterns of the PtGLK genes in the woody model plant Populus trichocarpa were analyzed in detail. In total, 55 putative PtGLKs (PtGLK1-PtGLK55) were identified and divided into 11 distinct subfamilies according to the gene structure, motif composition and phylogenetic analysis. Synteny analysis showed that 22 orthologous pairs and highly conservation between regions of GLK genes across P. trichocarpa and Arabidopsis were identified. Furthermore, analysis of the duplication events and divergence times provided insight into the evolutionary patterns of GLK genes. The previously published transcriptome data indicated that PtGLK genes exhibited distinct expression patterns in various tissues and different stages. Additionally, several PtGLKs were significantly upregulated under the responses of cold stress, osmotic stress, and methyl jasmonate (MeJA) and gibberellic acid (GA) treatments, implying that they might take part in abiotic stress and phytohormone responses. Overall, our results provide comprehensive information on the PtGLK gene family and elucidate the potential functional characterization of PtGLK genes in P. trichocarpa.

Genome-wide identification and bioinformatics analysis of the WRKY transcription factors and screening of candidate genes for anthocyanin biosynthesis in azalea (Rhododendron simsii)

WRKY transcription factors have been demonstrated to influence the anthocyanin biosynthesis in many plant species. However, there is limited knowledge about the structure and function of WRKY genes in the major ornamental plant azalea (Rhododendron simsii). In this study, we identified 57 RsWRKY genes in the R. simsii genome and classified them into three main groups and several subgroups based on their structural and phylogenetic characteristics. Comparative genomic analysis suggested WRKY gene family has significantly expanded during plant evolution from lower to higher species. Gene duplication analysis indicated that the expansion of the RsWRKY gene family was primarily due to whole-genome duplication (WGD). Additionally, selective pressure analysis (Ka/Ks) suggested that all RsWRKY duplication gene pairs underwent purifying selection. Synteny analysis indicated that 63 and 24 pairs of RsWRKY genes were orthologous to Arabidopsis thaliana and Oryza sativa, respectively. Furthermore, RNA-seq data was used to investigate the expression patterns of RsWRKYs, revealing that 17 and 9 candidate genes may be associated with anthocyanin synthesis at the bud and full bloom stages, respectively. These findings provide valuable insights into the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species and lay the foundation for future functional studies of WRKY genes.

Genome-wide transcriptional profiling provides clues to molecular mechanisms underlying cold tolerance in chickpea

Chickpea is an important food legume cultivated in several countries. A sudden drop in autumn temperature, freezing winter temperature, and late spring cold events result in significant losses in chickpea production. The current study used RNA sequencing of two cold tolerant (Saral) and sensitive (ILC533) Kabuli chickpea genotypes to identify cold tolerance-associated genes/pathways. A total of 200.85 million raw reads were acquired from the leaf samples by Illumina sequencing, and around 86% of the clean reads (199 million) were mapped to the chickpea reference genome. The results indicated that 3710 (1980 up- and 1730 down-regulated) and 3473 (1972 up- and 1501 down-regulated) genes were expressed differentially under cold stress in the tolerant and sensitive genotypes, respectively. According to the GO enrichment analysis of uniquely down-regulated genes under cold stress in ILC533, photosynthetic membrane, photosystem II, chloroplast part, and photosystem processes were enriched, revealing that the photosynthesis is severely sensitive to cold stress in this sensitive genotype. Many remarkable transcription factors (CaDREB1E, CaMYB4, CaNAC47, CaTCP4, and CaWRKY33), signaling/regulatory genes (CaCDPK4, CaPP2C6, CaMKK2, and CaHSFA3), and protective genes (CaCOR47, CaLEA3, and CaGST) were identified among the cold-responsive genes of the tolerant genotype. These findings would help improve cold tolerance across chickpea genotypes by molecular breeding or genetic engineering.

Comprehensive analysis of the WRKY gene family in Cucumis metuliferus and their expression profile in response to an early stage of root knot nematode infection

Root-knot nematode (RKN) is a major factor that limits the growth and productivity of important Cucumis crops, such as cucumber and melon, which lack RKN-resistance genes in their genome. Cucumis metuliferus is a wild Cucumis species that displays a high degree of RKN-resistance. WRKY transcription factors were involved in plant response to biotic stresses. However, little is known on the function of WRKY genes in response to RKN infection in Cucumis crops. In this study, Cucumis metuliferus 60 WRKY genes (CmWRKY) were identified in the C. metuliferus genome, and their conserved domains were classified into three main groups based on multiple sequence alignment and phylogenetic analysis. Synteny analysis indicated that the WRKY genes were highly conserved in Cucumis crops. Transcriptome data from of C. metuliferus roots inoculated with RKN revealed that 16 CmWRKY genes showed differential expression, of which 13 genes were upregulated and three genes were downregulated, indicating that these CmWRKY genes are important to C. metuliferus response to RKN infection. Two differentially expression CmWRKY genes (CmWRKY10 and CmWRKY28) were selected for further functional analysis. Both CmWRKY genes were localized in nucleus, indicating they may play roles in transcriptional regulation. This study provides a foundation for further research on the function of CmWRKY genes in RKN stress resistance and elucidation of the regulatory mechanism.

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.

Genetic and evolutionary dissection of melatonin response signaling facilitates the regulation of plant growth and stress responses

The expansion of gene families during evolution could generate functional diversity among their members to regulate plant growth and development. Melatonin, a phylogenetically ancient molecule, is vital for many aspects of a plant's life. Understanding the functional diversity of the molecular players involved in melatonin biosynthesis, signaling, and metabolism will facilitate the regulation of plant phenotypes. However, the molecular mechanism of melatonin response signaling elements in regulating this network still has many challenges. Here, we provide an in-depth analysis of the functional diversity and evolution of molecular components in melatonin signaling pathway. Genetic analysis of multiple mutants in plant species will shed light on the role of gene families in melatonin regulatory pathways. Phylogenetic analysis of these genes was performed, which will facilitate the identification of melatonin-related genes for future study. Based on the abovementioned signal networks, the mechanism of these genes was summarized to provide reference for studying the regulatory mechanism of melatonin in plant phenotypes. We hope that this work will facilitate melatonin research in higher plants and finely tuned spatio-temporal regulation of melatonin signaling.

Characterization of the WRKY Gene Family Related to Anthocyanin Biosynthesis and the Regulation Mechanism under Drought Stress and Methyl Jasmonate Treatment in Lycoris radiata

Lycoris radiata, belonging to the Amaryllidaceae family, is a well-known Chinese traditional medicinal plant and susceptible to many stresses. WRKY proteins are one of the largest families of transcription factors (TFs) in plants and play significant functions in regulating physiological metabolisms and abiotic stress responses. The WRKY TF family has been identified and investigated in many medicinal plants, but its members and functions are not identified in L. radiata. In this study, a total of 31 L. radiata WRKY (LrWRKY) genes were identified based on the transcriptome-sequencing data. Next, the LrWRKYs were divided into three major clades (Group I-III) based on the WRKY domains. A motif analysis showed the members within same group shared a similar motif component, indicating a conservational function. Furthermore, subcellular localization analysis exhibited that most LrWRKYs were localized in the nucleus. The expression pattern of the LrWRKY genes differed across tissues and might be important for Lycoris growth and flower development. There were large differences among the LrWRKYs based on the transcriptional levels under drought stress and MeJA treatments. Moreover, a total of 18 anthocyanin components were characterized using an ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis and pelargonidin-3-O-glucoside-5-O-arabinoside as well as cyanidin-3-O-sambubioside were identified as the major anthocyanin aglycones responsible for the coloration of the red petals in L. radiata. We further established a gene-to-metabolite correlation network and identified LrWRKY3 and LrWRKY27 significant association with the accumulation of pelargonidin-3-O-glucoside-5-O-arabinoside in the Lycoris red petals. These results provide an important theoretical basis for further exploring the molecular basis and regulatory mechanism of WRKY TFs in anthocyanin biosynthesis and in response to drought stress and MeJA treatment.

Combined Transcriptome and Metabolome Profiling Provide Insights into Cold Responses in Rapeseed (Brassica napus L.) Genotypes with Contrasting Cold-Stress Sensitivity

Low temperature is a major environmental factor, which limits rapeseed (Brassica napus L.) growth, development, and productivity. So far, the physiological and molecular mechanisms of rapeseed responses to cold stress are not fully understood. Here, we explored the transcriptome and metabolome profiles of two rapeseed genotypes with contrasting cold responses, i.e., XY15 (cold-sensitive) and GX74 (cold-tolerant). The global metabolome profiling detected 545 metabolites in siliques of both genotypes before (CK) and after cold-stress treatment (LW). The contents of several sugar metabolites were affected by cold stress with the most accumulated saccharides being 3-dehydro-L-threonic acid, D-xylonic acid, inositol, D-mannose, D-fructose, D-glucose, and L-glucose. A total of 1943 and 5239 differentially expressed genes were identified from the transcriptome sequencing in XY15CK_vs_XY15LW and GX74CK_vs_GX74LW, respectively. We observed that genes enriched in sugar metabolism and biosynthesis-related pathways, photosynthesis, reactive oxygen species scavenging, phytohormone, and MAPK signaling were highly expressed in GX74LW. In addition, several genes associated with cold-tolerance-related pathways, e.g., the CBF-COR pathway and MAPK signaling, were specifically expressed in GX74LW. Contrarily, genes in the above-mentioned pathways were mostly downregulated in XY15LW. Thus, our results indicate the involvement of these pathways in the differential cold-stress responses in XY15 and GX74.

The Banana MaWRKY18, MaWRKY45, MaWRKY60 and MaWRKY70 Genes Encode Functional Transcription Factors and Display Differential Expression in Response to Defense Phytohormones

WRKY transcription factors (TFs) play key roles in plant defense responses through phytohormone signaling pathways. However, their functions in tropical fruit crops, especially in banana, remain largely unknown. Several WRKY genes from the model plants rice (OsWRKY45) and Arabidopsis (AtWRKY18, AtWRKY60, AtWRKY70) have shown to be attractive TFs for engineering disease resistance. In this study, we isolated four banana cDNAs (MaWRKY18, MaWRKY45, MaWRKY60, and MaWRKY70) with homology to these rice and ArabidopsisWRKY genes. The MaWRKY cDNAs were isolated from the wild banana Musa acuminata ssp. malaccensis, which is resistant to several diseases of this crop and is a progenitor of most banana cultivars. The deduced amino acid sequences of the four MaWRKY cDNAs revealed the presence of the conserved WRKY domain of ~60 amino acids and a zinc-finger motif at the N-terminus. Based on the number of WRKY repeats and the structure of the zinc-finger motif, MaWRKY18 and MaWRKY60 belong to group II of WRKY TFs, while MaWRKY45 and MaWRKY70 are members of group III. Their corresponding proteins were located in the nuclei of onion epidermal cells and were shown to be functional TFs in yeast cells. Moreover, expression analyses revealed that the majority of these MaWRKY genes were upregulated by salicylic acid (SA) or methyl jasmonate (MeJA) phytohormones, although the expression levels were relatively higher with MeJA treatment. The fact that most of these banana WRKY genes were upregulated by SA or MeJA, which are involved in systemic acquired resistance (SAR) or induced systemic resistance (ISR), respectively, make them interesting candidates for bioengineering broad-spectrum resistance in this crop.

Genome-Wide Identification and Characterization of G2-Like Transcription Factor Genes in Moso Bamboo (Phyllostachys edulis)

G2-like (GLK) transcription factors contribute significantly and extensively in regulating chloroplast growth and development in plants. This study investigated the genome-wide identification, phylogenetic relationships, conserved motifs, promoter cis-elements, MCScanX, divergence times, and expression profile analysis of PeGLK genes in moso bamboo (Phyllostachys edulis). Overall, 78 putative PeGLKs (PeGLK1–PeGLK78) were identified and divided into 13 distinct subfamilies. Each subfamily contains members displaying similar gene structure and motif composition. By synteny analysis, 42 orthologous pairs and highly conserved microsynteny between regions of GLK genes across moso bamboo and maize were found. Furthermore, an analysis of the divergence times indicated that PeGLK genes had a duplication event around 15 million years ago (MYA) and a divergence happened around 38 MYA between PeGLK and ZmGLK. Tissue-specific expression analysis showed that PeGLK genes presented distinct expression profiles in various tissues, and many members were highly expressed in leaves. Additionally, several PeGLKs were significantly up-regulated under cold stress, osmotic stress, and MeJA and GA treatment, implying that they have a likelihood of affecting abiotic stress and phytohormone responses in plants. The results of this study provide a comprehensive understanding of the moso bamboo GLK gene family, as well as elucidating the potential functional characterization of PeGLK genes.

Identification of SNAT Family Genes Suggests GhSNAT3D Functional Reponse to Melatonin Synthesis Under Salinity Stress in Cotton

Serotonin N-acetyltransferase (SNAT) is a key enzyme in the biosynthesis of melatonin, and plays an important role in the regulation of melatonin synthesis. The study of SNAT is of great significance to understand the function of melatonin. In this study, we analyzed the structural characteristics, phylogenetic relationship, gene structure, expression pattern, evolutionary relationship and stress response of the members of the SNAT gene family in upland cotton through bioinformatics. A putative Serotonin n-acetyltransferase gene GhSNAT3D was identified, and preliminarily function of GhSNAT3D was verified by virus-induced gene silencing. We identified a total of 52 SNAT genes in the whole genome of G. hirsutum, and part of the GhSNATs were regulated by exogenous melatonin. The content of melatonin, antioxidant enzyme activity and Ca²⁺ content of GhSNAT3D gene silenced plants decreased, and the salt tolerance of GhSNAT3D gene silenced plants was reduced. Exogenous melatonin supplementation restored the salt tolerance of GhSNAT3D gene silenced plants. GhSNAT3D may interact with GhSNAT25D and ASMT to regulate melatonin synthesis. This study provided an important basis for further study on the regulation of melatonin in cotton against abiotic stress.

WRKY genes provide novel insights into their role against Ralstonia solanacearum infection in cultivated peanut (Arachis hypogaea L.)

As one of the most important and largest transcription factors, WRKY plays a critical role in plant disease resistance. However, little is known regarding the functions of the WRKY family in cultivated peanuts (Arachis hypogaea L.). In this study, a total of 174 WRKY genes (AhWRKY) were identified from the genome of cultivated peanuts. Phylogenetic analysis revealed that AhWRKY proteins could be divided into four groups, including 35 (20.12%) in group I, 107 (61.49%) in group II, 31 (17.82%) in group III, and 1 (0.57%) in group IV. This division is further supported by the conserved motif compositions and intron/exon structures. All AhWRKY genes were unevenly located on all 20 chromosomes, among which 132 pairs of fragment duplication and seven pairs of tandem duplications existed. Eighteen miRNAs were found to be targeting 50 AhWRKY genes. Most AhWRKY genes from some groups showed tissue-specific expression. AhWRKY46, AhWRKY94, AhWRKY156, AhWRKY68, AhWRKY41, AhWRKY128, AhWRKY104, AhWRKY19, AhWRKY62, AhWRKY155, AhWRKY170, AhWRKY78, AhWRKY34, AhWRKY12, AhWRKY95, and AhWRKY76 were upregulated in ganhua18 and kainong313 genotypes after Ralstonia solanacearum infection. Ten AhWRKY genes (AhWRKY34, AhWRKY76, AhWRKY78, AhWRKY120, AhWRKY153, AhWRKY155, AhWRKY159, AhWRKY160, AhWRKY161, and AhWRKY162) from group III displayed different expression patterns in R. solanacearum sensitive and resistant peanut genotypes infected with the R. solanacearum. Two AhWRKY genes (AhWRKY76 and AhWRKY77) from group III obtained the LRR domain. AhWRKY77 downregulated in both genotypes; AhWRKY76 showed lower-higher expression in ganhua18 and higher expression in kainong313. Both AhWRKY76 and AhWRKY77 are targeted by ahy-miR3512, which may have an important function in peanut disease resistance. This study identified candidate WRKY genes with possible roles in peanut resistance against R. solanacearum infection. These findings not only contribute to our understanding of the novel role of WRKY family genes but also provide valuable information for disease resistance in A. hypogaea.

Identification and Expression Profiling of WRKY Family Genes in Sugarcane in Response to Bacterial Pathogen Infection and Nitrogen Implantation Dosage

WRKY transcription factors (TFs) are essential players in different signaling cascades and regulatory networks involved in defense responses to various stressors. This study systematically analyzed and characterized WRKY family genes in the Saccharum spp. hybrid R570 and their expression in two sugarcane cultivars LCP85-384 (resistant to leaf scald) and ROC20 (susceptible to leaf scald) in response to bacterial pathogen infection and nitrogen implantation dosage. A total of 53 ShWRKY genes with 66 alleles were systematically identified in R570 based on the query sequence SsWRKY in S. spontaneum AP85-441. All ShRWKY alleles were further classified into four groups with 11 (16.7%) genes in group I, 36 (54.5%) genes in group II, 18 (27.3%) genes in group III, and 1 (1.5%) gene in group IV. Among them, 4 and 11 ShWRKY gene pairs displayed tandem and segmental duplication events, respectively. The ShWRKY genes exhibited conserved DNA-binding domains, which were accompanied by variations in introns, exons, and motifs. RT-qPCR analysis of two sugarcane cultivars triggered by Xanthomonas albilineans (Xa) revealed that four genes, ShWRKY13-2/39-1/49-3/125-3, exhibited significant upregulation in leaf scald-resistant LCP85-384. These WRKY genes were downregulated or unchanged in ROC20 at 24-72 h post-inoculation, suggesting that they play an important role in defense responses to Xa infection. Most of the 12 tested ShWRKYs, ShWRKY22-1/49-3/52-1 in particular, functioned as negative regulators in the two cultivars in response to a range of nitrogen (N) implantation doses. A total of 11 ShWRKY proteins were predicted to interact with each other. ShWRKY43 and ShWRKY49-3 are predicted to play core roles in the interaction network, as indicated by their interaction with six other ShWRKY proteins. Our results provide important candidate gene resources for the genetic improvement of sugarcane and lay the foundation for further functional characterization of ShWRKY genes in response to coupling effects of Xa infection and different N levels.

Genetic and Physiological Responses to Heat Stress in Brassica napus

Given the current rise in global temperatures, heat stress has become a major abiotic challenge affecting the growth and development of various crops and reducing their productivity. Brassica napus, the second largest source of vegetable oil worldwide, experiences a drastic reduction in seed yield and quality in response to heat. This review outlines the latest research that explores the genetic and physiological impact of heat stress on different developmental stages of B. napus with a special attention to the reproductive stages of floral progression, organogenesis, and post flowering. Several studies have shown that extreme temperature fluctuations during these crucial periods have detrimental effects on the plant and often leading to impaired growth and reduced seed production. The underlying mechanisms of heat stress adaptations and associated key regulatory genes are discussed. Furthermore, an overview and the implications of the polyploidy nature of B. napus and the regulatory role of alternative splicing in forming a priming-induced heat-stress memory are presented. New insights into the dynamics of epigenetic modifications during heat stress are discussed. Interestingly, while such studies are scarce in B. napus, opposite trends in expression of key genetic and epigenetic components have been identified in different species and in cultivars within the same species under various abiotic stresses, suggesting a complex role of these genes and their regulation in heat stress tolerance mechanisms. Additionally, omics-based studies are discussed with emphasis on the transcriptome, proteome and metabolome of B. napus, to gain a systems level understanding of how heat stress alters its yield and quality traits. The combination of omics approaches has revealed crucial interactions and regulatory networks taking part in the complex machinery of heat stress tolerance. We identify key knowledge gaps regarding the impact of heat stress on B. napus during its yield determining reproductive stages, where in-depth analysis of this subject is still needed. A deeper knowledge of heat stress response components and mechanisms in tissue specific models would serve as a stepping-stone to gaining insights into the regulation of thermotolerance that takes place in this important crop species and support future breeding of heat tolerant crops.

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 Identification and Analysis of the WRKY Gene Family and Cold Stress Response in Acer truncatum

WRKY transcription factors constitute one of the largest gene families in plants and are involved in many biological processes, including growth and development, physiological metabolism, and the stress response. In earlier studies, the WRKY gene family of proteins has been extensively studied and analyzed in many plant species. However, information on WRKY transcription factors in Acer truncatum has not been reported. In this study, we conducted genome-wide identification and analysis of the WRKY gene family in A. truncatum, 54 WRKY genes were unevenly located on all 13 chromosomes of A. truncatum, the highest number was found in chromosomes 5. Phylogenetic relationships, gene structure, and conserved motif identification were constructed, and the results affirmed 54 AtruWRKY genes were divided into nine subgroup groups. Tissue species analysis of AtruWRKY genes revealed which were differently exhibited upregulation in flower, leaf, root, seed and stem, and the upregulation number were 23, 14, 34, 18, and 8, respectively. In addition, the WRKY genes expression in leaf under cold stress showed that more genes were significantly expressed under 0, 6 and 12 h cold stress. The results of this study provide a new insight the regulatory function of WRKY genes under abiotic and biotic stresses.

Genome-wide identification and characteristic analysis of the downstream melatonin metabolism gene GhM2H in Gossypium hirsutum L

Background

Melatonin 2-hydroxylase (M2H) is the first enzyme in the catabolism pathway of melatonin, which catalyzes the production of 2-hydroxymelatonin (2-OHM) from melatonin. The content of 2-hydroxymelatonin in plants is much higher than that of melatonin. So M2H may be a key enzyme in the metabolic pathway of melatonin.

Method

We conducted a systematic analysis of the M2H gene family in Gossypium hirsutum based on the whole genome sequence by integrating the structural characteristics, phylogenetic relationships, expression profile, and biological stress of the members of the Gossypium hirsutum M2H gene family.

Result

We identified 265 M2H genes in the whole genome of Gossypium hirsutum, which were divided into 7 clades (clades I-VII) according to phylogenetic analysis. Most M2H members in each group had similar motif composition and gene structure characteristics. More than half of GhM2H members contain ABA-responsive elements and MeJA-responsive elements. Under different stress conditions, the expression levels of the gene changed, indicating that GhM2H members were involved in the regulation of abiotic stress. Some genes in the GhM2H family were involved in regulating melatonin levels in cotton under salt stress, and some genes were regulated by exogenous melatonin.

Conclusion

This study is helpful to explore the function of GhM2H, the downstream metabolism gene of melatonin in cotton, and lay the foundation for better exploring the molecular mechanism of melatonin improving cotton's response to abiotic stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-021-00358-y.

Genome-Wide Analysis of the YABBY Transcription Factor Family in Rapeseed (Brassica napus L.)

The YABBY family of plant-specific transcription factors play important regulatory roles during the development of leaves and floral organs, but their functions in Brassica species are incompletely understood. Here, we identified 79 YABBY genes from Arabidopsis thaliana and five Brassica species (B. rapa, B. nigra, B. oleracea, B. juncea, and B. napus). A phylogenetic analysis of YABBY proteins separated them into five clusters (YAB1–YAB5) with representatives from all five Brassica species, suggesting a high degree of conservation and similar functions within each subfamily. We determined the gene structure, chromosomal location, and expression patterns of the 21 BnaYAB genes identified, revealing extensive duplication events and gene loss following polyploidization. Changes in exon–intron structure during evolution may have driven differentiation in expression patterns and functions, combined with purifying selection, as evidenced by K_a/K_s values below 1. Based on transcriptome sequencing data, we selected nine genes with high expression at the flowering stage. qRT-PCR analysis further indicated that most BnaYAB family members are tissue-specific and exhibit different expression patterns in various tissues and organs of B. napus. This preliminary study of the characteristics of the YABBY gene family in the Brassica napus genome provides theoretical support and reference for the later functional identification of the family genes.

Root system architecture, physiological and transcriptional traits of soybean (Glycine max L.) in response to water deficit: A review

Drought stress is the main limiting factor for global soybean growth and production. Genetic improvement for water and nutrient uptake efficiency is critical to advance tolerance and enable more sustainable and resilient production, underpinning yield growth. The identification of quantitative traits and genes related to water and nutrient uptake will enhance our understanding of the mechanisms of drought tolerance in soybean. This review summarizes drought stress in the context of the physiological traits that enable effective acclimation, with a particular focus on roots. Genes controlling root system architecture play an important role in water and nutrient availability, and therefore important targets for breeding strategies to improve drought tolerance. This review highlights the candidate genes that have been identified as regulators of important root traits and responses to water stress. Progress in our understanding of the function of particular genes, including GmACX1, GmMS and GmPEPCK are discussed in the context of developing a system-based platform for genetic improvement of drought tolerance in soybean.

Omics: The way forward to enhance abiotic stress tolerance in Brassica napus L

Plant abiotic stresses negative affects growth and development, causing a massive reduction in global agricultural production. Rapeseed (Brassica napus L.) is a major oilseed crop because of its economic value and oilseed production. However, its productivity has been reduced by many environmental adversities. Therefore, it is a prime need to grow rapeseed cultivars, which can withstand numerous abiotic stresses. To understand the various molecular and cellular mechanisms underlying the abiotic stress tolerance and improvement in rapeseed, omics approaches have been extensively employed in recent years. This review summarized the recent advancement in genomics, transcriptomics, proteomics, metabolomics, and their imploration in abiotic stress regulation in rapeseed. Some persisting bottlenecks have been highlighted, demanding proper attention to fully explore the omics tools. Further, the potential prospects of the CRISPR/Cas9 system for genome editing to assist molecular breeding in developing abiotic stress-tolerant rapeseed genotypes have also been explained. In short, the combination of integrated omics, genome editing, and speed breeding can alter rapeseed production worldwide.

Genome-wide identification and functional characterization of the Camelina sativa WRKY gene family in response to abiotic stress

Background

WRKY transcription factors are a superfamily of regulators involved in diverse biological processes and stress responses in plants. However, there is limited knowledge about the WRKY family in camelina (Camelina sativa), an important Brassicaceae oil crop with strong tolerance for various stresses. Here, a genome-wide characterization of WRKY proteins is performed to examine their gene structures, phylogenetics, expression, conserved motif organizations, and functional annotation to identify candidate WRKYs that mediate stress resistance regulation in camelinas.

Results

A total of 242 CsWRKY proteins encoded by 224 gene loci distributed unevenly over the chromosomes were identified, and they were classified into three groups by phylogenetic analysis according to their WRKY domains and zinc finger motifs. The 15 CsWRKY gene loci generated 33 spliced variants. Orthologous WRKY gene pairs were identified, with 173 pairs in the C. sativa and Arabidopsis genomes as well as 282 pairs in the C. sativa and B. napus genomes, respectively. A total of 137 segmental duplication events were observed, but there was no tandem duplication in the camelina genome. Ten major conserved motifs were examined, with WRKYGQK being the most conserved, and several variants were present in many CsWRKYs. Expression analysis revealed that 50% more CsWRKY genes were expressed constitutively, and a set of them displayed tissue-specific expression. Notably, 11 CsWRKY genes exhibited significant expression changes in seedlings under cold, salt, and drought stresses, showing a preferentially inducible expression pattern in response to the stress.

Conclusions

The present article describes a detailed analysis of the CsWRKY gene family and its expression profiles in 12 tissues and under several stress conditions. Segmental duplication is the major force underlying the broad expansion of this gene family, and a strong purifying pressure occurred for CsWRKY proteins during their evolution. CsWRKY proteins play important roles in plant development, with differential functions in different tissues. Exceptionally, eleven CsWRKYs, particularly five alternative spliced isoforms, were found to be the possible key players in mediating plant responses to various stresses. Overall, our results provide a foundation for understanding the roles of CsWRKYs and the precise mechanism through which CsWRKYs regulate high stress resistance as well as the development of stress tolerance cultivars among Cruciferae crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07189-3.

Genome-wide analysis of the WRKY gene family in the cucumber genome and transcriptome-wide identification of WRKY transcription factors that respond to biotic and abiotic stresses

BACKGROUND

Cucumber (Cucumis sativus L.) is an economically important vegetable crop species. However, it is susceptible to various abiotic and biotic stresses. WRKY transcription factors play important roles in plant growth and development, particularly in the plant response to biotic and abiotic stresses. However, little is known about the expression pattern of WRKY genes under different stresses in cucumber.

RESULTS

In the present study, an analysis of the new assembly of the cucumber genome (v3.0) allowed the identification of 61 cucumber WRKY genes. Phylogenetic and synteny analyses were performed using related species to investigate the evolution of the cucumber WRKY genes. The 61 CsWRKYs were classified into three main groups, within which the gene structure and motif compositions were conserved. Tissue expression profiles of the WRKY genes demonstrated that 24 CsWRKY genes showed constitutive expression (FPKM > 1 in all samples), and some WRKY genes showed organ-specific expression, suggesting that these WRKYs might be important for plant growth and organ development in cucumber. Importantly, analysis of the CsWRKY gene expression patterns revealed that five CsWRKY genes strongly responded to both salt and heat stresses, 12 genes were observed to be expressed in response to infection from downy mildew and powdery mildew, and three CsWRKY genes simultaneously responded to all treatments analysed. Some CsWRKY genes were observed to be induced/repressed at different times after abiotic or biotic stress treatment, demonstrating that cucumber WRKY genes might play different roles during different stress responses and that their expression patterns vary in response to stresses.

CONCLUSIONS

Sixty-one WRKY genes were identified in cucumber, and insight into their classification, evolution, and expression patterns was gained in this study. Responses to different abiotic and biotic stresses in cucumber were also investigated. Our results provide a better understanding of the function of CsWRKY genes in improving abiotic and biotic stress resistance in cucumber.

Comprehensive Transcriptomic Analysis for Developing Seeds of a Synthetic Brassica Hexaploid

Polyploidization is a universal phenomenon in plants and plays a crucial role in evolution. In this study, the transcriptomes of developing seeds of a synthetic Brassica hexaploid and its parents (B. rapa and B. carinata) were analyzed to find the gene expression changes in hexaploid seeds. There were 3166 and 3893 DEGs between the Brassica hexaploid and its parents at the full-size stage and mature stage, respectively, most of which were upregulated in hexaploid seeds compared to its parents. At the mature stage, the hexaploid seeds showed a greater difference from its parents. These DEGs had a wide range of functions, which may account for the physiological and morphological differences between the Brassica hexaploid and its parents. The KEGG pathway analysis revealed that hexaploid seeds had higher levels of expression of genes involved in metabolic pathways, RNA transport and biosynthesis of secondary metabolites, and the expression levels in the photosynthesis-related pathways were significantly higher than those in B. rapa. Transgressive expression was the main non-additive expression pattern of the Brassica hexaploid. The gene expression difference between the Brassica hexaploid and its paternal parent was more significant than that with its maternal parent, which may be due in part to the cytoplasmic and maternal effects. Moreover, transcription factor genes, such as G2-like, MYB and mTERF, were highly expressed in hexaploid seeds, possibly promoting their resistance to stress. Our results may provide valuable insights into the adaptation mechanisms of polyploid plants.

Transcription factor BnaA9.WRKY47 contributes to the adaptation of Brassica napus to low boron stress by up-regulating the boric acid channel gene BnaA3.NIP5;1

Boron (B) deficiency is one of the major causes of growth inhibition and yield reduction in Brassica napus (B. napus). However, the molecular mechanisms of low B adaptation in B. napus are largely unknown. Here, fifty-one BnaWRKY transcription factors were identified as responsive to B deficiency in B. napus, in which BnaAn.WRKY26, BnaA9.WRKY47, BnaA1.WKRY53 and BnaCn.WRKY57 were tested in yeast one-hybrid assays and showed strong binding activity with conserved sequences containing a W box in the promoters of the B transport-related genes BnaNIP5;1s and BnaBOR1s. Green fluorescent protein fused to the target protein demonstrated the nuclear localization of BnaA9.WRKY47. CRISPR/Cas9-mediated knockout lines of BnaA9.WRKY47 in B. napus had increased sensitivity to low B and lower contents of B than wild-type plants. In contrast, overexpression of BnaA9.WRKY47 enhanced the adaptation to low B with higher B contents in tissues than in wild-type plants. Consistent with the phenotypic response and B accumulation in these transgenic lines, the transcription activity of BnaA3.NIP5;1, a B efficiency candidate gene, was decreased in the knockout lines but was significantly increased in the overexpressing lines under low B conditions. Electrophoretic mobility shift assays, transient expression experiments in tobacco and in situ hybridizations showed that BnaA9.WRKY47 directly activated BnaA3.NIP5;1 expression through binding to the specific cis-element. Taken together, our findings support BnaWRKYs as new participants in response to low B, and BnaA9.WRKY47 contributes to the adaptation of B. napus to B deficiency through up-regulating BnaA3.NIP5;1 expression to facilitate efficient B uptake.

Genome-wide investigation of WRKY transcription factors in sweet osmanthus and their potential regulation of aroma synthesis

WRKY transcription factors, one of the largest transcription factor families, play important roles in regulating the synthesis of secondary metabolites. In sweet osmanthus (Osmanthus fragrans), the monoterpenes have been demonstrated as the most important volatile compounds, and the W-box, which is the cognate binding site of WRKY transcription factors, could be identified in most of the terpene-synthesis-related genes' promoters. However, the role of the WRKY family in terpene synthesis in sweet osmanthus has rarely been examined. In this study, 154 WRKY genes with conserved WRKY domain were identified and classified into three groups. The group II was further divided into five subgroups, and almost all members of IId contained a plant zinc cluster domain. Eight OfWRKYs (OfWRKY7/19/36/38/42/84/95/139) were screened from 20 OfWRKYs for their flower-specific expression patterns in different tissues. Simultaneously, the expression patterns of OfWRKYs and emission patterns of volatile compounds during the flowering process were determined and gas chromatography-mass spectrometry results showed that monoterpenes, such as linalool and ocimene, accounted for the highest proportion, contributing to the floral scent of sweet osmanthus in two cultivars. In addition, correlation analysis revealed the expression patterns of OfWRKYs (OfWRKY7/19/36/139) were each correlated with distinct monoterpenes (linalool, linalool derivatives, ocimene and ocimene derivatives). Subcellular localization analysis showed that p35S::GFP-OfWRKY7/38/95/139 were localized in the nucleus and OfWRKY139 had very strong transactivation activity. Collectively, the results indicated potential roles of OfWRKY139 and OfWRKYs with plant zinc cluster domain in regulating synthesis of aromatic compounds in sweet osmanthus, laying the foundation for use of OfWRKYs to improve the aroma of ornamental plants.

Identification of WRKY gene family and characterization of cold stress-responsive WRKY genes in eggplant

BACKGROUND

WRKY proteins play a vital role in the plants response to different stresses, growth and development. Studies of WRKY proteins have been mainly focused on model plant Arabidopsis and a few other vegetable plants. However, the systematical study of eggplant WRKY transcription factor superfamily is scarce.

METHODS

Bioinformatics has been used to identify and characterize the eggplant WRKY gene family. For the exploration of the differentially expressed WRKY genes, two cultivars with different cold-tolerance were used. Finally, we performed a virus-induced gene silencing (VIGS) experiment to verify the functions of SmWRKY26 and SmWRKY32.

RESULTS

Fifty eight (58) genes encoding eggplant WRKY proteins were identified through searching the eggplant genome. Eggplant WRKY proteins could be classified into three groups or seven subgroups in accordance with other plants. WRKY variants were identified from the eggplant. Gene structure analysis showed that the number of intron in eggplant WRKY family was from 0 to 11, with an average of 4.4. Conserved motif analysis suggested that WRKY DNA-binding domain was conserved in eggplant WRKY proteins. Furthermore, RNA-seq data showed that WRKY genes were differentially expressed in eggplant response to cold stress. By using VIGS, the two differentially expressed genes-SmWRKY26 and SmWRKY32 were verified in response to cold stress.

DISCUSSIONS

This study provides a foundation for further exploring the functions of WRKY proteins in eggplant response to stresses and eggplant genetic improvement in stresses.

Genome-Wide Identification and Expression Profile Analysis of WRKY Family Genes in the Autopolyploid Saccharum spontaneum

WRKY is one of the largest transcription factor families in plants and plays important roles in the regulation of developmental and physiological processes. To date, the WRKY gene family has not been identified in Saccharum species because of its complex polyploid genome. In this study, a total of 294 sequences for 154 SsWRKY genes were identified in the polyploid Saccharum spontaneum genome and then named on the basis of their chromosome locations, including 13 (8.4%) genes with four alleles, 29 (18.8%) genes with three alleles and 41 (26.6%) genes with two alleles. Among them, 73.8% and 16.0% of the SsWRKY genes originated from segmental duplications and tandem duplications, respectively. The WRKY members exhibited conserved gene structures and amino acid sequences among the allelic haplotypes, which were accompanied by variations in intron sizes. Phylogenetic and collinearity analyses revealed that 27 SsWRKYs originated after the split of sorghum and Saccharum, resulting in a significantly higher number of WRKYs in sugarcane than in the proximal diploid species sorghum. The analysis of RNA-seq data revealed that SsWRKYs' expression profiles in 46 different samples including different developmental stages revealed distinct temporal and spatial patterns with 52 genes expressed in all tissues, four genes not expressed in any tissues and 21 SsWRKY genes likely to be involved in photosynthesis. The comprehensive analysis of SsWRKYs' expression will provide an important and valuable foundation for further investigation of the regulatory mechanisms of WRKYs in physiological roles in sugarcane S. spontaneum.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

Genome-wide identification of cold responsive transcription factors in Brassica napus L

BACKGROUND

Cold stress is one of the primary environmental factors that affect plant growth and productivity, especially for crops like Brassica napus that live through cold seasons. Till recently, although a number of genes and pathways involved in B. napus cold response have been revealed by independent studies, a genome-wide identification of the key regulators and the regulatory networks is still lack. In this study, we investigated the transcriptomes of cold stressed semi-winter and winter type rapeseeds in short day condition, mainly with the purpose to systematically identify the functional conserved transcription factors (TFs) in cold response of B. napus.

RESULTS

Global modulation of gene expression was observed in both the semi-winter type line (158A) and the winter type line (SGDH284) rapeseeds, in response to a seven-day chilling stress in short-day condition. Function analysis of differentially expressed genes (DEGs) revealed enhanced stresses response mechanisms and inhibited photosynthesis in both lines, as well as a more extensive inhibition of some primary biological processes in the semi-winter type line. Over 400 TFs were differentially expressed in response to cold stress, including 56 of them showed high similarity to the known cold response TFs and were consistently regulated in 158A and SGDH284, as well as 25 TFs which targets were over-represented in the total DEGs. A further investigation based on their interactions indicated the critical roles of several TFs in cold response of B. napus.

CONCLUSION

In summary, our results revealed the alteration of gene expression in cold stressed semi-winter and winter ecotype B. napus lines and provided a valuable collection of candidate key regulators involved in B. napus response to cold stress, which could expand our understanding of plant stress response and benefit the future improvement of the breed of rapeseeds.

A Comprehensive Transcriptome-Wide Identification and Screening of WRKY Gene Family Engaged in Abiotic Stress in Glycyrrhiza glabra

The study reports 147 full-length WRKY genes based on the transcriptome analysis of Glycyrrhiza genus (G. glabra and G. uralensis). Additional motifs in G. glabra included DivIVA (GgWRKY20) and SerS Superfamily (GgWRKY21) at the C-terminal, and Coat family motifs (GgWRKY55) at the N-terminal of the proteins, while Exo70 exo cyst complex subunit of 338 amino acid (GuWRKY9) was present at the N-terminal of G. uralensis only. Plant Zn cluster super-family domain (17 WRKYs) and bZIP domain (2 WRKYs) were common between the two species. Based on the number of WRKY domains, sequence alignment and phylogenesis, the study identified GuWRKY27 comprising of 3 WRKY domains in G. uralensis and a new subgroup-IIf (10 members), having novel zinc finger pattern (C-X4-C-X22-HXH) in G. glabra. Multiple WRKY binding domains (1-11) were identified in the promoter regions of the GgWRKY genes indicating strong interacting network between the WRKY proteins. Tissue-specific expression of 25 GgWRKYs, under normal and treated conditions, revealed 11 of the 18 induction factor triggered response corroborating to response observed in AtWRKYs. The study identified auxin-responsive GgWRKY 55 & GgWRKY38; GA3 responsive GgWRKYs15&59 in roots and GgWRKYs8, 20, 38, 57 &58 in the shoots of the treated plant. GgWRKYs induced under various stresses included GgWRKY33 (cold), GgWRKY4 (senescence), GgWRKYs2, 28 & 33 (salinity) and GgWRKY40 (wounding). Overall, 23 GgWRKYs responded to abiotic stress, and 17 WRKYs were induced by hormonal signals. Of them 13 WRKYs responded to both suggesting inter-connection between hormone signalling and stress response. The present study will help in understanding the transcriptional reprogramming, protein-protein interaction and cross-regulation during stress and other physiological processes in the plant.

Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L

BACKGROUND

Plant homeodomain (PHD) finger proteins are widely present in all eukaryotes and play important roles in chromatin remodeling and transcriptional regulation. The PHD finger can specifically bind a number of histone modifications as an "epigenome reader", and mediate the activation or repression of underlying genes. Many PHD finger genes have been characterized in animals, but only few studies were conducted on plant PHD finger genes to this day. Brassica rapa (AA, 2n = 20) is an economically important vegetal, oilseed and fodder crop, and also a good model crop for functional and evolutionary studies of important gene families among Brassica species due to its close relationship to Arabidopsis thaliana.

RESULTS

We identified a total of 145 putative PHD finger proteins containing 233 PHD domains from the current version of B. rapa genome database. Gene ontology analysis showed that 67.7% of them were predicted to be located in nucleus, and 91.3% were predicted to be involved in protein binding activity. Phylogenetic, gene structure, and additional domain analyses clustered them into different groups and subgroups, reflecting their diverse functional roles during plant growth and development. Chromosomal location analysis showed that they were unevenly distributed on the 10 B. rapa chromosomes. Expression analysis from RNA-Seq data showed that 55.7% of them were constitutively expressed in all the tested tissues or organs with relatively higher expression levels reflecting their important housekeeping roles in plant growth and development, while several other members were identified as preferentially expressed in specific tissues or organs. Expression analysis of a subset of 18 B. rapa PHD finger genes under drought and salt stresses showed that all these tested members were responsive to the two abiotic stress treatments.

CONCLUSIONS

Our results reveal that the PHD finger genes play diverse roles in plant growth and development, and can serve as a source of candidate genes for genetic engineering and improvement of Brassica crops against abiotic stresses. This study provides valuable information and lays the foundation for further functional determination of PHD finger genes across the Brassica species.

Genome-Wide Identification and Characterization of the UBP Gene Family in Moso Bamboo (Phyllostachys edulis)

The largest group of deubiquitinases-ubiquitin-specific proteases (UBPs)-perform extensive and significant roles in plants, including the regulation of development and stress responses. A comprehensive analysis of UBP genes has been performed in Arabidopsis thaliana, but no systematic study has been conducted in moso bamboo (Phyllostachys edulis). In this study, the genome-wide identification, classification, gene, protein, promoter region characterization, divergence time, and expression pattern analyses of the UBPs in moso bamboo were conducted. In total, 48 putative UBP genes were identified in moso bamboo, which were divided into 14 distinct subfamilies in accordance with a comparative phylogenetic analysis using 132 full-length protein sequences, including 48, 27, 25, and 32 sequences from moso bamboo, A. thaliana, rice (Oryza sativa), and purple false brome (Brachypodium distachyon), respectively. Analyses of the evolutionary patterns and divergence levels revealed that the PeUBP genes experienced a duplication event approximately 15 million years ago and that the divergence between PeUBP and OsUBP occurred approximately 27 million years ago. Additionally, several PeUBP members were significantly upregulated under abscisic acid, methyl jasmonate, and salicylic acid treatments, indicating their potential roles in abiotic stress responses in plants.

Transcriptome Analysis of Rlm2-Mediated Host Immunity in the Brassica napus-Leptosphaeria maculans Pathosystem

Our study investigated disease resistance in the Brassica napus-Leptosphaeria maculans pathosystem using a combination of laser microdissection, dual RNA sequencing, and physiological validations of large-scale gene sets. The use of laser microdissection improved pathogen detection and identified putative L. maculans effectors and lytic enzymes operative during host colonization. Within 24 h of inoculation, we detected large shifts in gene activity in resistant cotyledons associated with jasmonic acid and calcium signaling pathways that accelerated the plant defense response. Sequencing data were validated through the direct quantification of endogenous jasmonic acid levels. Additionally, resistance against L. maculans was abolished when the calcium chelator EGTA was applied to the inoculation site, providing physiological evidence of the role of calcium in B. napus immunity against L. maculans. We integrated gene expression data with all available information on cis-regulatory elements and transcription factor binding affinities to better understand the gene regulatory networks underpinning plant resistance to hemibiotrophic pathogens. These in silico analyses point to early cellular reprogramming during host immunity that are coordinated by CAMTA, BZIP, and bHLH transcription factors. Together, we provide compelling genetic and physiological evidence into the programming of plant resistance against fungal pathogens.

WRKYs, the Jack-of-various-Trades, Modulate Dehydration Stress in Populus davidiana-A Transcriptomic Approach

Populus davidiana, native to Korea and central Asian countries, is a major contributor to the Korean forest cover. In the current study, using high-throughput RNA-seq mediated transcriptome analysis, we identified about 87 P. davidiana WRKY transcription factors (PopdaWRKY TFs) that showed differential expression to dehydration stress in both sensitive and tolerant cultivars. Our results suggested that, on average, most of the WRKY genes were upregulated in tolerant cultivars but downregulated in sensitive cultivars. Based on protein sequence alignment, P. davidiana WRKYs were classified into three major groups, I, II, III, and further subgroups. Phylogenetic analysis showed that WRKY TFs and their orthologs in Arabidopsis and rice were clustered together in the same subgroups, suggesting similar functions across species. Significant correlation was found among qRT-PCR and RNA-seq analysis. In vivo analysis using model plant Arabidopsis showed that atwrky62 (orthologous to Potri.016G137900) knockout mutants were significantly sensitive to dehydration possibly due to an inability to close their stomata under dehydration conditions. In addition, a concomitant decrease in expression of ABA biosynthetic genes was observed. The AtHK1 that regulates stomatal movement was also downregulated in atwrky62 compared to the wild type. Taken together, our findings suggest a regulatory role of PopdaWRKYs under dehydration stress.

Genome-scale characterization of the vacuole nitrate transporter Chloride Channel (CLC) genes and their transcriptional responses to diverse nutrient stresses in allotetraploid rapeseed

The Chloride Channel (CLC) gene family is reported to be involved in vacuolar nitrate (NO₃^-) transport. Nitrate distribution to the cytoplasm is beneficial for enhancing NO₃^- assimilation and plays an important role in the regulation of nitrogen (N) use efficiency (NUE). In this study, genomic information, high-throughput transcriptional profiles, and gene co-expression analysis were integrated to identify the CLCs (BnaCLCs) in Brassica napus. The decreased NO₃^- concentration in the clca-2 mutant up-regulated the activities of nitrate reductase and glutamine synthetase, contributing to increase N assimilation and higher NUE in Arabidopsis thaliana. The genome-wide identification of 22BnaCLC genes experienced strong purifying selection. Segmental duplication was the major driving force in the expansion of the BnaCLC gene family. The most abundant cis-acting regulatory elements in the gene promoters, including DNA-binding One Zinc Finger, W-box, MYB, and GATA-box, might be involved in the transcriptional regulation of BnaCLCs expression. High-throughput transcriptional profiles and quantitative real-time PCR results showed that BnaCLCs responded differentially to distinct NO₃^- regimes. Transcriptomics-assisted gene co-expression network analysis identified BnaA7.CLCa-3 as the core member of the BnaCLC family, and this gene might play a central role in vacuolar NO₃^- transport in crops. The BnaCLC members also showed distinct expression patterns under phosphate depletion and cadmium toxicity. Taken together, our results provide comprehensive insights into the vacuolar BnaCLCs and establish baseline information for future studies on BnaCLCs-mediated vacuolar NO₃^- storage and its effect on NUE.

Exploring transcription factors reveals crucial members and regulatory networks involved in different abiotic stresses in Brassica napus L

BACKGROUND

Brassica napus (B. napus) encompasses diverse transcription factors (TFs), but thorough identification and characterization of TF families, as well as their transcriptional responsiveness to multifarious stresses are still not clear.

RESULTS

Totally 2167 TFs belonging to five families were genome-widely identified in B. napus, including 518 BnAP2/EREBPs, 252 BnbZIPs, 721 BnMYBs, 398 BnNACs and 278 BnWRKYs, which contained some novel members in comparison with existing results. Sub-genome distributions of BnAP2/EREBPs and BnMYBs indicated that the two families might have suffered from duplication and divergence during evolution. Synteny analysis revealed strong co-linearity between B. napus and its two ancestors, although chromosomal rearrangements have occurred and 85 TFs were lost. About 7.6% and 9.4% TFs of the five families in B. napus were novel genes and conserved genes, which both showed preference on the C sub-genome. RNA-Seq revealed that more than 80% TFs were abiotic stress inducible and 315 crucial differentially expressed genes (DEGs) were screened out. Network analysis revealed that the 315 DEGs are highly co-expressed. The homologous gene network in A. thaliana revealed that a considerable amount of TFs could trigger the differential expression of targeted genes, resulting in a complex clustered network with clusters of genes responsible for targeted stress responsiveness.

CONCLUSIONS

We identified and characterized five TF families in B. napus. Some crucial members and regulatory networks involved in different abiotic stresses have been explored. The investigations deepen our understanding of TFs for stress tolerance in B. napus.

Molecular dynamics simulations revealed structural differences among WRKY domain-DNA interaction in barley (Hordeum vulgare)

Background

The WRKY transcription factors are a class of DNA-binding proteins involved in diverse plant processes play critical roles in response to abiotic and biotic stresses. Genome-wide divergence analysis of WRKY gene family in Hordeum vulgare provided a framework for molecular evolution and functional roles. So far, the crystal structure of WRKY from barley has not been resolved; moreover, knowledge of the three-dimensional structure of WRKY domain is pre-requisites for exploring the protein-DNA recognition mechanisms. Homology modelling based approach was used to generate structures for WRKY DNA binding domain (DBD) and its variants using AtWRKY1 as a template. Finally, the stability and conformational changes of the generated model in unbound and bound form was examined through atomistic molecular dynamics (MD) simulations for 100 ns time period.

Results

In this study, we investigated the comparative binding pattern of WRKY domain and its variants with W-box cis-regulatory element using molecular docking and dynamics (MD) simulations assays. The atomic insight into WRKY domain exhibited significant variation in the intermolecular hydrogen bonding pattern, leading to the structural anomalies in the variant type and differences in the DNA-binding specificities. Based on the MD analysis, residual contribution and interaction contour, wild-type WRKY (HvWRKY46) were found to interact with DNA through highly conserved heptapeptide in the pre- and post-MD simulated complexes, whereas heptapeptide interaction with DNA was missing in variants (I and II) in post-MD complexes. Consequently, through principal component analysis, wild-type WRKY was also found to be more stable by obscuring a reduced conformational space than the variant I (HvWRKY34). Lastly, high binding free energy for wild-type and variant II allowed us to conclude that wild-type WRKY-DNA complex was more stable relative to variants I.

Conclusions

The results of our study revealed complete dynamic and structural information about WRKY domain-DNA interactions. However, no structure base information reported to date for WRKY variants and their mechanism of interaction with DNA. Our findings highlighted the importance of selecting a sequence to generate newer transgenic plants that would be increasingly tolerance to stress conditions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4506-3) contains supplementary material, which is available to authorized users.

Identification, Expression, and Functional Analysis of the Group IId WRKY Subfamily in Upland Cotton (Gossypium hirsutum L.)

WRKY transcription factors have diverse functions in regulating stress response, leaf senescence, and plant growth and development. However, knowledge of the group IId WRKY subfamily in cotton is largely absent. This study identified 34 group IId WRKY genes in the Gossypium hirsutum genome, and their genomic loci were investigated. Members clustered together in the phylogenetic tree had similar motif compositions and gene structural features, revealing similarity and conservation within group IId WRKY genes. During the evolutionary process, 14 duplicated genes appeared to undergo purification selection. Public RNA-seq data were used to examine the expression patterns of group IId WRKY genes in various tissues and under drought and salt stress conditions. Ten highly expressed genes were identified, and the ten candidate genes revealed distinct expression patterns under drought and salt treatments by qRT-PCR analysis. Among them, Gh_A11G1801 was used for functional characterization. GUS activity was differentially induced by various stresses in Gh_A11G1801p::GUS transgenic Arabidopsis plants. The virus-induced gene silencing (VIGS) of Gh_A11G1801 resulted in drought sensitivity in cotton plants, which was accompanied by elevated malondialdehyde (MDA) content and reduced catalase (CAT) content. Taken together, these findings obtained in this study provide valuable resources for further studying group IId WRKY genes in cotton. Our results also enrich the gene resources for the genetic improvements of cotton varieties that are suitable for growth in stressful conditions.

Drought Response in Wheat: Key Genes and Regulatory Mechanisms Controlling Root System Architecture and Transpiration Efficiency

Abiotic stresses such as, drought, heat, salinity, and flooding threaten global food security. Crop genetic improvement with increased resilience to abiotic stresses is a critical component of crop breeding strategies. Wheat is an important cereal crop and a staple food source globally. Enhanced drought tolerance in wheat is critical for sustainable food production and global food security. Recent advances in drought tolerance research have uncovered many key genes and transcription regulators governing morpho-physiological traits. Genes controlling root architecture and stomatal development play an important role in soil moisture extraction and its retention, and therefore have been targets of molecular breeding strategies for improving drought tolerance. In this systematic review, we have summarized evidence of beneficial contributions of root and stomatal traits to plant adaptation to drought stress. Specifically, we discuss a few key genes such as, DRO1 in rice and ERECTA in Arabidopsis and rice that were identified to be the enhancers of drought tolerance via regulation of root traits and transpiration efficiency. Additionally, we highlight several transcription factor families, such as, ERF (ethylene response factors), DREB (dehydration responsive element binding), ZFP (zinc finger proteins), WRKY, and MYB that were identified to be both positive and negative regulators of drought responses in wheat, rice, maize, and/or Arabidopsis. The overall aim of this review is to provide an overview of candidate genes that have been identified as regulators of drought response in plants. The lack of a reference genome sequence for wheat and non-transgenic approaches for manipulation of gene functions in wheat in the past had impeded high-resolution interrogation of functional elements, including genes and QTLs, and their application in cultivar improvement. The recent developments in wheat genomics and reverse genetics, including the availability of a gold-standard reference genome sequence and advent of genome editing technologies, are expected to aid in deciphering of the functional roles of genes and regulatory networks underlying adaptive phenological traits, and utilizing the outcomes of such studies in developing drought tolerant cultivars.

Genome-wide characterization of the WRKY gene family in radish (Raphanus sativus L.) reveals its critical functions under different abiotic stresses

The radish WRKY gene family was genome-widely identified and played critical roles in response to multiple abiotic stresses. The WRKY is among the largest transcription factors (TFs) associated with multiple biological activities for plant survival, including control response mechanisms against abiotic stresses such as heat, salinity, and heavy metals. Radish is an important root vegetable crop and therefore characterization and expression pattern investigation of WRKY transcription factors in radish is imperative. In the present study, 126 putative WRKY genes were retrieved from radish genome database. Protein sequence and annotation scrutiny confirmed that RsWRKY proteins possessed highly conserved domains and zinc finger motif. Based on phylogenetic analysis results, RsWRKYs candidate genes were divided into three groups (Group I, II and III) with the number 31, 74, and 20, respectively. Additionally, gene structure analysis revealed that intron-exon patterns of the WRKY genes are highly conserved in radish. Linkage map analysis indicated that RsWRKY genes were distributed with varying densities over nine linkage groups. Further, RT-qPCR analysis illustrated the significant variation of 36 RsWRKY genes under one or more abiotic stress treatments, implicating that they might be stress-responsive genes. In total, 126 WRKY TFs were identified from the R. sativus genome wherein, 35 of them showed abiotic stress-induced expression patterns. These results provide a genome-wide characterization of RsWRKY TFs and baseline for further functional dissection and molecular evolution investigation, specifically for improving abiotic stress resistances with an ultimate goal of increasing yield and quality of radish.