Soybean BARCSoySNP6K: An assay for soybean genetics and breeding research

A Novel Quantitative Trait Locus Reduces Fusarium graminearum Infection in Glycine max Seedlings

Over the past 5 years, seedling diseases have caused an average annual loss of $21.8 million worth of United States' soybean (Glycine max) production, with Fusarium graminearum (teleomorph Gibberella zeae) emerging as a significant threat within the seedling disease complex. Its cross-pathogenicity on wheat and maize, along with increasing reports of fungicide resistance, highlights the need for improved genetic resistance in soybean. In a previous germplasm screening and genome-wide association study, we identified a significantly resistant accession, PI 438500, from a panel of 208 diverse soybean accessions. This accession carried fewer marker-trait associations and lower predicted resistance than other significantly resistant accessions yet displayed a highly resistant phenotype with low standard deviations. In this study, we developed an F2:3 mapping population from a biparental cross between PI 438500 and PI 548631 (highly susceptible to F. graminearum) and identified a quantitative trait locus (QTL) that explains 20.29% of the variation in postemergent visual severity. This QTL region contains multiple candidate genes with predicted roles in plant defense mechanisms and helps elucidate PI 438500's significant resistance to F. graminearum. Molecular markers linked and within this QTL region can help facilitate marker-assisted selection for F. graminearum resistance in soybean breeding.

Characterization of a G. max × G. soja nested association mapping population and identification of loci controlling seed composition traits from wild soybean

Wild soybean (Glycine soja Siebold & Zucc.) has valuable genetic diversity for improved disease resistance, stress tolerance, seed protein content and seed sulfur-containing amino acid concentrations. Many studies have reported loci controlling seed composition traits based on cultivated soybean populations, but wild soybean has been largely overlooked. In this study, a nested association mapping (NAM) population consisting of 10 families and 1107 recombinant inbred lines was developed by crossing 10 wild accessions with the common cultivar NC-Raleigh. Seed composition of the F6 generation grown at two locations was phenotyped, and genetic markers were identified for each line. The average number of recombination events in the wild soybean-derived population was significantly higher than that in the cultivated soybean-derived population, which resulted in a higher resolution for QTL mapping. Segregation bias in almost all NAM families was significantly biased toward the alleles of the wild soybean parent. Through single-family linkage mapping and association analysis of the entire NAM population, new QTLs with positive allele effects were identified from wild parents, including 5, 6, 18, 9, 16, 17 and 20 for protein content, oil content, total protein and oil content, methionine content, cysteine content, lysine content and threonine content, respectively. Candidate genes associated with these traits were identified based on gene annotations and gene expression levels in different tissues. This is the first study to reveal the genetic characteristics of wild soybean-derived populations, landscapes and the extent of effects of QTLs and candidate genes controlling traits from different wild soybean parents.

Unlocking soybean potential: genetic resources and omics for breeding

Soybean (Glycine max) is a vital foundation of global food security, providing a primary source of high-quality protein and oil for human consumption and animal feed. The rising global population has significantly increased the demand for soybeans, emphasizing the urgency of developing high-yield, stress-tolerant, and nutritionally superior cultivars. The extensive collection of soybean germplasm resources-including wild relatives, landraces, and cultivars-represents a valuable reservoir of genetic diversity critical for breeding advancements. Recent breakthroughs in genomic technologies, particularly high-throughput sequencing and multi-omics approaches, have revolutionized the identification of key genes associated with essential agronomic traits within these resources. These innovations enable precise and strategic utilization of genetic diversity, empowering breeders to integrate traits that improve yield potential, resilience to biotic and abiotic stresses, and nutritional quality. This review highlights the critical role of genetic resources and omics-driven innovations in soybean breeding. It also offers insights into strategies for accelerating the development of elite soybean cultivars to meet the growing demands of global soybean production.

QTL mapping and whole-genome sequencing analysis for novel genetic resources associated with sucrose content in soybean [Glycine max (L.) Merr.]

KEY MESSAGE

A major QTL for sucrose content was mapped on chromosome 8 in PI 506593. The novel genetic variants and candidate genes were further identified within the major QTL. Sucrose in soybean [Glycine max (L.) Merr.] contribute to animal feed efficiency and natural sweetness of soy products. Thus, identifying novel genetic resources, such as quantitative trait loci (QTL), associated with sucrose content in soybean is essential for enhancing seed values. In this study, two recombinant inbred line populations derived from the same high sucrose donor parent, PI 506593, were used to identify significant QTLs. A total of 11 sucrose-related regions on chromosomes (Chrs.) 4, 5, 6, 8, 10, and 13 were identified using QTL analysis. Among them, four QTLs (qSUC_08.1, qSUC_08.2, qSUC_08.3, and qSUC_08.4) were clustered in the interval of 40,597,410-42,861,364 bp on Chr. 8, which was considered major QTL region. A desirable marker at 41,834,095 bp was tested in two populations, showing that two phenotypically extreme groups were efficiently differentiated. We further identified 44 and 54 candidate genes with non-synonymous mutations in the major QTL region based on the annotations of Wm82.a2.v1 and Wm82.a5.v1 assemblies, respectively. Among 54 candidate genes from Wm82.a5.v1, Protein Variation Effect Analyzer (PROVEAN) revealed that 18 genes contained 34 variants that had deleterious impacts on biological functions. RNA-seq analysis highlighted five candidate genes that were highly expressed in pod and seed tissues during reproductive stages and other plant parts. A gene, Gm_Wm82_23219 (Glyma.08G293800, Wm82.a2.v1) encoding proline-rich protein 4-like, was highlighted in both PROVEAN and RNA-seq analyses. Novel findings in this study will be valuable genetic resources in soybean breeding programs that aim to improve efficiency in animal feed and human food.

Integration of Genetic and Imaging Data to Detect QTL for Root Traits in Interspecific Soybean Populations

Wild soybean, which has many desirable traits, such as adaptability to climate change-related stresses, is a valuable resource for expanding the narrow genetic diversity of cultivated soybeans. Plants require roots to adapt to different environments and optimize water and nutrient uptake to support growth and facilitate the storage of metabolites; however, it is challenging and costly to evaluate root traits under field conditions. Previous studies of quantitative trait loci (QTL) have been mainly based on cultivated soybean populations. In this study, an interspecific mapping population from a cross between wild soybean 'PI483463' and cultivar 'Hutcheson' was used to investigate QTLs associated with root traits using image data. Our results showed that 39 putative QTLs were distributed across 10 chromosomes (chr.). Seventeen of these were clustered in regions on chr. 8, 14, 15, 16, and 17, accounting for 19.92% of the phenotypic variation. We identified five significant QTL clusters influencing root-related traits, such as total root length, surface area, lateral total length, and number of tips, across five chr., with favorable alleles from both wild and cultivated soybeans. Furthermore, we identified eight candidate genes controlling these traits based on functional annotation. These genes were highly expressed in root tissues and directly or indirectly affected soybean root growth, development, and stress responses. Our results provide valuable insights for breeders aiming to optimize soybean root traits and leveraging genetic diversity from wild soybean species to develop varieties with improved root morphological traits, ultimately enhancing overall plant growth, productivity, and resilience.

Allele-tagged TaqMan® PCR genotyping assays for high-throughput detection of soybean cyst nematode resistance

BACKGROUND

