A Population Structure and Genome-Wide Association Analysis on the USDA Soybean Germplasm Collection

A novel genome-wide association study method for detecting quantitative trait loci interacting with complex population structures in plant genetics

In plant genetics, most modern association analyses are performed on panels that bring together individuals from several populations, including admixed individuals whose genomes comprise chromosomal regions from different populations. These panels can identify quantitative trait loci (QTLs) with population-specific effects and epistatic interactions between QTLs and polygenic backgrounds. However, analyzing a diverse panel constitutes a challenge for statistical analysis. The statistical model must account for possible interactions between a QTL and the panel structure while strictly controlling the detection error rate. Although models to detect population-specific QTLs have already been developed, they rely on prior information about the population structure. In practice, this prior information may be missing as many genome-wide association study (GWAS) panels exhibit complex population structures. The present study introduces 2 new models for detecting QTLs interacting with complex population structures. Both incorporate an interaction term between single nucleotide polymorphism/haplotype block and genetic background into conventional GWAS models. The proposed models were compared with state-of-the-art models through simulation studies that considered QTLs with different levels of interaction with their genetic backgrounds. Results showed that models matching simulation settings were most effective for detecting corresponding QTLs while the proposed models outperformed classical models in detecting QTLs interacting with polygenes. Additionally, when applied to a soybean dataset, one of our models identified putative associated QTLs that conventional models failed to detect. The new models, implemented in the RAINBOWR package available on CRAN, are expected to help uncover complex trait genetic architectures.

Soybean genome‑wide association study of seed weight, protein, and oil content in the southeastern USA

Soybean is a globally significant legume crop, providing essential protein and oil for human and livestock nutrition. Improving oil and protein content simultaneously without compromising yield has been challenging due to the quantitative nature of these traits and their interrelationships. This study aims to deepen our understanding of the molecular basis soybean of seed weight, protein, and oil content to facilitate marker-assisted breeding to enhance these traits. In this research, a Genome-Wide Association Study (GWAS) was conducted utilizing 285 diverse soybean accessions from maturity group V, employing genotyping through the SoySNP50K platform. These accessions were tested in three environmental conditions of the southeast US for three traits: 100-seed weight, protein, and oil content. The study identified 18, 23, and 26 SNPs significantly associated with 100-seed weight, seed oil, and protein content. Colocalized protein and oil content regions were discovered on chromosomes 15, 16, and 20. Chromosomes 15 and 20 are well documented to have pleiotropic but opposite effects on oil and protein content, but both regions contain genes that affect individual traits, such as FAD2-1 and nodulin MtN21. A 1.92 Mb region on chromosome 11 exhibits a target region to improve oil and seed weight without affecting protein content. This study highlights key genomic regions and candidate genes influencing seed weight, protein, and oil content, with some regions affecting multiple traits. Hence, these findings provide a valuable foundation for marker-assisted selection to optimize seed weight and simultaneously enhance oil and protein content in soybean breeding programs.

A genome-wide association study prioritizes VRN1-2 as a candidate gene associated with plant height in soybean

Plant height is an important architectural trait that affects crop growth, yield, and stress resistance. Tremendous efforts have been dedicated to revealing the genetic basis or regulatory mechanism; however, the underlying molecular mechanism remains largely unknown primarily due to the lack of controlling genes. In this study, we conducted a single-nucleotide resolution genome-wide association study (GWAS) of plant height using a diverse soybean panel collected worldwide with 6.7 million genome-wide variants (SNPs and Indels). The GWAS of plant height identified three QTLs on chromosomes 10, 18, and 19, of which the one on chromosome 19 precisely co-localized with Dt1, known as a major stem growth habit-controlling gene. Other loci without reported genes for plant height were regarded to be new. A close investigation within QTL intervals proposed nine genes that were likely involved in the regulation of plant height according to the expression specificity in developing shoot tip meristems. VRN1-2 underlies the significant QTL on chromosome 10 was prioritized as the most promising candidate gene. VRN1-2 shows higher expression in Williams 82 with indeterminate growth habit than Dongnong50 with semi-determinate growth habit across vegetative (V2, V3) and reproductive (R1) growth stages. VRN1-2 carries non-synonymous variants in the coding region that were significantly associated with plant height variation. The GT allele conferring short plant height was likely subjected to artificial selection during domestication. These results provide a source of new loci and genes for further elaborating the regulatory mechanism of plant height and the key variants would facilitate soybean molecular breeding.

Detection of genes associated with soybean protein content using a genome-wide association study

The protein content in soybean seeds serves as a crucial measure of soybean quality. Breeding high-protein varieties remains the most cost-effective and efficient approach to increasing soybean protein levels. Nevertheless, limited research has focused on identifying the genes responsible for high protein content among the diverse soybean cultivars. To address this gap, a genome-wide association study (GWAS) was conducted on 455 soybean varieties with varying protein content to predict and validate novel genes involved in regulating protein levels in soybean seeds. Protein content data were obtained from three distinct environments, along with three environmental variables derived from oil content, which is closely related to protein levels. Genotyping was performed using the SoySNP180k BeadChip, yielding genotype data for 63,306 non-redundant single nucleotide polymorphisms (SNPs). Five multi-locus GWAS methods were employed, resulting in the identification of 81 significant quantitative trait nucleotides (QTNs), of which 37 QTNs detected across different methods and environments were further analyzed. Moreover, the simulation platform Blib was used to conduct single-crossing simulation breeding on 81 QTN loci for actual breeding prediction. Haplotype analysis based on re-sequencing data confirmed 2 genes closely linked to protein synthesis, providing a theoretical basis for breeding high-protein soybean varieties and developing molecular breeding strategies.

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.

Historic manioc genomes illuminate maintenance of diversity under long-lived clonal cultivation

Manioc-also called cassava and yuca-is among the world's most important crops, originating in South America in the early Holocene. Domestication for its starchy roots involved a near-total shift from sexual to clonal propagation, and almost all manioc worldwide is now grown from stem cuttings. In this work, we analyze 573 new and published genomes, focusing on traditional varieties from the Americas and wild relatives from herbaria, to reveal the effects of this shift to clonality. We observe kinship over large distances, maintenance of high genetic diversity, intergenerational heterozygosity enrichment, and genomic mosaics of identity-by-descent haploblocks that connect all manioc worldwide. Interviews with Indigenous traditional farmers in the Brazilian Cerrado illuminate how traditional management strategies for sustaining, diversifying, and sharing the gene pool have shaped manioc diversity.

Genomic and phenomic prediction for soybean seed yield, protein, and oil

Developments in genomics and phenomics have provided valuable tools for use in cultivar development. Genomic prediction (GP) has been used in commercial soybean [Glycine max L. (Merr.)] breeding programs to predict grain yield and seed composition traits. Phenomic prediction (PP) is a rapidly developing field that holds the potential to be used for the selection of genotypes early in the growing season. The objectives of this study were to compare the performance of GP and PP for predicting soybean seed yield, protein, and oil. We additionally conducted genome-wide association studies (GWAS) to identify significant single-nucleotide polymorphisms (SNPs) associated with the traits of interest. The GWAS panel of 292 diverse accessions was grown in six environments in replicated trials. Spectral data were collected at two time points during the growing season. A genomic best linear unbiased prediction (GBLUP) model was trained on 269 accessions, while three separate machine learning (ML) models were trained on vegetation indices (VIs) and canopy traits. We observed that PP had a higher correlation coefficient than GP for seed yield, while GP had higher correlation coefficients for seed protein and oil contents. VIs with high feature importance were used as covariates in a new GBLUP model, and a new random forest model was trained with the inclusion of selected SNPs. These models did not outperform the original GP and PP models. These results show the capability of using ML for in-season predictions for specific traits in soybean breeding and provide insights on PP and GP inclusions in breeding programs.

