Evolution and selection of Rhg1, a copy-number variant nematode-resistance locus

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.

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.

Involvement of cyclic nucleotide-gated channels in soybean cyst nematode chemotaxis and thermotaxis

The soybean cyst nematode (SCN) is one of the most damaging pests affecting soybean production. SCN displays important host recognition behaviors, such as hatching and infection, by recognizing several compounds produced by the host. Therefore, controlling SCN behaviors such as chemotaxis and thermotaxis is an attractive pest control strategy. In this study, we found that cyclic nucleotide-gated channels (CNG channels) regulate SCN chemotaxis and thermotaxis and Hg-tax-2, a gene encoding a CNG channel, is an important regulator of SCN behavior. Gene silencing of Hg-tax-2 and treatment with a CNG channel inhibitor reduced the attraction of second-stage juveniles to nitrate, an attractant with a different recognition mechanism from the host-derived chemoattractant(s), and to host soybean roots, as well as their avoidance behavior toward high temperatures. Co-treatment of ds Hg-tax-2 with the CNG channel inhibitor indicated that Hg-tax-2 is a major regulator of SCN chemotaxis and thermotaxis. These results suggest new avenues for research on control of SCN.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

Soybean Cyst Nematode Management Is Improved by Combining Native and Transgenic Resistance

Field trials were conducted to assess the benefit of combining a transgenic soybean cyst nematode (SCN) resistance trait, Cry14Ab-1 expressed by the event GMB151, with the native resistance allele rhg1b from PI 88788. The GMB151 event and rhg1b were crossed into common genetic backgrounds and segregated out to create four genetically related lines within each background. The lines created contained both native and transgenic resistance (rhg1b + GMB151), only native resistance (rhg1b alone), only transgenic resistance (GMB151 alone), or neither resistance type (susceptible). The benefit of GMB151 and rhg1b for SCN management was evaluated by measuring SCN control and yield protection. Soybean cyst nematode control was assessed by counting the number of females and cysts on roots early in the season and measuring the change in SCN egg population density over the entire season. The GMB151 transgenic event and the native resistance allele rhg1b both reduced early season SCN reproduction and contributed to significantly higher soybean yield. Compared to susceptible lines, the rhg1b allele improved yield by 33%, while GMB151 improved yield by 13%. Combining the GMB151 event and rhg1b allele resulted in greater SCN control and yield improvement than either provided alone. The combination of GMB151 and rhg1b reduced season-long SCN reproduction by 50% and resulted in 44% greater yield than the susceptible lines. Soybean cyst nematode virulence to rhg1b continues to increase due to the continuous planting of PI 88788-derived resistant cultivars. Pyramiding GMB151 with rhg1b provides a new management option to improve SCN control and soybean yield.

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.

Interactions Among Heterodera glycines, Macrophomina phaseolina, and Soybean Genotype

Heterodera glycines, the soybean cyst nematode (SCN), and fungal pathogen Macrophomina phaseolina are economically important soybean pathogens that may coinfest fields. Resistance remains the most effective management tactic for SCN, and the rhg1-b resistance allele derived from plant introduction 88788 is most commonly deployed in the northern United States. The concomitant effects of SCN and M. phaseolina on soybean performance, as well as the effect of the rhg1-b allele in two different genetic backgrounds, were evaluated in three environments (during 2013 to 2015) and a greenhouse bioassay. Within two soybean populations, half of the lines had the rhg1-b allele, and the other half had the susceptible allele in the backgrounds of the cultivars IA3023 and LD00-3309. Significant interactions between soybean rhg1-b allele and M. phaseolina-infested plots were observed in 2014. In all experiments, initial SCN populations (Pi) and M. phaseolina in roots were associated with reduced soybean yield. SCN reproduction factor (RF = final population/Pi) was affected by SCN Pi, rhg1-b, and genetic background. A background-by-genotype interaction on yield was observed only in 2015, with a stronger rhg1-b effect in the LD00-3309 background, which suggested that the susceptible parent 'IA3023' is tolerant to SCN. SCN female index from greenhouse experiments was compared with field RF, and Lin's concordance and Pearson's correlation coefficients decreased with increasing field SCN Pi in soil. In this study, both SCN and M. phaseolina reduced soybean yield asymptomatically, and the impact of SCN rhg1-b resistance was dependent on SCN virulence but also population density.

Genetic architecture of fresh-market tomato yield

BACKGROUND

The fresh-market tomato (Solanum lycopersicum) is bred for direct consumption and is selected for a high yield of large fruits. To understand the genetic variations (distinct types of DNA sequence polymorphism) that influence the yield, we collected the phenotypic variations in the yields of total fruit, extra-large-sized fruit, small-sized fruit, or red-colored fruit from 68 core inbred contemporary U.S. fresh-market tomatoes for three consecutive years and the genomic information in 8,289,741 single nucleotide polymorphism (SNP) positions from the whole-genome resequencing of these tomatoes.

RESULTS

Genome-wide association (GWA) mapping using the SNP data with or without SNP filtering steps using the regularization methods, validated with quantitative trait loci (QTL) linkage mapping, identified 18 significant association signals for traits evaluated. Among them, 10 of which were not located within genomic regions previously identified as being associated with fruit size/shape. When mapping-driven association signals [558 SNPs associated with 28 yield (component) traits] were used to calculate genomic estimated breeding values (GEBVs) of evaluated traits, the prediction accuracies of the extra-large-sized fruit and small-sized fruit yields were higher than those of the total and red-colored fruit yields, as we tested the generated breeding values in inbred tomatoes and F2 populations. Improved accuracy in GEBV calculation of evaluated traits was achieved by using 364 SNPs identified using the regularization methods.

