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

Heavy metal stress and mitogen activated kinase transcription factors in plants: Exploring heavy metal-ROS influences on plant signalling pathways

Due to their stationary nature, plants are exposed to a diverse range of biotic and abiotic stresses, of which heavy metal (HM) stress poses one of the most detrimental abiotic stresses, targeting diverse plant processes. HMs instigate the overproduction of reactive oxygen species (ROS), and to mitigate the adverse effects of ROS, plants induce multiple defence mechanisms. Besides the negative implications of overproduction of ROS, these molecules play a multitude of signalling roles in plants, acting as a central player in the complex signalling network of cells. One of the ROS-associated signalling mechanisms is the mitogen-activated protein kinase (MAPK) cascade, a signalling pathway which transduces extracellular stimuli into intracellular responses. Plant MAPKs have been implicated in signalling involved in stress response, phytohormone regulation, and cell cycle cues. However, the influence of various HMs on MAPK activation has not been well documented. In this review, we address and summarise several aspects related to various HM-induced ROS signalling. Additionally, we touch on how these signals activate the MAPK cascade and the downstream transcription factors that influence plant responses to HMs. Moreover, we propose a workflow that could characterise genes associated with MAPKs and their roles during plant HM stress responses.

Full-Length Transcriptome Sequencing Analysis and Characterization of WRKY Transcription Factors Responsive to Cadmium Stress in Arabis paniculata

Arabis paniculata is a newly discovered hyperaccumulator known for its ability to accumulate multiple metals. WRKY proteins play a significant role in plant responses to various stresses, including cadmium (Cd) stress. However, there is limited research on the molecular biology of Arabis paniculata, especially regarding the WRKY family. In this study, we conducted third-generation sequencing for functional annotation and structural analysis of Arabis paniculata. We obtained 41,196 high-quality isoforms from the full-length transcriptome, with an average length of 1043 bp. A total of 26,670 genes were predicted against NR, Swissprot, KOG, and KEGG databases. Functional comparison using the KOG database revealed excellent annotation in 25 functional categories, with general function prediction (1822 items) being the most predominant. MISA analysis identified 12,593 SSR loci, with single nucleotide repeats being the largest category (44.83% of the total). Moreover, our predictions provide insights into 20,022 coding sequences (CDS), 811 transcription factors, and 17,963 LncRNAs. In total, 34 WRKY gene sequences were identified in Arabis paniculata. Bioinformatics analysis revealed diverse numbers of amino acids in these WRKYs (113 to 545 aa), and a conserved WRKYGQK sequence within the N-terminus of the WRKY protein. Furthermore, all WRKYs were found to be localized in the nucleus. Phylogenetic analysis classified the WRKY genes into three categories: I (14 members), II (17 members), and III (3 members). Category II was subsequently divided into four sub-categories: II-a (8 members), II-b (1 member), II-c (1 member), and II-d (7 members). Our quantitative real-time polymerase chain reaction (qRT-PCR) experiments revealed that ApWRKY23 and ApWRKY34 exhibited the highest expression levels at the 24-h time point, suggesting their potential role as the candidate genes for Cd stress response. These findings contribute to our understanding of the genomic information of Arabis paniculata and provide a basis for the analysis of its genetic diversity. Additionally, this study paves the way for a comprehensive exploration of the molecular mechanisms underlying the WRKY genes in Arabis paniculata under Cd stress conditions.

The transcription factor BnaA9.WRKY47 coordinates leaf senescence and nitrogen remobilization in Brassica napus

Nitrogen (N) is an essential macronutrient for plants, and its remobilization is key for adaptation to deficiency stress. However, there is limited understanding of the regulatory mechanisms of N remobilization in the important crop species Brassica napus (oilseed rape). Here, we report the identification of a transcription factor, BnaA9.WRKY47, that is induced by N starvation in a canola variety. At the seedling stage, BnaA9.WRKY47-overexpressing (OE) lines displayed earlier senescence of older leaves and preferential growth of juvenile leaves compared to the wild type under N starvation. At the field scale, the seed yield was significantly increased in the BnaA9.WRKY47-OE lines compared with the wild type when grown under N deficiency conditions and, conversely, it was reduced in BnaA9.WRKY47-knockout mutants. Biochemical analyses demonstrated that BnaA9.WRKY47 directly activates BnaC7.SGR1 to accelerate senescence of older leaves. In line with leaf senescence, the concentration of amino acids in the older leaves of the OE lines was elevated, and the proportion of plant N that they contained was reduced. This was associated with BnaA9.WRKY47 activating the amino acid permease BnaA9.AAP1 and the nitrate transporter BnaA2.NRT1.7. Thus, the expression of BnaA9.WRKY47 efficiently facilitated N remobilization from older to younger leaves or to seeds. Taken together, our results demonstrate that BnaA9.WRKY47 up-regulates the expression of BnaC7.SGR1, BnaA2.NRT1.7, and BnaA9AAP1, thus promoting the remobilization of N in B. napus under starvation conditions.