Whole genome resequencing (WGRS) platforms provide exceptional fingerprinting of the entire genome but are expensive and less flexible to use as a routine genotyping tool for targeting causal polymorphisms within a germplasm collection or breeding program. Therefore, there has been a continuous effort to develop small-scale genotyping platforms that facilitate robust and quick assessments of the allelic status of causal variants for important traits within soybean breeding programs. The objective was to develop a comprehensive panel of soybean cyst nematode (SCN) resistance TaqMan® assays via selecting the causative genes and analyzing their associated alleles.

METHODS

The Soybean Allele Catalog was utilized to investigate WGRS-derived variants which are predicted to cause a change in the amino acid sequence of a gene product. This panel of TaqMan® assays reflects current knowledge about known SCN resistance-causing genes and their associated alleles: GmSNAP18-a and -b, GmSNAP11, GmSHMT08, GmSNAP15, GmNSFRAN07, and GmSNAP02-ins and -del. Developed assays were tested using elite breeding lines and segregating populations. TaqMan assays were compared to other currently available KASP and CAPS assays.

CONCLUSION

All assays showed excellent allele determination efficiencies. This SCN genotyping assay panel can be utilized as a simplified, accurate and reliable genotyping platform further equipping the updated soybean breeding toolbox.

Soybean genomics research community strategic plan: A vision for 2024-2028

This strategic plan summarizes the major accomplishments achieved in the last quinquennial by the soybean [Glycine max (L.) Merr.] genetics and genomics research community and outlines key priorities for the next 5 years (2024-2028). This work is the result of deliberations among over 50 soybean researchers during a 2-day workshop in St Louis, MO, USA, at the end of 2022. The plan is divided into seven traditional areas/disciplines: Breeding, Biotic Interactions, Physiology and Abiotic Stress, Functional Genomics, Biotechnology, Genomic Resources and Datasets, and Computational Resources. One additional section was added, Training the Next Generation of Soybean Researchers, when it was identified as a pressing issue during the workshop. This installment of the soybean genomics strategic plan provides a snapshot of recent progress while looking at future goals that will improve resources and enable innovation among the community of basic and applied soybean researchers. We hope that this work will inform our community and increase support for soybean research.

Glycine soja, PI424025, is a valuable genetic resource to improve soybean seed-protein content and composition

Soybean seed-protein content and composition is important because it contributes over half the value of this $61 billion crop. Historic negative correlation between seed-protein content and oil have been reported. Similarly, negative correlation between seed-protein content and yield have been reported but may be at least in part mitigated by increasing the genetic diversity of the elite soybean germplasm. Improvements in amino acid composition of seed-protein would increase the value of soybean meal and protein for animal and human consumption. We have identified a genetic resource in wild soybean germplasm, PI424025B, with elevated seed-protein, elevated cysteine in the seed and increased sulfur content in the seed. We have developed a population of Recombinant Inbred Lines derived from a cross between NC-Raleigh and PI424025B. We evaluated the carbon, nitrogen and sulfur (C, N and S) content of seeds from the progeny, the parents and other high-seed protein soybean lines. N content and C/N (ratio of C to N) were well correlated with protein content measured by NIR. PI424025B had a N content comparable to high-protein soybean lines and a superior N/S. These traits were inherited by some of the progeny. Quantitative Trait Loci (QTL) associated with high-seed protein were identified on chr2, chr20 and chr15 and colocalizes with well characterized loci on chr 15 and chr 20 known to affect seed-protein content and quality. A potentially unique QTL was identified on Chr15. Unlike the chr20 QTL, the chr15 QTL improved S content relative to N content and was superior to other high seed-protein phenotypes tested. These data indicate that PI424025B is a valuable resource to diversify the genetics of soybean while improving soybean seed-protein content and composition.

Development and implementation of nested single-nucleotide polymorphism (SNP) assays for breeding and genetic research applications

SoySNP50K and SoySNP6K are commonly used for soybean (Glycine max) genotyping. The SoySNP50K assay has been used to genetically analyze the entire USDA Soybean Germplasm Collection, while the SoySNP6K assay, containing a subset of 6000 single-nucleotide polymorphisms (SNPs) from SoySNP50K, has been used for quantitative trait loci mapping of different traits. To meet the needs for genomic selection, selection of parents for crosses, and characterization of breeding populations, especially early selection of ideal offspring from thousands of lines, we developed two assays, SoySNP3K and SoySNP1K, containing 3072 and 1252 SNPs, respectively, based on SoySNP50K and SoySNP6K mark sets. These two assays also contained the trait markers reported or contributed by soybean breeders. The SNPs in the SoySNP3K are a subset from SoySNP6K, while the SNPs in the SoySNP1K are a subset from SoySNP3K. These SNPs were chosen to reduce the SNP number in the large linkage blocks while capturing as much of the haplotype diversity as possible. They are highly polymorphic and of high quality. The mean minor allele frequencies of the SNPs in the southern and northern US elites were 0.25 and 0.27 for SoySNP3K, respectively, and 0.29 and 0.33 for SoySNP1K. The selected SNPs are a valuable source for developing targeted amplicon sequencing assay or beadchip assay in soybean. SoySNP3K and SoySNP1K assays are commercialized by Illumina Inc. and AgriPlex Genomics, respectively. Together with SoySNP50K and SoySNP6K, a series of nested assays with different marker densities will serve as additional low-cost genomic tools for genetic, genomic, and breeding research.

Comparison of Rps loci toward isolates, singly and combined inocula, of Phytophthora sojae in soybean PI 407985, PI 408029, PI 408097, and PI424477

For soybean, novel single dominant Resistance to Phytophthora sojae (Rps) genes are sought to manage Phytophthora root and stem rot. In this study, resistance to P. sojae was mapped individually in four recombinant inbred line (RIL) populations derived from crosses of the susceptible cultivar Williams with PI 407985, PI 408029, PI 408097, and PI424477 previously identified as putative novel sources of disease resistance. Each population was screened for resistance with five to seven isolates of P. sojae separately over multiple F7-F10 generations. Additionally, three of the populations were screened with inoculum from the combination of three P. sojae isolates (PPR), which comprised virulence to 14 Rps genes. Over 2,300 single-nucleotide polymorphism markers were used to construct genetic maps in each population to identify chromosomal regions associated with resistance to P. sojae. Resistance segregated as one or two genes to the individual isolates and one gene toward PPR in each population and mapped to chromosomes 3, 13, or 18 in one or more of the four RIL populations. Resistance to five isolates mapped to the same chromosome 3 region are as follows: OH7 (PI 424477 and PI408029), OH12168, OH7/8, PPR (PI 407985), and 1.S.1.1 (PI408029). The resistance regions on chromosome 13 also overlapped for OH1, OH25, OH-MIA (PI424477), PPR (PI 424477, PI 407985, and PI 408097), PPR and OH0217 (PI 408097), and OH4 (PI 408029), but were distinct for each population suggesting multiple genes confer resistance. Two regions were identified on chromosome 18 but all appear to map to known loci; notably, resistance to the combined inoculum (PPR) did not map at this locus. However, there are putative new alleles in three of four populations, three on chromosome 3 and two on chromosome 13 based on mapping location but also known virulence in the isolate used. This characterization of all the Rps genes segregating in these populations to these isolates will be informative for breeding, but the combined inoculum was able to map a novel loci. Furthermore, within each of these P. sojae isolates, there was virulence to more than the described Rps genes, and the effectiveness of the novel genes requires testing in larger populations.