Soybean2035: A decadal vision for soybean functional genomics and breeding

Soybean, the fourth most important crop in the world, uniquely serves as a source of both plant oil and plant protein for the world's food and animal feed. Although soybean production has increased approximately 13-fold over the past 60 years, the continually growing global population necessitates further increases in soybean production. In the past, especially in the last decade, significant progress has been made in both functional genomics and molecular breeding. However, many more challenges should be overcome to meet the anticipated future demand. Here, we summarize past achievements in the areas of soybean omics, functional genomics, and molecular breeding. Furthermore, we analyze trends in these areas, including shortages and challenges, and propose new directions, potential approaches, and possible outputs toward 2035. Our views and perspectives provide insight into accelerating the development of elite soybean varieties to meet the increasing demands of soybean production.

G2PDeep-v2: a web-based deep-learning framework for phenotype prediction and biomarker discovery for all organisms using multi-omics data

The G2PDeep-v2 server is a web-based platform powered by deep learning, for phenotype prediction and markers discovery from multi-omics data in any organisms including humans, plants, animals, and viruses. The server provides multiple services for researchers to create deep-learning models through an interactive interface and train these models using an automated hyperparameter tuning algorithm on high-performance computing resources. Users can visualize the results of phenotype and markers predictions and perform Gene Set Enrichment Analysis for the significant markers to provide insights into the molecular mechanisms underlying complex diseases, conditions and other biological phenotypes being studied. The G2PDeep-v2 server is publicly available at https://g2pdeep.org/ and can be utilized for all organisms.

Exploring SoySNP50K and USDA Germplasm Collection Data to Find New QTLs Associated with Protein and Oil Content in Brazilian Genotypes

Genetic diversity within a germplasm collection plays a vital role in the success of breeding programs. However, comprehending this diversity and identifying accessions with desirable traits pose significant challenges. This study utilized publicly available data to investigate SNP markers associated with protein and oil content in Brazilian soybeans. Through this research, twenty-two new QTLs (Quantitative Trait Loci) were identified, and we highlighted the substantial influence of Roanoke, Lee and Bragg ancestor on the genetic makeup of Brazilian soybean varieties. Our findings demonstrate that certain markers are being lost in modern cultivars, while others maintain or even increase their frequency. These observations indicate genomic regions that have undergone selection during soybean introduction in Brazil and could be valuable in breeding programs aimed at enhancing protein or oil content.

Identification of QTL-allele systems of seed size and oil content for simultaneous genomic improvement in Northeast China soybeans

Northeast China (NEC) is the major production area for soybeans in China, whereas its soybean germplasm has played key roles in world soybean production, especially in the Americas. For plant breeding, genomic selection involves two stages, cross design and progeny selection, with the former determining the latter's potential. In NEC, one of the major breeding purposes is for 100-seed weight (100SW) and seed oil content (SOC). A diverse sample with 361 NEC soybean germplasm accessions was evaluated for their 100SW and SOC in Tieling, Liaoning, China. Both traits exhibited significant phenotypic, genotypic, and G × E variation, with a trait heritability of 82.38% and 86.26%, respectively. A restricted two-stage multi-locus genome-wide association study (RTM-GWAS) with 15,501 SNPLDB (SNP linkage disequilibrium block) markers identified 80 and 92 QTLs with 230 and 299 alleles for 100SW and SOC, respectively. Corresponding to some increase of the two traits, almost all the alleles in the early maturity groups (MG 0 + 00 + 000) were inherited from the late MGs (MG I+II+III), indicating that genetic recombination was the major motivator in addition to a few allele emergence and some allele exclusion fluctuations among early MGs. Using the 95th percentile as indicator, the prediction of recombination potentials showed that 30.43 g 100SW and 27.73% SOC might be achieved, respectively. Three strategies of simultaneous genomic improvement of both traits in designing optimal crosses, namely, 100SW-first, SOC-first, and 100SW-SOC-balance, were proved to be efficient. Thus, the optimal cross design could be extended to multiple traits based on a relatively thorough identification of the QTL-alleles using RTM-GWAS.

Trait-customized sampling of core collections from a winter wheat genebank collection supports association studies

Subsampling a reduced number of accessions from ex situ genebank collections, known as core collections, is a widely applied method for the investigation of stored genetic diversity and for an exploitation by breeding and research. Optimizing core collections for genome-wide association studies could potentially maximize opportunities to discover relevant and rare variation. In the present study, eight strategies to sample core collections were implemented separately for two traits, namely susceptibility to yellow rust and stem lodging, on about 6,300 accessions of winter wheat (Triticum aestivum L.). Each strategy maximized different parameters or emphasized another aspect of the collection; the strategies relied on genomic data, phenotypic data or a combination thereof. The resulting trait-customized core collections of eight different sizes, covering the range between 100 and 800 accession samples, were analyzed based on characteristics such as population stratification, number of duplicate genotypes and genetic diversity. Furthermore, the statistical power for an association study was investigated as a key criterion for comparisons. While sampling extreme phenotypes boosts the power especially for smaller core collections of up to 500 accession samples, maximization of genetic diversity within the core collection minimizes population stratification and avoids the accumulation of less informative duplicate genotypes when increasing the size of a core collection. Advantages and limitations of different strategies to create trait-customized core collections are discussed for different scenarios of the availability of resources and data.

Deciphering the genetic basis of salinity tolerance in a diverse panel of cultivated and wild soybean accessions by genome-wide association mapping

KEY MESSAGE

In a genome-wide association study involving 269 cultivated and wild soybean accessions, potential salt tolerance donors were identified along with significant markers and candidate genes, such as GmKUP6 and GmWRKY33. Salt stress remains a significant challenge in agricultural systems, notably impacting soybean productivity worldwide. A comprehensive genome-wide association study (GWAS) was conducted to elucidate the genetic underpinnings of salt tolerance and identify novel source of salt tolerance among soybean genotypes. A diverse panel comprising 269 wild and cultivated soybean accessions was subjected to saline stress under controlled greenhouse conditions. Phenotypic data revealed that salt tolerance of soybean germplasm accessions was heavily compromised by the accumulation of sodium and chloride, as indicated by highly significant positive correlations of leaf scorching score with leaf sodium/chloride content. The GWAS analysis, leveraging a dataset of 32,832 SNPs, unveiled 32 significant marker-trait associations (MTAs) across seven traits associated with salt tolerance. These markers explained a substantial portion of the phenotypic variation, ranging from 14 to 52%. Notably, 11 markers surpassed Bonferroni's correction threshold, exhibiting highly significant associations with the respective traits. Gene Ontology enrichment analysis conducted within a 100 Kb range of the identified MTAs highlighted candidate genes such as potassium transporter 6 (GmKUP6), cation hydrogen exchanger (GmCHX15), and GmWRKY33. Expression levels of GmKUP6 and GmWRKY33 significantly varied between salt-tolerant and salt-susceptible soybean accessions under salt stress. The genetic markers and candidate genes identified in this study hold promise for developing soybean varieties resilient to salinity stress, thereby mitigating its adverse effects.

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 grain traits and genetic diversity between Chinese and Uruguayan soybeans (Glycine max L.)