Identification and analysis of isoflavone reductase gene family in Gossypium hirsutum L

Isoflavone reductase (IFR) is a key enzyme controlling isoflavone synthesis and widely involved in response to various stresses. In this study, the IFR genes in four Gossypium species and other 7 species were identified and analyzed in the whole genome, and the physicochemical properties, gene structures, cis-acting elements, chromosomal locations, collinearity relationships and expression patterns of IFR genes were systematically analyzed. 28, 28, 14 and 15 IFR genes were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii, respectively, which were divided into five clades according to the evolutionary tree and gene structure. Collinear analysis showed that segmental duplication and whole genome duplication were the main driving forces in the process of evolution, and most genes underwent pure selection. Gene structure analysis showed that IFR gene family was relatively conserved. Cis-element analysis of promoter showed that most GhIFR genes contain cis-elements related to abiotic stresses and plant hormones. Analysis of GhIFR gene expression under different stresses showed that GhIFR genes were involved in the response to drought, salt, heat and cold stresses through corresponding network mechanisms, especially GhIFR9A. Phenotypic analysis after silencing GhIFR9A gene by VIGS was shown that GhIFR9A gene was involved in the response to salt stress. This study laid a foundation for the subsequent functional study of cotton IFR genes.

Genome-wide investigation of the WRKY transcription factor gene family in weeping forsythia: expression profile and cold and drought stress responses

Weeping forsythia is a wide-spread shrub in China with important ornamental, medicinal and ecological values. It is widely distributed in China's warm temperate zone. In plants, WRKY transcription factors play important regulatory roles in seed germination, flower development, fruit ripening and coloring, and biotic and abiotic stress response. To date, WRKY transcription factors have not been systematically studied in weeping forsythia. In this study, we identified 79 WRKY genes in weeping forsythia and classified them according to their naming rules in Arabidopsis thaliana. Phylogenetic tree analysis showed that, except for IIe subfamily, whose clustering was inconsistent with A. thaliana clustering, other subfamily clustering groups were consistent. Cis-element analysis showed that WRKY genes related to pathogen resistance in weeping forsythia might be related to methyl jasmonate and salicylic acid-mediated signaling pathways. Combining cis-element and expression pattern analyses of WRKY genes showed that more than half of WRKY genes were involved in light-dependent development and morphogenesis in different tissues. The gene expression results showed that 13 WRKY genes were involved in drought response, most of which might be related to the abscisic acid signaling pathway, and a few of which might be regulated by MYB transcription factors. The gene expression results under cold stress showed that 17 WRKY genes were involved in low temperature response, and 9 of them had low temperature responsiveness cis-elements. Our study of WRKY family in weeping forsythia provided useful resources for molecular breeding and important clues for their functional verification.

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.

Genome-Wide Identification of the Hypericum perforatum WRKY Gene Family Implicates HpWRKY85 in Drought Resistance

WRKY, named for its special heptapeptide conserved sequence WRKYGOK, is one of the largest transcription factor families in plants and is widely involved in plant responses to biotic, abiotic, and hormonal stresses, especially the important regulatory function in response to drought stress. However, there is no complete comprehensive analysis of this family in H. perforatum, which is one of the most extensively studied plants and is probably the best-known herbal medicine on the market today, serving as an antidepressant, neuroprotective, an antineuralgic, and an antiviral. Here, we identified 86 HpWRKY genes according to the whole genome database of H. perforatum, and classified them into three groups through phylogenetic analysis. Gene structure, conserved domain, motif, cis-elements, gene ontology, and expression profiling were performed. Furthermore, it was found that HpWRKY85, a homologous gene of AtWRKY75, showed obvious responses to drought treatment. Subcellular localization analysis indicated that this protein was localized in the nucleus by the Arabidopsis protoplasts transient transfection. Meanwhile, HpWRKY85-overexpressing Arabidopsis plants showed a stronger ability of root growth and scavenging endogenous reactive oxygen species. The results provide a reference for further understanding the role of HpWRKY85 in the molecular mechanism of drought resistance of H. perforatum.

Selection and validation of reference genes for RT-qPCR analysis of different organs at various development stages in Caragana intermedia