Identification of Quantitative Trait Loci Controlling Root Morphological Traits in an Interspecific Soybean Population Using 2D Imagery Data

Roots are the hidden and most important part of plants. They serve as stabilizers and channels for uptaking water and nutrients and play a crucial role in the growth and development of plants. Here, two-dimensional image data were used to identify quantitative trait loci (QTL) controlling root traits in an interspecific mapping population derived from a cross between wild soybean 'PI366121' and cultivar 'Williams 82'. A total of 2830 single-nucleotide polymorphisms were used for genotyping, constructing genetic linkage maps, and analyzing QTLs. Forty-two QTLs were identified on twelve chromosomes, twelve of which were identified as major QTLs, with a phenotypic variation range of 36.12% to 39.11% and a logarithm of odds value range of 12.01 to 17.35. Two significant QTL regions for the average diameter, root volume, and link average diameter root traits were detected on chromosomes 3 and 13, and both wild and cultivated soybeans contributed positive alleles. Six candidate genes, Glyma.03G027500 (transketolase/glycoaldehyde transferase), Glyma.03G014500 (dehydrogenases), Glyma.13G341500 (leucine-rich repeat receptor-like protein kinase), Glyma.13G341400 (AGC kinase family protein), Glyma.13G331900 (60S ribosomal protein), and Glyma.13G333100 (aquaporin transporter) showed higher expression in root tissues based on publicly available transcriptome data. These results will help breeders improve soybean genetic components and enhance soybean root morphological traits using desirable alleles from wild soybeans.

Genetic mapping reveals the complex genetic architecture controlling slow canopy wilting in soybean

In soybean [Glycine max (L.) Merr.], drought stress is the leading cause of yield loss from abiotic stress in rain-fed US growing areas. Only 10% of the US soybean production is irrigated; therefore, plants must possess physiological mechanisms to tolerate drought stress. Slow canopy wilting is a physiological trait that is observed in a few exotic plant introductions (PIs) and may lead to yield improvement under drought stress. Canopy wilting of 130 recombinant inbred lines (RILs) derived from Hutcheson × PI 471938 grown under drought stress was visually evaluated and genotyped with the SoySNP6K BeadChip. Over four years, field evaluations of canopy wilting were conducted under rainfed conditions at three locations across the US (Georgia, Kansas, and North Carolina). Due to the variation in weather among locations and years, the phenotypic data were collected from seven environments. Substantial variation in canopy wilting was observed among the genotypes in the RIL population across environments. Three QTLs were identified for canopy wilting from the RIL population using composite interval mapping on chromosomes (Chrs) 2, 8, and 9 based on combined environmental analyses. These QTLs inherited the favorable alleles from PI 471938 and accounted for 11, 10, and 14% of phenotypic variation, respectively. A list of 106 candidate genes were narrowed down for these three QTLs based on the published information. The QTLs identified through this research can be used as targets for further investigation to understand the mechanisms of slow canopy wilting. These QTLs could be deployed to improve drought tolerance through a targeted selection of the genomic regions from PI 471938.

Soybean genetics, genomics, and breeding for improving nutritional value and reducing antinutritional traits in food and feed

Soybean [Glycine max (L.) Merr.] is a globally important crop due to its valuable seed composition, versatile feed, food, and industrial end-uses, and consistent genetic gain. Successful genetic gain in soybean has led to widespread adaptation and increased value for producers, processors, and consumers. Specific focus on the nutritional quality of soybean seed composition for food and feed has further elucidated genetic knowledge and bolstered breeding progress. Seed components are historical and current targets for soybean breeders seeking to improve nutritional quality of soybean. This article reviews genetic and genomic foundations for improvement of nutritionally important traits, such as protein and amino acids, oil and fatty acids, carbohydrates, and specific food-grade considerations; discusses the application of advanced breeding technology such as CRISPR/Cas9 in creating seed composition variations; and provides future directions and breeding recommendations regarding soybean seed composition traits.

Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Genomic selection of soybean (Glycine max) for genetic improvement of yield and seed composition in a breeding context

Genomic selection has been utilized for genetic improvement in both plant and animal breeding and is a favorable technique for quantitative trait development. Within this study, genomic selection was evaluated within a breeding program, using novel validation methods in addition to plant materials and data from a commercial soybean (Glycine max) breeding program. A total of 1501 inbred lines were used to test multiple genomic selection models for multiple traits. Validation included cross-validation, inter-environment, and empirical validation. The results indicated that the extended genomic best linear unbiased prediction (EGBLUP) model was the most effective model tested for yield, protein, and oil in cross-validation with accuracies of 0.50, 0.68, and 0.64, respectively. Increasing marker number from 1000 to 3000 to 6000 single nucleotide polymorphism markers leads to statistically significant increases in accuracy. Cross-environment predictions were statistically lower than cross-validation with accuracies of 0.24, 0.54, and 0.42 for yield, protein, and oil, respectively, using the extended genomic BLUP model. Empirical validation, predicting the yield of 510 soybean lines, had a prediction accuracy of 0.34, with the inclusion of a maturity covariate leading to a notable increase in accuracy. Genomic selection identified high-performance lines in inter-environment predictions: 34% of lines within the upper quartile of yield, and 51% and 48% of the highest quartile protein and oil lines, respectively. Statistically similar results occurred comparing rankings in empirical validation and selection for advancements in yield trials. These results indicate that genomic selection is a useful tool for selection decisions.

Identification and mapping of a recessive allele, dt3, specifying semideterminate stem growth habit in soybean

KEY MESSAGE

A locus, dt3, modulating semideterminancy in soybean, was discovered by a combination of genome-wide association studies and linkage mapping with multiple distinct biparental populations. Stem growth habit is a key architectural trait in many plants that contributes to plant productivity and environmental adaptation. In soybean, stem growth habit is classified as indeterminate, semideterminate, or determinate, of which semideterminacy is often considered as a counterpart of the "Green Revolution" trait in cereals that significantly increased grain yields. It has been demonstrated that semideterminacy in soybean is modulated by epistatic interaction between two loci, Dt1 on chromosome 19 and Dt2 on chromosome 18, with the latter as a negative regulator of the former. Here, we report the discovery of a third locus, Dt3, modulating soybean stem growth habit, which was delineated to a ~ 196-kb region on chromosome 10 by a combination of allelic and haplotypic analysis of the Dt1 and Dt2 loci in the USDA soybean Germplasm Collection, genome-wide association studies with three subsets of the collection, and linkage mapping with four biparental populations derived from crosses between one of two elite indeterminate cultivars and each of four semideterminate varieties possessing neither Dt2 nor dt1. These four semideterminate varieties are recessive mutants (i.e., dt3/dt3) in the Dt1/Dt1;dt2/dt2 background. As the semideterminacy modulated by the Dt2 allele has unfavorable pleotropic effects such as sensitivity to drought stress, dt3 may be an ideal alternative for use to develop semideterminate cultivars that are more resilient to such an environmental stress. This study enhances our understanding of the genetic factors underlying semideterminacy and enables more accurate marker-assisted selection for stem growth habit in soybean breeding.

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

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

Models to estimate genetic gain of soybean seed yield from annual multi-environment field trials

KEY MESSAGE