Soybeans (Glycine max L.), originating in China, were introduced to South America in the late 19th century after passing through North America. South America is now a major soybean-producing region, accounting for approximately 40% of the global soybean production. Crops like soybeans gradually adapt to the local climate and human-selected conditions, resulting in beneficial variations during cultivation in different regions. Comparing the phenotypic and genetic variations in soybeans across different regions is crucial to determining the variations that may enhance soybean productivity. This study identified seed-related traits and conducted a genetic diversity analysis using 46 breeding soybean varieties from China and Uruguay. Compared to the Chinese soybean germplasm, the Uruguayan equivalent had a lower 100-grain weight, higher oil content, lower protein content, and higher soluble sugar content. Using ZDX1 gene chips, genetic typing was performed on the 46 breeding varieties. Cluster analysis based on SNP sites revealed significant differences in the genetic basis of Sino-Uruguayan soybean germplasm. Selection analysis, including nucleotide polymorphism (π) and fixation indexes (Fst), identified several genomic regions under selection between Sino-Uruguayan soybean germplasm. The selected intervals significantly enriched gene ontology (GO) terms related to protein metabolism. Additionally, differentiation occurred in genes associated with the oil content, seed weight, and cyst nematodes between Sino-Uruguayan soybean germplasm, such as GmbZIP123 and GmSSS1. These findings highlight the differences in seed-related phenotypes between Sino-Uruguay soybean germplasm and provide genomic-level insights into the mechanisms behind phenotypic differences, offering valuable references for understanding soybean evolution and molecular breeding.

Across-environment seed protein stability and genetic architecture of seed components in soybean

The recent surge in the plant-based protein market has resulted in high demands for soybean genotypes with improved grain yield, seed protein and oil content, and essential amino acids (EAAs). Given the quantitative nature of these traits, complex interactions among seed components, as well as between seed components and environmental factors and management practices, add complexity to the development of desired genotypes. In this study, the across-environment seed protein stability of 449 genetically diverse plant introductions was assessed, revealing that genotypes may display varying sensitivities to such environmental stimuli. The EAAs valine, phenylalanine, and threonine showed the highest variable importance toward the variation in stability, while both seed protein and oil contents were among the explanatory variables with the lowest importance. In addition, 56 single nucleotide polymorphism (SNP) markers were significantly associated with various seed components. Despite the strong phenotypic Pearson's correlation observed among most seed components, many independent genomic regions associated with one or few seed components were identified. These findings provide insights for improving the seed concentration of specific EAAs and reducing the negative correlation between seed protein and oil contents.

Effective population size in field pea

BACKGROUND

Effective population size (Ne) is a pivotal parameter in population genetics as it can provide information on the rate of inbreeding and the contemporary status of genetic diversity in breeding populations. The population with smaller Ne can lead to faster inbreeding, with little potential for genetic gain making selections ineffective. The importance of Ne has become increasingly recognized in plant breeding, which can help breeders monitor and enhance the genetic variability or redesign their selection protocols. Here, we present the first Ne estimates based on linkage disequilibrium (LD) in the pea genome.

RESULTS

We calculated and compared Ne using SNP markers from North Dakota State University (NDSU) modern breeding lines and United States Department of Agriculture (USDA) diversity panel. The extent of LD was highly variable not only between populations but also among different regions and chromosomes of the genome. Overall, NDSU had a higher and longer-range LD than the USDA that could extend up to 500 Kb, with a genome-wide average r2 of 0.57 (vs 0.34), likely due to its lower recombination rates and the selection background. The estimated Ne for the USDA was nearly three-fold higher (Ne = 174) than NDSU (Ne = 64), which can be confounded by a high degree of population structure due to the selfing nature of pea.

CONCLUSIONS

Our results provided insights into the genetic diversity of the germplasm studied, which can guide plant breeders to actively monitor Ne in successive cycles of breeding to sustain viability of the breeding efforts in the long term.

Construction of Genetic Map and QTL Mapping for Seed Size and Quality Traits in Soybean (Glycine max L.)

Soybean (Glycine max L.) is the main source of vegetable protein and edible oil for humans, with an average content of about 40% crude protein and 20% crude fat. Soybean yield and quality traits are mostly quantitative traits controlled by multiple genes. The quantitative trait loci (QTL) mapping for yield and quality traits, as well as for the identification of mining-related candidate genes, is of great significance for the molecular breeding and understanding the genetic mechanism. In this study, 186 individual plants of the F2 generation derived from crosses between Changjiangchun 2 and Yushuxian 2 were selected as the mapping population to construct a molecular genetic linkage map. A genetic map containing 445 SSR markers with an average distance of 5.3 cM and a total length of 2375.6 cM was obtained. Based on constructed genetic map, 11 traits including hundred-seed weight (HSW), seed length (SL), seed width (SW), seed length-to-width ratio (SLW), oil content (OIL), protein content (PRO), oleic acid (OA), linoleic acid (LA), linolenic acid (LNA), palmitic acid (PA), stearic acid (SA) of yield and quality were detected by the multiple- d size traits and 113 QTLs related to quality were detected by the multiple QTL model (MQM) mapping method across generations F2, F2:3, F2:4, and F2:5. A total of 71 QTLs related to seed size traits and 113 QTLs related to quality traits were obtained in four generations. With those QTLs, 19 clusters for seed size traits and 20 QTL clusters for quality traits were summarized. Two promising clusters, one related to seed size traits and the other to quality traits, have been identified. The cluster associated with seed size traits spans from position 27876712 to 29009783 on Chromosome 16, while the cluster linked to quality traits spans from position 12575403 to 13875138 on Chromosome 6. Within these intervals, a reference genome of William82 was used for gene searching. A total of 36 candidate genes that may be involved in the regulation of soybean seed size and quality were screened by gene functional annotation and GO enrichment analysis. The results will lay the theoretical and technical foundation for molecularly assisted breeding in soybean.

Introgression of a Danbaekkong high-protein allele across different genetic backgrounds in soybean

Soybean meal is a major component of livestock feed due to its high content and quality of protein. Understanding the genetic control of protein is essential to develop new cultivars with improved meal protein. Previously, a genomic region on chromosome 20 significantly associated with elevated protein content was identified in the cultivar Danbaekkong. The present research aimed to introgress the Danbaekkong high-protein allele into elite lines with different genetic backgrounds by developing and deploying robust DNA markers. A multiparent population consisting of 10 F5-derived populations with a total of 1,115 recombinant inbred lines (RILs) was developed using "Benning HP" as the donor parent of the Danbaekkong high-protein allele. A new functional marker targeting the 321-bp insertion in the gene Glyma.20g085100 was developed and used to track the Danbaekkong high-protein allele across the different populations and enable assessment of its effect and stability. Across all populations, the high-protein allele consistently increased the content, with an increase of 3.3% in seed protein. A total of 103 RILs were selected from the multiparent population for yield testing in five environments to assess the impact of the high-protein allele on yield and to enable the selection of new breeding lines with high protein and high yield. The results indicated that the high-protein allele impacts yield negatively in general; however, it is possible to select high-yielding lines with high protein content. An analysis of inheritance of the Chr 20 high-protein allele in Danbaekkong indicated that it originated from a Glycine soja line (PI 163453) and is the same as other G. soja lines studied. A survey of the distribution of the allele across 79 G. soja accessions and 35 Glycine max ancestors of North American soybean cultivars showed that the high-protein allele is present in all G. soja lines evaluated but not in any of the 35 North American soybean ancestors. These results demonstrate that G. soja accessions are a valuable source of favorable alleles for improvement of protein composition.

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.