CONCLUSIONS

Together, these results provide an understanding of the genetic variations underlying the heritable phenotypic variability in yield in contemporary tomato breeding and the information necessary for improving such economically important and complex quantitative trait through breeding.

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.

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.

A mutant allele of the flowering promoting factor 1 gene at the tomato BRACHYTIC locus reduces plant height with high quality fruit

Reduced plant height due to shortened stems is beneficial for improving crop yield potential, better resilience to biotic/abiotic stresses, and rapid crop producer adoption of the agronomic and management practices. Breeding tomato plants with a reduced height, however, poses a particular challenge because this trait is often associated with a significant fruit size (weight) reduction. The tomato BRACHYTIC (BR) locus controls plant height. Genetic mapping and genome assembly revealed three flowering promoting factor 1 (FPF1) genes located within the BR mapping interval, and a complete coding sequence deletion of the telomere proximal FPF1 (Solyc01g066980) was found in the br allele but not in BR. The knock-out of Solyc01g066980 in BR large-fruited fresh-market tomato reduced the height and fruit yield, but the ability to produce large size fruits was retained. However, concurrent yield evaluation of a pair of sister lines with or without the br allele revealed that artificial selection contributes to commercially acceptable yield potential in br tomatoes. A network analysis of gene-expression patterns across genotypes, tissues, and the gibberellic acid (GA) treatment revealed that member(s) of the FPF1 family may play a role in the suppression of the GA biosynthesis in roots and provided a framework for identifying the responsible molecular signaling pathways in br-mediated phenotypic changes. Lastly, mutations of br homologs also resulted in reduced height. These results shed light on the genetic and physiological mechanisms by which the br allele alters tomato architecture.

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.

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.

Homoeologous chromosome exchange explains the creation of a QTL affecting soil-borne pathogen resistance in tobacco

Crop plant partial resistance to plant pathogens controlled by quantitative trait loci (QTL) is desirable in cultivar development programmes because of its increased durability. Mechanisms underlying such resistance are difficult to study. We performed RNA-seq analyses for tobacco (Nicotiana tabacum) nearly isogenic lines (NILs) with and without favourable allele(s) at Phn7.1, a major QTL influencing partial resistance to the soil-borne pathogens Phytophthora nicotianae and Ralstonia solanacearum. Based upon combined analyses of transcriptome-based sequence variation and gene expression profiles, we concluded that allelic variability at the Phn7.1 locus was likely generated from homoeologous exchange, which led to deletion of low-expressing members of the SAR8.2 gene family and duplication of high-expressing SAR8.2 genes from a different subgenome of allotetraploid tobacco. The high expression of endogenous Phn7.1-associated SAR8.2 genes was correlated with observed resistance to P. nicotianae. Our findings suggest a role for genomic rearrangements in the generation of favourable genetic variability affecting resistance to pathogens in plants.

A Bacillus thuringiensis Cry protein controls soybean cyst nematode in transgenic soybean plants

Plant-parasitic nematodes (PPNs) are economically important pests of agricultural crops, and soybean cyst nematode (SCN) in particular is responsible for a large amount of damage to soybean. The need for new solutions for controlling SCN is becoming increasingly urgent, due to the slow decline in effectiveness of the widely used native soybean resistance derived from genetic line PI 88788. Thus, developing transgenic traits for controlling SCN is of great interest. Here, we report a Bacillus thuringiensis delta-endotoxin, Cry14Ab, that controls SCN in transgenic soybean. Experiments in C. elegans suggest the mechanism by which the protein controls nematodes involves damaging the intestine, similar to the mechanism of Cry proteins used to control insects. Plants expressing Cry14Ab show a significant reduction in cyst numbers compared to control plants 30 days after infestation. Field trials also show a reduction in SCN egg counts compared with control plants, demonstrating that this protein has excellent potential to control PPNs in soybean.

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.

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.

Genomic Profiling of Virulence in the Soybean Cyst Nematode Using Single-Nematode Sequencing

Soybean cyst nematode (SCN) is one of the most important diseases in soybean. Currently, the main management strategy relies on planting resistant cultivars. However, the overuse of a single resistance source has led to the selection of virulent SCN populations, although the mechanisms by which the nematode overcomes the resistance genes remain unknown. In this study, we used a nematode-adapted single-cell RNA-seq approach to identify SCN genes potentially involved in resistance breakdown in Peking and PI 88788 parental soybean lines. We established for the first time the full transcriptome of single SCN individuals allowing us to identify a list of putative virulence genes against both major SCN resistance sources. Our analysis identified 48 differentially expressed putative effectors (secreted proteins required for infection) alongside 40 effectors showing evidence of novel structural variants, and 11 effector genes containing phenotype-specific sequence polymorphisms. Additionally, a differential expression analysis revealed an interesting phenomenon of coexpressed gene regions with some containing putative effectors. The selection of virulent SCN individuals on Peking resulted in a profoundly altered transcriptome, especially for genes known to be involved in parasitism. Several sequence polymorphisms were also specific to these virulent nematodes and could potentially play a role in the acquisition of nematode virulence. On the other hand, the transcriptome of virulent individuals on PI 88788 was very similar to avirulent ones with the exception of a few genes, which suggest a distinct virulence strategy to Peking.