Simulations demonstrated that estimates of realized genetic gain from linear mixed models using regional trials are biased to some degree. Thus, we recommend multiple selected models to obtain a range of reasonable estimates. Genetic improvements of discrete characteristics are obvious and easy to demonstrate, while quantitative traits require reliable and accurate methods to disentangle the confounding genetic and non-genetic components. Stochastic simulations of soybean [Glycine max (L.) Merr.] breeding programs were performed to evaluate linear mixed models to estimate the realized genetic gain (RGG) from annual multi-environment trials (MET). True breeding values were simulated under an infinitesimal model to represent the genetic contributions to soybean seed yield under various MET conditions. Estimators were evaluated using objective criteria of bias and linearity. Covariance modeling and direct versus indirect estimation-based models resulted in a substantial range of estimated values, all of which were biased to some degree. Although no models produced unbiased estimates, the three best-performing models resulted in an average bias of [Formula: see text] kg/ha[Formula: see text]/yr[Formula: see text] ([Formula: see text] bu/ac[Formula: see text]/yr[Formula: see text]). Rather than relying on a single model to estimate RGG, we recommend the application of several models with minimal and directional bias. Further, based on the parameters used in the simulations, we do not think it is appropriate to use any single model to compare breeding programs or quantify the efficiency of proposed new breeding strategies. Lastly, for public soybean programs breeding for maturity groups II and III in North America, the estimated RGG values ranged from 18.16 to 39.68 kg/ha[Formula: see text]/yr[Formula: see text] (0.27-0.59 bu/ac[Formula: see text]/yr[Formula: see text]) from 1989 to 2019. These results provide strong evidence that public breeders have significantly improved soybean germplasm for seed yield in the primary production areas of North America.

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.

Genetic Mapping for QTL Associated with Seed Nickel and Molybdenum Accumulation in the Soybean 'Forrest' by 'Williams 82' RIL Population

Understanding the genetic basis of seed Ni and Mo is essential. Since soybean is a major crop in the world and a major source for nutrients, including Ni and Mo, the objective of the current research was to map genetic regions (quantitative trait loci, QTL) linked to Ni and Mo concentrations in soybean seed. A recombinant inbred line (RIL) population was derived from a cross between 'Forrest' and 'Williams 82' (F × W82). A total of 306 lines was used for genotyping using 5405 single nucleotides polymorphism (SNP) markers using Infinium SNP6K BeadChips. A two-year experiment was conducted and included the parents and the RIL population. One experiment was conducted in 2018 in North Carolina (NC), and the second experiment was conducted in Illinois in 2020 (IL). Logarithm of the odds (LOD) of ≥2.5 was set as a threshold to report identified QTL using the composite interval mapping (CIM) method. A wide range of Ni and Mo concentrations among RILs was observed. A total of four QTL (qNi-01, qNi-02, and qNi-03 on Chr 2, 8, and 9, respectively, in 2018, and qNi-01 on Chr 20 in 2020) was identified for seed Ni. All these QTL were significantly (LOD threshold > 2.5) associated with seed Ni, with LOD scores ranging between 2.71-3.44, and with phenotypic variance ranging from 4.48-6.97%. A total of three QTL for Mo (qMo-01, qMo-02, and qMo-03 on Chr 1, 3, 17, respectively) was identified in 2018, and four QTL (qMo-01, qMo-02, qMo-03, and qMo-04, on Chr 5, 11, 14, and 16, respectively) were identified in 2020. Some of the current QTL had high LOD and significantly contributed to the phenotypic variance for the trait. For example, in 2018, Mo QTL qMo-01 on Chr 1 had LOD of 7.8, explaining a phenotypic variance of 41.17%, and qMo-03 on Chr 17 had LOD of 5.33, with phenotypic variance explained of 41.49%. In addition, one Mo QTL (qMo-03 on Chr 14) had LOD of 9.77, explaining 51.57% of phenotypic variance related to the trait, and another Mo QTL (qMo-04 on Chr 16) had LOD of 7.62 and explained 49.95% of phenotypic variance. None of the QTL identified here were identified twice across locations/years. Based on a search of the available literature and of SoyBase, the four QTL for Ni, identified on Chr 2, 8, 9, and 20, and the five QTL associated with Mo, identified on Chr 1, 17, 11, 14, and 16, are novel and not previously reported. This research contributes new insights into the genetic mapping of Ni and Mo, and provides valuable QTL and molecular markers that can potentially assist in selecting Ni and Mo levels in soybean seeds.

Genetic architecture of soybean tolerance to off-target dicamba

The adoption of dicamba-tolerant (DT) soybean in the United States resulted in extensive off-target dicamba damage to non-DT vegetation across soybean-producing states. Although soybeans are highly sensitive to dicamba, the intensity of observed symptoms and yield losses are affected by the genetic background of genotypes. Thus, the objective of this study was to detect novel marker-trait associations and expand on previously identified genomic regions related to soybean response to off-target dicamba. A total of 551 non-DT advanced breeding lines derived from 232 unique bi-parental populations were phenotyped for off-target dicamba across nine environments for three years. Breeding lines were genotyped using the Illumina Infinium BARCSoySNP6K BeadChip. Filtered SNPs were included as predictors in Random Forest (RF) and Support Vector Machine (SVM) models in a forward stepwise selection loop to identify the combination of SNPs yielding the highest classification accuracy. Both RF and SVM models yielded high classification accuracies (0.76 and 0.79, respectively) with minor extreme misclassifications (observed tolerant predicted as susceptible, and vice-versa). Eight genomic regions associated with off-target dicamba tolerance were identified on chromosomes 6 [Linkage Group (LG) C2], 8 (LG A2), 9 (LG K), 10 (LG O), and 19 (LG L). Although the genetic architecture of tolerance is complex, high classification accuracies were obtained when including the major effect SNP identified on chromosome 6 as the sole predictor. In addition, candidate genes with annotated functions associated with phases II (conjugation of hydroxylated herbicides to endogenous sugar molecules) and III (transportation of herbicide conjugates into the vacuole) of herbicide detoxification in plants were co-localized with significant markers within each genomic region. Genomic prediction models, as reported in this study, can greatly facilitate the identification of genotypes with superior tolerance to off-target dicamba.

Identification of quantitative trait loci controlling soybean seed protein and oil content

Soybean is a major source of seed protein and oil globally with an average composition of 40% protein and 20% oil in the seed. The goal of this study was to identify quantitative trait loci (QTL) conferring seed protein and oil content utilizing a population constructed by crossing an above average protein content line, PI 399084 to another line that had a low protein content value, PI 507429, both from the USDA soybean germplasm collection. The recombinant inbred line (RIL) population, PI 507429 x PI 399084, was evaluated in two replications over four years (2018-2021); the seeds were analyzed for seed protein and oil content using near-infrared reflectance spectroscopy. The recombinant inbred lines and the two parents were re-sequenced using genotyping by sequencing. A total of 12,761 molecular markers, which came from genotyping by sequencing, the SoySNP6k BeadChip and selected simple sequence repeat (SSR) markers from known protein QTL chromosomal regions were used for mapping. One QTL was identified on chromosome 2 explaining up to 56.8% of the variation for seed protein content and up to 43% for seed oil content. Another QTL identified on chromosome 15 explained up to 27.2% of the variation for seed protein and up to 41% of the variation for seed oil content. The protein and oil QTLs of this study and their associated molecular markers will be useful in breeding to improve nutritional quality in soybean.

Identification of candidate genes for soybean seed coat-related traits using QTL mapping and GWAS

