Genome-Wide Association Studies and QTL Mapping Reveal a New Locus Associated with Resistance to Bacterial Pustule Caused by Xanthomonas citri pv. glycines in Soybean.

Bacterial pustule (BP), caused by Xanthomonas citri pv. glycines, is an important disease that, under favorable conditions, can drastically affect soybean production. We performed a genome-wide association study (GWAS) with a panel containing Brazilian and American cultivars, which were screened qualitatively and quantitatively against two Brazilian X. citri isolates (IBS 333 and IBS 327). The panel was genotyped using a genotyping by sequencing (GBS) approach, and we identified two main new regions in soybeans associated with X. citri resistance on chromosomes 6 (IBS 333) and 18 (IBS 327), different from the traditional rxp gene located on chromosome 17. The region on chromosome 6 was also detected by QTL mapping using a biparental cross between Williams 82 (R) and PI 416937 (S), showing that Williams 82 has another recessive resistance gene besides rxp, which was also detected in nine BP-resistant ancestors of the Brazilian cultivars (including CNS, S-100), based on haplotype analysis. Furthermore, we identified additional SNPs in strong LD (0.8) with peak SNPs by exploring variation available in WGS (whole genome sequencing) data among 31 soybean accessions. In these regions in strong LD, two candidate resistance genes were identified (Glyma.06g311000 and Glyma.18g025100) for chromosomes 6 and 18, respectively. Therefore, our results allowed the identification of new chromosomal regions in soybeans associated with BP disease, which could be useful for marker-assisted selection and will enable a reduction in time and cost for the development of resistant cultivars.

Time Course RNA-seq Reveals Soybean Responses against Root-Lesion Nematode and Resistance Players.

Pratylenchus brachyurus causes serious damage to soybean production and other crops worldwide. Plant molecular responses to RLN infection remain largely unknown and no resistance genes have been identified in soybean. In this study, we analyzed molecular responses to RLN infection in moderately resistant BRSGO (Chapadões-BRS) and susceptible TMG115 RR (TMG) Glycine max genotypes. Differential expression analysis revealed two stages of response to RLN infection and a set of differentially expressed genes (DEGs) in the first stage suggested a pattern-triggered immunity (PTI) in both genotypes. The divergent time-point of DEGs between genotypes was observed four days post-infection, which included the activation of mitogen-activated protein kinase (MAPK) and plant-pathogen interaction genes in the BRS, suggesting the occurrence of an effector-triggered immunity response (ETI) in BRS. The co-expression analyses combined with single nucleotide polymorphism (SNP) uncovered a key element, a transcription factor phytochrome-interacting factor (PIF7) that is a potential regulator of moderate resistance to RLN infection. Two genes for resistance-related leucine-rich repeat (LRR) proteins were found as BRS-specific expressed genes. In addition, alternative splicing analysis revealed an intron retention in a myo-inositol oxygenase (MIOX) transcript, a gene related to susceptibility, may cause a loss of function in BRS.

TILLING-by-Sequencing+ Reveals the Role of Novel Fatty Acid Desaturases (GmFAD2-2s) in Increasing Soybean Seed Oleic Acid Content.

Soybean is the second largest source of oil worldwide. Developing soybean varieties with high levels of oleic acid is a primary goal of the soybean breeders and industry. Edible oils containing high level of oleic acid and low level of linoleic acid are considered with higher oxidative stability and can be used as a natural antioxidant in food stability. All developed high oleic acid soybeans carry two alleles; GmFAD2-1A and GmFAD2-1B. However, when planted in cold soil, a possible reduction in seed germination was reported when high seed oleic acid derived from GmFAD2-1 alleles were used. Besides the soybean fatty acid desaturase (GmFAD2-1) subfamily, the GmFAD2-2 subfamily is composed of five members, including GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E. Segmental duplication of GmFAD2-1A/GmFAD2-1B, GmFAD2-2A/GmFAD2-2C, GmFAD2-2A/GmFAD2-2D, and GmFAD2-2D/GmFAD2-2C have occurred about 10.65, 27.04, 100.81, and 106.55 Mya, respectively. Using TILLING-by-Sequencing+ technology, we successfully identified 12, 8, 10, 9, and 19 EMS mutants at the GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E genes, respectively. Functional analyses of newly identified mutants revealed unprecedented role of the five GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E members in controlling the seed oleic acid content. Most importantly, unlike GmFAD2-1 members, subcellular localization revealed that members of the GmFAD2-2 subfamily showed a cytoplasmic localization, which may suggest the presence of an alternative fatty acid desaturase pathway in soybean for converting oleic acid content without substantially altering the traditional plastidial/ER fatty acid production.

The soybean gene GmHsp22.4 is involved in the resistance response to Meloidogyne javanica in Arabidopsis thaliana.

BACKGROUND: Small heat shock proteins (sHSPs) belong to the class of molecular chaperones that respond to biotic and abiotic stresses in plants. A previous study has showed strong induction of the gene GmHsp22.4 in response to the nematode Meloidogyne javanica in a resistant soybean genotype, while repression in a susceptible one. This study aimed to investigate the functional involvement of this small chaperone in response to M. javanica in Arabidopsis thaliana. First, it was evaluated the activation of the promoter region after the nematode inoculation, and the occurrence of polymorphisms between resistant and susceptible re-sequenced soybean accessions. Then functional analysis using A. thaliana lines overexpressing the soybean GmHsp22.4 gene, and knocked-out mutants were challenged with M. javanica infestation. RESULTS: High expression levels of the GFP gene marker in transformed A. thaliana plants revealed that the promoter region of GmHsp22.4 was strongly activated after nematode inoculation. Moreover, the multiplication of the nematode was significantly reduced in plants overexpressing GmHsp22.4 gene in A. thaliana compared to the wild type. Additionally, the multiplication of M. javanica in the A. thaliana mutants was significantly increased mainly in the event athsp22.0-2. This increase was not that evident in the event athsp22.0-1, the one that preserved a portion of the promoter region, including the HSEs in the region around - 83 bp. However, structural analysis at sequence level among soybean resistant and susceptible genotypes did not detect any polymorphisms in the whole gene model. CONCLUSIONS: The soybean chaperone GmHsp22.4 is involved in the defense response to root-knot nematode M. javanica in A. thaliana. Specifically, the promoter region covering until - 191 from the transcriptional start site (TSS) is necessary to promoter activation after nematode infection in Arabidopsis. No polymorphisms that could explain these differences in the defense response were detected in the GmHsp22.4 gene between resistant and susceptible soybean genotypes. Therefore, further investigation is needed to elucidate the triggering factor of the plant's defense mechanism, both at the sequence level of the soybean genotypes presenting contrasting reaction to root-knot nematode and by detecting cis-elements that are essential for the activation of the GmHsp22.4 gene promoter.