Characterization of Soybean WRKY Gene Family and Identification of Soybean WRKY Genes that Promote Resistance to Soybean Cyst Nematode

Integrating Computational and Biological Hemodynamic Approaches to Improve Modeling of Atherosclerotic Arteries

Atherosclerosis is the primary cause of cardiovascular disease, resulting in mortality, elevated healthcare costs, diminished productivity, and reduced quality of life for individuals and their communities. This is exacerbated by the limited understanding of its underlying causes and limitations in current therapeutic interventions, highlighting the need for sophisticated models of atherosclerosis. This review critically evaluates the computational and biological models of atherosclerosis, focusing on the study of hemodynamics in atherosclerotic coronary arteries. Computational models account for the geometrical complexities and hemodynamics of the blood vessels and stenoses, but they fail to capture the complex biological processes involved in atherosclerosis. Different in vitro and in vivo biological models can capture aspects of the biological complexity of healthy and stenosed vessels, but rarely mimic the human anatomy and physiological hemodynamics, and require significantly more time, cost, and resources. Therefore, emerging strategies are examined that integrate computational and biological models, and the potential of advances in imaging, biofabrication, and machine learning is explored in developing more effective models of atherosclerosis.

A comprehensive analysis of the WRKY family in soybean and functional analysis of GmWRKY164-GmGSL7c in resistance to soybean mosaic virus

BACKGROUND

Soybean mosaic disease caused by soybean mosaic virus (SMV) is one of the most devastating and widespread diseases in soybean producing areas worldwide. The WRKY transcription factors (TFs) are widely involved in plant development and stress responses. However, the roles of the GmWRKY TFs in resistance to SMV are largely unclear.

RESULTS

Here, 185 GmWRKYs were characterized in soybean (Glycine max), among which 60 GmWRKY genes were differentially expressed during SMV infection according to the transcriptome data. The transcriptome data and RT-qPCR results showed that the expression of GmWRKY164 decreased after imidazole treatment and had higher expression levels in the incompatible combination between soybean cultivar variety Jidou 7 and SMV strain N3. Remarkably, the silencing of GmWRKY164 reduced callose deposition and enhanced virus spread during SMV infection. In addition, the transcript levels of the GmGSL7c were dramatically lower upon the silencing of GmWRKY164. Furthermore, EMSA and ChIP-qPCR revealed that GmWRKY164 can directly bind to the promoter of GmGSL7c, which contains the W-box element.

CONCLUSION

Our findings suggest that GmWRKY164 plays a positive role in resistance to SMV infection by regulating the expression of GmGSL7c, resulting in the deposition of callose and the inhibition of viral movement, which provides guidance for future studies in understanding virus-resistance mechanisms in soybean.

Genome-wide profiling of WRKY genes involved in flavonoid biosynthesis in Erigeron breviscapus

The transcription factors of WRKY genes play essential roles in plant growth, stress responses, and metabolite biosynthesis. Erigeron breviscapus, a traditional Chinese herb, is abundant in flavonoids and has been used for centuries to treat cardiovascular and cerebrovascular diseases. However, the WRKY transcription factors that regulate flavonoid biosynthesis in E. breviscapus remain unknown. In this study, a total of 75 EbWRKY transcription factors were predicted through comprehensive genome-wide characterization of E. breviscapus and the chromosomal localization of each EbWRKY gene was investigated. RNA sequencing revealed transient responses of 74 predicted EbWRKY genes to exogenous abscisic acid (ABA), salicylic acid (SA), and gibberellin 3 (GA3) after 4 h of treatment. In contrast, the expression of key structural genes involved in flavonoid biosynthesis increased after 4 h in GA3 treatment. However, the content of flavonoid metabolites in leaves significantly increased at 12 h. The qRT-PCR results showed that the expression patterns of EbWRKY11, EbWRKY30, EbWRKY31, EbWRKY36, and EbWRKY44 transcription factors exhibited a high degree of similarity to the 11 structural genes involved in flavonoid biosynthesis. Protein-DNA interactions were performed between the key genes involved in scutellarin biosynthesis and candidate WRKYs. The result showed that F7GAT interacts with EbWRKY11, EbWRKY36, and EbWRKY44, while EbF6H has a self-activation function. This study provides comprehensive information on the regulatory control network of flavonoid accumulation mechanisms, offering valuable insights for breeding E. breviscapus varieties with enhanced scutellarin content.

The soybean plasma membrane GmDR1 protein conferring broad-spectrum disease and pest resistance regulates several receptor kinases and NLR proteins

Overexpression of Glycine max disease resistant 1 (GmDR1) exhibits broad-spectrum resistance against Fusarium virguliforme, Heterodera glycines (soybean cyst nematode), Tetranychus urticae (Koch) (spider mites), and Aphis glycines Matsumura (soybean aphids) in soybean. To understand the mechanisms of broad-spectrum immunity mediated by GmDR1, the transcriptomes of a strong and a weak GmDR1-overexpressor following treatment with chitin, a pathogen- and pest-associated molecular pattern (PAMP) common to these organisms, were investigated. The strong and weak GmDR1-overexpressors exhibited altered expression of 6098 and 992 genes, respectively, as compared to the nontransgenic control following chitin treatment. However, only 192 chitin- and 115 buffer-responsive genes exhibited over two-fold changes in expression levels in both strong and weak GmDR1-overexpressors as compared to the control. MapMan analysis of the 192 chitin-responsive genes revealed 64 biotic stress-related genes, of which 53 were induced and 11 repressed as compared to the control. The 53 chitin-induced genes include nine genes that encode receptor kinases, 13 encode nucleotide-binding leucine-rich repeat (NLR) receptor proteins, seven encode WRKY transcription factors, four ethylene response factors, and three MYB-like transcription factors. Investigation of a subset of these genes revealed three receptor protein kinases, seven NLR proteins, and one WRKY transcription factor genes that are induced following F. virguliforme and H. glycines infection. The integral plasma membrane GmDR1 protein most likely recognizes PAMPs including chitin and activates transcription of genes encoding receptor kinases, NLR proteins and defense-related genes. GmDR1 could be a pattern recognition receptor that regulates the expression of several NLRs for expression of PAMP-triggered immunity and/or priming the effector triggered immunity.

GmSABP2-1 encodes methyl salicylate esterase and functions in soybean defense against soybean cyst nematode

KEY MESSAGE

The soybean gene GmSABP2-1 encodes methyl salicylate esterase and its overexpression led to significant reduction in development of pathogenic soybean cyst nematode. Soybean cyst nematode (SCN, Heterodera glycines) is one of the most devastating pests of soybean (Glycine max L. Merr.). In searching for SCN-defense genes, a soybean gene of the methylesterase (MES) family was found to be upregulated in an SCN-resistant soybean line and downregulated in an SCN-susceptible line upon SCN infection. This gene was designated as GmSABP2-1. Here, we report on biochemical and overexpression studies of GmSABP2-1 to examine its possible function in SCN resistance. The protein encoded by GmSABP2-1 is closely related to known methyl salicylate esterases. To determine the biochemical function of GmSABP2-1, a full-length cDNA of GmSABP2-1 was cloned into a protein expression vector and expressed in Escherichia coli. The resulting recombinant GmSABP2-1 was demonstrated to catalyze the demethylation of methyl salicylate. The biochemical properties of GmSABP2-1 were determined. Its apparent Km value was 46.2 ± 2.2 μM for methyl salicylate, comparable to those of the known methyl salicylate esterases. To explore the biological significance of GmSABP2-1 in soybean defense against SCN, we first overexpressed GmSABP2-1 in transgenic hairy roots of an SCN-susceptible soybean line. When infected with SCN, GmSABP2-1-overexpressing hairy roots showed 84.5% reduction in the development of SCN beyond J2 stage. To provide further genetic evidence for the role of GmSABP2-1 in SCN resistance, stable transgenic soybean plants overexpressing GmSABP2-1 were produced. Analysis of the GmSABP2-1-overexpressing lines showed a significant reduction in SCN development compared to non-transgenic plants. In conclusion, we demonstrated that GmSABP2-1 encodes methyl salicylate esterase and functions as a resistance-related gene against SCN.

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

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