t-SNAREs bind the Rhg1 α-SNAP and mediate soybean cyst nematode resistance

Soybean cyst nematode (SCN; Heterodera glycines) is the largest pathogenic cause of soybean yield loss. The Rhg1 locus is the most used and best characterized SCN resistance locus, and contains three genes including one encoding an α-SNAP protein. Although the Rhg1 α-SNAP is known to play an important role in vesicle trafficking and SCN resistance, the protein's binding partners and the molecular mechanisms underpinning SCN resistance remain unclear. In this report, we show that the Rhg1 α-SNAP strongly interacts with two syntaxins of the t-SNARE family (Glyma.12G194800 and Glyma.16G154200) in yeast and plants; importantly, the genes encoding these syntaxins co-localize with SCN resistance quantitative trait loci. Fluorescent visualization revealed that the α-SNAP and the two interacting syntaxins localize to the plasma membrane and perinuclear space in both tobacco epidermal and soybean root cells. The two syntaxins and their two homeologs were mutated, individually and in combination, using the CRISPR-Cas9 system in the SCN-resistant Peking and SCN-susceptible Essex soybean lines. Peking roots with deletions introduced into syntaxin genes exhibited significantly reduced resistance to SCN, confirming that t-SNAREs are critical to resisting SCN infection. The results presented here uncover a key step in the molecular mechanism of SCN resistance, and will be invaluable to soybean breeders aiming to develop highly SCN-resistant soybean varieties.

Dynamics of Population Density and Virulence Phenotype of the Soybean Cyst Nematode as Influenced by Resistance Source Sequence and Tillage

The soybean cyst nematode (SCN), Heterodera glycines, is the most damaging pathogen of soybean. Use of resistant cultivars is an effective strategy to manage SCN, but it also selects for virulent populations over time. A 12-year field experiment was initiated in 2003 to study how tillage and 11 different sequences of four cultivars impact SCN population dynamics and virulence. An SCN-susceptible cultivar and three resistant cultivars (R1, R2, and R3 derived from cultivars PI 88788, Peking, and PI 437654, respectively) were used. Tillage had minimal effect on SCN population density. Compared with no till, conventional tillage resulted in a faster increase of SCN virulence to Peking when the SCN was selected by R2 and virulence to PI 88788 by R3. Among the three SCN-resistant cultivars, R1 supported the greatest population density, R2 supported intermediate population density, and R3 supported the least SCN population density. The SCN populations selected by R1 overcame the resistance in PI 88788 but not in Peking and PI 437654. R2 selected SCN populations that overcame the resistance in Peking but not in PI 88788 and PI 437654. In contrast, R3 selected SCN populations that overcame both PI 88788 and Peking sources of resistance. There was no increase of virulence to PI 437654 in any cultivar sequence. R1 in rotation with R2 or R3 had a negative effect on female index on Peking. Susceptible soybean reduced SCN virulence to Peking, indicating that there was fitness cost of the Peking virulent SCN type. These results suggest that rotation of Peking with PI 88788 is a good strategy for managing the SCN, and susceptible cultivar and no till may reduce SCN virulence selection pressure in some rotations.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

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

KEY MESSAGE

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

Gynoecy instability in cucumber (Cucumis sativus L.) is due to unequal crossover at the copy number variation-dependent Femaleness (F) locus

Cucumber, Cucumis sativus is an important vegetable crop, and gynoecy has played a critical role in yield increase of hybrid cucumber production. Cucumber has a unique genetic system for gynoecious sex expression, which is determined by the copy number variation (CNV)-based, dominant, and dosage-dependent femaleness (F) locus. However, this gynoecy expression system seems unstable since monecious plants could often be found in F-dependent gynoecious cucumber inbreds. We hypothesized that gynoecy instability (gynoecy loss) may be due to unequal crossing over (UCO) during meiosis among repeat units of the CNV. In this study, using high throughput genome resequencing, fiber-FISH and genomic qPCR analyses, we first confirmed and refined the structure of the F locus, which was a CNV of a 30.2-kb tandem repeat. Gynoecious plants contained three genes: CsACS1, CsACS1G, and CsMYB, of which CsACS1G is a duplication of CsACS1 but with a recombinant distal promoter that may contribute to gynoecy sex expression. In two large populations from self-pollinated gynoecious inbred lines, 'gynoecy loss' mutants were identified with similar mutation rates (~0.12%). We show that these monecious mutants have lost CsACS1G. In addition, we identified gynoecious lines in natural populations that carry two copies of CSACS1G. We proposed a model to explain gynoecy instability in F-dependent cucumbers, which is caused by UCO among CSACS1/G units during meiosis. The findings present a convincing case that the phenotypic variation of an economically important trait is associated with the dynamic changes of copy numbers at the F locus. This work also has important implications in cucumber breeding.

Gene presence-absence variation associates with quantitative Verticillium longisporum disease resistance in Brassica napus

Although copy number variation (CNV) and presence-absence variation (PAV) have been discovered in selected gene families in most crop species, the global prevalence of these polymorphisms in most complex genomes is still unclear and their influence on quantitatively inherited agronomic traits is still largely unknown. Here we analyze the association of gene PAV with resistance of oilseed rape (Brassica napus) against the important fungal pathogen Verticillium longisporum, as an example for a complex, quantitative disease resistance in the strongly rearranged genome of a recent allopolyploid crop species. Using Single Nucleotide absence Polymorphism (SNaP) markers to efficiently trace PAV in breeding populations, we significantly increased the resolution of loci influencing V. longisporum resistance in biparental and multi-parental mapping populations. Gene PAV, assayed by resequencing mapping parents, was observed in 23-51% of the genes within confidence intervals of quantitative trait loci (QTL) for V. longisporum resistance, and high-priority candidate genes identified within QTL were all affected by PAV. The results demonstrate the prominent role of gene PAV in determining agronomic traits, suggesting that this important class of polymorphism should be exploited more systematically in future plant breeding.

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.

