Genetic architecture of protein and oil content in soybean seed and meal.

Soybean is grown primarily for the protein and oil extracted from its seed and its value is influenced by these components. The objective of this study was to map marker-trait associations (MTAs) for the concentration of seed protein, oil, and meal protein using the soybean nested association mapping (SoyNAM) population. The composition traits were evaluated on seed harvested from over 5000 inbred lines of the SoyNAM population grown in 10 field locations across 3 years. Estimated heritabilities were at least 0.85 for all three traits. The genotyping of lines with single nucleotide polymorphism markers resulted in the identification of 107 MTAs for the three traits. When MTAs for the three traits that mapped within 5 cM intervals were binned together, the MTAs were mapped to 64 intervals on 19 of the 20 soybean chromosomes. The majority of the MTA effects were small and of the 107 MTAs, 37 were for protein content, 39 for meal protein, and 31 for oil content. For cases where a protein and oil MTAs mapped to the same interval, most (94%) significant effects were opposite for the two traits, consistent with the negative correlation between these traits. A coexpression analysis identified candidate genes linked to MTAs and 18 candidate genes were identified. The large number of small effect MTAs for the composition traits suggest that genomic prediction would be more effective in improving these traits than marker-assisted selection.

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.

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.

QTL mapping and epistatic interaction analysis of field resistance to sudden death syndrome (Fusarium virguliforme) in soybean.

KEY MESSAGE: Two interactive quantitative trait loci (QTLs) controlled the field resistance to sudden death syndrome (SDS) in soybean. The interaction between them was confirmed. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is a major disease of soybean [Glycine max (L.) Merr.] in the United States. Breeding for soybean resistance to SDS is the most cost-effective method to manage the disease. The objective of this study was to identify and characterize quantitative trait loci (QTLs) underlying field resistance to SDS in a recombinant inbred line population from the cross GD2422 × LD01-5907. This population was genotyped with 1786 polymorphic single nucleotide polymorphisms (SNPs) using SoySNP6 K iSelect BeadChip and evaluated for SDS resistance in a naturally infested field. Four SDS resistance QTLs were mapped on Chromosomes 4, 8, 12 and 18. The resistant parent, LD01-5907, contributed the resistance alleles for the QTLs on Chromosomes 8 and 18 (qSDS-8 and qSDS-18), while the other parent, GD2422, provided the resistance alleles for the QTLs on Chromosomes 4 and 12 (qSDS-4 and qSDS-12). The minor QTL on Chromosome 12 (qSDS-12) is novel. The QTL on Chromosomes 8 and 18 (qSDS-8 and qSDS-18) overlapped with two soybean cyst nematode resistance-related loci, Rhg4 and Rhg1, respectively. A significant interaction between qSDS-8 and qSDS-18 was detected by disease incidence. Individual effects together with the interaction effect explained around 70% of the phenotypic variance. The epistatic interaction of qSDS-8 and qSDS-18 was confirmed by the field performance across multiple years. Furthermore, the resistance alleles at qSDS-8 and qSDS-18 were demonstrated to be recessive. The SNP markers linked to these QTLs will be useful for marker-assisted breeding to enhance the SDS resistance.

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.

A comparison of genotyping-by-sequencing analysis methods on low-coverage crop datasets shows advantages of a new workflow, GB-eaSy.