Regulation of seed traits in soybean

Soybean (Glycine max) is an essential economic crop that provides vegetative oil and protein for humans, worldwide. Increasing soybean yield as well as improving seed quality is of great importance. Seed weight/size, oil and protein content are the three major traits determining seed quality, and seed weight also influences soybean yield. In recent years, the availability of soybean omics data and the development of related techniques have paved the way for better research on soybean functional genomics, providing a comprehensive understanding of gene functions. This review summarizes the regulatory genes that influence seed size/weight, oil content and protein content in soybean. We also provided a general overview of the pleiotropic effect for the genes in controlling seed traits and environmental stresses. Ultimately, it is expected that this review will be beneficial in breeding improved traits in soybean.

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.

k-mer-based GWAS enhances the discovery of causal variants and candidate genes in soybean

Genome-wide association studies (GWAS) are powerful statistical methods that detect associations between genotype and phenotype at genome scale. Despite their power, GWAS frequently fail to pinpoint the causal variant or the gene controlling a given trait in crop species. Assessing genetic variants other than single-nucleotide polymorphisms (SNPs) could alleviate this problem. In this study, we tested the potential of structural variant (SV)- and k-mer-based GWAS in soybean by applying these methods as well as conventional SNP/indel-based GWAS to 13 traits. We assessed the performance of each GWAS approach based on loci for which the causal genes or variants were known from previous genetic studies. We found that k-mer-based GWAS was the most versatile approach and the best at pinpointing causal variants or candidate genes. Moreover, k-mer-based analyses identified promising candidate genes for loci related to pod color, pubescence form, and resistance to Phytophthora sojae. In our dataset, SV-based GWAS did not add value compared to k-mer-based GWAS and may not be worth the time and computational resources invested. Despite promising results, significant challenges remain regarding the downstream analysis of k-mer-based GWAS. Notably, better methods are needed to associate significant k-mers with sequence variation. Our results suggest that coupling k-mer- and SNP/indel-based GWAS is a powerful approach for discovering candidate genes in crop species.

Evaluating the Response of Glycine soja Accessions to Fungal Pathogen Macrophomina phaseolina during Seedling Growth

Charcoal rot caused by the fungal pathogen Macrophomina phaseolina (Tassi) Goid is one of various devastating soybean (Glycine max (L.) Merr.) diseases, which can severely reduce crop yield. The investigation into the genetic potential for charcoal rot resistance of wild soybean (Glycine soja) accessions will enrich our understanding of the impact of soybean domestication on disease resistance; moreover, the identified charcoal rot-resistant lines can be used to improve soybean resistance to charcoal rot. The objective of this study was to evaluate the resistance of wild soybean accessions to M. phaseolina at the seedling stage and thereby select the disease-resistant lines. The results show that the fungal pathogen infection reduced the growth of the root and hypocotyl in most G. soja accessions. The accession PI 507794 displayed the highest level of resistance response to M. phaseolina infection among the tested wild soybean accessions, while PI 487431 and PI 483660B were susceptible to charcoal rot in terms of the reduction in root and hypocotyl growth. The mean values of the root and hypocotyl parameters in PI 507794 were significantly higher (p < 0.05) than those of PI 487431 and PI 483460B. A analysis of the resistance of wild soybean accessions to M. phaseolina using the root and hypocotyl as the assessment parameters at the early seedling stage provides an alternative way to rapidly identify potential resistant genotypes and facilitate breeding for soybean resistance to charcoal rot.

Novel genetic resources associated with sucrose and stachyose content through genome-wide association study in soybean (Glycine max (L.) Merr.)

The nutritional value of soybean [Glycine max (L.) Merr.] for animals is influenced by soluble carbohydrates, such as sucrose and stachyose. Although sucrose is nutritionally desirable, stachyose is an antinutrient causing diarrhea and flatulence in non-ruminant animals. We conducted a genome-wide association study of 220 soybean accessions using 21,317 single nucleotide polymorphisms (SNPs) from the SoySNP50K iSelect Beadchip data to identify significant SNPs associated with sucrose and stachyose content. Seven significant SNPs were identified for sucrose content across chromosomes (Chrs.) 2, 8, 12, 17, and 20, while thirteen significant SNPs were identified for stachyose content across Chrs. 2, 5, 8, 9, 10, 13, 14, and 15. Among those significant SNPs, three sucrose-related SNPs on Chrs. 8 and 17 were novel, while twelve stachyose-related SNPs on Chrs. 2, 5, 8, 9, 10, 13, 14, and 15 were novel. Based on Phytozome, STRING, and GO annotation, 17 and 24 candidate genes for sucrose and stachyose content, respectively, were highly associated with the carbohydrate metabolic pathway. Among these, the publicly available RNA-seq Atlas database highlighted four candidate genes associated with sucrose (Glyma.08g361200 and Glyma.17g258100) and stachyose (Glyma.05g025300 and Glyma.13g077900) content, which had higher gene expression levels in developing seed and multiple parts of the soybean plant. The results of this study will extend knowledge of the molecular mechanism and genetic basis underlying sucrose and stachyose content in soybean seed. Furthermore, the novel candidate genes and SNPs can be valuable genetic resources that soybean breeders may utilize to modify carbohydrate profiles for animal and human usage.

QTLs and Candidate Genes for Seed Protein Content in Two Recombinant Inbred Line Populations of Soybean

This study aimed to discover the quantitative trait loci (QTL) associated with a high seed protein content in soybean and unravel the potential candidate genes. We developed two recombinant inbred line populations: YS and SI, by crossing Saedanbaek (high protein) with YS2035-B-91-1-B-1 (low protein) and Saedanbaek with Ilmi (low protein), respectively, and evaluated the protein content for three consecutive years. Using single-nucleotide polymorphism (SNP)-marker-based linkage maps, four QTLs were located on chromosomes 15, 18, and 20 with high logarithm of odds values (5.9-55.0), contributing 5.5-66.0% phenotypic variance. In all three experimental years, qPSD20-1 and qPSD20-2 were stable and identified in overlapping positions in the YS and SI populations, respectively. Additionally, novel QTLs were identified on chromosomes 15 and 18. Considering the allelic sequence variation between parental lines, 28 annotated genes related to soybean seed protein-including starch, lipid, and fatty acid biosynthesis-related genes-were identified within the QTL regions. These genes could potentially affect protein accumulation during seed development, as well as sucrose and oil metabolism. Overall, this study offers insights into the genetic mechanisms underlying a high soybean protein content. The identified potential candidate genes can aid marker-assisted selection for developing soybean lines with an increased protein content.

Population Structure and Genetic Diversity of the 175 Soybean Breeding Lines and Varieties Cultivated in West Siberia and Other Regions of Russia