Identification of WRKY gene family in Dioscorea opposita Thunb. reveals that DoWRKY71 enhanced the tolerance to cold and ABA stress

WRKY transcription factors constitute one of the largest plant-specific gene families, regulating various aspects of plant growth, development, physiological processes, and responses to abiotic stresses. This study aimed to comprehensively analyze the WRKY gene family of yam (Dioscorea opposita Thunb.), to understand their expression patterns during the growth and development process and their response to different treatments of yam and analyze the function of DoWRKY71 in detail. A total of 25 DoWRKY genes were identified from the transcriptome of yam, which were divided into six clades (I, IIa, IIc, IId, IIe, III) based on phylogenetic analysis. The analysis of conserved motifs revealed 10 motifs, varying in length from 16 to 50 amino acids. Based on real-time quantitative PCR (qRT-PCR) analysis, DoWRKY genes were expressed at different stages of growth and development and responded differentially to various abiotic stresses. The expression level of DoWRKY71 genes was up-regulated in the early stage and then down-regulated in tuber enlargement. This gene showed responsiveness to cold and abiotic stresses, such as abscisic acid (ABA) and methyl jasmonate (MeJA). Therefore, further study was conducted on this gene. Subcellular localization analysis revealed that the DoWRKY71 protein was localized in the nucleus. Moreover, the overexpression of DoWRKY71 enhanced the cold tolerance of transgenic tobacco and promoted ABA mediated stomatal closure. This study presents the first systematic analysis of the WRKY gene family in yam, offering new insights for studying WRKY transcription factors in yam. The functional study of DoWRKY71 lays theoretical foundation for further exploring the regulatory function of the DoWRKY71 gene in the growth and development related signaling pathway of yam.

In Silico Comparison of WRKY Transcription Factors in Wild and Cultivated Soybean and Their Co-expression Network Arbitrating Disease Resistance

WRKY Transcription factors (TFs) play critical roles in plant defence mechanisms that are activated in response to biotic and abiotic stresses. However, information on the Glycine soja WRKYs (GsoWRKYs) is scarce. Owing to its importance in soybean breeding, here we identified putative WRKY TFs in wild soybean, and compared the results with Glycine max WRKYs (GmaWRKYs) by phylogenetic, conserved motif, and duplication analyses. Moreover, we explored the expression trends of WRKYs in G. max (oomycete, fungi, virus, bacteria, and soybean cyst nematode) and G. soja (soybean cyst nematode), and identified commonly expressed WRKYs and their co-expressed genes. We identified, 181 and 180 putative WRKYs in G. max and G. soja, respectively. Though the number of WRKYs in both studied species is almost the same, they differ in many ways, i.e., the number of WRKYs on corresponding chromosomes, conserved domain structures, WRKYGQK motif variants, and zinc-finger motifs. WRKYs in both species grouped in three major clads, i.e., I-III, where group-II had sub-clads IIa-IIe. We found that GsoWRKYs expanded mostly through segmental duplication. A large number of WRKYs were expressed in response to biotic stresses, i.e., Phakospora pachyrhizi, Phytoplasma, Heterodera glycines, Macrophomina phaseolina, and Soybean mosaic virus; 56 GmaWRKYs were commonly expressed in soybean plants infected with these diseases. Finally, 30 and 63 GmaWRKYs and GsoWRKYs co-expressed with 205 and 123 non-WRKY genes, respectively, indicating that WRKYs play essential roles in biotic stress tolerance in Glycine species.

Genome-wide identification, structural characterization and gene expression analysis of the WRKY transcription factor family in pea (Pisum sativum L.)

BACKGROUND

The WRKY gene family is one of the largest families of transcription factors in higher plants, and WRKY transcription factors play important roles in plant growth and development as well as in response to abiotic stresses; however, the WRKY gene family in pea has not been systematically reported.

RESULTS

In this study, 89 pea WRKY genes were identified and named according to the random distribution of PsWRKY genes on seven chromosomes. The gene family was found to have nine pairs of tandem duplicates and 19 pairs of segment duplicates. Phylogenetic analyses of the PsWRKY and 60 Arabidopsis WRKY proteins were performed to determine their homology, and the PsWRKYs were classified into seven subfamilies. Analysis of the physicochemical properties, motif composition, and gene structure of pea WRKYs revealed significant differences in the physicochemical properties within the PsWRKY family; however, their gene structure and protein-conserved motifs were highly conserved among the subfamilies. To further investigate the evolutionary relationships of the PsWRKY family, we constructed comparative syntenic maps of pea with representative monocotyledonous and dicotyledonous plants and found that it was most recently homologous to the dicotyledonous WRKY gene families. Cis-acting element analysis of PsWRKY genes revealed that this gene family can respond to hormones, such as abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin (GA), methyl jasmonate (MeJA), and salicylic acid (SA). Further analysis of the expression of 14 PsWRKY genes from different subfamilies in different tissues and fruit developmental stages, as well as under five different hormone treatments, revealed differences in their expression patterns in the different tissues and fruit developmental stages, as well as under hormone treatments, suggesting that PsWRKY genes may have different physiological functions and respond to hormones.

CONCLUSIONS

In this study, we systematically identified WRKY genes in pea for the first time and further investigated their physicochemical properties, evolution, and expression patterns, providing a theoretical basis for future studies on the functional characterization of pea WRKY genes during plant growth and development.

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

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

Unraveling the involvement of WRKY TFs in regulating plant disease defense signaling

MAIN CONCLUSION

This review article explores the intricate role, regulation, and signaling mechanisms of WRKY TFs in response to biotic stress, particularly emphasizing their pivotal role in the trophism of plant-pathogen interactions. Transcription factors (TFs) play a vital role in governing both plant defense and development by controlling the expression of various downstream target genes. Early studies have shown the differential expression of certain WRKY transcription factors by microbial infections. Several transcriptome-wide studies later demonstrated that diverse sets of WRKYs are significantly activated in the early stages of viral, bacterial, and fungal infections. Furthermore, functional investigations indicated that overexpression or silencing of certain WRKY genes in plants can drastically alter disease symptoms as well as pathogen multiplication rates. Hence the new aspects of pathogen-triggered WRKY TFs mediated regulation of plant defense can be explored. The already recognized roles of WRKYs include transcriptional regulation of defense-related genes, modulation of hormonal signaling, and participation in signal transduction pathways. Some WRKYs have been shown to directly bind to pathogen effectors, acting as decoys or resistance proteins. Notably, the signaling molecules like salicylic acid, jasmonic acid, and ethylene which are associated with plant defense significantly increase the expression of several WRKYs. Moreover, induction of WRKY genes or heightened WRKY activities is also observed during ISR triggered by the beneficial microbes which protect the plants from subsequent pathogen infection. To understand the contribution of WRKY TFs towards disease resistance and their exact metabolic functions in infected plants, further studies are required. This review article explores the intrinsic transcriptional regulation, signaling mechanisms, and hormonal crosstalk governed by WRKY TFs in plant disease defense response, particularly emphasizing their specific role against different biotrophic, hemibiotrophic, and necrotrophic pathogen infections.