Reverse transcription quantitative PCR (RT-qPCR) is a technique widely used to investigate the expression of genes. An appropriate reference gene (RG) is essential for RT-qPCR analysis to obtain accurate and reliable results. Caragana intermedia plays an important role in afforestation as a bush. However, due to the lack of appropriate RGs, the research on development-related genes is limited. In this study, the selection for suitable RGs of different organs at various development stages to normalize the results of RT-qPCR about development-related genes was performed. To test the expression stability across all samples, we used the software algorithms such as geNorm, NormFinder, BestKeeper, and RefFinder to evaluate all the candidate RGs. Our results showed that CiEF1α was the most stable RG with little fluctuation among all samples. In addition, CiGAPDH in roots, CiSKIP1 in stems and leaves, and CiEF1α in different organs were selected as the most stable RGs. To confirm the applicability of the most stable RGs, the relative expression of CiWRKY17 was normalized using different candidate RGs. Taken together, our research laid a foundation for the study of development-related genes in C. intermedia.

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.

MAPK Cascades and Transcriptional Factors: Regulation of Heavy Metal Tolerance in Plants

In nature, heavy metal (HM) stress is one of the most destructive abiotic stresses for plants. Heavy metals produce toxicity by targeting key molecules and important processes in plant cells. The mitogen-activated protein kinase (MAPK) cascade transfers the signals perceived by cell membrane surface receptors to cells through phosphorylation and dephosphorylation and targets various effector proteins or transcriptional factors so as to result in the stress response. Signal molecules such as plant hormones, reactive oxygen species (ROS), and nitric oxide (NO) can activate the MAPK cascade through differentially expressed genes, the activation of the antioxidant system and synergistic crosstalk between different signal molecules in order to regulate plant responses to HMs. Transcriptional factors, located downstream of MAPK, are key factors in regulating plant responses to heavy metals and improving plant heavy metal tolerance and accumulation. Thus, understanding how HMs activate the expression of the genes related to the MAPK cascade pathway and then phosphorylate those transcriptional factors may allow us to develop a regulation network to increase our knowledge of HMs tolerance and accumulation. This review highlighted MAPK pathway activation and responses under HMs and mainly focused on the specificity of MAPK activation mediated by ROS, NO and plant hormones. Here, we also described the signaling pathways and their interactions under heavy metal stresses. Moreover, the process of MAPK phosphorylation and the response of downstream transcriptional factors exhibited the importance of regulating targets. It was conducive to analyzing the molecular mechanisms underlying heavy metal accumulation and tolerance.

Approaches Involved in the Vegetable Crops Salt Stress Tolerance Improvement: Present Status and Way Ahead

Salt stress is one of the most important abiotic stresses as it persists throughout the plant life cycle. The productivity of crops is prominently affected by soil salinization due to faulty agricultural practices, increasing human activities, and natural processes. Approximately 10% of the total land area (950 Mha) and 50% of the total irrigated area (230 Mha) in the world are under salt stress. As a consequence, an annual loss of 12 billion US$ is estimated because of reduction in agriculture production inflicted by salt stress. The severity of salt stress will increase in the upcoming years with the increasing world population, and hence the forced use of poor-quality soil and irrigation water. Unfortunately, majority of the vegetable crops, such as bean, carrot, celery, eggplant, lettuce, muskmelon, okra, pea, pepper, potato, spinach, and tomato, have very low salinity threshold (ECt, which ranged from 1 to 2.5 dS m-1 in saturated soil). These crops used almost every part of the world and lakes' novel salt tolerance gene within their gene pool. Salt stress severely affects the yield and quality of these crops. To resolve this issue, novel genes governing salt tolerance under extreme salt stress were identified and transferred to the vegetable crops. The vegetable improvement for salt tolerance will require not only the yield influencing trait but also target those characters or traits that directly influence the salt stress to the crop developmental stage. Genetic engineering and grafting is the potential tool which can improve salt tolerance in vegetable crop regardless of species barriers. In the present review, an updated detail of the various physio-biochemical and molecular aspects involved in salt stress have been explored.

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.

Genome-Wide Analysis of the Auxin/Indoleacetic Acid Gene Family and Response to Indole-3-Acetic Acid Stress in Tartary Buckwheat (Fagopyrum tataricum)

Auxin/indoleacetic acid (Aux/IAA) family genes respond to the hormone auxin, which have been implicated in the regulation of multiple biological processes. In this study, all 25 Aux/IAA family genes were identified in Tartary buckwheat (Fagopyrum tataricum) by a reiterative database search and manual annotation. Our study provided comprehensive information of Aux/IAA family genes in buckwheat, including gene structures, chromosome locations, phylogenetic relationships, and expression patterns. Aux/IAA family genes were nonuniformly distributed in the buckwheat chromosomes and divided into seven groups by phylogenetic analysis. Aux/IAA family genes maintained a certain correlation and a certain species-specificity through evolutionary analysis with Arabidopsis and other grain crops. In addition, all Aux/IAA genes showed a complex response pattern under treatment of indole-3-acetic acid (IAA). These results provide valuable reference information for dissecting function and molecular mechanism of Aux/IAA family genes in buckwheat.

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.