Soybean is a leguminous plant cultivated in many countries and is considered important in the food industry due to the high levels of oil and protein content in the beans. The high demand for soybeans and its products in the industry requires the expansion of cultivation areas. Despite climatic restrictions, West Siberia is gradually expanding its area of soybean cultivation. In this study, we present the first analysis of the population structure and genetic diversity of the 175 soybean Glycine max breeding lines and varieties cultivated in West Siberia (103 accessions) and other regions of Russia (72 accessions), and we compare them with the cultivated soybean varieties from other geographical locations. Principal component analysis revealed several genetic clusters with different levels of genetic heterogeneity. Studied accessions are genetically similar to varieties from China, Japan, and the USA and are genetically distant to varieties from South Korea. Admixture analysis revealed four ancestry groups based on genetic ancestry and geographical origin, which are consistent with the regions of cultivation and origin of accessions and correspond to the principal component analysis result. Population statistics, including nucleotide diversity, Tajima's D, and linkage disequilibrium, are comparatively similar to those observed for studied accessions of a different origin. This study provides essential population and genetic information about the unique collection of breeding lines and varieties cultivated in West Siberia and other Russian regions to foster further evolutionary, genome-wide associations and functional breeding studies.

Haplotype blocks for genomic prediction: a comparative evaluation in multiple crop datasets

In modern plant breeding, genomic selection is becoming the gold standard for selection of superior genotypes. The basis for genomic prediction models is a set of phenotyped lines along with their genotypic profile. With high marker density and linkage disequilibrium (LD) between markers, genotype data in breeding populations tends to exhibit considerable redundancy. Therefore, interest is growing in the use of haplotype blocks to overcome redundancy by summarizing co-inherited features. Moreover, haplotype blocks can help to capture local epistasis caused by interacting loci. Here, we compared genomic prediction methods that either used single SNPs or haplotype blocks with regards to their prediction accuracy for important traits in crop datasets. We used four published datasets from canola, maize, wheat and soybean. Different approaches to construct haplotype blocks were compared, including blocks based on LD, physical distance, number of adjacent markers and the algorithms implemented in the software "Haploview" and "HaploBlocker". The tested prediction methods included Genomic Best Linear Unbiased Prediction (GBLUP), Extended GBLUP to account for additive by additive epistasis (EGBLUP), Bayesian LASSO and Reproducing Kernel Hilbert Space (RKHS) regression. We found improved prediction accuracy in some traits when using haplotype blocks compared to SNP-based predictions, however the magnitude of improvement was very trait- and model-specific. Especially in settings with low marker density, haplotype blocks can improve genomic prediction accuracy. In most cases, physically large haplotype blocks yielded a strong decrease in prediction accuracy. Especially when prediction accuracy varies greatly across different prediction models, prediction based on haplotype blocks can improve prediction accuracy of underperforming models. However, there is no "best" method to build haplotype blocks, since prediction accuracy varied considerably across methods and traits. Hence, criteria used to define haplotype blocks should not be viewed as fixed biological parameters, but rather as hyperparameters that need to be adjusted for every dataset.

Genome-wide association study of soybean (Glycine max [L.] Merr.) germplasm for dissecting the quantitative trait nucleotides and candidate genes underlying yield-related traits

Soybean (Glycine max [L.] Merr.) is one of the most significant crops in the world in terms of oil and protein. Owing to the rising demand for soybean products, there is an increasing need for improved varieties for more productive farming. However, complex correlation patterns among quantitative traits along with genetic interactions pose a challenge for soybean breeding. Association studies play an important role in the identification of accession with useful alleles by locating genomic sites associated with the phenotype in germplasm collections. In the present study, a genome-wide association study was carried out for seven agronomic and yield-related traits. A field experiment was conducted in 2015/2016 at two locations that include 155 diverse soybean germplasm. These germplasms were genotyped using SoySNP50K Illumina Infinium Bead-Chip. A total of 51 markers were identified for node number, plant height, pods per plant, seeds per plant, seed weight per plant, hundred-grain weight, and total yield using a multi-locus linear mixed model (MLMM) in FarmCPU. Among these significant SNPs, 18 were putative novel QTNs, while 33 co-localized with previously reported QTLs. A total of 2,356 genes were found in 250 kb upstream and downstream of significant SNPs, of which 17 genes were functional and the rest were hypothetical proteins. These 17 candidate genes were located in the region of 14 QTNs, of which ss715580365, ss715608427, ss715632502, and ss715620131 are novel QTNs for PH, PPP, SDPP, and TY respectively. Four candidate genes, Glyma.01g199200, Glyma.10g065700, Glyma.18g297900, and Glyma.14g009900, were identified in the vicinity of these novel QTNs, which encode lsd one like 1, Ergosterol biosynthesis ERG4/ERG24 family, HEAT repeat-containing protein, and RbcX2, respectively. Although further experimental validation of these candidate genes is required, several appear to be involved in growth and developmental processes related to the respective agronomic traits when compared with their homologs in Arabidopsis thaliana. This study supports the usefulness of association studies and provides valuable data for functional markers and investigating candidate genes within a diverse germplasm collection in future breeding programs.

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.

The genetic architecture of soybean photothermal adaptation to high latitudes

Soybean is a major plant protein source for both human food and animal feed, but to meet global demands as well as a trend towards regional production, soybean cultivation needs to be expanded to higher latitudes. In this study, we developed a large diversity panel consisting of 1503 early-maturing soybean lines and used genome-wide association mapping to dissect the genetic architecture underlying two crucial adaptation traits, flowering time and maturity. This revealed several known maturity loci, E1, E2, E3, and E4, and the growth habit locus Dt2 as causal candidate loci, and also a novel putative causal locus, GmFRL1, encoding a homolog of the vernalization pathway gene FRIGIDA-like 1. In addition, the scan for quantitative trait locus (QTL)-by-environment interactions identified GmAPETALA1d as a candidate gene for a QTL with environment-dependent reversed allelic effects. The polymorphisms of these candidate genes were identified using whole-genome resequencing data of 338 soybeans, which also revealed a novel E4 variant, e4-par, carried by 11 lines, with nine of them originating from Central Europe. Collectively, our results illustrate how combinations of QTL and their interactions with the environment facilitate the photothermal adaptation of soybean to regions far beyond its center of origin.

Genetic Diversity and Population Structure of European Soybean Germplasm Revealed by Single Nucleotide Polymorphism

Soybean is the most grown high-protein crop in the world. Despite the rapid increase of acreage and production volume, European soybean production accounts for only 34% of its consumption in Europe. This study aims to support the optimal exploitation of genetic resources by European breeding programs by investigating the genetic diversity and the genetic structure of 207 European cultivars or American introductions registered in Europe, which were genotyped by the SoySNP50K array. The expected heterozygosity (H_e) was 0.34 for the entire collection and ranged among countries from 0.24 for Swiss cultivars to 0.32 for American cultivars (partly reflecting differences in sample size between countries). Cluster analysis grouped all genotypes into two main clusters with eight subgroups that corresponded to the country of origin of cultivars and their maturity group. Pairwise F_st values between countries of origin showed the highest differentiation of Swiss cultivars from the rest of the European gene pool, while the lowest mean differentiation was found between American introductions and all other European countries. On the other hand, F_st values between maturity groups were much lower compared to those observed between countries. In analysis of molecular variance, the total genetic variation was partitioned either by country of origin or by maturity group, explaining 9.1% and 3.5% of the total genetic variance, respectively. On the whole, our results suggest that the European soybean gene pool still has sufficient diversity due to the different historical breeding practices in western and eastern countries and the relatively short period of breeding in Europe.

Genetic mapping and functional genomics of soybean seed protein