BACKGROUND: Genotyping-by-sequencing (GBS), a method to identify genetic variants and quickly genotype samples, reduces genome complexity by using restriction enzymes to divide the genome into fragments whose ends are sequenced on short-read sequencing platforms. While cost-effective, this method produces extensive missing data and requires complex bioinformatics analysis. GBS is most commonly used on crop plant genomes, and because crop plants have highly variable ploidy and repeat content, the performance of GBS analysis software can vary by target organism. Here we focus our analysis on soybean, a polyploid crop with a highly duplicated genome, relatively little public GBS data and few dedicated tools. RESULTS: We compared the performance of five GBS pipelines using low-coverage Illumina sequence data from three soybean populations. To address issues identified with existing methods, we developed GB-eaSy, a GBS bioinformatics workflow that incorporates widely used genomics tools, parallelization and automation to increase the accuracy and accessibility of GBS data analysis. Compared to other GBS pipelines, GB-eaSy rapidly and accurately identified the greatest number of SNPs, with SNP calls closely concordant with whole-genome sequencing of selected lines. Across all five GBS analysis platforms, SNP calls showed unexpectedly low convergence but generally high accuracy, indicating that the workflows arrived at largely complementary sets of valid SNP calls on the low-coverage data analyzed. CONCLUSIONS: We show that GB-eaSy is approximately as good as, or better than, other leading software solutions in the accuracy, yield and missing data fraction of variant calling, as tested on low-coverage genomic data from soybean. It also performs well relative to other solutions in terms of the run time and disk space required. In addition, GB-eaSy is built from existing open-source, modular software packages that are regularly updated and commonly used, making it straightforward to install and maintain. While GB-eaSy outperformed other individual methods on the datasets analyzed, our findings suggest that a comprehensive approach integrating the results from multiple GBS bioinformatics pipelines may be the optimal strategy to obtain the largest, most highly accurate SNP yield possible from low-coverage polyploid sequence data.

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.

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.

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

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2-kb unit contains four genes. One allele of Rhg1, Rhg1-b, is responsible for protecting most US soybean production from SCN. Whole-genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2-kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high-density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non-neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.

Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean.

The rhg1-b allele of soybean is widely used for resistance against soybean cyst nematode (SCN), the most economically damaging pathogen of soybeans in the United States. Gene silencing showed that genes in a 31-kilobase segment at rhg1-b, encoding an amino acid transporter, an α-SNAP protein, and a WI12 (wound-inducible domain) protein, each contribute to resistance. There is one copy of the 31-kilobase segment per haploid genome in susceptible varieties, but 10 tandem copies are present in an rhg1-b haplotype. Overexpression of the individual genes in roots was ineffective, but overexpression of the genes together conferred enhanced SCN resistance. Hence, SCN resistance mediated by the soybean quantitative trait locus Rhg1 is conferred by copy number variation that increases the expression of a set of dissimilar genes in a repeated multigene segment.

A nematode demographics assay in transgenic roots reveals no significant impacts of the Rhg1 locus LRR-Kinase on soybean cyst nematode resistance.

BACKGROUND: Soybean cyst nematode (Heterodera glycines, SCN) is the most economically damaging pathogen of soybean (Glycine max) in the U.S. The Rhg1 locus is repeatedly observed as the quantitative trait locus with the greatest impact on SCN resistance. The Glyma18g02680.1 gene at the Rhg1 locus that encodes an apparent leucine-rich repeat transmembrane receptor-kinase (LRR-kinase) has been proposed to be the SCN resistance gene, but its function has not been confirmed. Generation of fertile transgenic soybean lines is difficult but methods have been published that test SCN resistance in transgenic roots generated with Agrobacterium rhizogenes. RESULTS: We report use of artificial microRNA (amiRNA) for gene silencing in soybean, refinements to transgenic root SCN resistance assays, and functional tests of the Rhg1 locus LRR-kinase gene. A nematode demographics assay monitored infecting nematode populations for their progress through developmental stages two weeks after inoculation, as a metric for SCN resistance. Significant differences were observed between resistant and susceptible control genotypes. Introduction of the Rhg1 locus LRR-kinase gene (genomic promoter/coding region/terminator; Peking/PI 437654-derived SCN-resistant source), into rhg1- SCN-susceptible plant lines carrying the resistant-source Rhg4+ locus, provided no significant increases in SCN resistance. Use of amiRNA to reduce expression of the LRR-kinase gene from the Rhg1 locus of Fayette (PI 88788 source of Rhg1) also did not detectably alter resistance to SCN. However, silencing of the LRR-kinase gene did have impacts on root development. CONCLUSION: The nematode demographics assay can expedite testing of transgenic roots for SCN resistance. amiRNAs and the pSM103 vector that drives interchangeable amiRNA constructs through a soybean polyubiqutin promoter (Gmubi), with an intron-GFP marker for detection of transgenic roots, may have widespread use in legume biology. Studies in which expression of the Rhg1 locus LRR-kinase gene from different resistance sources was either reduced or complemented did not reveal significant impacts on SCN resistance.