Molecular Basis of Soybean Resistance to Soybean Aphids and Soybean Cyst Nematodes

Soybean aphid (SBA; Aphis glycines Matsumura) and soybean cyst nematode (SCN; Heterodera glycines Ichninohe) are major pests of the soybean (Glycine max [L.] Merr.). Substantial progress has been made in identifying the genetic basis of limiting these pests in both model and non-model plant systems. Classical linkage mapping and genome-wide association studies (GWAS) have identified major and minor quantitative trait loci (QTLs) in soybean. Studies on interactions of SBA and SCN effectors with host proteins have identified molecular cues in various signaling pathways, including those involved in plant disease resistance and phytohormone regulations. In this paper, we review the molecular basis of soybean resistance to SBA and SCN, and we provide a synthesis of recent studies of soybean QTLs/genes that could mitigate the effects of virulent SBA and SCN populations. We also review relevant studies of aphid–nematode interactions, particularly in the soybean–SBA–SCN system.

Effective identification of soybean candidate genes involved in resistance to soybean cyst nematode via direct whole genome re-sequencing of two segregating mutants

KEY MESSAGE

Three soybean candidate genes involved in resistance to soybean cyst nematode race 4 were identified via direct whole genome re-sequencing of two segregating mutants. The genes conferring resistance to soybean cyst nematode (SCN) race 4 (Hg type 1.2.3.5.7) in soybean (Glycine max L. Merr.) remains unknown. Next generation sequencing-based methods identify a wide range of targets, it is difficult to identify genes underlying traits. Use of the MutMap and QTL-seq methods to identify trait candidate genes needs backcrossing and is very time-consuming. Here we report a simple method to effectively identify candidate genes involved in resistance to SCN race 4. Two ethane methylsulfonate mutagenized mutants of soybean 'PI 437654', whose SCN race 4-infection phenotype altered, were selected. Six relevant whole genomes were re-sequenced, and then calling of genomic variants (SNPs and InDels) was conducted and compared to 'Williams 82'. The comparison eliminated many genomic variants from the mutant lines that overlapped two non-phenotypic but mutant progeny plants, wild-type PI 437654 and 'Zhonghuang 13'. Finally, only 27 mutations were found among 10 genes. Of these 10 genes, 3 genes, Glyma.09g054000, Glyma.16g065700 and Glyma.18g192200 were overlapped between two phenotypic mutant progeny plants. Therefore, the three genes may be the candidate genes involved in resistance of PI 437654 to soybean cyst nematode race 4. This method simplifies the effective identification of candidate genes.

The rhg1-a (Rhg1 low-copy) nematode resistance source harbors a copia-family retrotransposon within the Rhg1-encoded α-SNAP gene

Soybean growers widely use the Resistance to Heterodera glycines 1 (Rhg1) locus to reduce yield losses caused by soybean cyst nematode (SCN). Rhg1 is a tandemly repeated four gene block. Two classes of SCN resistance-conferring Rhg1 haplotypes are recognized: rhg1-a ("Peking-type," low-copy number, three or fewer Rhg1 repeats) and rhg1-b ("PI 88788-type," high-copy number, four or more Rhg1 repeats). The rhg1-a and rhg1-b haplotypes encode α-SNAP (alpha-Soluble NSF Attachment Protein) variants α-SNAP Rhg1 LC and α-SNAP Rhg1 HC, respectively, with differing atypical C-terminal domains, that contribute to SCN resistance. Here we report that rhg1-a soybean accessions harbor a copia retrotransposon within their Rhg1 Glyma.18G022500 (α-SNAP-encoding) gene. We termed this retrotransposon "RAC," for Rhg1 alpha-SNAP copia. Soybean carries multiple RAC-like retrotransposon sequences. The Rhg1 RAC insertion is in the Glyma.18G022500 genes of all true rhg1-a haplotypes we tested and was not detected in any examined rhg1-b or Rhg1WT (single-copy) soybeans. RAC is an intact element residing within intron 1, anti-sense to the rhg1-a α-SNAP open reading frame. RAC has intrinsic promoter activities, but overt impacts of RAC on transgenic α-SNAP Rhg1 LC mRNA and protein abundance were not detected. From the native rhg1-a RAC+ genomic context, elevated α-SNAP Rhg1 LC protein abundance was observed in syncytium cells, as was previously observed for α-SNAP Rhg1 HC (whose rhg1-b does not carry RAC). Using a SoySNP50K SNP corresponding with RAC presence, just ~42% of USDA accessions bearing previously identified rhg1-a SoySNP50K SNP signatures harbor the RAC insertion. Subsequent analysis of several of these putative rhg1-a accessions lacking RAC revealed that none encoded α-SNAPRhg1LC, and thus, they are not rhg1-a. rhg1-a haplotypes are of rising interest, with Rhg4, for combating SCN populations that exhibit increased virulence against the widely used rhg1-b resistance. The present study reveals another unexpected structural feature of many Rhg1 loci, and a selectable feature that is predictive of rhg1-a haplotypes.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.