Binding of stress-responsive OsWRKY proteins through WRKYGQK heptapeptide residue with the promoter region of two rice blast disease resistance genes Pi2 and Pi54 is important for development of blast resistance

Molecular docking was done to investigate the interactions between five differentially expressed rice WRKY proteins when challenged with the rice blast disease caused by Magnaporthe oryzae and drought stresses applied either individually or overlapped, with the promoter region of two blast resistance genes (Pi2 and Pi54). Molecular docking was performed using the HDOCK server. Initially, the homology models for each of the five rice WRKY proteins were prepared using I-TASSER server, and then the secondary structure as well as the DNA-binding pockets were predicted using PSIPRED and BindUP servers, respectively. The molecular docking study revealed a differential binding pattern of the rice WRKYs with the two blast resistance genes. The WRKY proteins (OsWRKY88 and OsWRKY102), whose transcript levels decrease when drought and blast stresses are overlapped, interact with the two resistance genes mostly involving the residues of the zinc finger structure. On the other hand, the WRKY proteins (OsWRKY53-1 and OsWRKY113), whose transcript levels did not reduce significantly when challenged by drought and blast overlapped condition compared to individual treatment of blast, interact mostly involving the residues of the conserved WRKYGQK heptapeptide sequence. Interestingly, the protein OsWRKY74 whose transcript levels are unaffected in both individual and overlapped stresses, interacts with both the blast resistance genes involving few residues of both WRKYGQK heptapeptide and the zinc finger structure. The findings thus indicate that the interaction of OsWRKY proteins involving the conserved WRKYGQK heptapeptide sequence with the blast resistance genes Pi2 and Pi54 is important to mitigate the blast challenge in rice even during overlapping challenges of drought.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03711-y.

Integrated transcriptome and physiological analysis revealed core transcription factors that promote flavonoid biosynthesis in apricot in response to pathogenic fungal infection

MAIN CONCLUSION

Integrated transcriptome and physiological analysis of apricot leaves after Fusarium solani treatment. In addition, we identified core transcription factors and flavonoid-related synthase genes which may function in apricot disease resistance. Apricot (Prunus armeniaca) is an important economic fruit species, whose yield and quality of fruit are limited owing to its susceptibility to diseases. However, the molecular mechanisms underlying the response of P. armeniaca to diseases is still unknown. In this study, we used physiology and transcriptome analysis to characterize responses of P. armeniaca subjected to Fusarium solani. The results showed increasing malondialdehyde (MDA) content, enhanced peroxidase (POD) and catalase (CAT) activity during F. solani infestation. A large number of differentially expressed genes (DEGs), which included 4281 upregulated DEGs and 3305 downregulated DEGs, were detected in P. armeniaca leaves exposed to F. solani infestation. Changes in expression of transcription factors (TFs), including bHLH, AP2/ERF, and WRKY indicated their role in triggering pathogen-responsive genes in P. armeniaca. During the P. armeniaca response to F. solani infestation, the content of total flavonoid was changed, and we identified enzyme genes associated with flavonoid biosynthesis. Ectopic overexpression of PabHLH15 and PabHLH102 in Nicotiana benthamiana conferred elevated resistance to Fspa_1. Moreover, PabHLH15 and PabHLH102 positively interact with the promoter of flavonoid biosynthesis-related genes. A regulatory network of TFs regulating enzyme genes related to flavonoid synthesis affecting apricot disease resistance was constructed. These results reveal the potential underlying mechanisms of the F. solani response of P. armeniaca, which would help improve the disease resistance of P. armeniaca and may cultivate high-quality disease-resistant varieties in the future.

Multi-algorithm cooperation research of WRKY genes under nitrogen stress in Panax notoginseng

WRKY transcription factors play an important role in the immune system and the innate defense response of plants. WRKY transcription factors have great feedback on nitrogen stress. In this study, bioinformatics was used to detect the WRKYs of Panax notoginseng (PnWRKYs). The response of PnWRKYs under nitrogen stress was also well studied. PnWRKYs were distributed on 11 chromosomes. According to PnWRKY and Arabidopsis thaliana WRKY (AtWRKY) domains, these PnWRKY proteins were divided into three groups by phylogenetic analysis. MEME analysis showed that almost every member contained motif 1 and motif 2. PlantCARE online predicted the cis-acting elements of the promoter. PnWRKY gene family members obtained 22 pairs of repeat fragments by collinearity analysis. The expression levels of PnWRKYs in different parts (roots, flowers, and leafs) were analyzed by the gene expression pattern. They reflected tissue-specific expressions. The qRT-PCR experiments were used to detect 74 PnWRKYs under nitrogen stress. The results showed that the expression levels of 8 PnWRKYs were significantly induced. The PnWRKY gene family may be involved in biotic/abiotic stresses and hormone induction. This study will not only lay the foundation to explore the functions of PnWRKYs but also provide candidate genes for the future improvement of P. notoginseng.

Biotechnological Tools to Elucidate the Mechanism of Plant and Nematode Interactions

Plant-parasitic nematodes (PPNs) pose a threat to global food security in both the developed and developing worlds. PPNs cause crop losses worth a total of more than USD 150 billion worldwide. The sedentary root-knot nematodes (RKNs) also cause severe damage to various agricultural crops and establish compatible relationships with a broad range of host plants. This review aims to provide a broad overview of the strategies used to identify the morpho-physiological and molecular events that occur during RKN parasitism. It describes the most current developments in the transcriptomic, proteomic, and metabolomic strategies of nematodes, which are important for understanding compatible interactions of plants and nematodes, and several strategies for enhancing plant resistance against RKNs. We will highlight recent rapid advances in molecular strategies, such as gene-silencing technologies, RNA interference (RNAi), and small interfering RNA (siRNA) effector proteins, that are leading to considerable progress in understanding the mechanism of plant-nematode interactions. We also take into account genetic engineering strategies, such as targeted genome editing techniques, the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) (CRISPR/Cas-9) system, and quantitative trait loci (QTL), to enhance the resistance of plants against nematodes.

Genome-Wide Association Analysis for Resistance to Coniothyrium glycines Causing Red Leaf Blotch Disease in Soybean

