The soybean GmSNAP18 gene underlies two types of resistance to soybean cyst nematode

Genome scans for selection signatures identify candidate virulence genes for adaptation of the soybean cyst nematode to host resistance

Plant pathogens are constantly under selection pressure for host resistance adaptation. Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean primarily managed through resistant cultivars; however, SCN populations have evolved virulence in response to selection pressures driven by repeated monoculture of the same genetic resistance. Resistance to SCN is mediated by multiple epistatic interactions between Rhg (for resistance to H. glycines) genes. However, the identity of SCN virulence genes that confer the ability to overcome resistance remains unknown. To identify candidate genomic regions showing signatures of selection for increased virulence, we conducted whole genome resequencing of pooled individuals (Pool-Seq) from two pairs of SCN populations adapted on soybeans with Peking-type (rhg1-a, rhg2, and Rhg4) resistance. Population differentiation and principal component analysis-based approaches identified approximately 0.72-0.79 million SNPs, the frequency of which showed potential selection signatures across multiple genomic regions. Chromosomes 3 and 6 between population pairs showed the greatest density of outlier SNPs with high population differentiation. Conducting multiple outlier detection tests to identify overlapping SNPs resulted in a total of 966 significantly differentiated SNPs, of which 285 exon SNPs were mapped to 97 genes. Of these, six genes encoded members of known stylet-secreted effector protein families potentially involved in host defence modulation including venom-allergen-like, annexin, glutathione synthetase, SPRYSEC, chitinase, and CLE effector proteins. Further functional analysis of identified candidate genes will provide new insights into the genetic mechanisms by which SCN overcomes soybean resistance and inform the development of molecular markers for rapidly screening the virulence profile of an SCN-infested field.

Structural insights into binding of polyglutamylated tetrahydrofolate by serine hydroxymethyltransferase 8 from soybean

Tetrahydrofolate and its derivatives participate in one-carbon transfer reactions in all organisms. The cellular form of tetrahydrofolate (THF) is modified by multiple glutamate residues and polyglutamylation plays a key role in organellar and cellular folate homeostasis. In addition, polyglutamylation of THF is known to increase the binding affinity to enzymes in the folate cycle, many of which can utilize polyglutamylated THF as a substrate. Here, we use X-ray crystallography to provide a high-resolution view of interactions between the enzyme serine hydroxymethyltransferase (SHMT), which provides one carbon precursors for the folate cycle, and a polyglutamylated form of THF. Our 1.7 Å crystal structure of soybean SHMT8 in complex with diglutamylated 5-formyl-THF reveals, for the first time, a structural rearrangement of a loop at the entrance to the folate binding site accompanied by the formation of novel specific interactions between the enzyme and the diglutamyl tail of the ligand. Biochemical assays show that additional glutamate moieties on the folate ligand increase both enzyme stability and binding affinity. Together these studies provide new information on SHMT structure and function and inform the design of anti-folate agents.

No Pairwise Interactions of GmSNAP18, GmSHMT08 and AtPR1 with Suppressed AtPR1 Expression Enhance the Susceptibility of Arabidopsis to Beet Cyst Nematode

GmSNAP18 and GmSHMT08 are two major genes conferring soybean cyst nematode (SCN) resistance in soybean. Overexpression of either of these two soybean genes would enhance the susceptibility of Arabidopsis to beet cyst nematode (BCN), while overexpression of either of their corresponding orthologs in Arabidopsis, AtSNAP2 and AtSHMT4, would suppress it. However, the mechanism by which these two pairs of orthologous genes boost or inhibit BCN susceptibility of Arabidopsis still remains elusive. In this study, Arabidopsis with simultaneously overexpressed GmSNAP18 and GmSHMT0 suppressed the growth of underground as well as above-ground parts of plants. Furthermore, Arabidopsis that simultaneously overexpressed GmSNAP18 and GmSHMT08 substantially stimulated BCN susceptibility and remarkably suppressed expression of AtPR1 in the salicylic acid signaling pathway. However, simultaneous overexpression of GmSNAP18 and GmSHMT08 did not impact the expression of AtJAR1 and AtHEL1 in the jasmonic acid and ethylene signaling pathways. GmSNAP18, GmSHMT08, and a pathogenesis-related (PR) protein, GmPR08-Bet VI, in soybean, and AtSNAP2, AtSHMT4, and AtPR1 in Arabidopsis could interact pair-wisely for mediating SCN and BCN resistance in soybean and Arabidopsis, respectively. Both AtSNAP2 and AtPR1 were localized on the plasma membrane, and AtSHMT4 was localized both on the plasma membrane and in the nucleus of cells. Nevertheless, after interactions, AtSNAP2 and AtPR1 could partially translocate into the cell nucleus. GmSNAP18 interacted with AtSHMT4, and GmSHMT4 interacted with AtSNAP2. However, neither GmSNAP18 nor GmSHMT08 interacted with AtPR1. Thus, no pairwise interactions among α-SNAPs, SHMTs, and AtPR1 occurred in Arabidopsis overexpressing either GmSNAP18 or GmSHMT08, or both of them. Transgenic Arabidopsis overexpressing either GmSNAP18 or GmSHMT08 substantially suppressed AtPR1 expression, while transgenic Arabidopsis overexpressing either AtSNAP2 or AtSHMT4 remarkably enhanced it. Taken together, no pairwise interactions of GmSNAP18, GmSHMT08, and AtPR1 with suppressed expression of AtPR1 enhanced BCN susceptibility in Arabidopsis. This study may provide a clue that nematode-resistant or -susceptible functions of plant genes likely depend on both hosts and nematode species.

Loss-of-function of an α-SNAP gene confers resistance to soybean cyst nematode

Plant-parasitic nematodes are one of the most economically impactful pests in agriculture resulting in billions of dollars in realized annual losses worldwide. Soybean cyst nematode (SCN) is the number one biotic constraint on soybean production making it a priority for the discovery, validation and functional characterization of native plant resistance genes and genetic modes of action that can be deployed to improve soybean yield across the globe. Here, we present the discovery and functional characterization of a soybean resistance gene, GmSNAP02. We use unique bi-parental populations to fine-map the precise genomic location, and a combination of whole genome resequencing and gene fragment PCR amplifications to identify and confirm causal haplotypes. Lastly, we validate our candidate gene using CRISPR-Cas9 genome editing and observe a gain of resistance in edited plants. This demonstrates that the GmSNAP02 gene confers a unique mode of resistance to SCN through loss-of-function mutations that implicate GmSNAP02 as a nematode virulence target. We highlight the immediate impact of utilizing GmSNAP02 as a genome-editing-amenable target to diversify nematode resistance in commercially available cultivars.

Cataloging SCN resistance loci in North American public soybean breeding programs

Soybean cyst nematode (SCN) is a destructive pathogen of soybeans responsible for annual yield loss exceeding $1.5 billion in the United States. Here, we conducted a series of genome-wide association studies (GWASs) to understand the genetic landscape of SCN resistance in the University of Missouri soybean breeding programs (Missouri panel), as well as germplasm and cultivars within the United States Department of Agriculture (USDA) Uniform Soybean Tests-Northern Region (NUST). For the Missouri panel, we evaluated the resistance of breeding lines to SCN populations HG 2.5.7 (Race 1), HG 1.2.5.7 (Race 2), HG 0 (Race 3), HG 2.5.7 (Race 5), and HG 1.3.6.7 (Race 14) and identified seven quantitative trait nucleotides (QTNs) associated with SCN resistance on chromosomes 2, 8, 11, 14, 17, and 18. Additionally, we evaluated breeding lines in the NUST panel for resistance to SCN populations HG 2.5.7 (Race 1) and HG 0 (Race 3), and we found three SCN resistance-associated QTNs on chromosomes 7 and 18. Through these analyses, we were able to decipher the impact of seven major genetic loci, including three novel loci, on resistance to several SCN populations and identified candidate genes within each locus. Further, we identified favorable allelic combinations for resistance to individual SCN HG types and provided a list of available germplasm for integration of these unique alleles into soybean breeding programs. Overall, this study offers valuable insight into the landscape of SCN resistance loci in U.S. public soybean breeding programs and provides a framework to develop new and improved soybean cultivars with diverse plant genetic modes of SCN resistance.