BASIC PENTACYSTEINE2 negatively regulates osmotic stress tolerance by modulating LEA4-5 expression in Arabidopsis thaliana

Osmotic stress substantially affects plant growth and development. Study of plant transcription factors involved in osmotic stress can enhance our understanding of the mechanisms of plant osmotic stress tolerance and how the tolerance of plants to osmotic stress can be improved. Here, we identified the specific function of Arabidopsis thaliana BARLEY B RECOMBINANT/BASIC PENTACYSTEINE transcription factor, BPC2, in the osmotic stress response. Phenotypic analysis showed that loss-of-function of BPC2 led to an increase in osmotic stress tolerance in the seedling growth stage. Physiological analysis showed that mutation of BPC2 in Arabidopsis alleviated osmotic-induced increases in H₂O₂ accumulation, the malondialdehyde (MDA) content, and percent electrolyte leakage. BPC2 was localized in the nucleus. RNA-seq and qRT-PCR analysis showed that BPC2 could negatively regulate the expression of late embryogenesis abundant (LEA) genes (LEA3, LEA4-2, and LEA4-5). Further analysis showed that BPC2 could directly bind to the promoter of LEA4-5 in vitro and in vivo. Overexpression of BPC2 enhanced hypersensitivity to osmotic stress in the seedling growth stage. Overexpression of BPC2 led to decreases in LEA4-5 expression and aggravated osmotic-induced increases in H₂O₂ accumulation, the MDA content, and percent electrolyte leakage. Overall, our results indicate that BPC2 negatively regulates LEA4-5 expression to modulate osmotic-induced H₂O₂ accumulation, the MDA content, and percent electrolyte leakage, all of which affect the osmotic stress response in Arabidopsis thaliana.

Genome-wide identification and functional prediction of salt- stress related long non-coding RNAs (lncRNAs) in chickpea (Cicer arietinum L.)

LncRNAs (long noncoding RNAs) are 200 bp length crucial RNA molecules, lacking coding potential and having important roles in regulating gene expression, particularly in response to abiotic stresses. In this study, we identified salt stress-induced lncRNAs in chickpea roots and predicted their intricate regulatory roles. A total of 3452 novel lncRNAs were identified to be distributed across all 08 chickpea chromosomes. On comparing salt-tolerant (ICCV 10, JG 11) and salt-sensitive cultivars (DCP 92-3, Pusa 256), 4446 differentially expressed lncRNAs were detected under various salt  treatments. We predicted 3373 lncRNAs to be regulating their target genes in cis regulating manner and 80 unique lncRNAs were observed as interacting with 136 different miRNAs, as eTMs (endogenous target mimic) targets of miRNAs and implicated them in the regulatory network of salt stress response. Functional analysis of these lncRNA revealed their association in targeting salt stress response-related genes like potassium transporter, transporter family genes, serine/threonine-protein kinase, aquaporins like TIP1-2, PIP2-5 and transcription factors like, AP2, NAC, bZIP, ERF, MYB and WRKY. Furthermore, about 614 lncRNA-SSRs (simple sequence repeats) were identified as a new generation of molecular markers with higher efficiency and specificity in chickpea. Overall, these findings will pave the understanding of comprehensive functional role of potential lncRNAs, which can help in providing insight into the molecular mechanism of salt tolerance in chickpea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01093-0.

Genome-Wide Identification and Transcriptional Expression Profiles of Transcription Factor WRKY in Common Walnut (Juglans regia L.)

The transcription factor WRKY is widely distributed in the plant kingdom, playing a significant role in plant growth, development and response to stresses. Walnut is an economically important temperate tree species valued for both its edible nuts and high-quality wood, and its response to various stresses is an important factor that determines the quality of its fruit. However, in walnut trees themselves, information about the WRKY gene family remains scarce. In this paper, we perform a comprehensive study of the WRKY gene family in walnut. In total, we identified 103 WRKY genes in the common walnut that are clustered into 4 groups and distributed on 14 chromosomes. The conserved domains all contained a WRKY domain, and motif 2 was observed in most WRKYs, suggesting a high degree of conservation and similar functions within each subfamily. However, gene structure was significantly differentiated between different subfamilies. Synteny analysis indicates that there were 56 gene pairs in J. regia and A. thaliana, 76 in J. regia and J. mandshurica, 75 in J. regia and J. microcarpa, 76 in J. regia and P. trichocarpa, and 33 in J. regia and Q. robur, indicating that the WRKY gene family may come from a common ancestor. GO and KEGG enrichment analysis showed that the WRKY gene family was involved in resistance traits and the plant-pathogen interaction pathway. In anthracnose-resistant F26 fruits (AR) and anthracnose-susceptible F423 fruits (AS), transcriptome and qPCR analysis results showed that JrWRKY83, JrWRKY73 and JrWRKY74 were expressed significantly more highly in resistant cultivars, indicating that these three genes may be important contributors to stress resistance in walnut trees. Furthermore, we investigate how these three genes potentially target miRNAs and interact with proteins. JrWRKY73 was target by the miR156 family, including 12 miRNAs; this miRNA family targets WRKY genes to enhance plant defense. JrWRKY73 also interacted with the resistance gene AtMPK6, showing that it may play a crucial role in walnut defense.