Soybean is a high oil and protein-rich legume with several production constraints. Globally, several fungi, viruses, nematodes, and bacteria cause significant yield losses in soybean. Coniothyrium glycines (CG), the causal pathogen for red leaf blotch disease, is the least researched and causes severe damage to soybean. The identification of resistant soybean genotypes and mapping of genomic regions associated with resistance to CG is critical for developing improved cultivars for sustainable soybean production. This study used single nucleotide polymorphism (SNP) markers generated from a Diversity Arrays Technology (DArT) platform to conduct a genome-wide association (GWAS) analysis of resistance to CG using 279 soybean genotypes grown in three environments. A total of 6395 SNPs was used to perform the GWAS applying a multilocus model Fixed and random model Circulating Probability Unification (FarmCPU) with correction of the population structure and a statistical test p-value threshold of 5%. A total of 19 significant marker-trait associations for resistance to CG were identified on chromosomes 1, 5, 6, 9, 10, 12, 13, 15, 16, 17, 19, and 20. Approximately 113 putative genes associated with significant markers for resistance to red leaf blotch disease were identified across soybean genome. Positional candidate genes associated with significant SNP loci-encoding proteins involved in plant defense responses and that could be associated with soybean defenses against CG infection were identified. The results of this study provide valuable insight for further dissection of the genetic architecture of resistance to CG in soybean. They also highlight SNP variants and genes useful for genomics-informed selection decisions in the breeding process for improving resistance traits in soybean.

Genome-Wide Analysis of WRKY Transcription Factors Involved in Abiotic Stress and ABA Response in Caragana korshinskii

The WRKY transcription factor family plays a vital role in plant development and environmental response. However, the information of WRKY genes at the genome-wide level is rarely reported in Caragana korshinskii. In this study, we identified and renamed 86 CkWRKY genes, which were further classified into three groups through phylogenetic analysis. Most of these WRKY genes were clustered and distributed on eight chromosomes. Multiple sequence alignment revealed that the conserved domain (WRKYGQK) of the CkWRKYs was basically consistent, but there were also six variation types (WRKYGKK, GRKYGQK, WRMYGQK, WRKYGHK, WKKYEEK and RRKYGQK) that appeared. The motif composition of the CkWRKYs was quite conservative in each group. In general, the number of WRKY genes gradually increased from lower to higher plant species in the evolutionary analysis of 28 species, with some exceptions. Transcriptomics data and RT-qPCR analysis showed that the CkWRKYs in different groups were involved in abiotic stresses and ABA response. Our results provided a basis for the functional characterization of the CkWRKYs involved in stress resistance in C. korshinskii.

Genome-wide identification of the WRKY gene family in Camellia oleifera and expression analysis under phosphorus deficiency

Camellia oleifera Abel. is an economically important woody edible-oil species that is mainly cultivated in hilly areas of South China. The phosphorus (P) deficiency in the acidic soils poses severe challenges for the growth and productivity of C. oleifera. WRKY transcription factors (TFs) have been proven to play important roles in biological processes and plant responses to various biotic/abiotic stresses, including P deficiency tolerance. In this study, 89 WRKY proteins with conserved domain were identified from the C. oleifera diploid genome and divided into three groups, with group II further classified into five subgroups based on the phylogenetic relationships. WRKY variants and mutations were detected in the gene structure and conserved motifs of CoWRKYs. Segmental duplication events were considered as the primary driver in the expanding process of WRKY gene family in C. oleifera. Based on transcriptomic analysis of two C. oleifera varieties characterized with different P deficiency tolerances, 32 CoWRKY genes exhibited divergent expression patterns in response to P deficiency stress. qRT-PCR analysis demonstrated that CoWRKY11, -14, -20, -29 and -56 had higher positive impact on P-efficient CL40 variety compared with P-inefficient CL3 variety. Similar expression trends of these CoWRKY genes were further observed under P deficiency with longer treatment period of 120d. The result indicated the expression sensitivity of CoWRKYs on the P-efficient variety and the C. oleifera cultivar specificity on the P deficiency tolerance. Tissue expression difference showed CoWRKYs may play a crucial role in the transportation and recycling P in leaves by affecting diverse metabolic pathways. The available evidences in the study conclusively shed light on the evolution of the CoWRKY genes in C. oleifera genome and provided a valuable resource for further investigation of functional characterization of WRKY genes involved to enhance the P deficiency tolerance in C. oleifera.

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

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

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.

Recent genome-wide replication promoted expansion and functional differentiation of the JAZs in soybeans

Jasmonate Zim-domain (JAZ) protein is an inhibitor of the jasmonate (JA) signal transduction pathway, and plays an important role in regulating plant growth, development, and defense. However, there have been few studies on its function under environmental stress in soybeans. In this study, a total of 275 JAZs protein-coding genes were identified in 29 soybean genomes. SoyC13 contained the least JAZ family members (26 JAZs), which was twice as high as AtJAZs. The genes are mainly generated by recent genome-wide replication (WGD), which replicated during the Late Cenozoic Ice Age. In addition, transcriptome analysis showed that the differences in gene expression patterns in the roots, stems, and leaves of the 29 cultivars at the V1 stage were not significant, but there was a significant difference among the three seed development stages. Finally, qRT-PCR results showed that GmJAZs responded the most strongly to heat stress, followed by drought and cold stress. This is consistent with the reason for their expansion and promoter analysis results. Therefore, we explored the significant role of conserved, duplicated, and neofunctionalized JAZs in the evolution of soybeans, which will contribute to the functional characterization of GmJAZ and the improvement of crops.

Effect of Pochonia chlamydosporia Endophytism and Meloidogyne incognita Parasitism on WRKYs and Defense Gene Expression in Tomato Roots

The expression of WRKY transcription factors and plant defense-related genes was studied in the roots of Apulian tomato (Solanum lycopersicum) cv Regina di Fasano (accessions MRT and PLZ) endophytically colonized by Pochonia chlamydosporia and parasitized or not by the root-knot nematode (RKN) Meloidogyne incognita. The effect on plant growth, nematode parasitism and histological aspect of the interaction were considered. The association of P. chlamydosporia with RKN-parasitized MRT plants increased the total biomass and shoot fresh weight in comparison with healthy plants and with those only parasitized by RKN, without the endophyte. However, the PLZ accession showed no significant difference in the observed biometric parameters. The number of RKN-induced galls per plant was not affected by endophytism eight days after inoculation. No histological changes were observed in the nematode feeding sites in the presence of the fungus. Gene expression analysis showed an accession-specific response to P. chlamydosporia with differential activation of WRKY-related genes. No significant change was found for WRKY76 expression in nematode-parasitized plants compared with control roots, confirming cultivar susceptibility. Data indicate genotype-specific responses of the WRKY genes to parasitism examined in roots with nematodes and/or endophytic P. chlamydosporia. At 25 days post-inoculation with P. chlamydosporia, no significant difference was observed in the expression of defense-related genes in both accessions, suggesting that salicylic acid (SA) (PAL and PR1) and jasmonate (JA) related genes (Pin II) are not active during endophytism.

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.

Identifications of QTLs and Candidate Genes Associated with Pseudomonas syringae Responses in Cultivated Soybean (Glycine max) and Wild Soybean (Glycine soja)

Soybeans (Glycine max) are a key food crop, serving as a valuable source of both oil and plant-derived protein. Pseudomonas syringae pv. glycinea (Psg) is among the most aggressive and prevalent pathogens affecting soybean production, causing a form of bacterial spot disease that impacts soybean leaves and thereby reduces crop yields. In this study, 310 natural soybean varieties were screened for Psg resistance and susceptibility. The identified susceptible and resistant varieties were then used for linkage mapping, BSA-seq, and whole genome sequencing (WGS) analyses aimed at identifying key QTLs associated with Psg responses. Candidate Psg-related genes were further confirmed through WGS and qPCR analyses. Candidate gene haplotype analyses were used to explore the associations between haplotypes and soybean Psg resistance. In addition, landrace and wild soybean plants were found to exhibit a higher degree of Psg resistance as compared to cultivated soybean varieties. In total, 10 QTLs were identified using chromosome segment substitution lines derived from Suinong14 (cultivated soybean) and ZYD00006 (wild soybean). Glyma.10g230200 was found to be induced in response to Psg, with the Glyma.10g230200 haplotype corresponding to soybean disease resistance. The QTLs identified herein can be leveraged to guide the marker-assisted breeding of soybean cultivars that exhibit partial resistance to Psg. Moreover, further functional and molecular studies of Glyma.10g230200 have the potential to offer insight into the mechanistic basis for soybean Psg resistance.