Large-scale data mining pipeline for identifying novel soybean genes involved in resistance against the soybean cyst nematode

The soybean cyst nematode (SCN) [Heterodera glycines Ichinohe] is a devastating pathogen of soybean [Glycine max (L.) Merr.] that is rapidly becoming a global economic issue. Two loci conferring SCN resistance have been identified in soybean, Rhg1 and Rhg4; however, they offer declining protection. Therefore, it is imperative that we identify additional mechanisms for SCN resistance. In this paper, we develop a bioinformatics pipeline to identify protein-protein interactions related to SCN resistance by data mining massive-scale datasets. The pipeline combines two leading sequence-based protein-protein interaction predictors, the Protein-protein Interaction Prediction Engine (PIPE), PIPE4, and Scoring PRotein INTeractions (SPRINT) to predict high-confidence interactomes. First, we predicted the top soy interacting protein partners of the Rhg1 and Rhg4 proteins. Both PIPE4 and SPRINT overlap in their predictions with 58 soybean interacting partners, 19 of which had GO terms related to defense. Beginning with the top predicted interactors of Rhg1 and Rhg4, we implement a "guilt by association" in silico proteome-wide approach to identify novel soybean genes that may be involved in SCN resistance. This pipeline identified 1,082 candidate genes whose local interactomes overlap significantly with the Rhg1 and Rhg4 interactomes. Using GO enrichment tools, we highlighted many important genes including five genes with GO terms related to response to the nematode (GO:0009624), namely, Glyma.18G029000, Glyma.11G228300, Glyma.08G120500, Glyma.17G152300, and Glyma.08G265700. This study is the first of its kind to predict interacting partners of known resistance proteins Rhg1 and Rhg4, forming an analysis pipeline that enables researchers to focus their search on high-confidence targets to identify novel SCN resistance genes in soybean.

Twenty years of mining salt tolerance genes in soybean

Current combined challenges of rising food demand, climate change and farmland degradation exert enormous pressure on agricultural production. Worldwide soil salinization, in particular, necessitates the development of salt-tolerant crops. Soybean, being a globally important produce, has its genetic resources increasingly examined to facilitate crop improvement based on functional genomics. In response to the multifaceted physiological challenge that salt stress imposes, soybean has evolved an array of defences against salinity. These include maintaining cell homeostasis by ion transportation, osmoregulation, and restoring oxidative balance. Other adaptations include cell wall alterations, transcriptomic reprogramming, and efficient signal transduction for detecting and responding to salt stress. Here, we reviewed functionally verified genes that underly different salt tolerance mechanisms employed by soybean in the past two decades, and discussed the strategy in selecting salt tolerance genes for crop improvement. Future studies could adopt an integrated multi-omic approach in characterizing soybean salt tolerance adaptations and put our existing knowledge into practice via omic-assisted breeding and gene editing. This review serves as a guide and inspiration for crop developers in enhancing soybean tolerance against abiotic stresses, thereby fulfilling the role of science in solving real-life problems.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01383-3.

Single Nematode Transcriptomic Analysis, Using Long-Read Technology, Reveals Two Novel Virulence Gene Candidates in the Soybean Cyst Nematode, Heterodera glycines

The soybean cyst nematode (Heterodera glycines, SCN), is the most damaging disease of soybean in North America. While management of this pest using resistant soybean is generally still effective, prolonged exposure to cultivars derived from the same source of resistance (PI 88788) has led to the emergence of virulence. Currently, the underlying mechanisms responsible for resistance breakdown remain unknown. In this study, we combined a single nematode transcriptomic profiling approach with long-read sequencing to reannotate the SCN genome. This resulted in the annotation of 1932 novel transcripts and 281 novel gene features. Using a transcript-level quantification approach, we identified eight novel effector candidates overexpressed in PI 88788 virulent nematodes in the late infection stage. Among these were the novel gene Hg-CPZ-1 and a pioneer effector transcript generated through the alternative splicing of the non-effector gene Hetgly21698. While our results demonstrate that alternative splicing in effectors does occur, we found limited evidence of direct involvement in the breakdown of resistance. However, our analysis highlighted a distinct pattern of effector upregulation in response to PI 88788 resistance indicative of a possible adaptation process by SCN to host resistance.

Soybean transporter AAT Rhg1 abundance increases along the nematode migration path and impacts vesiculation and ROS

Rhg1 (Resistance to Heterodera glycines 1) mediates soybean (Glycine max) resistance to soybean cyst nematode (SCN; H. glycines). Rhg1 is a 4-gene, ∼30-kb block that exhibits copy number variation, and the common PI 88788-type rhg1-b haplotype carries 9 to 10 tandem Rhg1 repeats. Glyma.18G022400 (Rhg1-GmAAT), 1 of 3 resistance-conferring genes at the complex Rhg1 locus, encodes the putative amino acid transporter AATRhg1 whose mode of action is largely unknown. We discovered that AATRhg1 protein abundance increases 7- to 15-fold throughout root cells along the migration path of SCN. These root cells develop an increased abundance of vesicles and large vesicle-like bodies (VLB) as well as multivesicular and paramural bodies. AATRhg1 protein is often present in these structures. AATRhg1 abundance remained low in syncytia (plant cells reprogrammed by SCN for feeding), unlike the Rhg1 α-SNAP protein, whose abundance has previously been shown to increase in syncytia. In Nicotiana benthamiana, if soybean AATRhg1 was present, oxidative stress promoted the formation of large VLB, many of which contained AATRhg1. AATRhg1 interacted with the soybean NADPH oxidase GmRBOHG, the ortholog of Arabidopsis thaliana RBOHD previously found to exhibit upregulated expression upon SCN infection. AATRhg1 stimulated reactive oxygen species (ROS) generation when AATRhg1 and GmRBOHG were co-expressed. These findings suggest that AATRhg1 contributes to SCN resistance along the migration path as SCN invades the plant and does so, at least in part, by increasing ROS production. In light of previous findings about α-SNAPRhg1, this study also shows that different Rhg1 resistance proteins function via at least 2 spatially and temporally separate modes of action.

Deciphering the genetic basis of resistance to soybean cyst nematode combining IBD and association mapping

IBD analysis clarified the dynamics of chromosomal recombination during the ZP pedigree breeding process and identified ten genomic regions resistant to SCN race3 combining association mapping. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pathogens for soybean production worldwide. The cultivar Zhongpin03-5373 (ZP), derived from SCN-resistant progenitor parents, Peking, PI 437654 and Huipizhi Heidou, is an elite line with high resistance to SCN race3. In the current study, a pedigree variation map was generated for ZP and its ten progenitors using 3,025,264 high-quality SNPs identified from an average of 16.2 × re-sequencing for each genome. Through identity by decent (IBD) tracking, we showed the dynamic change of genome and detected important IBD fragments, which revealed the comprehensively artificial selection of important traits during ZP breeding process. A total of 2,353 IBD fragments related to SCN resistance including SCN-resistant genes rhg1, rhg4 and NSF_RAN07 were identified based on the resistant-related genetic paths. Moreover, 23 genomic regions underlying resistance to SCN race3 were identified by genome-wide association study (GWAS) in 481 re-sequenced cultivated soybeans. Ten common loci were found by both IBD tracking and GWAS analysis. Haplotype analysis of 16 potential candidate genes suggested a causative SNP (C/T, - 1065) located in the promoter of Glyma.08G096500 and encoding a predicted TIFY5b-related protein on chr8 was highly correlated with SCN race3 resistance. Our results more thoroughly elucidated the dynamics of genomic fragments during ZP pedigree breeding and the genetic basis of SCN resistance, which will provide useful information for gene cloning and the development of resistant soybean cultivars using a marker-assisted selection approach.

Identification of resistant sources from Glycine max against soybean cyst nematode

Soybean cyst nematode (SCN, Heterodera glycines, HG) is one of the severe pests in plant-parasitic nematodes, which impairs root development and causes severe losses in soybean production worldwide. Breeding SCN-resistant cultivars is an important measure for securing harvests without affecting the environment, as can be done with pesticides. The majority of genetic resources for plant pest resistances are found in wild or closely related species which are often difficult to use in breeding due to crossing barriers or close linkage with unfavorable agronomic traits. In this study, 12 soybean cultivars were evaluated for their marker haplotype at the rhg1 and Rhg4 SCN resistance loci and their SCN resistance tested against multiple races in environmentally controlled bioassays. The results showed that all 12 cultivars displayed Peking-type resistance marker haplotypes and all of them proved to be resistant to multiple SCN races. These cultivars provide potential for improving H. glycines resistance of soybean as donor parent in breeding and can contribute to reduce SCN field populations as part of a sustainable agriculture management.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Quantitative disease resistance: Multifaceted players in plant defense