WRKY transcription factors and plant defense responses: latest discoveries and future prospects

WRKY transcription factors are among the largest families of transcriptional regulators. In this review, their pivotal role in modulating various signal transduction pathways during biotic and abiotic stresses is discussed. Transcription factors (TFs) are important constituents of plant signaling pathways that define plant responses against biotic and abiotic stimuli besides playing a role in response to internal signals which coordinate different interacting partners during developmental processes. WRKY TFs, deriving their nomenclature from their signature DNA-binding sequence, represent one of the largest families of transcriptional regulators found exclusively in plants. By modulating different signal transduction pathways, these TFs contribute to various plant processes including nutrient deprivation, embryogenesis, seed and trichome development, senescence as well as other developmental and hormone-regulated processes. A growing body of research suggests transcriptional regulation of WRKY TFs in adapting plant to a variety of stressed environments. WRKY TFs can regulate diverse biological functions from receptors for pathogen triggered immunity, modulator of chromatin for specific interaction and signal transfer through a complicated network of genes. Latest discoveries illustrate the interaction of WRKY proteins with other TFs to form an integral part of signaling webs that regulate several seemingly disparate processes and defense-related genes, thus establishing their significant contributions to plant immune response. The present review starts with a brief description on the structural characteristics of WRKY TFs followed by the sections that present recent evidence on their roles in diverse biological processes in plants. We provide a comprehensive overview on regulatory crosstalks involving WRKY TFs during multiple stress responses in plants and future prospects of WRKY TFs as promising molecular diagnostics for enhancing crop improvement.

Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants

Plants produce thousands of chemically diverse secondary metabolites, many of which have valuable pharmaceutical properties. There is much interest in the synthesis of these pharmaceuticallyvaluable compounds, including the key enzymes and the transcription factors involved. The function and regulatory mechanism of transcription factors in biotic and abiotic stresses have been studied in depth. However, their regulatory roles in the biosynthesis of bioactive compounds, especially in medicinal plants, have only begun. Here, we review what is currently known about how transcription factors contribute to the synthesis of bioactive compounds (alkaloids, terpenoids, flavonoids, and phenolic acids) in medicinal plants. Recent progress has been made in the cloning and characterization of transcription factors in medicinal plants on the genome scale. So far, several large transcription factors have been identified in MYB, WRKY, bHLH, ZIP, AP2/ERF transcription factors. These transcription factors have been predicted to regulate bioactive compound production. These transcription factors positively or negatively regulate the expression of multiple genes encoding key enzymes, and thereby control the metabolic flow through the biosynthetic pathway. Although the research addressing this niche topic is in its infancy, significant progress has been made, and advances in high-throughput sequencing technology are expected to accelerate the discovery of key regulatory transcription factors in medicinal plants. This review is likely to be useful for those interested in the synthesis of pharmaceutically- valuable plant compounds, especially those aiming to breed or engineer plants that produce greater yields of these compounds.

Transcriptome-Wide Identification of WRKY Transcription Factors and Their Expression Profiles under Different Types of Biological and Abiotic Stress in Pinus massoniana Lamb

Pinus massoniana Lamb, an economically important conifer tree, is widely distributed in China. WRKY transcription factors (TFs) play important roles in plant growth and development, biological and abiotic stress. Nevertheless, there is little information about the WRKY genes in P. massoniana. By searching for conserved WRKY motifs in transcriptomic RNA sequencing data for P. massoniana, 31 sequences were identified as WRKY TFs. Then, phylogenetic and conserved motif analyses of the WRKY family in P. massoniana, Pinus taeda and Arabidopsis thaliana were used to classify WRKY genes. The expression patterns of six PmWRKY genes from different groups were determined using real-time quantitative PCR for 2-year-old P. massoniana seedings grown in their natural environment and challenged by phytohormones (salicylic acid, methyl jasmonate, or ethephon), abiotic stress (H₂O₂) and mechanical damage stress. As a result, the 31 PmWRKY genes identified were divided into three major groups and several subgroups based on structural and phylogenetic features. PmWRKY genes are regulated in response to abiotic stress and phytohormone treatment and may participate in signaling to improve plant stress resistance. Some PmWRKY genes behaved as predicted based on their homology with A. thaliana WRKY genes, but others showed divergent behavior. This systematic analysis lays the foundation for further identification of WRKY gene functions to aid further exploration of the functions and regulatory mechanisms of PmWRKY genes in biological and abiotic stress in P. massoniana.