Genome wide association study discovers genomic regions involved in resistance to soybean cyst nematode (Heterodera glycines) in common bean

Common bean (Phaseolus vulgaris L.) is an important high protein crop grown worldwide. North Dakota and Minnesota are the largest producers of common beans in the USA, but crop production is threatened by soybean cyst nematode (SCN; Heterodera glycines Ichinohe) because most current cultivars are susceptible. Greenhouse screening data using SCN HG type 0 from 317 plant introductions (PI's) from the USDA core collection was used to conduct a genome wide association study (GWAS). These lines were divided into two subpopulations based on principal component analysis (Middle American vs. Andean). Phenotypic results based on the female index showed that accessions could be classified as highly resistant (21% and 27%), moderately resistant (51% and 48%), moderately susceptible (27% and 22%) and highly susceptible (1% and 3%) for Middle American and Andean gene pools, respectively. Mixed models with two principal components (PCs) and kinship matrix for Middle American genotypes and Andean genotypes were used in the GWAS analysis using 3,985 and 4,811 single nucleotide polymorphic (SNP) markers, respectively which were evenly distributed across all 11 chromosomes. Significant peaks on Pv07, and Pv11 in Middle American and on Pv07, Pv08, Pv09 and Pv11 in Andean group were found to be associated with SCN resistance. Homologs of soybean rhg1, a locus which confers resistance to SCN in soybean, were identified on chromosomes Pv01 and Pv08 in the Middle American and Andean gene pools, respectively. These genomic regions may be the key to develop SCN-resistant common bean cultivars.

The soybean Rhg1 amino acid transporter gene alters glutamate homeostasis and jasmonic acid-induced resistance to soybean cyst nematode

Rhg1 (resistance to Heterodera glycines 1) is an important locus that contributes to resistance against soybean cyst nematode (SCN; Heterodera glycines Ichinohe), which is the most economically damaging disease of soybean worldwide. Simultaneous overexpression of three genes encoding a predicted amino acid transporter, an α-soluble N-ethylmaleimide-sensitive factor attachment protein (α-SNAP) and a predicted wound-induced protein resulted in resistance to SCN provided by this locus. However, the roles of two of these genes (excluding α-SNAP) remain unknown. Here, we report the functional characterization of Glyma.18G022400, a gene at the Rhg1 locus that encodes the predicted amino acid transporter Rhg1-GmAAT. Although the direct role of Rhg1-GmAAT in glutamate transport was not demonstrated, multiple lines of evidence showed that Rhg1-GmAAT impacts glutamic acid tolerance and glutamate transportation in soybean. Transcriptomic and metabolite profiling indicated that overexpression of Rhg1-GmAAT activated the jasmonic acid (JA) pathway. Treatment with a JA biosynthesis inhibitor reduced the resistance provided by the Rhg1-containing PI88788 to SCN, which suggested that the JA pathway might play a role in Rhg1-mediated resistance to SCN. Our results could be helpful for the clarification of the mechanism of resistance to SCN provided by Rhg1 in soybean.

Survey of Heterodera glycines Population Densities and Virulence Phenotypes During 2015-2016 in Missouri

The soybean cyst nematode (SCN), Heterodera glycines, is one of the most important pathogens of soybean. Periodic monitoring of SCN population densities and virulence phenotypes is necessary for developing management strategies utilizing resistant cultivars, the primary strategy used to combat this pest. Therefore, we conducted a statewide survey of Missouri to determine SCN population densities and virulence phenotypes during 2015-2016 and compared these results with a similar survey conducted in 2005. SCN population densities were determined for 393 soil samples representing 74 soybean-producing counties across eight geographical regions of Missouri. Eighty-eight percent of samples tested positive for SCN, up from 50% in 2005, and population densities ranged from 125 to 99,000 eggs per 250 cm3 of soil. The virulence phenotypes of 48 SCN populations also were determined. For this, female indices (FI) were calculated by dividing the mean number of females that develop on the roots of a set of resistant soybean indicator lines by the mean number of females that develop on the roots of susceptible cultivar Lee74 after 30 days in the greenhouse then multiplying by 100 to obtain a percentage. Notably, all SCN populations evaluated during 2015-2016 had a FI > 10 on PI 88788, the most widely used source of resistance in Missouri, in contrast to 78% in 2005. Moreover, 50% of these populations had a FI > 50 on PI 88788, up from 16% in 2005. Forty-three percent of populations tested also had a FI > 10 on Peking, the second most common source of resistance by farmers. Our results show that over the last decade, SCN has become more prevalent in Missouri fields. Additionally, the percentage of individuals within SCN field populations that are virulent on PI 88788 and Peking has markedly increased. The results stress the importance of rotating cultivars with different types of resistance when using resistant cultivars to manage SCN.

An atypical N-ethylmaleimide sensitive factor enables the viability of nematode-resistant Rhg1 soybeans

N-ethylmaleimide sensitive factor (NSF) and α-soluble NSF attachment protein (α-SNAP) are essential eukaryotic housekeeping proteins that cooperatively function to sustain vesicular trafficking. The "resistance to Heterodera glycines 1" (Rhg1) locus of soybean (Glycine max) confers resistance to soybean cyst nematode, a highly damaging soybean pest. Rhg1 loci encode repeat copies of atypical α-SNAP proteins that are defective in promoting NSF function and are cytotoxic in certain contexts. Here, we discovered an unusual NSF allele (Rhg1-associated NSF on chromosome 07; NSFRAN07 ) in Rhg1+ germplasm. NSFRAN07 protein modeling to mammalian NSF/α-SNAP complex structures indicated that at least three of the five NSFRAN07 polymorphisms reside adjacent to the α-SNAP binding interface. NSFRAN07 exhibited stronger in vitro binding with Rhg1 resistance-type α-SNAPs. NSFRAN07 coexpression in planta was more protective against Rhg1 α-SNAP cytotoxicity, relative to WT NSFCh07 Investigation of a previously reported segregation distortion between chromosome 18 Rhg1 and a chromosome 07 interval now known to contain the Glyma.07G195900 NSF gene revealed 100% coinheritance of the NSFRAN07 allele with disease resistance Rhg1 alleles, across 855 soybean accessions and in all examined Rhg1+ progeny from biparental crosses. Additionally, we show that some Rhg1-mediated resistance is associated with depletion of WT α-SNAP abundance via selective loss of WT α-SNAP loci. Hence atypical coevolution of the soybean SNARE-recycling machinery has balanced the acquisition of an otherwise disruptive housekeeping protein, enabling a valuable disease resistance trait. Our findings further indicate that successful engineering of Rhg1-related resistance in plants will require a compatible NSF partner for the resistance-conferring α-SNAP.

Copy number variation and disease resistance in plants

Plant genome diversity varies from single nucleotide polymorphisms to large-scale deletions, insertions, duplications, or re-arrangements. These re-arrangements of sequences resulting from duplication, gains or losses of DNA segments are termed copy number variations (CNVs). During the last decade, numerous studies have emphasized the importance of CNVs as a factor affecting human phenotype; in particular, CNVs have been associated with risks for several severe diseases. In plants, the exploration of the extent and role of CNVs in resistance against pathogens and pests is just beginning. Since CNVs are likely to be associated with disease resistance in plants, an understanding of the distribution of CNVs could assist in the identification of novel plant disease-resistance genes. In this paper, we review existing information about CNVs; their importance, role and function, as well as their association with disease resistance in plants.

Systematic Mutagenesis of Serine Hydroxymethyltransferase Reveals an Essential Role in Nematode Resistance

Rhg4 is a major genetic locus that contributes to soybean cyst nematode (SCN) resistance in the Peking-type resistance of soybean (Glycine max), which also requires the rhg1 gene. By map-based cloning and functional genomic approaches, we previously showed that the Rhg4 gene encodes a predicted cytosolic serine hydroxymethyltransferase (GmSHMT08); however, the novel gain of function of GmSHMT08 in SCN resistance remains to be characterized. Using a forward genetic screen, we identified an allelic series of GmSHMT08 mutants that shed new light on the mechanistic aspects of GmSHMT08-mediated resistance. The new mutants provide compelling genetic evidence that Peking-type rhg1 resistance in cv Forrest is fully dependent on the GmSHMT08 gene and demonstrates that this resistance is mechanistically different from the PI 88788-type of resistance that only requires rhg1 We also demonstrated that rhg1-a from cv Forrest, although required, does not exert selection pressure on the nematode to shift from HG type 7, which further validates the bigenic nature of this resistance. Mapping of the identified mutations onto the SHMT structural model uncovered key residues for structural stability, ligand binding, enzyme activity, and protein interactions, suggesting that GmSHMT08 has additional functions aside from its main enzymatic role in SCN resistance. Lastly, we demonstrate the functionality of the GmSHMT08 SCN resistance gene in a transgenic soybean plant.

Soybean domestication: the origin, genetic architecture and molecular bases

Domestication provides an important model for the study of evolution, and information learned from domestication research aids in the continued improvement of crop species. Recent progress in de novo assembly and whole-genome resequencing of wild and cultivated soybean genomes, in addition to new archeological discoveries, sheds light on the origin of this important crop and provides a clearer view on the modes of artificial selection that drove soybean domestication and diversification. This novel genomic information enables the search for polymorphisms that underlie variation in agronomic traits and highlights genes that exhibit a signature of selection, leading to the identification of a number of candidate genes that may have played important roles in soybean domestication, diversification and improvement. These discoveries provide a novel point of comparison on the evolutionary bases of important agronomic traits among different crop species.

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

Two types of resistant soybean (Glycine max (L.) Merr.) sources are widely used against soybean cyst nematode (SCN, Heterodera glycines Ichinohe). These include Peking-type soybean, whose resistance requires both the rhg1-a and Rhg4 alleles, and PI 88788-type soybean, whose resistance requires only the rhg1-b allele. Multiple copy number of PI 88788-type GmSNAP18, GmAAT, and GmWI12 in one genomic segment simultaneously contribute to rhg1-b resistance. Using an integrated set of genetic and genomic approaches, we demonstrate that the rhg1-a Peking-type GmSNAP18 is sufficient for resistance to SCN in combination with Rhg4. The two SNAPs (soluble NSF attachment proteins) differ by only five amino acids. Our findings suggest that Peking-type GmSNAP18 is performing a different role in SCN resistance than PI 88788-type GmSNAP18. As such, this is an example of a pathogen resistance gene that has evolved to underlie two types of resistance, yet ensure the same function within a single plant species.

Genetics and Adaptation of Soybean Cyst Nematode to Broad Spectrum Soybean Resistance

The soybean cyst nematode (SCN) Heterodera glycines is a major threat to soybean production, made more challenging by the current limitations of natural resistance for managing this pathogen. The use of resistant host cultivars is effective, but, over time, results in the generation of virulent nematode populations able to robustly parasitize the resistant host. In order to understand how virulence develops in SCN, we utilized a single backcross BC1F2 strategy to mate a highly virulent inbred population (TN20), capable of reproducing on all current sources of resistance, with an avirulent one (PA3), unable to reproduce on any of the resistant soybean lines. The offspring were then investigated to determine how virulence is inherited on the main sources of SCN resistance, derived from soybean lines Peking, PI 88788, PI 90763, and the broad spectrum resistance source PI 437654. Significantly, our results suggest virulence on PI 437654 is a multigenic recessive trait that allows the nematode to reproduce on all current sources of resistance. In addition, we examined how virulence on different sources of resistance interact by placing virulent SCN populations under secondary selection, and identified a strong counter-selection between virulence on PI 88788- and PI 90763-derived resistances, while no such counter-selection existed between virulence on Peking and PI 88788 resistance sources. Our results suggest that the genes responsible for virulence on PI 88788 and PI 90763 may be different alleles at a common locus. If so, rotation of cultivars with resistance from these two sources may be an effective management protocol.

Advancements in breeding, genetics, and genomics for resistance to three nematode species in soybean

KEY MESSAGE

Integration of genetic analysis, molecular biology, and genomic approaches drastically enhanced our understanding of genetic control of nematode resistance and provided effective breeding strategies in soybeans. Three nematode species, including soybean cyst (SCN, Heterodera glycine), root-knot (RKN, Meloidogyne incognita), and reniform (RN, Rotylenchulus reniformis), are the most destructive pests and have spread to soybean growing areas worldwide. Host plant resistance has played an important role in their control. This review focuses on genetic, genomic studies, and breeding efforts over the past two decades to identify and improve host resistance to these three nematode species. Advancements in genetics, genomics, and bioinformatics have improved our understanding of the molecular and genetic mechanisms of nematode resistance and enabled researchers to generate large-scale genomic resources and marker-trait associations. Whole-genome resequencing, genotyping-by-sequencing, genome-wide association studies, and haplotype analyses have been employed to map and dissect genomic locations for nematode resistance. Recently, two major SCN-resistant loci, Rhg1 and Rhg4, were cloned and other novel resistance quantitative trait loci (QTL) have been discovered. Based on these discoveries, gene-specific DNA markers have been developed for both Rhg1 and Rhg4 loci, which were useful for marker-assisted selection. With RKN resistance QTL being mapped, candidate genes responsible for RKN resistance were identified, leading to the development of functional single nucleotide polymorphism markers. So far, three resistances QTL have been genetically mapped for RN resistance. With nematode species overcoming the host plant resistance, continuous efforts in the identification and deployment of new resistance genes are required to support the development of soybean cultivars with multiple and durable resistance to these pests.

Soybean Resistance to the Soybean Cyst Nematode Heterodera glycines: An Update

The soybean cyst nematode (SCN), Heterodera glycines, remains a serious threat to soybean production throughout the world. A lack of genetic diversity in resistant soybean cultivars has led to a widespread shift toward virulence in SCN populations, leaving farmers with few proven options other than nonhost rotation to manage this nematode. Recent advances in our understanding of the genes controlling resistance to the nematode have led to improved molecular markers, which are, in turn, increasing the efficiency and precision of the breeding pipeline. A better understanding of the molecular and biochemical basis of SCN resistance and nematode virulence will provide information useful for the development of a long-term strategic plan for diversification and the deployment of cultivars that protect current sources of natural resistance while identifying new targets for engineering novel resistance.

Impact of Rhg1 copy number, type, and interaction with Rhg4 on resistance to Heterodera glycines in soybean

KEY MESSAGE

Evaluations of soybean populations showed that both Rhg1 copy number and type were important in determining soybean cyst nematode resistance with higher copy number within Rhg1 type conferring greater resistance. Rhg1 and Rhg4 are important loci conferring resistance to soybean cyst nematode (SCN; Heterodera glycines). Alleles at Rhg1 have been shown to vary for copy number and type and the importance of this variation in conferring resistance is not well defined. The repeat number ranges from one to 10 and there are three variant repeat sequence types [plant introduction (PI) 88788-'Fayette' type (F), 'Peking' type (P) and Williams 82 type (W)] across diverse soybean germplasm. We developed populations segregating for Rhg1 copy number and type and Rhg4 allele type to investigate the effect of these factors and their interaction on SCN resistance. F2 plants from each cross were evaluated for the segregation of Rhg1 and Rhg4 alleles and for SCN reproduction after infesting plants with HG type 2.5.7 and HG type 7 populations. Within repeat types, an increase in repeat number was associated with greater resistance. The P type Rhg1 showed an advantage over F + W type for SCN population HG type 2.5.7 but this was not observed for SCN HG type 7. While plants with P type Rhg1 required Rhg4 to achieve full resistance, Rhg4 did not increase resistance in the background of F + W type Rhg1 repeat. This study demonstrates the importance of both Rhg1 copy number and type in determining resistance and can assist soybean breeders in determining what alleles would best fit their breeding goals.

Disease resistance through impairment of α-SNAP-NSF interaction and vesicular trafficking by soybean Rhg1

α-SNAP [soluble NSF (N-ethylmaleimide-sensitive factor) attachment protein] and NSF proteins are conserved across eukaryotes and sustain cellular vesicle trafficking by mediating disassembly and reuse of SNARE protein complexes, which facilitate fusion of vesicles to target membranes. However, certain haplotypes of the Rhg1 (resistance to Heterodera glycines 1) locus of soybean possess multiple repeat copies of an α-SNAP gene (Glyma.18G022500) that encodes atypical amino acids at a highly conserved functional site. These Rhg1 loci mediate resistance to soybean cyst nematode (SCN; H. glycines), the most economically damaging pathogen of soybeans worldwide. Rhg1 is widely used in agriculture, but the mechanisms of Rhg1 disease resistance have remained unclear. In the present study, we found that the resistance-type Rhg1 α-SNAP is defective in interaction with NSF. Elevated in planta expression of resistance-type Rhg1 α-SNAPs depleted the abundance of SNARE-recycling 20S complexes, disrupted vesicle trafficking, induced elevated abundance of NSF, and caused cytotoxicity. Soybean, due to ancient genome duplication events, carries other loci that encode canonical (wild-type) α-SNAPs. Expression of these α-SNAPs counteracted the cytotoxicity of resistance-type Rhg1 α-SNAPs. For successful growth and reproduction, SCN dramatically reprograms a set of plant root cells and must sustain this sedentary feeding site for 2-4 weeks. Immunoblots and electron microscopy immunolocalization revealed that resistance-type α-SNAPs specifically hyperaccumulate relative to wild-type α-SNAPs at the nematode feeding site, promoting the demise of this biotrophic interface. The paradigm of disease resistance through a dysfunctional variant of an essential gene may be applicable to other plant-pathogen interactions.

An efficient method for measuring copy number variation applied to improvement of nematode resistance in soybean

Copy number variation (CNV) is implicated in important traits in multiple crop plants, but can be challenging to genotype using conventional methods. The Rhg1 locus of soybean, which confers resistance to soybean cyst nematode (SCN), is a CNV of multiple 31.2-kb genomic units each containing four genes. Reliable, high-throughput methods to quantify Rhg1 and other CNVs for selective breeding were developed. The CNV genotyping assay described here uses a homeologous gene copy within the paleopolyploid soybean genome to provide the internal control for a single-tube TaqMan copy number assay. Using this assay, CNV in breeding populations can be tracked with high precision. We also show that extensive CNV exists within Fayette, a released, inbred SCN-resistant soybean cultivar with a high copy number at Rhg1 derived from a single donor parent. Copy number at Rhg1 is therefore unstable within a released variety over a relatively small number of generations. Using this assay to select for individuals with altered copy number, plants were obtained with both increased copy number and increased SCN resistance relative to control plants. Thus, CNV genotyping technologies can be used as a new type of marker-assisted selection to select for desirable traits in breeding populations, and to control for undesirable variation within cultivars.

Identification of haplotypes at the Rsv4 genomic region in soybean associated with durable resistance to soybean mosaic virus

KEY MESSAGE

Discovery of new germplasm sources and identification of haplotypes for the durable Soybean mosaic virus resistance gene, Rsv 4, provide novel resources for map-based cloning and genetic improvement efforts in soybean. The Soybean mosaic virus (SMV) resistance locus Rsv4 is of interest because it provides a durable type of resistance in soybean [Glycine max (L.) Merr.]. To better understand its molecular basis, we used a population of 309 BC3F2 individuals to fine-map Rsv4 to a ~120 kb interval and leveraged this genetic information in a second study to identify accessions 'Haman' and 'Ilpumgeomjeong' as new sources of Rsv4. These two accessions along with three other Rsv4 and 14 rsv4 accessions were used to examine the patterns of nucleotide diversity at the Rsv4 region based on high-depth resequencing data. Through a targeted association analysis of these 19 accessions within the ~120 kb interval, a cluster of four intergenic single-nucleotide polymorphisms (SNPs) was found to perfectly associate with SMV resistance. Interestingly, this ~120 kb interval did not contain any genes similar to previously characterized dominant disease resistance genes. Therefore, a haplotype analysis was used to further resolve the association signal to a ~94 kb region, which also resulted in the identification of at least two Rsv4 haplotypes. A haplotype phylogenetic analysis of this region suggests that the Rsv4 locus in G. max is recently introgressed from G. soja. This integrated study provides a strong foundation for efforts focused on the cloning of this durable virus resistance gene and marker-assisted selection of Rsv4-mediated SMV resistance in soybean breeding programs.

Genomic-assisted phylogenetic analysis and marker development for next generation soybean cyst nematode resistance breeding

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is a serious soybean pest. The use of resistant cultivars is an effective approach for preventing yield loss. In this study, 19,652 publicly available soybean accessions that were previously genotyped with the SoySNP50K iSelect BeadChip were used to evaluate the phylogenetic diversity of SCN resistance genes Rhg1 and Rhg4 in an attempt to identify novel sources of resistance. The sequence information of soybean lines was utilized to develop KASPar (KBioscience Competitive Allele-Specific PCR) assays from single nucleotide polymorphisms (SNPs) of Rhg1, Rhg4, and other novel quantitative trait loci (QTL). These markers were used to genotype a diverse set of 95 soybean germplasm lines and three recombinant inbred line (RIL) populations. SNP markers from the Rhg1 gene were able to differentiate copy number variation (CNV), such as resistant-high copy (PI 88788-type), low copy (Peking-type), and susceptible-single copy (Williams 82) numbers. Similarly, markers for the Rhg4 gene were able to detect Peking-type (resistance) genotypes. The phylogenetic information of SCN resistance loci from a large set of soybean accessions and the gene/QTL specific markers that were developed in this study will accelerate SCN resistance breeding programs.