In contrast to large-effect qualitative disease resistance, quantitative disease resistance (QDR) exhibits partial and generally durable resistance and has been extensively utilized in crop breeding. The molecular mechanisms underlying QDR remain largely unknown but considerable progress has been made in this area in recent years. In this review, we summarize the genes that have been associated with plant QDR and their biological functions. Many QDR genes belong to the canonical resistance gene categories with predicted functions in pathogen perception, signal transduction, phytohormone homeostasis, metabolite transport and biosynthesis, and epigenetic regulation. However, other "atypical" QDR genes are predicted to be involved in processes that are not commonly associated with disease resistance, such as vesicle trafficking, molecular chaperones, and others. This diversity of function for QDR genes contrasts with qualitative resistance, which is often based on the actions of nucleotide-binding leucine-rich repeat (NLR) resistance proteins. An understanding of the diversity of QDR mechanisms and of which mechanisms are effective against which classes of pathogens will enable the more effective deployment of QDR to produce more durably resistant, resilient crops.

Characterization of microRNAs in the cyst nematode Heterodera glycines identifies possible candidates involved in cross-kingdom interactions with its host Glycine max

The soybean cyst nematode (SCN - Heterodera glycines) is one of the most damaging pests to the cultivated soybean worldwide. Using a wide array of stylet-secreted effector proteins, this nematode can restructure its host cells into a complex and highly active feeding structure called the syncytium. Tight regulation of these proteins is thought to be essential to the successful formation of this syncytium. To date, multiple mechanisms have been proposed to regulate the expression of these proteins including through post-transcriptional regulation. MicroRNAs (miRNAs) are a class of small, roughly 22-nucleotide-long, non-coding RNA shown to regulate gene expression through its interaction with the 3' untranslated region of genes. These same small RNAs have also been hypothesized to be able to cross over kingdom barriers and regulate genes in other species in a process called cross-kingdom interactions. In this study, we characterized the miRNome of the SCN via sequencing of small-RNAs isolated from whole nematodes and exosomes representing all developmental stages. We identified 121 miRNA loci encoding 96 distinct miRNA families including multiple lineage- and species-specific candidates. Using a combination of plant- and animal-specific miRNA target predictors, we generated a unique repertoire of miRNA:mRNA interacting partners in the nematode and its host plant leading to the identification of a set of nine probable cross-kingdom miRNA candidates.

Effector-Dependent and -Independent Molecular Mechanisms of Soybean-Microbe Interaction

Soybean is a pivotal staple crop worldwide, supplying the main food and feed plant proteins in some countries. In addition to interacting with mutualistic microbes, soybean also needs to protect itself against pathogens. However, to grow inside plant tissues, plant defense mechanisms ranging from passive barriers to induced defense reactions have to be overcome. Pathogenic but also symbiotic micro-organisms effectors can be delivered into the host cell by secretion systems and can interfere with the immunity system and disrupt cellular processes. This review summarizes the latest advances in our understanding of the interaction between secreted effectors and soybean feedback mechanism and uncovers the conserved and special signaling pathway induced by pathogenic soybean cyst nematode, Pseudomonas, Xanthomonas as well as by symbiotic rhizobium.

Comparisons of constitutive resistances to soybean cyst nematode between PI 88788- and Peking-type sources of resistance in soybean by transcriptomic and metabolomic profilings

Soybean cyst nematode (SCN) is a serious damaging disease in soybean worldwide. Peking- and PI 88788-type sources of resistance are two most important germplasm used in breeding resistant soybean cultivars against this disease. However, until now, no comparisons of constitutive resistances to soybean cyst nematode between these two types of sources had been conducted, probably due to the influences of different backgrounds. In this study, we used pooled-sample analysis strategy to minimize the influence of different backgrounds and directly compared the molecular mechanisms underlying constitutive resistance to soybean cyst nematode between these two types of sources via transcriptomic and metabolomic profilings. Six resistant soybean accessions that have identical haplotypes as Peking at Rgh1 and Rhg4 loci were pooled to represent Peking-type sources. The PI88788-type and control pools were also constructed in a same way. Through transcriptomic and metabolomics anaylses, differentially expressed genes and metabolites were identified. The molecular pathways involved in the metabolism of toxic metabolites were predicted to play important roles in conferring soybean cyst nematode resistance to soybean. Functions of two resistant candidate genes were confirmed by hairy roots transformation methods in soybean. Our studies can be helpful for soybean scientists to further learn about the molecular mechanism of resistance to soybean cyst nematode in soybean.

Genome-wide identification of SHMT family genes in cucumber (Cucumis sativus L.) and functional analyses of CsSHMTs in response to hormones and abiotic stresses

Serine hydroxymethyltransferase (SHMT) is a pyridoxal phosphate-dependent enzyme that plays crucial roles in the photorespiration and one-carbon metabolism of plants. In the present research, we conducted a systematic analysis of the SHMT gene family in cucumber (Cucumis sativus L). Results show that a total of 6 SHMT members were identified from the cucumber genome database. CsSHMT1 and CsSHMT2 participate in a fragment duplication event, indicating that CsSHMTs may complete the expansion of family members through fragment duplication. Gene structure analysis found that the number of exons of CsSHMTs ranges from 4 to 15. Members with the same number of exons are classified into the same class in the phylogenetic analysis. Each class reflects its subcellular distribution. Expression and function analysis reveals that CsSHMTs express in a variety of plant tissues, indicating that SHMT gene expression pattern is not organ-specific. qRT-PCR analysis found that CsSHMT3 and CsSHMT5 positively respond to abscisic acid (ABA), and CsSHMT2-6 are induced by indole-3-acetic acid (IAA) and methyl jasmonate (MeJA). Abiotic stress analysis shows that CsSHMT3 is significantly induced by drought and salt stress. These results may provide useful information for further function and evolution analysis of cucumber SHMT genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03378-x.

Opposite Beet Cyst Nematode Infection Phenotypes of Transgenic Arabidopsis Between Overexpressing GmSNAP18 and AtSNAP2 and Between Overexpressing GmSHMT08 and AtSHMT4

The rhg1-a GmSNAP18 (an α-SNAP) and Rhg4 GmSHMT08 are two major cloned genes conferring soybean cyst nematode resistance in Peking-type soybeans, but the application of α-SNAPs and SHMTs in cyst nematode management remains elusive. In this study, GmSNAP18 and GmSHMT08, together with their orthologs in Arabidopsis, AtSNAP2 (an α-SNAP) and AtSHMT4, were individually transformed into Arabidopsis Col-0 to generate the transgenic lines, and the growth of transgenic plants, beet cyst nematode (BCN) infection phenotypes, and AtSNAP2, AtSHMT4, and AtPR1 expression patterns were analyzed using Arabidopsis-BCN compatible interaction system, in addition with protein-protein interaction assay. Pulldown and BiFC assays revealed that GmSNAP18 and GmSHMT08 interacted with AtSHMT4 and AtSNAP2, respectively. Plant root growth was not impacted by overexpression of GmSNAP18 and AtSNAP2. However, overexpression of GmSHMT08 and AtSHMT4 both increased plant height, additionally, overexpression of GmSHMT08 decreased rosette leaf size. Overexpression of GmSNAP18 and GmSHMT08 both suppressed AtPR1 expression and significantly enhanced BCN susceptibility, while overexpression of AtSNAP2 and AtSHMT4 both substantially boosted AtPR1 expression and remarkably enhanced BCN resistance, in transgenic Arabidopsis. Overexpression of GmSNAP18 reduced, while overexpression of AtSNAP2 unaltered AtSHMT4 expression. Overexpression of GmSHMT08 and AtSHMT4 both suppressed AtSNAP2 expression in transgenic Arabidopsis. Thus, different expression patterns of AtPR1 and AtSHMT4 are likely associated with opposite BCN infection phenotypes of Arabidopsis between overexpressing GmSNAP18 and AtSNAP2, and between overexpressing GmSHMT08 and AtSHMT4; and boosted AtPR1 expression are required for enhanced BCN resistance in Arabidopsis. All these results establish a basis for extension of α-SNAPs and SHMTs in cyst nematode management.