Seed coat color is a typical morphological trait that can be used to reveal the evolution of soybean. The study of seed coat color-related traits in soybeans is of great significance for both evolutionary theory and breeding practices. In this study, 180 F₁₀ recombinant inbred lines (RILs) derived from the cross between the yellow-seed coat cultivar Jidou12 (ZDD23040, JD12) and the wild black-seed coat accession Y9 (ZYD02739) were used as materials. Three methods, single-marker analysis (SMA), interval mapping (IM), and inclusive composite interval mapping (ICIM), were used to identify quantitative trait loci (QTLs) controlling seed coat color and seed hilum color. Simultaneously, two genome-wide association study (GWAS) models, the generalized linear model (GLM) and mixed linear model (MLM), were used to jointly identify seed coat color and seed hilum color QTLs in 250 natural populations. By integrating the results from QTL mapping and GWAS analysis, we identified two stable QTLs (qSCC02 and qSCC08) associated with seed coat color and one stable QTL (qSHC08) related to seed hilum color. By combining the results of linkage analysis and association analysis, two stable QTLs (qSCC02, qSCC08) for seed coat color and one stable QTL (qSHC08) for seed hilum color were identified. Upon further investigation using Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, we validated the previous findings that two candidate genes (CHS3C and CHS4A) reside within the qSCC08 region and identified a new QTL, qSCC02. There were a total of 28 candidate genes in the interval, among which Glyma.02G024600, Glyma.02G024700, and Glyma.02G024800 were mapped to the glutathione metabolic pathway, which is related to the transport or accumulation of anthocyanin. We considered the three genes as potential candidate genes for soybean seed coat-related traits. The QTLs and candidate genes detected in this study provide a foundation for further understanding the genetic mechanisms underlying soybean seed coat color and seed hilum color and are of significant value in marker-assisted breeding.

Genomic prediction of optimal cross combinations to accelerate genetic improvement of soybean (Glycine max)

Improving yield is a primary soybean breeding goal, as yield is the main determinant of soybean's profitability. Within the breeding process, selection of cross combinations is one of most important elements. Cross prediction will assist soybean breeders in identifying the best cross combinations among parental genotypes prior to crossing, increasing genetic gain and breeding efficiency. In this study optimal cross selection methods were created and applied in soybean and validated using historical data from the University of Georgia soybean breeding program, under multiple training set compositions and marker densities utilizing multiple genomic selection models for marker evaluation. Plant materials consisted of 702 advanced breeding lines evaluated in multiple environments and genotyped using SoySNP6k BeadChips. An additional marker set, the SoySNP3k marker set, was tested in this study as well. Optimal cross selection methods were used to predict the yield of 42 previously made crosses and compared to the performance of the cross's offspring in replicated field trials. The best prediction accuracy was obtained when using Extended Genomic BLUP with the SoySNP6k marker set, consisting of 3,762 polymorphic markers, with an accuracy of 0.56 with a training set maximally related to the crosses predicted and 0.4 in a training set with minimized relatedness to predicted crosses. Prediction accuracy was most significantly impacted by training set relatedness to the predicted crosses, marker density, and the genomic model used to predict marker effects. The usefulness criterion selected had an impact on prediction accuracy within training sets with low relatedness to the crosses predicted. Optimal cross prediction provides a useful method that assists plant breeders in selecting crosses in soybean breeding.

Development of SNP marker panels for genotyping by target sequencing (GBTS) and its application in soybean

A high-throughput genotyping platform with customized flexibility, high genotyping accuracy, and low cost is critical for marker-assisted selection and genetic mapping in soybean. Three assay panels were selected from the SoySNP50K, 40K, 20K, and 10K arrays, containing 41,541, 20,748, and 9670 SNP markers, respectively, for genotyping by target sequencing (GBTS). Fifteen representative accessions were used to assess the accuracy and consistency of the SNP alleles identified by the SNP panels and sequencing platform. The SNP alleles were 99.87% identical between technical replicates and 98.86% identical between the 40K SNP GBTS panel and 10× resequencing analysis. The GBTS method was also accurate in the sense that the genotypic dataset of the 15 representative accessions correctly revealed the pedigree of the accessions, and the biparental progeny datasets correctly constructed the linkage maps of the SNPs. The 10K panel was also used to genotype two parent-derived populations and analyze QTLs controlling 100-seed weight, resulting in the identification of the stable associated genetic locus Locus_OSW_06 on chromosome 06. The markers flanking the QTL explained 7.05% and 9.83% of the phenotypic variation, respectively. Compared with GBS and DNA chips, the 40K, 20K, and 10K panels reduced costs by 5.07% and 58.28%, 21.44% and 65.48%, and 35.74% and 71.76%, respectively. Low-cost genotyping panels could facilitate soybean germplasm assessment, genetic linkage map construction, QTL identification, and genomic selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01372-6.

Multi-omics assisted breeding for biotic stress resistance in soybean

Biotic stress is a critical factor limiting soybean growth and development. Soybean responses to biotic stresses such as insects, nematodes, fungal, bacterial, and viral pathogens are governed by complex regulatory and defense mechanisms. Next-generation sequencing has availed research techniques and strategies in genomics and post-genomics. This review summarizes the available information on marker resources, quantitative trait loci, and marker-trait associations involved in regulating biotic stress responses in soybean. We discuss the differential expression of related genes and proteins reported in different transcriptomics and proteomics studies and the role of signaling pathways and metabolites reported in metabolomic studies. Recent advances in omics technologies offer opportunities to reshape and improve biotic stress resistance in soybean by altering gene regulation and/or other regulatory networks. We suggest using 'integrated omics' to precisely understand how soybean responds to different biotic stresses. We also discuss the potential challenges of integrating multi-omics for the functional analysis of genes and their regulatory networks and the development of biotic stress-resistant cultivars. This review will help direct soybean breeding programs to develop resistance against different biotic stresses.

Development of molecular inversion probes for soybean progeny genomic selection genotyping

Increasing rate of genetic gain for key agronomic traits through genomic selection requires the development of new molecular methods to run genome-wide single-nucleotide polymorphisms (SNPs). The main limitation of current methods is the cost is too high to screen breeding populations. Molecular inversion probes (MIPs) are a targeted genotyping-by-sequencing (GBS) method that could be used for soybean [Glycine max (L.) Merr.] that is both cost-effective, high-throughput, and provides high data quality to screen breeder's germplasm for genomic selection. A 1K MIP SNP set was developed for soybean with uniformly distributed markers across the genome. The SNPs were selected to maximize the number of informative markers in germplasm being tested in soybean breeding programs located in the northern-central and middle-southern regions of the United States. The 1K SNP MIP set was tested on diverse germplasm and a recombinant inbred line (RIL) population. Targeted sequencing with MIPs obtained an 85% enrichment for the targeted SNPs. The MIP genotyping accuracy was 93% overall, whereas homozygous call accuracy was 98% with <10% missing data. The accuracy of MIPs combined with its low per-sample cost makes it a powerful tool to enable genomic selection within soybean breeding programs.

3D-GBS: a universal genotyping-by-sequencing approach for genomic selection and other high-throughput low-cost applications in species with small to medium-sized genomes

