Genome-wide annotation of the soybean WRKY family and functional characterization of genes involved in response to Phakopsora pachyrhizi infection.

BACKGROUND: Many previous studies have shown that soybean WRKY transcription factors are involved in the plant response to biotic and abiotic stresses. Phakopsora pachyrhizi is the causal agent of Asian Soybean Rust, one of the most important soybean diseases. There are evidences that WRKYs are involved in the resistance of some soybean genotypes against that fungus. The number of WRKY genes already annotated in soybean genome was underrepresented. In the present study, a genome-wide annotation of the soybean WRKY family was carried out and members involved in the response to P. pachyrhizi were identified. RESULTS: As a result of a soybean genomic databases search, 182 WRKY-encoding genes were annotated and 33 putative pseudogenes identified. Genes involved in the response to P. pachyrhizi infection were identified using superSAGE, RNA-Seq of microdissected lesions and microarray experiments. Seventy-five genes were differentially expressed during fungal infection. The expression of eight WRKY genes was validated by RT-qPCR. The expression of these genes in a resistant genotype was earlier and/or stronger compared with a susceptible genotype in response to P. pachyrhizi infection. Soybean somatic embryos were transformed in order to overexpress or silence WRKY genes. Embryos overexpressing a WRKY gene were obtained, but they were unable to convert into plants. When infected with P. pachyrhizi, the leaves of the silenced transgenic line showed a higher number of lesions than the wild-type plants. CONCLUSIONS: The present study reports a genome-wide annotation of soybean WRKY family. The participation of some members in response to P. pachyrhizi infection was demonstrated. The results contribute to the elucidation of gene function and suggest the manipulation of WRKYs as a strategy to increase fungal resistance in soybean plants.

Genome-wide analysis of the Hsp20 gene family in soybean: comprehensive sequence, genomic organization and expression profile analysis under abiotic and biotic stresses.

BACKGROUND: The Hsp20 genes are associated with stress caused by HS and other abiotic factors, but have recently been found to be associated with the response to biotic stresses. These genes represent the most abundant class among the HSPs in plants, but little is known about this gene family in soybean. Because of their apparent multifunctionality, these proteins are promising targets for developing crop varieties that are better adapted to biotic and abiotic stresses. Thus, in the present study an in silico identification of GmHsp20 gene family members was performed, and the genes were characterized and subjected to in vivo expression analysis under biotic and abiotic stresses. RESULTS: A search of the available soybean genome databases revealed 51 gene models as potential GmHsp20 candidates. The 51 GmHsp20 genes were distributed across a total of 15 subfamilies where a specific predicted secondary structure was identified. Based on in vivo analysis, only 47 soybean Hsp20 genes were responsive to heat shock stress. Among the GmHsp20 genes that were potentials HSR, five were also cold-induced, and another five, in addition to one GmAcd gene, were responsive to Meloidogyne javanica infection. Furthermore, one predicted GmHsp20 was shown to be responsive only to nematode infection; no expression change was detected under other stress conditions. Some of the biotic stress-responsive GmHsp20 genes exhibited a divergent expression pattern between resistant and susceptible soybean genotypes under M. javanica infection. The putative regulatory elements presenting some conservation level in the GmHsp20 promoters included HSE, W-box, CAAT box, and TA-rich elements. Some of these putative elements showed a unique occurrence pattern among genes responsive to nematode infection. CONCLUSIONS: The evolution of Hsp20 family in soybean genome has most likely involved a total of 23 gene duplications. The obtained expression profiles revealed that the majority of the 51 GmHsp20 candidates are induced under HT, but other members of this family could also be involved in normal cellular functions, unrelated to HT. Some of the GmHsp20 genes might be specialized to respond to nematode stress, and the predicted promoter structure of these genes seems to have a particular conserved pattern related to their biological function.

Ubiquitous urease affects soybean susceptibility to fungi.

The soybean ubiquitous urease (encoded by GmEu4) is responsible for recycling metabolically derived urea. Additional biological roles have been demonstrated for plant ureases, notably in toxicity to other organisms. However, urease enzymatic activity is not related to its toxicity. The role of GmEu4 in soybean susceptibility to fungi was investigated in this study. A differential expression pattern of GmEu4 was observed in susceptible and resistant genotypes of soybeans over the course of a Phakopsora pachyrhizi infection, especially 24 h after infection. Twenty-nine adult, transgenic soybean plants, representing six independently transformed lines, were obtained. Although the initial aim of this study was to overexpress GmEu4, the transgenic plants exhibited GmEu4 co-suppression and decreased ureolytic activity. The growth of Rhizoctonia solani, Phomopsis sp., and Penicillium herguei in media containing a crude protein extract from either transgenic or non-transgenic leaves was evaluated. The fungal growth was higher in the protein extracts from transgenic urease-deprived plants than in extracts from non-transgenic controls. When infected by P. pachyrhizi uredospores, detached leaves of urease-deprived plants developed a significantly higher number of lesions, pustules and erupted pustules than leaves of non-transgenic plants containing normal levels of the enzyme. The results of the present work show that the soybean plants were more susceptible to fungi in the absence of urease. It was not possible to overexpress active GmEu4. For future work, overexpression of urease fungitoxic peptides could be attempted as an alternative approach.

Genomic and transcriptomic characterization of the transcription factor family R2R3-MYB in soybean and its involvement in the resistance responses to Phakopsora pachyrhizi.

Myb genes constitute one of the largest transcription factor families in the plant kingdom. Soybean MYB transcription factors have been related to the plant response to biotic stresses. Their involvement in response to Phakopsora pachyrhizi infection has been reported by several transcriptional studies. Due to their apparently highly diverse functions, these genes are promising targets for developing crop varieties resistant to diseases. In the present study, the identification and phylogenetic analysis of the soybean R2R3-MYB (GmMYB) transcription factor family was performed and the expression profiles of these genes under biotic stress were determined. GmMYBs were identified from the soybean genome using bioinformatic tools, and their putative functions were determined based on the phylogenetic tree and classified into subfamilies using guides AtMYBs describing known functions. The transcriptional profiles of GmMYBs upon infection with different pathogen were revealed by in vivo and in silico analyses. Selected target genes potentially involved in disease responses were assessed by RT-qPCR after different times of inoculation with P. pachyrhizi using different genetic backgrounds related to resistance genes (Rpp2 and Rpp5). R2R3-MYB transcription factors related to lignin synthesis and genes responsive to chitin were significantly induced in the resistant genotypes.