Abstract:The WRKY transcription factor gene family is a plant-specific transcription factor that plays important roles defense responses. Studies in model plant Arabidopsis demonstrated that WRKYs function downstream of mitogen activated-protein kinase (MAPK) signaling cascade and participate in defense responses through activating the expression of defense-related genes. However, the roles of WRKYs in defense responses have not been previously investigated in paleopolyploidy soybean. Bioinfomatic analysis revealed that there are three pair of GmWRKY33 genes in the soybean genome. The identity of first two pair of GmWRKY33 genes is greater than 84% (named as GmWRKY33A). The identity of genes within the same pair is greater than 95%. A 300 bp fragment highly homologous to these four GmWRKY33A was chosen to clone into bean pod mosaic virus (BPMV)-based silencing vector (BPMV-VIGS) to achieve the goal of silencing four GmWRKY33A genes simultaneously. In this study, we simultaneously silenced four homologous genes of GmWRKY33A using a bean pod mottle virus (BPMV) vector carrying a single fragment of GmWRKY33A. Comparing the silenced plants with the vector control plants, no evident morphological phenotypes were observed. However, the GmWRKY33A-silenced plants exhibited significantly reduced resistance to Pseudomonas syringae pv. glycinea (Psg), Xanthomonas axonopodis pv. glycine (Xag), as well as to soybean mosaic virus (SMV). Furthermore, we demonstrated that silencing these GmWRKY33A genes significantly inhibited the activation of GmMPK3/GmMPK6 induced by Psg infection. Collectively, our results suggest that GmWRKY33As are involved in soybean immunity through regulating the transcription of GmMPK3/6 genes or activating the kinase activities of GmMPK3/6. Taken together, our results demonstrated that GmWRKY33As are positive regulators of soybean immune responses.