Soybean is an utterly important crop for high-quality meal protein and vegetative oil. Soybean seed protein content has become a key factor in nutrients for livestock feed as well as human dietary consumption. Genetic improvement of soybean seed protein is highly desired to meet the demands of rapidly growing world population. Molecular mapping and genomic analysis in soybean have identified many quantitative trait loci (QTL) underlying seed protein content control. Exploring the mechanisms of seed storage protein regulation will be helpful to achieve the improvement of protein content. However, the practice of breeding higher protein soybean is challenging because soybean seed protein is negatively correlated with seed oil content and yield. To overcome the limitation of such inverse relationship, deeper insights into the property and genetic control of seed protein are required. Recent advances of soybean genomics have strongly enhanced the understandings for molecular mechanisms of soybean with better seed quality. Here, we review the research progress in the genetic characteristics of soybean storage protein, and up-to-date advances of molecular mappings and genomics of soybean protein. The key factors underlying the mechanisms of the negative correlation between protein and oil in soybean seeds are elaborated. We also briefly discuss the future prospects of breaking the bottleneck of the negative correlation to develop high protein soybean without penalty of oil and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01373-5.

Branch angle and leaflet shape are associated with canopy coverage in soybean

Early canopy coverage is a desirable trait that is a major determinant of yield in soybean (Glycine max). Variation in traits comprising shoot architecture can influence canopy coverage, canopy light interception, canopy-level photosynthesis, and source-sink partitioning efficiency. However, little is known about the extent of phenotypic diversity of shoot architecture traits and their genetic control in soybean. Thus, we sought to understand the contribution of shoot architecture traits to canopy coverage and to determine the genetic control of these traits. We examined the natural variation for shoot architecture traits in a set of 399 diverse maturity group I soybean (SoyMGI) accessions to identify relationships between traits, and to identify loci that are associated with canopy coverage and shoot architecture traits. Canopy coverage was correlated with branch angle, number of branches, plant height, and leaf shape. Using previously collected 50K single nucleotide polymorphism data, we identified quantitative trait locus (QTL) associated with branch angle, number of branches, branch density, leaflet shape, days to flowering, maturity, plant height, number of nodes, and stem termination. In many cases, QTL intervals overlapped with previously described genes or QTL. We also found QTL associated with branch angle and leaflet shape located on chromosomes 19 and 4, respectively, and these QTL overlapped with QTL associated with canopy coverage, suggesting the importance of branch angle and leaflet shape in determining canopy coverage. Our results highlight the role individual architecture traits play in canopy coverage and contribute information on their genetic control that could help facilitate future efforts in their genetic manipulation.

Genome-Wide Detection of Quantitative Trait Loci and Prediction of Candidate Genes for Seed Sugar Composition in Early Mature Soybean

Seed sugar composition, mainly including fructose, glucose, sucrose, raffinose, and stachyose, is an important indicator of soybean [Glycine max (L.) Merr.] seed quality. However, research on soybean sugar composition is limited. To better understand the genetic architecture underlying the sugar composition in soybean seeds, we conducted a genome-wide association study (GWAS) using a population of 323 soybean germplasm accessions which were grown and evaluated under three different environments. A total of 31,245 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) ≥ 5% and missing data ≤ 10% were selected and used in the GWAS. The analysis identified 72 quantitative trait loci (QTLs) associated with individual sugars and 14 with total sugar. Ten candidate genes within the 100 Kb flanking regions of the lead SNPs across six chromosomes were significantly associated with sugar contents. According to GO and KEGG classification, eight genes were involved in the sugar metabolism in soybean and showed similar functions in Arabidopsis. The other two, located in known QTL regions associated with sugar composition, may play a role in sugar metabolism in soybean. This study advances our understanding of the genetic basis of soybean sugar composition and facilitates the identification of genes controlling this trait. The identified candidate genes will help improve seed sugar composition in soybean.

Genome-wide association analysis of sucrose and alanine contents in edamame beans

The sucrose and Alanine (Ala) content in edamame beans significantly impacts the sweetness flavor of edamame-derived products as an important attribute to consumers' acceptance. Unlike grain-type soybeans, edamame beans are harvested as fresh beans at the R6 to R7 growth stages when beans are filled 80-90% of the pod capacity. The genetic basis of sucrose and Ala contents in fresh edamame beans may differ from those in dry seeds. To date, there is no report on the genetic basis of sucrose and Ala contents in the edamame beans. In this study, a genome-wide association study was conducted to identify single nucleotide polymorphisms (SNPs) related to sucrose and Ala levels in edamame beans using an association mapping panel of 189 edamame accessions genotyped with a SoySNP50K BeadChip. A total of 43 and 25 SNPs was associated with sucrose content and Ala content in the edamame beans, respectively. Four genes (Glyma.10g270800, Glyma.08g137500, Glyma.10g268500, and Glyma.18g193600) with known effects on the process of sucrose biosynthesis and 37 novel sucrose-related genes were characterized. Three genes (Gm17g070500, Glyma.14g201100 and Glyma.18g269600) with likely relevant effects in regulating Ala content and 22 novel Ala-related genes were identified. In addition, by summarizing the phenotypic data of edamame beans from three locations in two years, three PI accessions (PI 532469, PI 243551, and PI 407748) were selected as the high sucrose and high Ala parental lines for the perspective breeding of sweet edamame varieties. Thus, the beneficial alleles, candidate genes, and selected PI accessions identified in this study will be fundamental to develop edamame varieties with improved consumers' acceptance, and eventually promote edamame production as a specialty crop in the United States.

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.

Soybean Protein and Oil Variants Identified through a Forward Genetic Screen for Seed Composition

Mutagenesis remains an important tool in soybean biology. In classical plant mutation breeding, mutagenesis has been a trusted approach for decades, creating stable non-transgenic variation, and many mutations have been incorporated into germplasm for several crops, especially to introduce favorable seed composition traits. We performed a genetic screen for aberrant oil or protein composition of soybean seeds, and as a result isolated over 100 mutant lines for seed composition phenotypes, with particular interest in high protein or high oil phenotypes. These lines were followed for multiple seasons and generations to select the most stable traits for further characterization. Through backcrossing and outcrossing experiments, we determined that a subset of the lines showed recessive inheritance, while others showed a dominant inheritance pattern that suggests the involvement of multiple loci and genetic mechanisms. These lines can be used as a resource for future studies of the genetic control of seed protein and oil content in soybean.

Soybean genetic resources contributing to sustainable protein production

KEY MESSAGE

Genetic resources contributes to the sustainable protein production in soybean. Soybean is an important crop for food, oil, and forage and is the main source of edible vegetable oil and vegetable protein. It plays an important role in maintaining balanced dietary nutrients for human health. The soybean protein content is a quantitative trait mainly controlled by gene additive effects and is usually negatively correlated with agronomic traits such as the oil content and yield. The selection of soybean varieties with high protein content and high yield to secure sustainable protein production is one of the difficulties in soybean breeding. The abundant genetic variation of soybean germplasm resources is the basis for overcoming the obstacles in breeding for soybean varieties with high yield and high protein content. Soybean has been cultivated for more than 5000 years and has spread from China to other parts of the world. The rich genetic resources play an important role in promoting the sustainable production of soybean protein worldwide. In this paper, the origin and spread of soybean and the current status of soybean production are reviewed; the genetic characteristics of soybean protein and the distribution of resources are expounded based on phenotypes; the discovery of soybean seed protein-related genes as well as transcriptomic, metabolomic, and proteomic studies in soybean are elaborated; the creation and utilization of high-protein germplasm resources are introduced; and the prospect of high-protein soybean breeding is described.