Phylogenetic and Transcription Analysis of Hibiscus hamabo Sieb. et Zucc. WRKY Transcription Factors

WRKY transcription factors are known to play important roles in the regulation of various aspects of plant growth and development, including germination, stress resistance, and senescence. Nevertheless, there is little information about the WRKY genes in Hibiscus hamabo Sieb. et Zucc., an important semimangrove plant. In this study, HhWRKY genes in H. hamabo were identificated based on Illumina RNA-sequencing and isoform sequencing from salt-treated roots. Then phylogenetic analysis and conserved motif analysis of the WRKY family in H. hamabo and Arabidopsis thaliana were used to classify WRKY genes. Sixteen HhWRKY genes were selected from different groups to detect their expression patterns using real-time quantitative PCR in different organ (root, old leaf, tender leaf, receptacle, petal, or stamen) from 10-year-old H. hamabo plants grown in their natural environment and in seedlings with 8 to 10 true leaves challenged by phytohormone (salicylic acid, methyl jasmonate, or abscisic acid) and abiotic stress (drought, salt, or high temperature). As a result, the identified 78 HhWRKY genes were divided into two major groups and several subgroups based on their structural and phylogenetic features. Most transcripts of the selected 16 HhWRKY genes were more abundant in old than in tender leaves of H. hamabo. HhWRKY genes were regulated in reaction to abiotic stresses and phytohormone treatments and may participate in signaling networks to improve plant stress resistance. Some of HhWRKY genes behaved as would be predicted based on their homology with A. thaliana WRKY genes, but others showed divergent behavior. This systematic analysis lays the foundation for further identification of WRKY gene functions, with the aim of improving woody plants.

WRKY genes family study reveals tissue-specific and stress-responsive TFs in wild potato species

Wild potatoes, as dynamic resource adapted to various environmental conditions, represent a powerful and informative reservoir of genes useful for breeding efforts. WRKY transcription factors (TFs) are encoded by one of the largest families in plants and are involved in several biological processes such as growth and development, signal transduction, and plant defence against stress. In this study, 79 and 84 genes encoding putative WRKY TFs have been identified in two wild potato relatives, Solanum commersonii and S. chacoense. Phylogenetic analysis of WRKY proteins divided ScWRKYs and SchWRKYs into three Groups and seven subGroups. Structural and phylogenetic comparative analyses suggested an interspecific variability of WRKYs. Analysis of gene expression profiles in different tissues and under various stresses allowed to select ScWRKY045 as a good candidate in wounding-response, ScWRKY055 as a bacterial infection triggered WRKY and ScWRKY023 as a multiple stress-responsive WRKY gene. Those WRKYs were further studied through interactome analysis allowing the identification of potential co-expression relationships between ScWRKYs/SchWRKYs and genes of various pathways. Overall, this study enabled the discrimination of WRKY genes that could be considered as potential candidates in both breeding programs and functional studies.

A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

Genome-wide identification, characterisation and functional evaluation of WRKY genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.) G. Don. under abiotic stresses

BACKGROUND

WRKY DNA-binding protein (WRKY) is a large gene family involved in plant responses and adaptation to salt, drought, cold and heat stresses. Sweet potato from the genus Ipomoea is a staple food crop, but the WRKY genes in Ipomoea species remain unknown to date. Hence, we carried out a genome-wide analysis of WRKYs in Ipomoea trifida (H.B.K.) G. Don., the wild ancestor of sweet potato.

RESULTS

A total of 83 WRKY genes encoding 96 proteins were identified in I. trifida, and their gene distribution, duplication, structure, phylogeny and expression patterns were studied. ItfWRKYs were distributed on 15 chromosomes of I. trifida. Gene duplication analysis showed that segmental duplication played an important role in the WRKY gene family expansion in I. trifida. Gene structure analysis showed that the intron-exon model of the ItfWRKY gene was highly conserved. Meanwhile, the ItfWRKYs were divided into five groups (I, IIa + IIb, IIc, IId + IIe and III) on the basis of the phylogenetic analysis on I. trifida and Arabidopsis thaliana WRKY proteins. In addition, gene expression profiles confirmed by quantitative polymerase chain reaction showed that ItfWRKYs were highly up-regulated or down-regulated under salt, drought, cold and heat stress conditions, implying that these genes play important roles in response and adaptation to abiotic stresses.