Breeding for disease resistance in soybean: a global perspective

KEY MESSAGE

This review provides a comprehensive atlas of QTLs, genes, and alleles conferring resistance to 28 important diseases in all major soybean production regions in the world. Breeding disease-resistant soybean [Glycine max (L.) Merr.] varieties is a common goal for soybean breeding programs to ensure the sustainability and growth of soybean production worldwide. However, due to global climate change, soybean breeders are facing strong challenges to defeat diseases. Marker-assisted selection and genomic selection have been demonstrated to be successful methods in quickly integrating vertical resistance or horizontal resistance into improved soybean varieties, where vertical resistance refers to R genes and major effect QTLs, and horizontal resistance is a combination of major and minor effect genes or QTLs. This review summarized more than 800 resistant loci/alleles and their tightly linked markers for 28 soybean diseases worldwide, caused by nematodes, oomycetes, fungi, bacteria, and viruses. The major breakthroughs in the discovery of disease resistance gene atlas of soybean were also emphasized which include: (1) identification and characterization of vertical resistance genes reside rhg1 and Rhg4 for soybean cyst nematode, and exploration of the underlying regulation mechanisms through copy number variation and (2) map-based cloning and characterization of Rps11 conferring resistance to 80% isolates of Phytophthora sojae across the USA. In this review, we also highlight the validated QTLs in overlapping genomic regions from at least two studies and applied a consistent naming nomenclature for these QTLs. Our review provides a comprehensive summary of important resistant genes/QTLs and can be used as a toolbox for soybean improvement. Finally, the summarized genetic knowledge sheds light on future directions of accelerated soybean breeding and translational genomics studies.

Quantitative trait loci and gene-specific markers associated with resistance to soybean cyst nematode HG type 2.5.7

Soybean cyst nematode (SCN) is one of the most damaging soybean (Glycine max) pests worldwide. More than 95% of SCN-resistant commercial cultivars in North America are derived from a single source of resistance named PI 88788, and the widespread use of this source in the past three decades has led to the selection of virulent biotypes of SCN, such as HG (Heterodera glycines) type 2.5.7 that can overcome the PI 88788-type resistance. The objectives of this study were to identify quantitative trait loci (QTL) and candidate genes underlying the resistance to HG type 2.5.7 isolate and to measure the impact of the resistance factors on seed yield. To achieve the goals, a recombinant inbreed line (RIL) population was established from a cross between an SCN-susceptible high-yielding elite soybean cultivar, OAC Calypso, and the cultivar LD07-3419, resistant to SCN HG type 2.5.7. The RILs resistant to HG type 2.5.7 were identified using greenhouse bioassay tested for differentiation of resistant sources using Kompetitive Allele-Specific PCR (KASP) assay at rhg1 and Rhg4 loci and also for rhg1 copy number variation using TaqMan assay. The RILs were also genotyped using genotype-by-sequencing and three SCN-related QTL were identified on chromosomes 9, 12, and 18 using composite interval mapping. In addition, 31 genes involved in protein kinase activity were identified within QTL regions as potential causal candidate genes underlying the resistance. No significant correlation was found between seed yield and the resistance to SCN in the RIL population evaluated under non-SCN-infested environments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01330-8.

Exploring Soybean Resistance to Soybean Cyst Nematode

Resistance to the soybean cyst nematode (SCN) is a topic incorporating multiple mechanisms and multiple types of science. It is also a topic of substantial agricultural importance, as SCN is estimated to cause more yield damage than any other pathogen of soybean, one of the world's main food crops. Both soybean and SCN have experienced jumps in experimental tractability in the past decade, and significant advances have been made. The rhg1-b locus, deployed on millions of farm acres, has been durable and will remain important, but local SCN populations are gradually evolving to overcome rhg1-b. Multiple other SCN resistance quantitative trait loci (QTL) of proven value are now in play with soybean breeders. QTL causal gene discovery and mechanistic insights into SCN resistance are contributing to both basic and applied disciplines. Additional understanding of SCN and other cyst nematodes will also grow in importance and lead to novel disease control strategies.

Characterization of Heterodera sojae Virulence Phenotypes in Korea

The white soybean cyst nematode Heterodera sojae, isolated from the roots of soybean in Korea, is widespread in most provinces of the country and has the potential to be as harmful to soybean as H. glycines. Determining the virulence phenotypes of H. sojae is essential to devising management strategies that use resistant cultivars. Consequently, virulence phenotypes of 15 H. sojae populations from Korea were determined on seven soybean lines and one susceptible check variety. Two different HS types were found to be present in Korea; the more common HS type 2.5.7, comprising 73.3% of the H. sojae populations and the less common HS type 0, constituting only 26.7% of the tested populations. Considering the high frequency of H. sojae adaptation to soybean indicator lines, the PI 88788 group may not be a possible source of resistance while PI 548402, PI 90763, PI 437654, and PI 89772 can be used as resistance sources for soybean breeding programs aimed at developing H. sojae-resistant soybean cultivars in Korea.

The GmSNAP11 Contributes to Resistance to Soybean Cyst Nematode Race 4 in Glycine max

Soybean cyst nematode (SCN) has devastating effects on soybean production, making it crucial to identify genes conferring SCN resistance. Here we employed next-generation sequencing-based bulked segregant analysis (BSA) to discover genomic regions, candidate genes, and diagnostic markers for resistance to SCN race 4 (SCN4) in soybean. Phenotypic analysis revealed highly significant differences among the reactions of 145 recombinant inbred lines (RILs) to SCN4. In combination with euclidean distance (ED) and Δsingle-nucleotide polymorphism (SNP)-index analyses, we identified a genomic region on Gm11 (designated as rhg1-paralog) associated with SCN4 resistance. Overexpression and RNA interference analyzes of the two candidate genes identified in this region (GmPLAC8 and GmSNAP11) revealed that only GmSNAP11 significantly contributes to SCN4 resistance. We developed a diagnostic marker for GmSNAP11. Using this marker, together with previously developed markers for SCN-resistant loci, rhg1 and Rhg4, we evaluated the relationship between genotypes and SCN4 resistance in 145 RILs and 30 soybean accessions. The results showed that all the SCN4-resistant lines harbored all the three loci, however, some lines harboring the three loci were still susceptible to SCN4. This suggests that these three loci are necessary for the resistance to SCN4, but they alone cannot confer full resistance. The GmSNAP11 and the diagnostic markers developed could be used in genomic-assisted breeding to develop soybean varieties with increased resistance to SCN4.

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.

Breeding a Soybean Cultivar Heinong 531 with Peking-Type Cyst Nematode Resistance, Enhanced Yield, and High Seed-Oil Contents

Soybean cyst nematode (SCN) is a destructive threat to soybean production. It is economically important to develop a new SCN-resistant soybean cultivar with high yield and other good agronomic traits. In this study, a yellow-seed-coated and yellow-hilum-pigmented cultivar Heinong 531 belonging to maturity group I was developed by a pedigree breeding method through a test-cross between a female parental SCN-resistant soybean cultivar Pengdou 158 and a male parental line F1 (high-yield but SCN-susceptible Hefeng 55 × SCN-resistant Kangxian 12). Heinong 531 was evaluated for SCN resistance in both SCN-infested field and autoclaved soil inoculated with hatched second-stage juveniles of SCN HG Type 0. The results indicated that SCN development at all stages in Heinong 531 was suppressed and the female index was only 1.6 to 5.6%. Heinong 531 as well as Pengdou 158 and Kangxian 12 were identified as carrying the Peking-type resistance with both rhg1-a GmSNAP18 and Rhg4 GmSHMT08 genes. In the 2-year regional trials, the average yield of Heinong 531 reached 2805.0 kg/ha, and the 1-year production trial demonstrated an average yield of 2,751.5 kg/ha with yield increase of >12.0% when compared with the local cultivars. The average seed-fat (oil) contents of Heinong 531 reached up to 22.3%. The Peking-type SCN-resistant Heilong 531 with enhanced yield and high seed-oil contents was released in China in June 2021 with the certified number of 'Heishendou 20210004'. These agronomic traits make Heinong 531 a good prospect in a wide attempt to control SCN in the main soybean-producing areas of Northeast China.