Identification and Functional Characterization of WRKY, PHD and MYB Three Salt Stress Responsive Gene Families in Mungbean (Vigna radiata L.)

WRKY-, PHD-, and MYB-like proteins are three important types of transcription factors in mungbeans, and play an important role in development and stress resistance. The genes' structures and characteristics were clearly reported and were shown to contain the conservative WRKYGQK heptapeptide sequence, Cys4-His-cys3 zinc binding motif, and HTH (helix) tryptophan cluster W structure, respectively. Knowledge on the response of these genes to salt stress is largely unknown. To address this issue, 83 VrWRKYs, 47 VrPHDs, and 149 VrMYBs were identified by using comparative genomics, transcriptomics, and molecular biology methods in mungbeans. An intraspecific synteny analysis revealed that the three gene families had strong co-linearity and an interspecies synteny analysis showed that mungbean and Arabidopsis were relatively close in genetic relationship. Moreover, 20, 10, and 20 genes showed significantly different expression levels after 15 days of salt treatment (p < 0.05; Log2 FC > 0.5), respectively. Additionally, in the qRT-PCR analysis, VrPHD14 had varying degrees of response to NaCl and PEG treatments after 12 h. VrWRKY49 was upregulated by ABA treatment, especially in the beginning (within 24 h). VrMYB96 was significantly upregulated in the early stages of ABA, NaCl, and PEG stress treatments (during the first 4 h). VrWRKY38 was significantly upregulated by ABA and NaCl treatments, but downregulated by PEG treatment. We also constructed a gene network centered on the seven DEGs under NaCl treatment; the results showed that VrWRKY38 was in the center of the PPI network and most of the homologous Arabidopsis genes of the interacted genes were reported to have response to biological stress. Candidate genes identified in this study provide abundant gene resources for the study of salt tolerance in mungbeans.

Genome-wide analysis of WRKY transcription factor genes in Toona sinensis: An insight into evolutionary characteristics and terpene synthesis

WRKY transcription factors (TFs), one of the largest TF families, serve critical roles in the regulation of secondary metabolite production. However, little is known about the expression pattern of WRKY genes during the germination and maturation processes of Toona sinensis buds. In the present study, the new assembly of the T. sinensis genome was used for the identification of 78 TsWRKY genes, including gene structures, phylogenetic features, chromosomal locations, conserved protein domains, cis-regulatory elements, synteny, and expression profiles. Gene duplication analysis revealed that gene tandem and segmental duplication events drove the expansion of the TsWRKYs family, with the latter playing a key role in the creation of new TsWRKY genes. The synteny and evolutionary constraint analyses of the WRKY proteins among T. sinensis and several distinct species provided more detailed evidence of gene evolution for TsWRKYs. Besides, the expression patterns and co-expression network analysis show TsWRKYs may multi-genes co-participate in regulating terpenoid biosynthesis. The findings revealed that TsWRKYs potentially play a regulatory role in secondary metabolite synthesis, forming the basis for further functional characterization of WRKY genes with the intention of improving T. sinensis.

Gene fusions, micro-exons and splice variants define stress signaling by AP2/ERF and WRKY transcription factors in the sesame pan-genome

Evolutionary dynamics of AP2/ERF and WRKY genes, the major components of defense response were studied extensively in the sesame pan-genome. Massive variation was observed for gene copy numbers, genome location, domain structure, exon-intron structure and protein parameters. In the pan-genome, 63% of AP2/ERF members were devoid of introns whereas >99% of WRKY genes contained multiple introns. AP2 subfamily was found to be micro-exon rich with the adjoining intronic sequences sharing sequence similarity to many stress-responsive and fatty acid metabolism genes. WRKY family included extensive multi-domain gene fusions where the additional domains significantly enhanced gene and exonic sizes as well as gene copy numbers. The fusion genes were found to have roles in acquired immunity, stress response, cell and membrane integrity as well as ROS signaling. The individual genomes shared extensive synteny and collinearity although ecological adaptation was evident among the Chinese and Indian accessions. Significant positive selection effects were noticed for both micro-exon and multi-domain genes. Splice variants with changes in acceptor, donor and branch sites were common and 6-7 splice variants were detected per gene. The study ascertained vital roles of lipid metabolism and chlorophyll biosynthesis in the defense response and stress signaling pathways. 60% of the studied genes localized in the nucleus while 20% preferred chloroplast. Unique cis-element distribution was noticed in the upstream promoter region with MYB and STRE in WRKY genes while MYC was present in the AP2/ERF genes. Intron-less genes exhibited great diversity in the promoter sequences wherein the predominance of dosage effect indicated variable gene expression levels. Mimicking the NBS-LRR genes, a chloroplast localized WRKY gene, Swetha_24868, with additional domains of chorismate mutase, cAMP and voltage-dependent potassium channel was found to act as a master regulator of defense signaling, triggering immunity and reducing ROS levels.

Melatonin and nitric oxide: Dual players inhibiting hazardous metal toxicity in soybean plants via molecular and antioxidant signaling cascades

Melatonin (MT), a ubiquitous signaling molecule, is known to improve plant growth. Its regulatory function alongside nitric oxide (NO) is known to induce heavy metal (Cd and Pb) stress tolerance, although the underlying mechanisms remain unknown. Here, we observed that the combined application of MT and NO remarkably enhanced plant biomass by reducing oxidative stress. Both MT and NO minimized metal toxicity by significantly lowering the levels of endogenous abscisic acid and jasmonic acid via downregulating NCED3 and upregulating catabolic genes (CYP707A1 and CYP707A2). MT/NO-induced mitigation of Cd and Pb stress was associated with increased endo-melatonin and variable endo-S-nitrosothiol levels caused by enhanced expression of gmNR and gmGSNOR mRNAs. Remarkably, the combined application of MT/NO reduced soil Cd and Pb mobilization by increasing the uptake of Ca²⁺ and K⁺ and increasing the exudation of organic acids into the rhizosphere. These results correlated with the upregulation of MTF-1 and WARKY27 during metal translocation. MT/NO regulates the MAPK and CDPK cascades to promote plant cell survival and Ca²⁺ signaling, thereby imparting resistance to heavy metal toxicity. In conclusion, MT/NO modulates the stress-resistance machinery to mitigate Cd and Pb toxicity by regulating the activation of antioxidant and molecular transcription factors.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Akebiatrifoliata: A Bioinformatics Study