CONCLUSIONS

In summary, genome-wide identification, gene structure, phylogeny and expression analysis of WRKY gene in I. trifida provide basic information for further functional studies of ItfWRKYs and for the molecular breeding of sweet potato.

Genome-Wide Identification of WRKY Transcription Factors in the Asteranae

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

Genome-Wide Identification of WRKY Genes and Their Response to Cold Stress in Coffea canephora

WRKY transcription factors are known to play roles in diverse stress responses in plants. Low temperatures limit the geographic distribution of Coffea canephora Pierre ex A.Froehner. The WRKYs of C. canephora are still not well characterized, and the response of C. canephora WRKYs (CcWRKYs) under cold stress is still largely unknown. We identified 49 CcWRKYs from the C. canephora genome to gain insight into these mechanisms. These CcWRKYs were divided into three groups that were based on the conserved WRKY domains and zinc-finger structure. Gene expression analysis demonstrated that 14 CcWRKYs were induced during the cold acclimation stage, 17 CcWRKYs were preferentially upregulated by 4 °C treatment, and 12 CcWRKYs were downregulated by cold stress. Subsequently, we carried out a genome-wide analysis to predict 14,513 potential CcWRKY target genes in C. canephora. These isolated genes were involved in multiple biological processes, and most of them could be grouped by the response to stimulus. Among the putative CcWRKY target genes, 235 genes were categorized into response to the cold process, including carbohydrate metabolic, lipid metabolic, and photosynthesis process-related genes. Furthermore, the qRT-PCR and correlation analysis indicated that CcWRKY might control their putative targets that respond to cold stress. These results provide a basis for understanding the molecular mechanism for CcWRKY-mediated cold responses.

Genome-wide Analysis of the WRKY Gene Family and its Response to Abiotic Stress in Buckwheat (Fagopyrum Tataricum)

The WRKY gene family is an ancient plant transcription factor (TF) family with a vital role in plant growth and development, especially in response to biotic and abiotic stresses. Although many researchers have studied WRKY TFs in numerous plant species, little is known of them in Tartary buckwheat (Fagopyrum tataricum). Based on the recently reported genome sequence of Tartary buckwheat, we identified 78 FtWRKY proteins that could be classified into three major groups. All 77 WRKY genes were distributed unevenly across all eight chromosomes. Exon-intron analysis and motif composition prediction revealed the complexity and diversity of FtWRKYs, indicating that WRKY TFs may be of significance in plant growth regulation and stress response. Two separate pairs of tandem duplication genes were found, but no segmental duplications were identified. Overall, most orthologous gene-pairs between Tartary and common buckwheat evolved under strong purifying selection. qRT-PCR was used to analyze differences in expression among four FtWRKYs (FtWRKY6, 74, 31, and 7) under salt, drought, cold, and heat treatments. The results revealed that all four proteins are related to abiotic stress responses, although they exhibited various expression patterns. In particular, the relative expression levels of FtWRKY6, 74, and 31 were significantly upregulated under salt stress, while the highest expression of FtWRKY7 was observed from heat treatment. This study provides comprehensive insights into the WRKY gene family in Tartary buckwheat, and can support the screening of additional candidate genes for further functional characterization of WRKYs under various stresses.

CaWRKY27 Negatively Regulates H2O2-Mediated Thermotolerance in Pepper (Capsicum annuum)

Heat stress, an important and damaging abiotic stress, regulates numerous WRKY transcription factors, but their roles in heat stress responses remain largely unexplored. Here, we show that pepper (Capsicum annuum) CaWRKY27 negatively regulates basal thermotolerance mediated by H2O2 signaling. CaWRKY27 expression increased during heat stress and persisted during recovery. CaWRKY27 overexpression impaired basal thermotolerance in tobacco (Nicotiana tabacum) and Arabidopsis thaliana, CaWRKY27-overexpressing plants had a lower survival rate under heat stress, accompanied by decreased expression of multiple thermotolerance-associated genes. Accordingly, silencing of CaWRKY27 increased basal thermotolerance in pepper plants. Exogenously applied H2O2 induced CaWRKY27 expression, and CaWRKY27 overexpression repressed the scavenging of H2O2 in Arabidopsis, indicating a positive feedback loop between H2O2 accumulation and CaWRKY27 expression. Consistent with this, CaWRKY27 expression was repressed under heat stress in the presence H2O2 scavengers and CaWRKY27 silencing decreased H2O2 accumulation in pepper leaves. These changes may result from changes in levels of reactive oxygen species (ROS)-scavenging enzymes, since the heat stress-challenged CaWRKY27-silenced pepper plants had significantly higher expression of multiple genes encoding ROS-scavenging enzymes, such as CaCAT1, CaAPX1, CaAPX2, CaCSD2, and CaSOD1. Therefore, CaWRKY27 acts as a downstream negative regulator of H2O2-mediated heat stress responses, preventing inappropriate responses during heat stress and recovery.