Epistatic interaction between Rhg1-a and Rhg2 in PI 90763 confers resistance to virulent soybean cyst nematode populations

KEY MESSAGE

An epistatic interaction between SCN resistance loci rhg1-a and rhg2 in PI 90763 imparts resistance against virulent SCN populations which can be employed to diversify SCN resistance in soybean cultivars. With more than 95% of the $46.1B soybean market dominated by a single type of genetic resistance, breeding for soybean cyst nematode (SCN)-resistant soybean that can effectively combat the widespread increase in virulent SCN populations presents a significant challenge. Rhg genes (for Resistance to Heterodera glycines) play a key role in resistance to SCN; however, their deployment beyond the use of the rhg1-b allele has been limited. In this study, quantitative trait loci (QTL) were mapped using PI 90763 through two biparental F3:4 recombinant inbred line (RIL) populations segregating for rhg1-a and rhg1-b alleles against a SCN HG type 1.2.5.7 (Race 2) population. QTL located on chromosome 18 (rhg1-a) and chromosome 11 (rhg2) were determined to confer SCN resistance in PI 90763. The rhg2 gene was fine-mapped to a 169-Kbp region pinpointing GmSNAP11 as the strongest candidate gene. We demonstrated a unique epistatic interaction between rhg1-a and rhg2 loci that not only confers resistance to multiple virulent SCN populations. Further, we showed that pyramiding rhg2 with the conventional mode of resistance, rhg1-b, is ineffective against these virulent SCN populations. This highlights the importance of pyramiding rhg1-a and rhg2 to maximize the impact of gene pyramiding strategies toward management of SCN populations virulent on rhg1-b sources of resistance. Our results lay the foundation for the next generation of soybean resistance breeding to combat the number one pathogen of soybean.

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.

Identification of Candidate Genes Controlling Soybean Cyst Nematode Resistance in "Handou 10" Based on Genome and Transcriptome Analyzes

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) is a highly destructive pathogen for soybean production worldwide. The use of resistant varieties is the most effective way of preventing yield loss. Handou 10 is a commercial soybean variety with desirable agronomic traits and SCN resistance, however genes underlying the SCN resistance in the variety are unknown. An F_2:8 recombinant inbred line (RIL) population derived from a cross between Zheng 9525 (susceptible) and Handou 10 was developed and its resistance to SCN HG type 2.5.7 (race 1) and 1.2.5.7 (race 2) was identified. We identified seven quantitative trait loci (QTLs) with additive effects. Among these, three QTLs on Chromosomes 7, 8, and 18 were resistant to both races. These QTLs could explain 1.91-7.73% of the phenotypic variation of SCN's female index. The QTLs on chromosomes 8 and 18 have already been reported and were most likely overlapped with rhg1 and Rhg4 loci, respectively. However, the QTL on chromosome 7 was novel. Candidate genes for the three QTLs were predicted through genes functional analysis and transcriptome analysis of infected roots of Handou 10 vs. Zheng 9525. Transcriptome analysis performed also indicated that the plant-pathogen interaction played an important role in the SCN resistance for Handou 10. The information will facilitate SCN-resistant gene cloning, and the novel resistant gene will be a source for improving soybeans' resistance to SCN.

Detection of rare nematode resistance Rhg1 haplotypes in Glycine soja and a novel Rhg1 α-SNAP

This study pursued the hypothesis that wild plant germplasm accessions carrying alleles of interest can be identified using available single nucleotide polymorphism (SNP) genotypes for particular alleles of other (unlinked) genes that contribute to the trait of interest. The soybean cyst nematode (SCN, Heterodera glycines [HG]) resistance locus Rhg1 is widely used in farmed soybean [Glycine max (L.) Merr.]. The two known resistance-conferring haplotypes, rhg1-a and rhg1-b, typically contain three or seven to 10 tandemly duplicated Rhg1 segments, respectively. Each Rhg1 repeat carries four genes, including Glyma.18G022500, which encodes unusual isoforms of the vesicle-trafficking chaperone α-SNAP. Using SoySNP50K data for NSFRAN07 allele presence, we discovered a new Rhg1 haplotype, rhg1-ds, in six accessions of wild soybean, Glycine soja Siebold & Zucc. (0.5% of the ∼1,100 G. soja accessions in the USDA collection). The α-SNAP encoded by rhg1-ds is unique at an important site of amino acid variation and shares with the rhg1-a and rhg1-b α-SNAP proteins the traits of cytotoxicity and altered N-ethylmaleimide sensitive factor (NSF) protein interaction. Copy number assays indicate three repeats of rhg1-ds. G. soja PI 507613 and PI 507623 exhibit resistance to HG type 2.5.7 SCN populations, in part because of contributions from other loci. In a segregating F2 population, rhg1-b and rhg1-ds made statistically indistinguishable contributions to resistance to a partially virulent HG type 2.5.7 SCN population. Hence, the unusual multigene copy number variation Rhg1 haplotype was present but rare in ancestral G. soja and was present in accessions that offer multiple traits for SCN resistance breeding. The accessions were initially identified for study based on an unlinked SNP.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

WI12Rhg1 interacts with DELLAs and mediates soybean cyst nematode resistance through hormone pathways

The soybean cyst nematode (SCN) is one of the most important causes of soybean yield loss. The major source of genetic resistance to SCN is the Rhg1 repeat, a tandem copy number polymorphism of three genes. The roles of these genes are only partially understood. Moreover, nematode populations virulent on Rhg1-carrying soybeans are becoming more common, increasing the need to understand the most successful genetic resistance mechanism. Here, we show that a Rhg1-locus gene (Glyma.18G02270) encoding a wound-inducible protein (WI12Rhg1 ) is needed for SCN resistance. Furthermore, knockout of WI12Rhg1 reduces the expression of DELLA18, and the expression of WI12Rhg1 is itself induced by either JA, SA or GA. The content of the defence hormone SA is significantly lower whilst GA12 and GA53 are increased in WI12Rhg1 knockout roots compared with unedited hairy roots. We find that WI12Rhg1 directly interacts with DELLA18 (Glyma.18G040000) in yeast and plants and that double knockout of DELLA18 and its homeolog DELLA11 (Glyma.11G216500) significantly reduces SCN resistance and alters the root morphology. As DELLA proteins are implicated in hormone signalling, we explored the content of defence hormones (JA and SA) in DELLA knockout and unedited roots, finding reduced levels of JA and SA after the knockout of DELLA. Additionally, the treatment of DELLA-knockout roots with JA or SA rescues SCN resistance lost by the knockout. Meanwhile, the SCN resistance of unedited roots decreases after the treatment with GA, but increases with JA or SA. Our findings highlight the critical roles of WI12Rhg1 and DELLA proteins in SCN resistance through interconnection with hormone signalling.

A Broad Review of Soybean Research on the Ongoing Race to Overcome Soybean Cyst Nematode

Plant pathogens greatly impact food security of the ever-growing human population. Breeding resistant crops is one of the most sustainable strategies to overcome the negative effects of these biotic stressors. In order to efficiently breed for resistant plants, the specific plant-pathogen interactions should be understood. Soybean is a short-day legume that is a staple in human food and animal feed due to its high nutritional content. Soybean cyst nematode (SCN) is a major soybean stressor infecting soybean worldwide including in China, Brazil, Argentina, USA and Canada. There are many Quantitative Trait Loci (QTLs) conferring resistance to SCN that have been identified; however, only two are widely used: rhg1 and Rhg4. Overuse of cultivars containing these QTLs/genes can lead to SCN resistance breakdown, necessitating the use of additional strategies. In this manuscript, a literature review is conducted on research related to soybean resistance to SCN. The main goal is to provide a current understanding of the mechanisms of SCN resistance and list the areas of research that could be further explored.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

A Qa-SNARE complex contributes to soybean cyst nematode resistance via regulation of mitochondria-mediated cell death