WRKY transcription factors have been found in most plants and play an important role in regulating organ growth and disease response. Outlining the profile of WRKY genes is a very useful project for studying morphogenesis and resistance formation. In the present study, a total of 63 WRKY genes consisting of 13 class I, 41 class II, and 9 class III genes were identified from the newly published A. trifoliata genome, of which 62 were physically distributed on all 16 chromosomes. Structurally, two AkWRKY genes (AkWRKY6 and AkWRKY52) contained four domains, and AkWRKY17 lacked the typical heptapeptide structure. Evolutionarily, 42, 16, and 5 AkWRKY genes experienced whole genome duplication (WGD) or fragmentation, dispersed duplication, and tandem duplication, respectively; 28 Ka/Ks values of 30 pairs of homologous genes were far lower than 1, while those of orthologous gene pairs between AkWRKY41 and AkWRKY52 reached up to 2.07. Transcriptome analysis showed that many of the genes were generally expressed at a low level in 12 fruit samples consisting of three tissues, including rind, flesh, and seeds, at four developmental stages, and interaction analysis between AkWRKY and AkNBS genes containing W-boxes suggested that AkWRKY24 could play a role in plant disease resistance by positively regulating AkNBS18. In summary, the WRKY gene family of A. trifoliata was systemically characterized for the first time, and the data and information obtained regarding AkWRKY could be very useful in further theoretically elucidating the molecular mechanisms of plant development and response to pathogens and practically improving favorable traits such as disease resistance.

Exogenous genistein enhances soybean resistance to Xanthomonas axonopodis pv. glycines

BACKGROUND

Xanthomonas axonopodis pv. glycines (Xag) is the causal agent of bacterial pustule disease and results in enormous losses in soybean production. Although isoflavones are known to be involved in soybean resistance against pathogen infection, the effects of exogenous isoflavones on soybean plants remain unexplored.

RESULTS

Irrigation of soybean plants with isoflavone genistein inhibited plant growth for short periods, probably by inhibiting the tyrosine (brassinosteroids) kinase pathway, and increased disease resistance against Xag. The number of lesions was reduced by 59%-63% when applying 50 μg ml-1 genistein. The effects on disease resistance were observed for 15 days after treatment. Genistein also enhanced the disease resistance of soybean against the fungal pathogen Sclerotinia sclerotiorum. Exogenous genistein increased antioxidant capacity, decreased H2 O2 level and promoted the accumulation of phenolics in Xag-infected soybean leaves. Exogenous genistein reduced the amounts of endogenous daidzein, genistein and glycitein and increased the concentration of genistin, which was found to show strong antibacterial activity against the pathogen and to reduce the expression of virulence factor yapH, and flagella formation gene flgK. The expression of several soybean defense genes, such as chalcone isomerase, glutathione S-transferase and 1-aminocyclopropane-1-carboxylate oxidase 1, was upregulated after genistein treatment.

CONCLUSIONS

The effects of exogenous genistein on soybean plants were examined for the first time, revealing new insights into the roles of isoflavones in soybean defense and demonstrating that irrigation with genistein can be a suitable method to induce disease resistance in soybean plants. © 2022 Society of Chemical Industry.

Essential Roles of Cupredoxin Family Proteins in Soybean Cyst Nematode Resistance

Soybean cyst nematode (SCN, Heterodera glycines), one of the most devastating soybean pathogens, causes a significant yield loss in soybean production. One of the most effective ways to manage SCN is to grow resistant cultivars. Therefore, comparative study using resistant and susceptible soybean cultivars provides a powerful tool to identify new genes involved in soybean SCN resistance. In the present study, a transcriptome analysis was carried out using both the resistant (PI88788) and susceptible (Williams 82) soybean cultivars to characterize the responses to nematode infection. Various defense-related genes and different pathways involved in nematode resistance were recognized as being highly expressed in resistant cultivar. Promoter-GUS analysis was conducted to monitor the spatial expression pattern of the genes highly induced by nematode infection. Two nematode-inducible promoters for Glyma.05g147000 (encoding caffeoyl-CoA O-methyltransferase) and Glyma.06g036700 (encoding cupredoxin superfamily protein) were characterized, and the promoters could efficiently drive the expression of known nematode resistance genes (α-SNAPRhg1HC or GmSHMT) to affect soybean SCN resistance. Interestingly, expression of the cupredoxin family genes was upregulated not only by SCN, but also by jasmonic acid treatment. DNA sequence analysis identified that a conserved motif (GGTGCATG) with high similarity to SCNbox1 and GC-rich element is enriched in their promoter regions, suggesting its potential to serve as a nematode-responsive regulatory element. Overexpression of Glyma.06g036700 significantly enhanced soybean resistance to cyst nematode. Overall, our findings not only highlight the essential role of cupredoxin family genes in SCN resistance, but also offer potential functional tools to develop nematode resistance in crops.

RNA-Seq and Gene Ontology Analysis Reveal Differences Associated With Low R/FR-Induced Shade Responses in Cultivated Lentil and a Wild Relative

Lentil is an important pulse crop not only because of its high nutrient value but also because of its ecological advantage in a sustainable agricultural system. Our previous work showed that the cultivated lentil and wild lentil germplasm respond differently to light environments, especially to low R/FR-induced shade conditions. Little is known about how cultivated and wild lentils respond to shade at the level of gene expression and function. In this study, transcriptomic profiling of a cultivated lentil (Lupa, L. culinaris) and a wild lentil (BGE 016880, L. orientalis) at several growth stages is presented. De novo transcriptomes were assembled for both genotypes, and differential gene expression analysis and gene ontology enrichment analysis were performed. The transcriptomic resources generated in this study provide fundamental information regarding biological processes and genes associated with shade responses in lentils. BGE 016880 and Lupa shared a high similarity in their transcriptomes; however, differential gene expression profiles were not consistent between these two genotypes. The wild lentil BGE 016880 had more differentially expressed genes than the cultivated lentil Lupa. Upregulation of genes involved in gibberellin, brassinosteroid, and auxin synthesis and signaling pathways, as well as cell wall modification, in both genotypes explains their similarity in stem elongation response under the shade. Genes involved in jasmonic acid and flavonoid biosynthesis pathways were downregulated in BGE 016880 only, and biological processes involved in defense responses were significantly enriched in the wild lentil BGE 016880 only. Downregulation of WRKY and MYB transcription factors could contribute to the reduced defense response in BGE 016880 but not in Lupa under shade conditions. A better understanding of shade responses of pulse crop species and their wild relatives will play an important role in developing genetic strategies for crop improvement in response to changes in light environments.

A Comprehensive Study of the WRKY Transcription Factor Family in Strawberry

WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.

GmWRKY81 Encoding a WRKY Transcription Factor Enhances Aluminum Tolerance in Soybean

Aluminum (Al) toxicity is an essential factor that adversely limits soybean (Glycine max (L.) Merr.) growth in acid soils. WRKY transcription factors play important roles in soybean responses to abiotic stresses. Here, GmWRKY81 was screened from genes that were differentially expressed under Al treatment in Al-tolerant soybean Baxi10 and Al-sensitive soybean Bendi2. We found that GmWRKY81 was significantly induced by 20 μM AlCl3 and upregulated by AlCl3 treatment for 2 h. In different tissues, the expression of GmWRKY81 was differentially induced. In 0-1 cm root tips, the expression of GmWRKY81 was induced to the highest level. The overexpression of GmWRKY81 in soybean resulted in higher relative root elongation, root weight, depth, root length, volume, number of root tips and peroxidase activity but lower root average diameter, malonaldehyde and H2O2 contents, indicating enhanced Al tolerance. Moreover, RNA-seq identified 205 upregulated and 108 downregulated genes in GmWRKY81 transgenic lines. Fifteen of these genes that were differentially expressed in both AlCl3-treated and GmWRKY81-overexpressing soybean had the W-box element, which can bind to the upstream-conserved WRKY domain. Overall, the combined functional analysis indicates that GmWRKY81 may improve soybean Al tolerance by regulating downstream genes participating in Al3+ transport, organic acid secretion and antioxidant reactions.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