Despite the increased efficiency of sequencing technologies and the development of reduced-representation sequencing (RRS) approaches allowing high-throughput sequencing (HTS) of multiplexed samples, the per-sample genotyping cost remains the most limiting factor in the context of large-scale studies. For example, in the context of genomic selection (GS), breeders need genome-wide markers to predict the breeding value of large cohorts of progenies, requiring the genotyping of thousands candidates. Here, we introduce 3D-GBS, an optimized GBS procedure, to provide an ultra-high-throughput and ultra-low-cost genotyping solution for species with small to medium-sized genome and illustrate its use in soybean. Using a combination of three restriction enzymes (PstI/NsiI/MspI), the portion of the genome that is captured was reduced fourfold (compared to a "standard" ApeKI-based protocol) while reducing the number of markers by only 40%. By better focusing the sequencing effort on limited set of restriction fragments, fourfold more samples can be genotyped at the same minimal depth of coverage. This GBS protocol also resulted in a lower proportion of missing data and provided a more uniform distribution of SNPs across the genome. Moreover, we investigated the optimal number of reads per sample needed to obtain an adequate number of markers for GS and QTL mapping (500-1000 markers per biparental cross). This optimization allows sequencing costs to be decreased by ~ 92% and ~ 86% for GS and QTL mapping studies, respectively, compared to previously published work. Overall, 3D-GBS represents a unique and affordable solution for applications requiring extremely high-throughput genotyping where cost remains the most limiting factor.

Identification of genomic regions associated with soybean responses to off-target dicamba exposure

The widespread adoption of genetically modified (GM) dicamba-tolerant (DT) soybean was followed by numerous reports of off-target dicamba damage and yield losses across most soybean-producing states. In this study, a subset of the USDA Soybean Germplasm Collection consisting of 382 genetically diverse soybean accessions originating from 15 countries was used to identify genomic regions associated with soybean response to off-target dicamba exposure. Accessions were genotyped with the SoySNP50K BeadChip and visually screened for damage in environments with prolonged exposure to off-target dicamba. Two models were implemented to detect significant marker-trait associations: the Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) and a model that allows the inclusion of population structure in interaction with the environment (G×E) to account for variable patterns of genotype responses in different environments. Most accessions (84%) showed a moderate response, either moderately tolerant or moderately susceptible, with approximately 8% showing tolerance and susceptibility. No differences in off-target dicamba damage were observed across maturity groups and centers of origin. Both models identified significant associations in regions of chromosomes 10 and 19. The BLINK model identified additional significant marker-trait associations on chromosomes 11, 14, and 18, while the G×E model identified another significant marker-trait association on chromosome 15. The significant SNPs identified by both models are located within candidate genes possessing annotated functions involving different phases of herbicide detoxification in plants. These results entertain the possibility of developing non-GM soybean cultivars with improved tolerance to off-target dicamba exposure and potentially other synthetic auxin herbicides. Identification of genetic sources of tolerance and genomic regions conferring higher tolerance to off-target dicamba may sustain and improve the production of other non-DT herbicide soybean production systems, including the growing niche markets of organic and conventional soybean.

Breeding for disease resistance in soybean: a global perspective

KEY MESSAGE

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

Incorporation of Soil-Derived Covariates in Progeny Testing and Line Selection to Enhance Genomic Prediction Accuracy in Soybean Breeding

The availability of high-dimensional molecular markers has allowed plant breeding programs to maximize their efficiency through the genomic prediction of a phenotype of interest. Yield is a complex quantitative trait whose expression is sensitive to environmental stimuli. In this research, we investigated the potential of incorporating soil texture information and its interaction with molecular markers via covariance structures for enhancing predictive ability across breeding scenarios. A total of 797 soybean lines derived from 367 unique bi-parental populations were genotyped using the Illumina BARCSoySNP6K and tested for yield during 5 years in Tiptonville silt loam, Sharkey clay, and Malden fine sand environments. Four statistical models were considered, including the GBLUP model (M1), the reaction norm model (M2) including the interaction between molecular markers and the environment (G×E), an extended version of M2 that also includes soil type (S), and the interaction between soil type and molecular markers (G×S) (M3), and a parsimonious version of M3 which discards the G×E term (M4). Four cross-validation scenarios simulating progeny testing and line selection of tested–untested genotypes (TG, UG) in observed–unobserved environments [OE, UE] were implemented (CV2 [TG, OE], CV1 [UG, OE], CV0 [TG, UE], and CV00 [UG, UE]). Across environments, the addition of G×S interaction in M3 decreased the amount of variability captured by the environment (−30.4%) and residual (−39.2%) terms as compared to M1. Within environments, the G×S term in M3 reduced the variability captured by the residual term by 60 and 30% when compared to M1 and M2, respectively. M3 outperformed all the other models in CV2 (0.577), CV1 (0.480), and CV0 (0.488). In addition to the Pearson correlation, other measures were considered to assess predictive ability and these showed that the addition of soil texture seems to structure/dissect the environmental term revealing its components that could enhance or hinder the predictability of a model, especially in the most complex prediction scenario (CV00). Hence, the availability of soil texture information before the growing season could be used to optimize the efficiency of a breeding program by allowing the reconsideration of field experimental design, allocation of resources, reduction of preliminary trials, and shortening of the breeding cycle.

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.

Linkage analysis and residual heterozygotes derived near isogenic lines reveals a novel protein quantitative trait loci from a Glycine soja accession

Modern soybean [Glycine max (L.) Merr] cultivars have low overall genetic variation due to repeated bottleneck events that arose during domestication and from selection strategies typical of many soybean breeding programs. In both public and private soybean breeding programs, the introgression of wild soybean (Glycine soja Siebold and Zucc.) alleles is a viable option to increase genetic diversity and identify new sources for traits of value. The objectives of our study were to examine the genetic architecture responsible for seed protein and oil using a recombinant inbred line (RIL) population derived from hybridizing a G. max line ('Osage') with a G. soja accession (PI 593983). Linkage mapping identified a total of seven significant quantitative trait loci on chromosomes 14 and 20 for seed protein and on chromosome 8 for seed oil with LOD scores ranging from 5.3 to 31.7 for seed protein content and from 9.8 to 25.9 for seed oil content. We analyzed 3,015 single F4:9 soybean plants to develop two residual heterozygotes derived near isogenic lines (RHD-NIL) populations by targeting nine SNP markers from genotype-by-sequencing, which corresponded to two novel quantitative trait loci (QTL) derived from G. soja: one for a novel seed oil QTL on chromosome 8 and another for a novel protein QTL on chromosome 14. Single marker analysis and linkage analysis using 50 RHD-NILs validated the chromosome 14 protein QTL, and whole genome sequencing of RHD-NILs allowed us to reduce the QTL interval from ∼16.5 to ∼4.6 Mbp. We identified two genomic regions based on recombination events which had significant increases of 0.65 and 0.72% in seed protein content without a significant decrease in seed oil content. A new Kompetitive allele-specific polymerase chain reaction (KASP) assay, which will be useful for introgression of this trait into modern elite G. max cultivars, was developed in one region. Within the significantly associated genomic regions, a total of eight genes are considered as candidate genes, based on the presence of gene annotations associated with the protein or amino acid metabolism/movement. Our results provide better insights into utilizing wild soybean as a source of genetic diversity for soybean cultivar improvement utilizing native traits.

Genetic Augmentation of Legume Crops Using Genomic Resources and Genotyping Platforms for Nutritional Food Security