The resistance to Heterodera glycines 1 (Rhg1) locus is widely used by soybean breeders to reduce yield loss caused by soybean cyst nematode (SCN). α-SNAP (α-soluble NSF attachment protein) within Rhg1 locus contributes to SCN resistance by modulation of cell status at the SCN feeding site; however, the underlying mechanism is largely unclear. Here, we identified an α-SNAP-interacting protein, GmSYP31A, a Qa-SNARE (soluble NSF attachment protein receptor) protein from soybean. Expression of GmSYP31A significantly induced cell death in Nicotiana benthamiana leaves, and co-expression of α-SNAP and GmSYP31A could accelerate cell death. Overexpression of GmSYP31A increased SCN resistance, while silencing or overexpression of a dominant-negative form of GmSYP31A increased SCN sensitivity. GmSYP31A expression also disrupted endoplasmic reticulum-Golgi trafficking, and the exocytosis pathway. Moreover, α-SNAP was also found to interact with GmVDAC1D (voltage-dependent anion channel). The cytotoxicity induced by the expression of GmSYP31A could be relieved either with the addition of an inhibitor of VDAC protein, or by silencing the VDAC gene. Taken together, our data not only demonstrate that α-SNAP works together with GmSYP31A to increase SCN resistance through triggering cell death, but also highlight the unexplored link between the mitochondrial apoptosis pathway and vesicle trafficking.

Resistance Gene Pyramiding and Rotation to Combat Widespread Soybean Cyst Nematode Virulence

Soybean cyst nematode (SCN) is an important pathogen of soybean causing >$1 billion in yield losses annually in the United States. Planting SCN-resistant soybean cultivars is the primary management strategy. Resistance genes derived from the plant introduction (PI) 88788 (rhg1-b) and PI 548402 (Peking; rhg1-a and Rhg4) are the main types of resistance available in commercial cultivars. The PI 88788 rhg1-b resistance allele is found in the majority of SCN-resistant cultivars in the north central United States. The widespread use of PI 88788 rhg1-b has led to limited options for farmers to rotate resistance sources to manage SCN. Consequently, overreliance on a single type of resistance has resulted in the selection of SCN populations that have adapted to reproduce on these resistant cultivars. Here we evaluated the effectiveness of rotating soybean lines with different combinations of resistance genes to determine the best strategy for combating the widespread increase in virulent SCN and limit future nematode adaptation to resistant cultivars. Eight SCN populations were developed by continuous selection of a virulent SCN field population (Heterodera glycines [HG] type 1.2.5.7) on a single resistance source or in rotation with soybean pyramiding different resistance gene alleles derived from PI 88788 (rhg1-b), PI 437654 (rhg1-a and Rhg4), PI 468916 (cqSCN-006 and cqSCN-007), and PI 567516C (Chr10). SCN population densities were determined for eight generations. HG type tests were conducted after the eighth generation to evaluate population shifts. The continued use of rhg1-b or 006/007 had limited effectiveness for reducing SCN type 1.2.5.7 population density, whereas rotation to the use of rhg1-a/Rhg4 resistance significantly reduced SCN population density but selected for broader SCN virulence (HG type 1.2.3.5.6.7). A rotation of rhg1-a/Rhg4 with a pyramid of rhg1-b/006/007/Chr10 was the most effective combination at both reducing population density and minimizing selection pressure. Our results provide guidance for implementation of a strategic SCN resistance rotation plan to manage the widespread virulence on PI 88788 and sustain the future durability of SCN resistance genes.

Identification of genomic loci conferring broad-spectrum resistance to multiple nematode species in exotic soybean accession PI 567305

KEY MESSAGE

Genetic analysis identified a unique combination of major QTL for resistance to important soybean nematodes concurrently present in a single soybean accession, which has not been reported earlier. An exotic soybean [Glycine max (L.) Merr.] accession, PI 567305, was reported to be highly resistant to three important nematode species, soybean cyst (SCN), root-knot (RKN), and reniform (RN) nematodes. However, genetic basis controlling broad-spectrum resistance in this germplasm has not been investigated. We report results of genetic analysis to identify genomic loci conferring resistance to these nematode species. A bi-parental population consisting of 242 F8-derived recombinant inbred lines (RILs) was developed from a cross of a nematode susceptible cultivar, Magellan, and resistant accession, PI 567305. The RILs were phenotyped for nematode resistance to three SCN HG types. They were genotyped using the Infinium SoySNP6K BeadChips and genotype-by-sequencing (GBS) methods in an attempt to evaluate the cost-effectiveness and efficiency of these two genotyping platforms. Genetic analysis confirmed the major QTL on chromosomes (Chrs) 10 and 18 with broad-spectrum resistance to the three nematodes present in this germplasm. Haplotype and copy number variation analyses of SCN resistance QTL indicated that PI 567305 has a different haplotype, which is associated with likely a unique SCN resistance mechanism different from Peking- or PI 88788-type resistance. The evaluations of both Infinium Beadchip- and GBS-based genotyping technologies provided comprehensive insights for researchers to choose a cost-effective and efficient platform for QTL mapping and for other genomic studies in soybeans.

Dissecting nematode resistance regions in soybean revealed pleiotropic effect of soybean cyst and reniform nematode resistance genes

Reniform nematode (RN, Rotylenchulus reniformis Linford & Oliveira) has emerged as one of the most important plant parasitic nematodes of soybean [Glycine max (L.) Merr.]. Planting resistant varieties is the most effective strategy for nematode management. The objective of this study was to identify quantitative trait loci (QTL) for RN resistance in an exotic soybean line, PI 438489B, using two linkage maps constructed from the Universal Soybean Linkage Panel (USLP 1.0) and next-generation whole-genome resequencing (WGRS) technology. Two QTL controlling RN resistance were identified-the soybean cyst nematode (SCN, Heterodera glycines) resistance gene GmSNAP18 at the rhg1 locus and its paralog GmSNAP11. Strong association between resistant phenotype and haplotypes of the GmSNAP11 and GmSNAP18 was observed. The results indicated that GmSNAP11 possibly could have epistatic effect on GmSNAP18, or vice versa, with the presence of a significant correlation in RN resistance of rhg1-a GmSNAP18 vs. rhg1-b GmSNAP18. Most importantly, our preliminary data suggested that GmSNAP18 and GmSNAP11 proteins physically interact in planta, suggesting that they belong to the same pathway for resistance. Unlike GmSNAP18, no indication of GmSNAP11 copy number variation was found. Moreover, gene-based single nucleotide polymorphism (SNP) markers were developed for rapid detection of RN or SCN resistance at these loci. Our analysis substantiates synergic interaction between GmSNAP11 and GmSNAP18 genes and confirms their roles in RN as well as SCN resistance. These results could contribute to a better understanding of evolution and subfunctionalization of genes conferring resistance to multiple nematode species and provide a framework for further investigations.

Transgressive resistance to Heterodera glycines in chromosome segment substitution lines derived from susceptible soybean parents

Chromosome segment substitution lines (CSSLs) are valuable genetic resources for quantitative trait loci (QTL) mapping of complex agronomic traits especially suitable for minor effect QTL. Here, 162 BC3 F7 -BC7 F3 CSSLs derived from crossing two susceptible parent lines, soybean [Glycine max (L.) Merr.] 'Suinong14' (recurrent parent) × wild soybean (G. soja Siebold & Zucc.) ZYD00006, were used for QTL mapping of soybean cyst nematode (SCN, Heterodera glycine Ichinohe) resistance based on female index (FI) and cysts per gram root (CGR) through phenotypic screening and whole-genome resequencing of CSSLs. Phenotypic results displayed a wide range of distribution and transgressive lines in both HG Type 2.5.7 FI and CGR and demonstrated a higher correlation between CGR and root weight (R2 = .5424) compared with than between FI and CGR (R2 = .0018). Using the single-marker analysis nonparametric mapping test, 33 significant QTL were detected on 18 chromosomes contributing resistance to FI and CGR. Fourteen QTL contributing 5.6-15.5% phenotypic variance (PVE) to FI were revealed on 11 chromosomes, and 16 QTL accounting for 6.1-36.2% PVE in CGR were detected on 14 chromosomes with strong additive effect by multiple-QTL model (MQM) mapping. Twenty-five and 13 out of all 38 QTL identified for FI and CGR on 20 chromosomes were from ZYD00006 and Suinong14, respectively. The CSSLs with the combination of positive alleles for FI, CGR, and root weight exhibited low nematode reproduction. For the first time, QTL associated with CGR have been detected, and both FI and CGR should be considered for breeding purposes in the absence of strong resistance genes such as rhg1 and Rhg4.