Genomic regions associated with resistance to soybean rust (Phakopsora pachyrhizi) under field conditions in soybean germplasm accessions from Japan, Indonesia and Vietnam

Eight soybean genomic regions, including six never before reported, were found to be associated with resistance to soybean rust (Phakopsora pachyrhizi) in the southeastern USA. Soybean rust caused by Phakopsora pachyrhizi is one of the most important foliar diseases of soybean [Glycine max (L.) Merr.]. Although seven Rpp resistance gene loci have been reported, extensive pathotype variation in and among fungal populations increases the importance of identifying additional genes and loci associated with rust resistance. One hundred and ninety-one soybean plant introductions from Japan, Indonesia and Vietnam, and 65 plant introductions from other countries were screened for resistance to P. pachyrhizi under field conditions in the southeastern USA between 2008 and 2015. The results indicated that 84, 69, and 49% of the accessions from southern Japan, Vietnam or central Indonesia, respectively, had negative BLUP values, indicating less disease than the panel mean. A genome-wide association analysis using SoySNP50K Infinium BeadChip data identified eight genomic regions on seven chromosomes associated with SBR resistance, including previously unreported regions of Chromosomes 1, 4, 6, 9, 13, and 15, in addition to the locations of the Rpp3 and Rpp6 loci. The six unreported genomic regions might contain novel Rpp loci. The identification of additional sources of rust resistance and associated genomic regions will further efforts to develop soybean cultivars with broad and durable resistance to soybean rust in the southern USA.

Genome-Wide Identification of Candidate Genes Underlying Soluble Sugar Content in Vegetable Soybean (Glycine max L.) via Association and Expression Analysis

Soluble sugar is a major indicator of the intrinsic quality of vegetable soybean [Glycine max (L.) Merr. ]. The improvement of soluble sugar content in soybean is very important due to its healthcare functions for humans. The genetic mechanism of soluble sugar in soybean is unclear. In this study, 278 diverse soybean accessions were utilized to identify the quantitative trait nucleotides (QTNs) for total soluble sugar content in soybean seeds based on a genome-wide association study (GWAS). A total of 25,921 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) ≥ 5% and missing data ≤ 10% were selected for GWAS. Totally, thirteen QTNs associated with total soluble sugar content were identified, which were distributed on ten chromosomes. One hundred and fifteen genes near the 200-kb flanking region of these identified QTNs were considered candidate genes associated with total soluble sugar content in soybean seed. Gene-based association analysis and haplotype analysis were utilized to further identify the effect of candidate genes on total soluble sugar content. Totally, 84 SNPs from seventeen genes across four chromosomes were significantly associated with the total soluble sugar content. Among them, three SNPs from Glyma.02G292900 were identified at two locations, and other eighty-one SNPs from sixteen genes were detected at three locations. Furthermore, expression level analysis of candidate genes revealed that Glyma.02G293200 and Glyma.02G294900 were significantly positively associated with soluble sugar content and Glyma.02G294000 was significantly negatively associated with soluble sugar content. Six genes (i.e., Glyma.02G292600, Glyma.02G292700, Glyma.02G294000, Glyma.02G294300, Glyma.02G294400, and Glyma.15G264200) identified by GWAS were also detected by the analysis of differential expression genes based on soybean germplasms with higher and lower soluble sugar content. Among them, Glyma.02G294000 is the only gene that was identified by gene-based association analysis with total soluble sugar content and was considered an important candidate gene for soluble sugar content. These candidate genes and beneficial alleles would be useful for improving the soluble sugar content of soybean.

Transgressive Potential Prediction and Optimal Cross Design of Seed Protein Content in the Northeast China Soybean Population Based on Full Exploration of the QTL-Allele System

Northeast China is a major soybean production region in China. A representative sample of the Northeast China soybean germplasm population (NECSGP) composed of 361 accessions was evaluated for their seed protein content (SPC) in Tieling, Northeast China. This SPC varied greatly, with a mean SPC of 40.77%, ranging from 36.60 to 46.07%, but it was lower than that of the Chinese soybean landrace population (43.10%, ranging from 37.51 to 50.46%). The SPC increased slightly from 40.32-40.97% in the old maturity groups (MG, MGIII + II + I) to 40.93-41.58% in the new MGs (MG0 + 00 + 000). The restricted two-stage multi-locus genome-wide association study (RTM-GWAS) with 15,501 SNP linkage-disequilibrium block (SNPLDB) markers identified 73 SPC quantitative trait loci (QTLs) with 273 alleles, explaining 71.70% of the phenotypic variation, wherein 28 QTLs were new ones. The evolutionary changes of QTL-allele structures from old MGs to new MGs were analyzed, and 97.79% of the alleles in new MGs were inherited from the old MGs and 2.21% were new. The small amount of new positive allele emergence and possible recombination between alleles might explain the slight SPC increase in the new MGs. The prediction of recombination potentials in the SPC of all the possible crosses indicated that the mean of SPC overall crosses was 43.29% (+2.52%) and the maximum was 50.00% (+9.23%) in the SPC, and the maximum transgressive potential was 3.93%, suggesting that SPC breeding potentials do exist in the NECSGP. A total of 120 candidate genes were annotated and functionally classified into 13 categories, indicating that SPC is a complex trait conferred by a gene network.

Ensuring Global Food Security by Improving Protein Content in Major Grain Legumes Using Breeding and 'Omics' Tools

Grain legumes are a rich source of dietary protein for millions of people globally and thus a key driver for securing global food security. Legume plant-based 'dietary protein' biofortification is an economic strategy for alleviating the menace of rising malnutrition-related problems and hidden hunger. Malnutrition from protein deficiency is predominant in human populations with an insufficient daily intake of animal protein/dietary protein due to economic limitations, especially in developing countries. Therefore, enhancing grain legume protein content will help eradicate protein-related malnutrition problems in low-income and underprivileged countries. Here, we review the exploitable genetic variability for grain protein content in various major grain legumes for improving the protein content of high-yielding, low-protein genotypes. We highlight classical genetics-based inheritance of protein content in various legumes and discuss advances in molecular marker technology that have enabled us to underpin various quantitative trait loci controlling seed protein content (SPC) in biparental-based mapping populations and genome-wide association studies. We also review the progress of functional genomics in deciphering the underlying candidate gene(s) controlling SPC in various grain legumes and the role of proteomics and metabolomics in shedding light on the accumulation of various novel proteins and metabolites in high-protein legume genotypes. Lastly, we detail the scope of genomic selection, high-throughput phenotyping, emerging genome editing tools, and speed breeding protocols for enhancing SPC in grain legumes to achieve legume-based dietary protein security and thus reduce the global hunger risk.

The Population Structure of a Globe Artichoke Worldwide Collection, as Revealed by Molecular and Phenotypic Analyzes