Genome-wide identification and characterization of CONSTANS-like gene family in radish (Raphanus sativus)

Floral induction that initiates bolting and flowering is crucial for reproductive fitness in radishes. CONSTANS-like (CO-like, COL) genes play an important role in the circadian clock, which ensures regular development through complicated time-keeping mechanisms. However, the specific biological and functional roles of each COL transcription factor gene in the radish remain unknown. In this study, we performed a genome-wide identification of COL genes in the radish genome of three cultivars including ‘Aokubi’, ‘kazusa’ and ‘WK10039’, and we analyzed their exon-intron structure, gene phylogeny and synteny, and expression levels in different tissues. The bioinformatics analysis identified 20 COL transcription factors in the radish genome, which were divided into three subgroups (Group I to Group III). RsaCOL-09 and RsaCOL-12 might be tandem duplicated genes, whereas the others may have resulted from segmental duplication. The Ka/Ks ratio indicated that all the COL genes in radish, Arabidopsis, Brassica rapa, Brassica oleracea, Capsella rubella and rice were under purifying selection. We identified 6 orthologous and 19 co-orthologous COL gene pairs between the radish and Arabidopsis, and we constructed an interaction network among these gene pairs. The expression values for each COL gene during vegetable and flower development showed that the majority of Group I members had similar expression patterns. In general, the expression of radish COL genes in Groups I and III decreased during development, whereas the expression of radish COL genes in Group II first increased and then decreased. Substantial numbers of radish COL genes were differentially expressed after vernalization treatment. The expression levels of RsaCOL-02 and RsaCOL-04 were significantly increased during vernalization treatment, while the expression of RsaCOL-10 was significantly decreased. These outcomes provide insights for improving the genetic control of bolting and flowering in radish and other root vegetable crops, and they facilitate genetic improvements to radish yields and quality.

Identification and expression analyses of WRKY genes reveal their involvement in growth and abiotic stress response in watermelon (Citrullus lanatus)

Despite identification of WRKY family genes in numerous plant species, a little is known about WRKY genes in watermelon, one of the most economically important fruit crops around the world. Here, we identified a total of 63 putative WRKY genes in watermelon and classified them into three major groups (I-III) and five subgroups (IIa-IIe) in group II. The structure analysis indicated that ClWRKYs with different WRKY domains or motifs may play different roles by regulating respective target genes. The expressions of ClWRKYs in different tissues indicate that they are involved in various tissue growth and development. Furthermore, the diverse responses of ClWRKYs to drought, salt, or cold stress suggest that they positively or negatively affect plant tolerance to various abiotic stresses. In addition, the altered expression patterns of ClWRKYs in response to phytohormones such as, ABA, SA, MeJA, and ETH, imply the occurrence of complex cross-talks between ClWRKYs and plant hormone signals in regulating plant physiological and biological processes. Taken together, our findings provide valuable clues to further explore the function and regulatory mechanisms of ClWRKY genes in watermelon growth, development, and adaption to environmental stresses.

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.

Genome-wide characterization and expression profiling of NAC transcription factor genes under abiotic stresses in radish (Raphanus sativus L.)

NAC (NAM, no apical meristem; ATAF, Arabidopsis transcription activation factor and CUC, cup-shaped cotyledon) proteins are among the largest transcription factor (TF) families playing fundamental biological processes, including cell expansion and differentiation, and hormone signaling in response to biotic and abiotic stresses. In this study, 172 RsNACs comprising 17 membrane-bound members were identified from the whole radish genome. In total, 98 RsNAC genes were non-uniformly distributed across the nine radish chromosomes. In silico analysis revealed that expression patterns of several NAC genes were tissue-specific such as a preferential expression in roots and leaves. In addition, 21 representative NAC genes were selected to investigate their responses to heavy metals (HMs), salt, heat, drought and abscisic acid (ABA) stresses using real-time polymerase chain reaction (RT-qPCR). As a result, differential expressions among these genes were identified where RsNAC023 and RsNAC080 genes responded positively to all stresses except ABA, while RsNAC145 responded more actively to salt, heat and drought stresses compared with other genes. The results provides more valuable information and robust candidate genes for future functional analysis for improving abiotic stress tolerances in radish.