Analysis of WRKY Resistance Gene Family in Boehmeria nivea (L.) Gaudich: Crosstalk Mechanisms of Secondary Cell Wall Thickening and Cadmium Stress

A total of 60 WRKY family genes of ramie were identified in the ramie. The genes were unevenly distributed across 14 chromosomes in the specie and highly concentrated (72%) in the distal telomeric region. Phylogenetic analysis placed these genes into seven distinct subfamilies groups: I, II (a, b, c, d, e), and III, with group IIc containing only the variant of heptapetide sequence (WRKYGKK). Segmental duplication events (41.7%) was found to be the main driver of BnGWRKY evolution. Thirty eight from among the genes showed collinear relationships with WRKY genes from Arabidopsis thaliana, Cannabis sativa, Oryza sativa, and Zea mays. The number and density of stress and hormone responsives cis-acting elements were comparably higher than other elements, with abundant ARE and rare LTR cis-acting elements indicating the long-standing adaptability of ramie to its natural environment. GO and KEGG enrichment analysis of the WRKY target genes revealed their involvement in response to stimuli, immune system processes, transporter protein activity and antioxidant activity. Expression analysis show that most WRKYs were activated by the cadmium stress, more especially the BnGWRKY2, BnGWRKY15, BnGWRKY20, BnGWRKY50 and BnGWRKY58. Combining transcriptome, orthologous gene relationships and qPCR result, we established the possible involvement of BnGWRKY50 and BnGWRKY58 in crosstalk mechanism between secondary cell wall thickening and Cd²⁺ stress. This provided information into the role of BnGWRKY proteins in ramie secondary wall development and cadmium stress response to, and could serve as basis for improvement of the ramie.

Large-Scale Integrative Analysis of Soybean Transcriptome Using an Unsupervised Autoencoder Model

Plant tissues are distinguished by their gene expression patterns, which can help identify tissue-specific highly expressed genes and their differential functional modules. For this purpose, large-scale soybean transcriptome samples were collected and processed starting from raw sequencing reads in a uniform analysis pipeline. To address the gene expression heterogeneity in different tissues, we utilized an adversarial deconfounding autoencoder (AD-AE) model to map gene expressions into a latent space and adapted a standard unsupervised autoencoder (AE) model to help effectively extract meaningful biological signals from the noisy data. As a result, four groups of 1,743, 914, 2,107, and 1,451 genes were found highly expressed specifically in leaf, root, seed and nodule tissues, respectively. To obtain key transcription factors (TFs), hub genes and their functional modules in each tissue, we constructed tissue-specific gene regulatory networks (GRNs), and differential correlation networks by using corrected and compressed gene expression data. We validated our results from the literature and gene enrichment analysis, which confirmed many identified tissue-specific genes. Our study represents the largest gene expression analysis in soybean tissues to date. It provides valuable targets for tissue-specific research and helps uncover broader biological patterns. Code is publicly available with open source at https://github.com/LingtaoSu/SoyMeta.

Interaction between Rag genes results in a unique synergistic transcriptional response that enhances soybean resistance to soybean aphids

BACKGROUND

Pyramiding different resistance genes into one plant genotype confers enhanced resistance at the phenotypic level, but the molecular mechanisms underlying this effect are not well-understood. In soybean, aphid resistance is conferred by Rag genes. We compared the transcriptional response of four soybean genotypes to aphid feeding to assess how the combination of Rag genes enhanced the soybean resistance to aphid infestation.

RESULTS

A strong synergistic interaction between Rag1 and Rag2, defined as genes differentially expressed only in the pyramid genotype, was identified. This synergistic effect in the Rag1/2 phenotype was very evident early (6 h after infestation) and involved unique biological processes. However, the response of susceptible and resistant genotypes had a large overlap 12 h after aphid infestation. Transcription factor (TF) analyses identified a network of interacting TF that potentially integrates signaling from Rag1 and Rag2 to produce the unique Rag1/2 response. Pyramiding resulted in rapid induction of phytochemicals production and deposition of lignin to strengthen the secondary cell wall, while repressing photosynthesis. We also identified Glyma.07G063700 as a novel, strong candidate for the Rag1 gene.

CONCLUSIONS

The synergistic interaction between Rag1 and Rag2 in the Rag1/2 genotype can explain its enhanced resistance phenotype. Understanding molecular mechanisms that support enhanced resistance in pyramid genotypes could facilitate more directed approaches for crop improvement.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Molecular interaction of charcoal rot pathogenesis in soybean: a complex interaction

Charcoal rot (CR) is a major disease of soybean, which is caused by Macrophomina phaseolina (Mp). Increasing temperatures and low rainfall in recent years have immensely benefitted the pathogen. Hence, the search for genetically acquired resistance to this pathogen is essential. The pathogen is a hemibiotroph, which germinates on the root surface and colonizes epidermal tissue. Several surface receptors initiate pathogenesis, followed by the secretion of various enzymes that provide entry to host tissue. Several enzymes and other converging cascades in the pathogen participate against host defensive responses. β-glucan of the fungal cell wall is recognized as MAMPs (microbe-associated molecular patterns) in plants, which trigger host immune responses. Kinase receptors, resistance, and pathogenesis-related genes correspond to host defense response. They work in conjunction with hormone-mediated defense pathway especially, the systemic acquired resistance, calcium-signaling, and production of phytoalexins. Due to its quantitative nature, limited QTLs have been identified in soybean for CR resistance. The present review attempts to provide a functional link between M. phaseolina pathogenicity and soybean responses. Elucidation of CR resistance responses would facilitate improved designing of breeding programs, and may help in the selection of corresponding genes to introgress CR resistant traits.

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.

Based on the whole genome clarified the evolution and expression process of fatty acid desaturase genes in three soybeans

Soybean is an important oil crop cultivated worldwide. With the increasing global population crossed with growing challenging cultivation conditions, improving soybean breeding by selecting important traits is urgent needed. Genes coding for plant fatty acid desaturases (FADs) genes are major candidates for that, because they are involving in controlling fatty acid composition and holding membrane fluidity under abiotic stress. Here, 75 FADs were found in three soybean genomes, which were further classified into four sub-groups. Phylogenetic tree, gene structure, motif and promoter analysis showed that the FAD gene family was conserved in the three soybeans. In addition, the numbers of omega desaturase from Chinese cultivated varieties were significantly higher than those in Chinese wild soybean and ancient polyploid soybean, respectively. However, it was the opposite for the sphingolipid subfamily. These results indicated that each subfamily was subjected to different selection pressures during cultivation and domestication. As the extra genes of the subfamily were very close to other family members' positions on chromosomes, they should be produced by duplication. The cis-element analysis of FAD promoter sequences revealed that upstream sequences of FAD contained abundant light, hormone and abiotic stress responsive cis-elements, suggesting that the quality of soybean could be improved by regulating these stresses. Expression analysis of Chinese wild soybean under salt stress showed that GsDES1.1, GsDES1.2, GsFAD2.1 and GsSLD1 in leaves and GsSLD2, GsSLD5 and GsSLD6 in roots were not closely related to salt stress response. Therefore, we explored the significant role of conserved, duplicated and neofunctionalized FAD in the domestication of soybean, which contributes to the importance of soybean as a global oil crop.