Recent advances in next generation sequencing (NGS) technologies have led the surge of genomic resources for the improvement legume crops. Advances in high throughput genotyping (HTG) and high throughput phenotyping (HTP) enable legume breeders to improve legume crops more precisely and efficiently. Now, the legume breeder can reshuffle the natural gene combinations of their choice to enhance the genetic potential of crops. These genomic resources are efficiently deployed through molecular breeding approaches for genetic augmentation of important legume crops, such as chickpea, cowpea, pigeonpea, groundnut, common bean, lentil, pea, as well as other underutilized legume crops. In the future, advances in NGS, HTG, and HTP technologies will help in the identification and assembly of superior haplotypes to tailor the legume crop varieties through haplotype-based breeding. This review article focuses on the recent development of genomic resource databases and their deployment in legume molecular breeding programmes to secure global food security.

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.

Exploring Machine Learning Algorithms to Unveil Genomic Regions Associated With Resistance to Southern Root-Knot Nematode in Soybeans

Southern root-knot nematode [SRKN, Meloidogyne incognita (Kofold & White) Chitwood] is a plant-parasitic nematode challenging to control due to its short life cycle, a wide range of hosts, and limited management options, of which genetic resistance is the main option to efficiently control the damage caused by SRKN. To date, a major quantitative trait locus (QTL) mapped on chromosome (Chr.) 10 plays an essential role in resistance to SRKN in soybean varieties. The confidence of discovered trait-loci associations by traditional methods is often limited by the assumptions of individual single nucleotide polymorphisms (SNPs) always acting independently as well as the phenotype following a Gaussian distribution. Therefore, the objective of this study was to conduct machine learning (ML)-based genome-wide association studies (GWAS) utilizing Random Forest (RF) and Support Vector Machine (SVM) algorithms to unveil novel regions of the soybean genome associated with resistance to SRKN. A total of 717 breeding lines derived from 330 unique bi-parental populations were genotyped with the Illumina Infinium BARCSoySNP6K BeadChip and phenotyped for SRKN resistance in a greenhouse. A GWAS pipeline involving a supervised feature dimension reduction based on Variable Importance in Projection (VIP) and SNP detection based on classification accuracy was proposed. Minor effect SNPs were detected by the proposed ML-GWAS methodology but not identified using Bayesian-information and linkage-disequilibrium Iteratively Nested Keyway (BLINK), Fixed and Random Model Circulating Probability Unification (FarmCPU), and Enriched Compressed Mixed Linear Model (ECMLM) models. Besides the genomic region on Chr. 10 that can explain most of SRKN resistance variance, additional minor effects SNPs were also identified on Chrs. 10 and 11. The findings in this study demonstrated that overfitting in GWAS may lead to lower prediction accuracy, and the detection of significant SNPs based on classification accuracy limited false-positive associations. The expansion of the basis of the genetic resistance to SRKN can potentially reduce the selection pressure over the major QTL on Chr. 10 and achieve higher levels of resistance.

Genetic and Genomic Resources for Soybean Breeding Research

Soybean (Glycine max) is a legume species of significant economic and nutritional value. The yield of soybean continues to increase with the breeding of improved varieties, and this is likely to continue with the application of advanced genetic and genomic approaches for breeding. Genome technologies continue to advance rapidly, with an increasing number of high-quality genome assemblies becoming available. With accumulating data from marker arrays and whole-genome resequencing, studying variations between individuals and populations is becoming increasingly accessible. Furthermore, the recent development of soybean pangenomes has highlighted the significant structural variation between individuals, together with knowledge of what has been selected for or lost during domestication and breeding, information that can be applied for the breeding of improved cultivars. Because of this, resources such as genome assemblies, SNP datasets, pangenomes and associated databases are becoming increasingly important for research underlying soybean crop improvement.

The SoyaGen Project: Putting Genomics to Work for Soybean Breeders

The SoyaGen project was a collaborative endeavor involving Canadian soybean researchers and breeders from academia and the private sector as well as international collaborators. Its aims were to develop genomics-derived solutions to real-world challenges faced by breeders. Based on the needs expressed by the stakeholders, the research efforts were focused on maximizing realized yield through optimization of maturity and improved disease resistance. The main deliverables related to molecular breeding in soybean will be reviewed here. These include: (1) SNP datasets capturing the genetic diversity within cultivated soybean (both within a worldwide collection of > 1,000 soybean accessions and a subset of 102 short-season accessions (MG0 and earlier) directly relevant to this group); (2) SNP markers for selecting favorable alleles at key maturity genes as well as loci associated with increased resistance to key pathogens and pests (Phytophthora sojae, Heterodera glycines, Sclerotinia sclerotiorum); (3) diagnostic tools to facilitate the identification and mapping of specific pathotypes of P. sojae; and (4) a genomic prediction approach to identify the most promising combinations of parents. As a result of this fruitful collaboration, breeders have gained new tools and approaches to implement molecular, genomics-informed breeding strategies. We believe these tools and approaches are broadly applicable to soybean breeding efforts around the world.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array

KEY MESSAGE

We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Quality Evaluation Indices for Soybean Oil in Relation to Cultivar, Application of N Fertiliser and Seed Inoculation with Bradyrhizobium japonicum

Soybean ranks second in production and consumption of vegetable oils worldwide and these are expected to continue to increase. The suitability of soybean oil for specific uses is determined by the fatty acid composition from which a number of indices and indicators can be calculated. The aim of this study was to evaluate the indices of nutritional and health-promoting fat in seeds of soybean cultivars grown in 2016–2019 under the influence of varying doses of N and inoculation with Bradyrhizobium japonicum. Omega 3 and Omega 6, unsaturated fatty acids (UFA), saturated fatty acids (SFA), polyunsaturated fatty acids (PUFA), index of desirable fatty acids (DFA), sum of hypercholesterolemic fatty acids (OFA), index of atherogenicity (AI), index of thrombogenicity (TI), oleic desaturation ratio (ODR), linoleic desaturation ratio (LDR), calculated oxidizability value (COX) and the hypocholesterolemic/hypercholesterolemic ratio (HH), saturation fat index (S/P) and ALA/LA, OL/(LA+ALA) ratios and the consumer index (CI) were included. Fat quality indices for soybean seeds were strongly determined by weather conditions. Seeds of the cv. Aldana contained higher amounts of Omega 6 and featured more favourable MUFA/PUFA and OL/(LA+ALA) ratios, while the seeds of the cv. Annushka had more favourable CI and higher ODR, COX and S/P indices. No important differences were observed regarding the effect of nitrogen dose and seed inoculation on the formation of the DFA, OFA, HH, AI, TI and CI indices. The value of the S/P index suggests that higher nitrogen rates (60 kg∙ha−1) and the lack of inoculation treatment produce seeds with a more favourable dietary fatty acid balance.

Progress in soybean functional genomics over the past decade

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

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

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

Phenotyping and Quantitative Trait Locus Analysis for the Limited Transpiration Trait in an Upper-Mid South Soybean Recombinant Inbred Line Population ("Jackson" × "KS4895"): High Throughput Aquaporin Inhibitor Screening