Novel resistance strategies to soybean cyst nematode (SCN) in wild soybean

Soybean cyst nematode (SCN, Heterodera glycine Ichinohe) is the most damaging soybean pest worldwide and management of SCN remains challenging. The current SCN resistant soybean cultivars, mainly developed from the cultivated soybean gene pool, are losing resistance due to SCN race shifts. The domestication process and modern breeding practices of soybean cultivars often involve strong selection for desired agronomic traits, and thus, decreased genetic variation in modern cultivars, which consequently resulted in limited sources of SCN resistance. Wild soybean (Glycine soja) is the wild ancestor of cultivated soybean (Glycine max) and it's gene pool is indisputably more diverse than G. max. Our aim is to identify novel resistant genetic resources from wild soybean for the development of new SCN resistant cultivars. In this study, resistance response to HG type 2.5.7 (race 5) of SCN was investigated in a newly identified SCN resistant ecotype, NRS100. To understand the resistance mechanism in this ecotype, we compared RNA seq-based transcriptomes of NRS100 with two SCN-susceptible accessions of G. soja and G. max, as well as an extensively studied SCN resistant cultivar, Peking, under both control and nematode J2-treated conditions. The proposed mechanisms of resistance in NRS100 includes the suppression of the jasmonic acid (JA) signaling pathway in order to allow for salicylic acid (SA) signaling-activated resistance response and polyamine synthesis to promote structural integrity of root cell walls. Our study identifies a set of novel candidate genes and associated pathways involved in SCN resistance and the finding provides insight into the mechanism of SCN resistance in wild soybean, advancing the understanding of resistance and the use of wild soybean-sourced resistance for soybean improvement.

Genetic characterization of qSCN10 from an exotic soybean accession PI 567516C reveals a novel source conferring broad-spectrum resistance to soybean cyst nematode

KEY MESSAGE

The qSCN10 locus with broad-spectrum SCN resistance was fine-mapped to a 379-kb region on chromosome 10 in soybean accession PI 567516C. Candidate genes and potential application benefits of this locus were discussed. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pests of soybean, causing significant yield losses worldwide every year. Genetic resistance has been the major strategy to control this pest. However, the overuse of the same genetic resistance derived primarily from PI 88788 has led to the genetic shifts in nematode populations and resulted in the reduced effectiveness in soybean resistance to SCN. Therefore, novel genetic resistance resources, especially those with broad-spectrum resistance, are needed to develop new resistant cultivars to cope with the genetic shifts in nematode populations. In this study, a quantitative trait locus (QTL) qSCN10 previously identified from a soybean landrace PI 567516C was confirmed to confer resistance to multiple SCN HG Types. This QTL was further fine-mapped to a 379-kb region. There are 51 genes in this region. Four of them are defense-related and were regulated by SCN infection, suggesting their potential role in mediating resistance to SCN. The phylogenetic and haplotype analyses of qSCN10 as well as other information indicate that this locus is different from other reported resistance QTL or genes. There was no yield drag or other unfavorable traits associated with this QTL when near-isogenic lines with and without qSCN10 were tested in a SCN-free field. Therefore, our study not only provides further insight into the genetic basis of soybean resistance to SCN, but also identifies a novel genetic resistance resource for breeding soybean for durable, broad-spectrum resistance to this pest.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Fine-mapping and characterization of qSCN18, a novel QTL controlling soybean cyst nematode resistance in PI 567516C

KEY MESSAGE

The qSCN18 QTL from PI 56756C was confirmed and fine-mapped to improve soybean resistance to the SCN population HG Type 2.5.7 using near-isogenic lines carrying recombination crossovers within the QTL region. The QTL underlying resistance was fine-mapped to a 166-Kbp region on chromosome 18, and the candidate genes were selected based on genomic analyses. Soybean cyst nematode (SCN, Heterodera glycines, Ichinohe) is the most devastating pathogen of soybean. Understanding the genetic basis of SCN resistance is crucial for managing this parasite in the field. Two major loci, rhg1 and Rhg4, were previously characterized as valuable resources for SCN resistance. However, their continuous use has caused shifts in the virulence of SCN populations, which can overcome the resistance conferred by these two major loci. Reduced effectiveness became a major concern in the soybean industry due to continuous use of rhg1 for decades. Thus, it is imperative to identify sources of SCN resistance for durable SCN management. A novel QTL qSCN18 was identified in PI567516C. To fine-map qSCN18 and identify resistance genes, a large backcross population was developed. Nineteen near-isogenic lines (NILs) carrying recombination crossovers within the QTL region were identified. The first phase of fine-mapping narrowed the QTL region to 549-Kbp, whereas the second phase confined the region to 166-Kbp containing 23 genes. Two flanking markers, MK-1 and MK-6, were developed and validated to detect the presence of the qSCN18 resistance allele. Haplotype analysis clustered the fine-mapped qSCN18 region from PI 567516C with the cqSCN-007 locus previously mapped in the wild soybean accession PI 468916. The NILs were developed to further characterize the causal gene(s) harbored in this QTL. This study also confirmed the previously identified qSCN18. The results will facilitate marker-assisted selection (MAS) introducing the qSCN18 locus from PI 567516C into high-yielding soybean cultivars with durable resistance to SCN.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Mapping soybean rhg2 locus, which confers resistance to soybean cyst nematode race 1 in combination with rhg1 and Rhg4 derived from PI 84751

The soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is a devastating pest of soybean (Glycine max (L.) Merr.) in the world. Three soybean QTLs for resistance to SCN race 1 were detected through QTL analyses using recombinant inbred lines (RILs) derived from a cross between 'Tokei 758' (susceptible) and 'To-8E' (resistant to races 1 and 3, derived from 'PI 84751' and 'Gedenshirazu'). Two of the three QTLs appear to be rhg1 and Rhg4 from their locations on the linkage map. The third QTL, detected around Satt359 on chromosome 11, was tentatively identified as rhg2. All RILs resistant to race 1 had all three QTLs. We developed lines carrying the three loci in various combinations, including all and none, from descendants of a cross between 'NIL-SCN' (with resistance derived from 'PI 84751' in the 'Natto-shoryu' background) and 'Natto-shoryu'. Evaluating these lines in a race 1-infected field in Mito, Ibaraki, showed that resistance to race 1 required all three loci. Through field evaluation of 10 recombinant fixed pairs that we developed, we located the rhg2 locus to an 821 kb-region between SSR markers Sat_123 (=WGSP11_0140) and BARCSOYSSR11_1420 on chromosome 11.

Proteome-Wide Analyses Provide New Insights into the Compatible Interaction of Rice with the Root-Knot Nematode Meloidogyne graminicola

The root-knot nematode Meloidogyne graminicola is an important pathogen in rice, causing huge yield losses annually worldwide. Details of the interaction between rice and M. graminicola and the resistance genes in rice still remain unclear. In this study, proteome-wide analyses of the compatible interaction of the japonica rice cultivar “Nipponbare” (NPB) with M. graminicola were performed. In total, 6072 proteins were identified in NPB roots with and without infection of M. graminicola by label-free quantitative mass spectrometry. Of these, 513 specifically or significantly differentially expressed proteins were identified to be uniquely caused by nematode infection. Among these unique proteins, 99 proteins were enriched on seven Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. By comparison of protein expression and gene transcription, LOC_Os01g06600 (ACX, a glutaryl-CoA dehydrogenase), LOC_Os09g23560 (CAD, a cinnamyl-alcohol dehydrogenase), LOC_Os03g39850 (GST, a glutathione S-transferase) and LOC_Os11g11960 (RPM1, a disease resistance protein) on the alpha-linolenic acid metabolism, phenylpropanoid biosynthesis, glutathione metabolism and plant–pathogen interaction pathways, respectively, were all associated with disease defense and identified to be significantly down-regulated in the compatible interaction of NPB with nematodes, while the corresponding genes were remarkably up-regulated in the roots of a resistant rice accession “Khao Pahk Maw” with infection of nematodes. These four genes likely played important roles in the compatible interaction of rice with M. graminicola. Conversely, these disease defense-related genes were hypothesized to be likely involved in the resistance of resistant rice lines to this nematode. The proteome-wide analyses provided many new insights into the interaction of rice with M. graminicola.