Harnessing the potential of plant transcription factors in developing climate resilient crops to improve global food security: Current and future perspectives

Crop plants should be resilient to climatic factors in order to feed ever-increasing populations. Plants have developed stress-responsive mechanisms by changing their metabolic pathways and switching the stress-responsive genes. The discovery of plant transcriptional factors (TFs), as key regulators of different biotic and abiotic stresses, has opened up new horizons for plant scientists. TFs perceive the signal and switch certain stress-responsive genes on and off by binding to different cis-regulatory elements. More than 50 families of plant TFs have been reported in nature. Among them, DREB, bZIP, MYB, NAC, Zinc-finger, HSF, Dof, WRKY, and NF-Y are important with respect to biotic and abiotic stresses, but the potential of many TFs in the improvement of crops is untapped. In this review, we summarize the role of different stress-responsive TFs with respect to biotic and abiotic stresses. Further, challenges and future opportunities linked with TFs for developing climate-resilient crops are also elaborated.

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

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

Genome-Wide Identification of WRKY Genes in Artemisia annua: Characterization of a Putative Ortholog of AtWRKY40

Artemisia annua L. is well-known as the plant source of artemisinin, a sesquiterpene lactone with effective antimalarial activity. Here, a putative ortholog of the Arabidopsis thaliana WRKY40 transcription factor (TF) was isolated via reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends in A. annua and named AaWRKY40. A putative nuclear localization domain was identified in silico and experimentally confirmed by using protoplasts of A. annua transiently transformed with AaWRKY40-GFP. A genome-wide analysis identified 122 WRKY genes in A. annua, and a manually curated database was obtained. The deduced proteins were categorized into the major WRKY groups, with group IIa containing eight WRKY members including AaWRKY40. Protein motifs, gene structure, and promoter regions of group IIa WRKY TFs of A. annua were characterized. The promoter region of AaWRKY group IIa genes contained several abiotic stress cis-acting regulatory elements, among which a highly conserved W-box motif was identified. Expression analysis of AaWRKY40 compared to AaWRKY1 in A. annua cell cultures treated with methyl jasmonate known to enhance artemisinin production, suggested a possible involvement of AaWRKY40 in terpenoid metabolism. Further investigation is necessary to study the role of AaWRKY40 and possible interactions with other TFs in A. annua.

The Effect of Virulence and Resistance Mechanisms on the Interactions between Parasitic Plants and Their Hosts

Parasitic plants have a unique heterotrophic lifestyle based on the extraction of water and nutrients from host plants. Some parasitic plant species, particularly those of the family Orobanchaceae, attack crops and cause substantial yield losses. The breeding of resistant crop varieties is an inexpensive way to control parasitic weeds, but often does not provide a long-lasting solution because the parasites rapidly evolve to overcome resistance. Understanding mechanisms underlying naturally occurring parasitic plant resistance is of great interest and could help to develop methods to control parasitic plants. In this review, we describe the virulence mechanisms of parasitic plants and resistance mechanisms in their hosts, focusing on obligate root parasites of the genera Orobanche and Striga. We noticed that the resistance (R) genes in the host genome often encode proteins with nucleotide-binding and leucine-rich repeat domains (NLR proteins), hence we proposed a mechanism by which host plants use NLR proteins to activate downstream resistance gene expression. We speculated how parasitic plants and their hosts co-evolved and discussed what drives the evolution of virulence effectors in parasitic plants by considering concepts from similar studies of plant-microbe interaction. Most previous studies have focused on the host rather than the parasite, so we also provided an updated summary of genomic resources for parasitic plants and parasitic genes for further research to test our hypotheses. Finally, we discussed new approaches such as CRISPR/Cas9-mediated genome editing and RNAi silencing that can provide deeper insight into the intriguing life cycle of parasitic plants and could potentially contribute to the development of novel strategies for controlling parasitic weeds, thereby enhancing crop productivity and food security globally.

Transcription Factor GmWRKY142 Confers Cadmium Resistance by Up-Regulating the Cadmium Tolerance 1-Like Genes

Cadmium (Cd) is a widespread pollutant that is toxic to living organisms. Previous studies have identified certain WRKY transcription factors, which confer Cd tolerance in different plant species. In the present study, we have identified 29 Cd-responsive WRKY genes in Soybean [Glycine max (L.) Merr.], and confirmed that 26 of those GmWRKY genes were up-regulated, while 3 were down-regulated. We have also cloned the novel, positively regulated GmWRKY142 gene from soybean and investigated its regulatory mechanism in Cd tolerance. GmWRKY142 was highly expressed in the root, drastically up-regulated by Cd, localized in the nucleus, and displayed transcriptional activity. The overexpression of GmWRKY142 in Arabidopsis thaliana and soybean hairy roots significantly enhanced Cd tolerance and lead to extensive transcriptional reprogramming of stress-responsive genes. ATCDT1, GmCDT1-1, and GmCDT1-2 encoding cadmium tolerance 1 were induced in overexpression lines. Further analysis showed that GmWRKY142 activated the transcription of ATCDT1, GmCDT1-1, and GmCDT1-2 by directly binding to the W-box element in their promoters. In addition, the functions of GmCDT1-1 and GmCDT1-2, responsible for decreasing Cd uptake, were validated by heterologous expression in A. thaliana. Our combined results have determined GmWRKYs to be newly discovered participants in response to Cd stress, and have confirmed that GmWRKY142 directly targets ATCDT1, GmCDT1-1, and GmCDT1-2 to decrease Cd uptake and positively regulate Cd tolerance. The GmWRKY142-GmCDT1-1/2 cascade module provides a potential strategy to lower Cd accumulation in soybean.

Characterization of Pingliang xiaoheidou (ZDD 11047), a soybean variety with resistance to soybean cyst nematode Heterodera glycines

KEY MESSAGE

A novel QTL (qSCN-PL10) for SCN resistance and related candidate genes were identified in the soybean variety Pingliang xiaoheidou, and plant basal immunity seems to contribute to the SCN resistance. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating soybean pests worldwide. The development of host plant resistance represents an effective strategy to control SCN. However, owing to the lack of diversity of resistance genes in soybean varieties, further investigation is necessary to identify new SCN resistance genes. By analyzing the resistance phenotypes of soybean variety Pingliang xiaoheidou (Pingliang, ZDD 11047), we found that it exhibited the different resistance phenotypes from PI 88788 and Peking varieties. Because Pingliang variety contains the Rhg1-a (low copy) haplotype and lacks the resistant Rhg4 haplotype, novel quantitative trait locus might account for their SCN resistance. After sequencing parental lines (Magellan and Pingliang) and 200 F2:3 progenies, a high-density genetic map was constructed using the specific length amplified fragment sequencing method and qSCN-PL10 was identified as a novel locus for SCN resistance. Candidate genes were predicted by RNA sequencing (RNA-seq) in the qSCN-PL10 locus region. The RNA-seq analysis performed also indicated that plant basal immunity plays an important role in the resistance of Pingliang to SCN. These results lay a foundation for the use of marker-assisted breeding to enhance the resistance to SCN.