The knowledge of the organization of the domesticated gene pool of crop species is an essential requirement to understand crop evolution, to rationalize conservation programs, and to support practical decisions in plant breeding. Here, we integrate simple sequence repeat (SSR) analysis and phenotypic characterization to investigate a globe artichoke collection that comprises most of the varieties cultivated worldwide. We show that the cultivated gene pool of globe artichoke includes five distinct genetic groups associated with the major phenotypic typologies: Catanesi (which based on our analysis corresponds to Violetti di Provenza), Spinosi, Violetti di Toscana, Romaneschi, and Macau. We observed that 17 and 11% of the molecular and phenotypic variance, respectively, is between these groups, while within groups, strong linkage disequilibrium and heterozygote excess are evident. The divergence between groups for quantitative traits correlates with the average broad-sense heritability within the groups. The phenotypic divergence between groups for both qualitative and quantitative traits is strongly and positively correlated with SSR divergence (FST) between groups. All this implies a low population size and strong bottleneck effects, and indicates a long history of clonal propagation and selection during the evolution of the domesticated gene pool of globe artichoke. Moreover, the comparison between molecular and phenotypic population structures suggests that harvest time, plant architecture (i.e., plant height, stem length), leaf spininess, head morphology (i.e., head shape, bract shape, spininess) together with the number of heads per plant were the main targets of selection during the evolution of the cultivated germplasm. We emphasize our findings in light of the potential exploitation of this collection for association mapping studies.

Identification of quantitative trait loci associated with seed quality traits between Canadian and Ukrainian mega-environments using genome-wide association study

KEY MESSAGE

Identifying QTL associated with soybean seed quality traits from a diverse GWAS panel cultivated in Canadian and Ukrainian mega-environments may facilitate future cultivar development for foreign markets. Understanding the complex genetic basis of seed quality traits for soybean in the mega-environments (MEs) is critical for developing a marker-assisted selection program that will lead to breeding superior cultivars adapted to specific regions. This study aimed to analyze the accumulation of 14 soybean seed quality traits in Canadian ME and two seed quality traits in Ukrainian ME and identify associated ME specific quantitative trait loci (QTLSP) and ME universal QTL (QTLU) for protein and oil using a genome-wide association study (GWAS) panel consisting of 184 soybean genotypes. The panel was planted in three locations in Canada and two locations in Ukraine in 2018 and 2019. Genotype plus genotype-by-environment biplot analysis was conducted to assess the accumulation of individual seed compounds across different locations. The protein accumulation was high in the Canadian ME and low in the Ukrainian ME, whereas the oil concentration showed the opposite trends between the two MEs. No QTLU were identified across the MEs for protein and oil concentrations. In contrast, nine Canadian QTLSP for protein were identified on various chromosomes, which were co-located with QTL controlling other traits identified in the Canadian ME. The lack of common QTLU for protein and oil suggests that it may be necessary to use QTLSP associated with these traits separately for the Canadian and Ukrainian ME. Additional Ukrainian data for seed compounds other than oil and protein are required to identify novel QTLSP and QTLU for such traits for the individual or combined Canadian and Ukrainian MEs.

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.

A novel Synthetic phenotype association study approach reveals the landscape of association for genomic variants and phenotypes

INTRODUCTION

Genome-Wide Association Studies (GWAS) identify tagging variants in the genome that are statistically associated with the phenotype because of their linkage disequilibrium (LD) relationship with the causative mutation (CM). When both low-density genotyped accession panels with phenotypes and resequenced data accession panels are available, tagging variants can assist with post-GWAS challenges in CM discovery.

OBJECTIVES

Our objective was to identify additional GWAS evaluation criteria to assess correspondence between genomic variants and phenotypes, as well as enable deeper analysis of the localized landscape of association.

METHODS

We used genomic variant positions as Synthetic phenotypes in GWAS that we named "Synthetic phenotype association study" (SPAS). The extreme case of SPAS is what we call an "Inverse GWAS" where we used CM positions of cloned soybean genes. We developed and validated the Accuracy concept as a measure of the correspondence between variant positions and phenotypes.

RESULTS

The SPAS approach demonstrated that the genotype status of an associated variant used as a Synthetic phenotype enabled us to explore the relationships between tagging variants and CMs, and further, that utilizing CMs as Synthetic phenotypes in Inverse GWAS illuminated the landscape of association. We implemented the Accuracy calculation for a curated accession panel to an online Accuracy calculation tool (AccuTool) as a resource for gene identification in soybean. We demonstrated our concepts on three examples of soybean cloned genes. As a result of our findings, we devised an enhanced "GWAS to Genes" analysis (Synthetic phenotype to CM strategy, SP2CM). Using SP2CM, we identified a CM for a novel gene.

CONCLUSION

The SP2CM strategy utilizing Synthetic phenotypes and the Accuracy calculation of correspondence provides crucial information to assist researchers in CM discovery. The impact of this work is a more effective evaluation of landscapes of GWAS associations.

Haplotype mapping uncovers unexplored variation in wild and domesticated soybean at the major protein locus cqProt-003

KEY MESSAGE

The major soy protein QTL, cqProt-003, was analysed for haplotype diversity and global distribution, and results indicate 304 bp deletion and variable tandem repeats in protein coding regions are likely causal candidates. Here, we present association and linkage analysis of 985 wild, landrace and cultivar soybean accessions in a pan genomic dataset to characterize the major high-protein/low-oil associated locus cqProt-003 located on chromosome 20. A significant trait-associated region within a 173 kb linkage block was identified, and variants in the region were characterized, identifying 34 high confidence SNPs, 4 insertions, 1 deletion and a larger 304 bp structural variant in the high-protein haplotype. Trinucleotide tandem repeats of variable length present in the second exon of gene Glyma.20G085100 are strongly correlated with the high-protein phenotype and likely represent causal variation. Structural variation has previously been found in the same gene, for which we report the global distribution of the 304 bp deletion and have identified additional nested variation present in high-protein individuals. Mapping variation at the cqProt-003 locus across demographic groups suggests that the high-protein haplotype is common in wild accessions (94.7%), rare in landraces (10.6%) and near absent in cultivated breeding pools (4.1%), suggesting its decrease in frequency primarily correlates with domestication and continued during subsequent improvement. However, the variation that has persisted in under-utilized wild and landrace populations holds high breeding potential for breeders willing to forego seed oil to maximize protein content. The results of this study include the identification of distinct haplotype structures within the high-protein population, and a broad characterization of the genomic context and linkage patterns of cqProt-003 across global populations, supporting future functional characterization and modification.

Fine mapping and cloning of the major seed protein quantitative trait loci on soybean chromosome 20

Soybean is the most important source of protein meal worldwide and the quantitative trait loci (QTL) cqSeed protein‐003 on chromosome 20 exerts the greatest additive effect of any protein QTL mapped in the crop. Through genetic mapping and candidate gene downregulation, we identified that an insertion/deletion variant in Glyma.20G85100 is the likely gene that underlies this important QTL.

Bioinformatics for plant and agricultural discoveries in the age of multiomics: A review and case study of maize nodal root growth under water deficit

Advances in next-generation sequencing and other high-throughput technologies have facilitated multiomics research, such as genomics, epigenomics, transcriptomics, proteomics, metabolomics, and phenomics. The resultant emerging multiomics data have brought new challenges as well as opportunities, as seen in the plant and agriculture science domains. We reviewed several bioinformatic and computational methods, models, and platforms, and we have highlighted some of our in-house developed efforts aimed at multiomics data analysis, integration, and management issues faced by the research community. A case study using multiomics datasets generated from our studies of maize nodal root growth under water deficit stress demonstrates the power of these datasets and some other publicly available tools. This analysis also sheds light on the landscape of such applied bioinformatic tools currently available for plant and crop science studies and introduces emerging trends and how they may affect the future.