A pathogenesis-related protein GmPR08-Bet VI promotes a molecular interaction between the GmSHMT08 and GmSNAP18 in resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is the most devastating pest affecting soybean production worldwide. SCN resistance requires both the GmSHMT08 and the GmSNAP18 in 'Peking'-type resistance. Here, we describe the molecular interaction between GmSHMT08 and GmSNAP18, which is potentiated by a pathogenesis-related protein GmPR08-Bet VI. Like GmSNAP18 and GmSHMT08, GmPR08-Bet VI expression was induced in response to SCN and its overexpression decreased SCN cysts by 65% in infected transgenic soybean roots. Overexpression of GmPR08-Bet VI did not have an effect on SCN resistance when the two cytokinin-binding sites in GmPR08-Bet VI were mutated, indicating a new role of GmPR08-Bet VI in SCN resistance. GmPR08-Bet VI was mapped to a QTL for resistance to SCN using different mapping populations. GmSHMT08, GmSNAP18 and GmPR08-Bet VI localize to the cytosol and plasma membrane. GmSNAP18 expression and localization hyper-accumulated at the plasma membrane and was specific to the root cells surrounding the nematode in SCN-resistant soybeans. Genes encoding key components of the salicylic acid signalling pathway were induced under SCN infection. GmSNAP18 and GmPR08-Bet VI were also induced under salicylic acid and cytokinin exogenous treatments, while GmSHMT08 was induced only when the resistant GmSNAP18 was present, pointing to the presence of a molecular crosstalk between SCN-resistant genes and defence genes. Expression analysis of GmSHMT08 and GmSNAP18 identified the need of a minimum expression requirement to trigger the SCN resistance reaction. These results provide insight into a new response mechanism towards plant nematode resistance involving haplotype compatibility, gene dosage and hormone signalling.

Mutations at the Serine Hydroxymethyltransferase Impact its Interaction with a Soluble NSF Attachment Protein and a Pathogenesis-Related Protein in Soybean

Resistance to soybean cyst nematodes (SCN) in “Peking-type” resistance is bigenic, requiring Rhg4-a and rhg1-a. Rhg4-a encodes a serine hydroxymethyltransferase (GmSHMT08) and rhg1-a encodes a soluble NSF attachment protein (GmSNAP18). Recently, it has been shown that a pathogenesis-related protein, GmPR08-Bet VI, potentiates the interaction between GmSHMT08 and GmSNAP18. Mutational analysis using spontaneously occurring and ethyl methanesulfonate (EMS)-induced mutations was carried out to increase our knowledge of the interacting GmSHMT08/GmSNAP18/GmPR08-Bet VI multi-protein complex. Mutations affecting the GmSHMT08 protein structure (dimerization and tetramerization) and interaction sites with GmSNAP18 and GmPR08-Bet VI proteins were found to impact the multi-protein complex. Interestingly, mutations affecting the PLP/THF substrate binding and catalysis did not affect the multi-protein complex, although they resulted in increased susceptibility to SCN. Most importantly, GmSHMT08 and GmSNAP18 from PI88788 were shown to interact within the cell, being potentiated in the presence of GmPR08-Bet VI. In addition, we have shown the presence of incompatibility between the GmSNAP18 (rhg1-b) of PI88788 and GmSHMT08 (Rhg4-a) from Peking. Components of the reactive oxygen species (ROS) pathway were shown to be induced in the SCN incompatible reaction and were mapped to QTLs for resistance to SCN using different mapping populations.

Haplotype diversity underlying quantitative traits in Canadian soybean breeding germplasm

KEY MESSAGE

Identification of marker-trait associations and trait-associated haplotypes in breeding germplasm identifies regions under selection and highlights changes in haplotype diversity over decades of soybean improvement in Canada. Understanding marker-trait associations using genome-wide association in soybean is typically carried out in diverse germplasm groups where identified loci are often not applicable to soybean breeding efforts. To address this challenge, this study focuses on defining marker-trait associations in breeding germplasm and studying the underlying haplotypes in these regions to assess genetic change through decades of selection. Phenotype data were generated for 175 accessions across multiple environments in Ontario, Canada. A set of 76,549 SNPs were used in the association analysis. A total of 23 genomic regions were identified as significantly associated with yield (5), days to maturity (5), seed oil (3), seed protein (5) and 100-seed weight (5), of which 14 are novel. Each significant region was haplotyped to assess haplotype diversity of the underlying genomic region, identifying ten regions with trait-associated haplotypes in the breeding germplasm. The range of genomic length for these regions (7.2 kb to 6.8 Mb) indicates variation in regional LD for the trait-associated regions. Six of these regions showed changes between eras of breeding, from historical to modern and experimental soybean accessions. Continued selection on these regions may necessitate introgression of novel parental genetic diversity as some haplotypes were fixed within the breeding germplasm. This finding highlights the importance of studying associations and haplotype diversity at a breeding program scale to understand breeders' selections and trends in soybean improvement over time. The haplotypes may also be used as a tool for selection of parental germplasm to inform breeder's decisions on further soybean improvement.

Soybean Resistance Locus Rhg1 Confers Resistance to Multiple Cyst Nematodes in Diverse Plant Species

Cyst nematodes consistently threaten agricultural production, causing billions of dollars in losses globally. The Rhg1 (resistance to Heterodera glycines 1) locus of soybean (Glycine max) is the most popular resistance source used against soybean cyst nematodes (H. glycines). Rhg1 is a complex locus that has multiple repeats of an ≈30-kilobase segment carrying three genes that contribute to resistance. We investigated whether soybean Rhg1 could function in different plant families, conferring resistance to their respective cyst nematode parasites. Transgenic Arabidopsis thaliana and potato (Solanum tuberosum) plants expressing the three soybean Rhg1 genes were generated. The recipient Brassicaceae and Solanaceae plant species exhibited elevated resistance to H. schachtii and Globodera rostochiensis and to G. pallida, respectively. However, some negative consequences including reduced root growth and tuber biomass were observed upon Rhg1 expression in heterologous species. One of the genes at Rhg1 encodes a toxic version of an alpha-SNAP protein that has been demonstrated to interfere with vesicle trafficking. Using a transient expression assay for Nicotiana benthamiana, native Arabidopsis and potato alpha-SNAPs (soluble NSF [N-ethylamine sensitive factor] attachment protein) were found to compensate for the toxicity of soybean Rhg1 alpha-SNAP proteins. Hence, future manipulation of the balance between Rhg1 alpha-SNAP and the endogenous wild-type alpha-SNAPs (as well as the recently discovered soybean NSF-RAN07) may mitigate impacts of Rhg1 on plant productivity. The multispecies efficacy of soybean Rhg1 demonstrates that the encoded mechanisms can function across plant and cyst nematode species and offers a possible avenue for engineered resistance in diverse crop species.

Comparative genomic analyses of two segregating mutants reveal seven genes likely involved in resistance to Fusarium equiseti in soybean via whole genome re-sequencing

KEY MESSAGE

The candidate genes involved in resistance to Fusarium equiseti in soybean PI 437654 were identified through comparative genomic analyses of mutants via whole genome re-sequencing. The fungus Fusarium infects each stage of the growth and development of soybean and causes soybean (Glycine max (L.)) seed and root rot and seedling damping-off and wilt with a large quantity of annual yield loss worldwide. It is very important to identify the resistant genes in soybean to prevent and control this pathogen. One Fusarium equiseti isolate was previously identified to be incompatible with 'PI 437654' but compatible with a Chinese soybean cultivar 'Zhonghuang 13'. In this study, with the infection of this isolate on the seedling roots of developed PI 437654 mutants, 6 mutants were identified from 500 mutants to significantly alter their phenotypes to F. equiseti compared to wild-type PI 437654. Then, two identified segregating mutants were selected to directly perform whole genome re-sequencing. Finally, through comparative genomic analyses 7 genes including one cluster of 4 nucleotide binding site-leucine-rich repeat genes on one genomic region of chromosome 7, a 60S ribosomal protein L12 gene and 2 uncharacterized genes were identified to be likely involved in the resistance to F. equiseti. These genes will facilitate the breeding of resistant germplasm resources and the identification of resistance of soybean to Fusarium spp.