Soybean is most often grown under rainfed conditions and negatively impacted by drought stress in the upper mid-south of the United States. Therefore, identification of drought-tolerance traits and their corresponding genetic components are required to minimize drought impacts on productivity. Limited transpiration (TRlim) under high vapor pressure deficit (VPD) is one trait that can help conserve soybean water-use during late-season drought. The main research objective was to evaluate a recombinant inbred line (RIL) population, from crossing two mid-south soybean lines ("Jackson" × "KS4895"), using a high-throughput technique with an aquaporin inhibitor, AgNO3, for the TRlim trait. A secondary objective was to undertake a genetic marker/quantitative trait locus (QTL) genetic analysis using the AgNO3 phenotyping results. A set of 122 soybean genotypes (120-RILs and parents) were grown in controlled environments (32/25-d/n °C). The transpiration rate (TR) responses of derooted soybean shoots before and after application of AgNO3 were measured under 37°C and >3.0 kPa VPD. Then, the decrease in transpiration rate (DTR) for each genotype was determined. Based on DTR rate, a diverse group (slow, moderate, and high wilting) of 26 RILs were selected and tested for the whole plant TRs under varying levels of VPD (0.0-4.0 kPa) at 32 and 37°C. The phenotyping results showed that 88% of slow, 50% of moderate, and 11% of high wilting genotypes expressed the TRlim trait at 32°C and 43, 10, and 0% at 37°C, respectively. Genetic mapping with the phenotypic data we collected revealed three QTL across two chromosomes, two associated with TRlim traits and one associated with leaf temperature. Analysis of Gene Ontologies of genes within QTL regions identified several intriguing candidate genes, including one gene that when overexpressed had previously been shown to confer enhanced tolerance to abiotic stress. Collectively these results will inform and guide ongoing efforts to understand how to deploy genetic tolerance for drought stress.

Mapping Major Disease Resistance Genes in Soybean by Genome-Wide Association Studies

Soybean is one of the most valuable agricultural crops in the world. Besides, this legume is constantly attacked by a wide range of pathogens (fungi, bacteria, viruses, and nematodes) compromising yield and increasing production costs. One of the major disease management strategies is the genetic resistance provided by single genes and quantitative trait loci (QTL). Identifying the genomic regions underlying the resistance against these pathogens on soybean is one of the first steps performed by molecular breeders. In the past, genetic mapping studies have been widely used to discover these genomic regions. However, over the last decade, advances in next-generation sequencing technologies and their subsequent cost decreasing led to the development of cost-effective approaches to high-throughput genotyping. Thus, genome-wide association studies applying thousands of SNPs in large sets composed of diverse soybean accessions have been successfully done. In this chapter, a comprehensive review of the majority of GWAS for soybean diseases published since this approach was developed is provided. Important diseases caused by Heterodera glycines, Phytophthora sojae, and Sclerotinia sclerotiorum have been the focus of the several GWAS. However, other bacterial and fungi diseases also have been targets of GWAS. As such, this GWAS summary can serve as a guide for future studies of these diseases. The protocol begins by describing several considerations about the pathogens and bringing different procedures of molecular characterization of them. Advice to choose the best isolate/race to maximize the discovery of multiple R genes or to directly map an effective R gene is provided. A summary of protocols, methods, and tools to phenotyping the soybean panel is given to several diseases. We also give details of options of DNA extraction protocols and genotyping methods, and we describe parameters of SNP quality to soybean data. Websites and their online tools to obtain genotypic and phenotypic data for thousands of soybean accessions are highlighted. Finally, we report several tricks and tips in Subheading 4, especially related to composing the soybean panel as well as generating and analyzing the phenotype data. We hope this protocol will be helpful to achieve GWAS success in identifying resistance genes on soybean.

Identification of Candidate Genes and Genomic Selection for Seed Protein in Soybean Breeding Pipeline

Soybean is a primary meal protein for human consumption, poultry, and livestock feed. In this study, quantitative trait locus (QTL) controlling protein content was explored via genome-wide association studies (GWAS) and linkage mapping approaches based on 284 soybean accessions and 180 recombinant inbred lines (RILs), respectively, which were evaluated for protein content for 4 years. A total of 22 single nucleotide polymorphisms (SNPs) associated with protein content were detected using mixed linear model (MLM) and general linear model (GLM) methods in Tassel and 5 QTLs using Bayesian interval mapping (IM), single-trait multiple interval mapping (SMIM), single-trait composite interval mapping maximum likelihood estimation (SMLE), and single marker regression (SMR) models in Q-Gene and IciMapping. Major QTLs were detected on chromosomes 6 and 20 in both populations. The new QTL genomic region on chromosome 6 (Chr6_18844283-19315351) included 7 candidate genes and the Hap.X ^AA at the Chr6_19172961 position was associated with high protein content. Genomic selection (GS) of protein content was performed using Bayesian Lasso (BL) and ridge regression best linear unbiased prediction (rrBULP) based on all the SNPs and the SNPs significantly associated with protein content resulted from GWAS. The results showed that BL and rrBLUP performed similarly; GS accuracy was dependent on the SNP set and training population size. GS efficiency was higher for the SNPs derived from GWAS than random SNPs and reached a plateau when the number of markers was >2,000. The SNP markers identified in this study and other information were essential in establishing an efficient marker-assisted selection (MAS) and GS pipelines for improving soybean protein content.

Genotyping Platforms for Genome-Wide Association Studies: Options and Practical Considerations

Genome-wide association studies (GWAS) in crops requires genotyping platforms that are capable of producing accurate high density genotyping data on hundreds of plants in a cost-effective manner. Currently there are multiple commercial platforms available that are being effectively used across crops. These platforms include genotyping arrays such as the Illumina Infinium arrays and the Applied Biosystems Axiom Arrays along with a variety of resequencing methods. These methods are being used to genotype tens of thousands of markers up to millions of markers on GWAS panels. They are being used on crops with simple genomes to crops with very complex, large, polyploid genomes. Depending on the crop and the goal of the GWAS, there are several options and practical considerations to take into account when selecting a genotyping technology to ensure that the right coverage, accuracy, and cost for the study is achieved.

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

KEY MESSAGE

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

The Soybean High Density 'Forrest' by 'Williams 82' SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306); map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020); identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

Genome-wide association study and genomic selection for yield and related traits in soybean

Soybean [Glycine max (L.) Merr.] is a crop of great interest worldwide. Exploring molecular approaches to increase yield genetic gain has been one of the main challenges for soybean breeders and geneticists. Agronomic traits such as maturity, plant height, and seed weight have been found to contribute to yield. In this study, a total of 250 soybean accessions were genotyped with 10,259 high-quality SNPs postulated from genotyping by sequencing (GBS) and evaluated for grain yield, maturity, plant height, and seed weight over three years. A genome-wide association study (GWAS) was performed using a Bayesian Information and Linkage Disequilibrium Iteratively Nested Keyway (BLINK) model. Genomic selection (GS) was evaluated using a ridge regression best linear unbiased predictor (rrBLUP) model. The results revealed that 20, 31, 37, and 23 SNPs were significantly associated with maturity, plant height, seed weight, and yield, respectively; Many SNPs were mapped to previously described maturity and plant height loci (E2, E4, and Dt1) and a new plant height locus was mapped to chromosome 20. Candidate genes were found in the vicinity of the two SNPs with the highest significant levels associated with yield, maturity, plant height, seed weight, respectively. A 11.5-Mb region of chromosome 10 was associated with both yield and seed weight. Overall, the accuracy of GS was dependent on the trait, year, and population structure, and high accuracy indicates that these agronomic traits can be selected in molecular breeding through GS. The SNP markers identified in this study can be used to improve yield and agronomic traits through the marker-assisted selection and GS in breeding programs.