Genetic separation of southern and northern soybean breeding programs in North America and their associated allelic variation at four maturity loci

Genome-wide association study reveals GmFulb as candidate gene for maturity time and reproductive length in soybeans (Glycine max)

The ability of soybean [Glycine max (L.) Merr.] to adapt to different latitudes is attributed to genetic variation in major E genes and quantitative trait loci (QTLs) determining flowering time (R1), maturity (R8), and reproductive length (RL). Fully revealing the genetic basis of R1, R8, and RL in soybeans is necessary to enhance genetic gains in soybean yield improvement. Here, we performed a genome-wide association analysis (GWA) with 31,689 single nucleotide polymorphisms (SNPs) to detect novel loci for R1, R8, and RL using a soybean panel of 329 accessions with the same genotype for three major E genes (e1-as/E2/E3). The studied accessions were grown in nine environments and observed for R1, R8 and RL in all environments. This study identified two stable peaks on Chr 4, simultaneously controlling R8 and RL. In addition, we identified a third peak on Chr 10 controlling R1. Association peaks overlap with previously reported QTLs for R1, R8, and RL. Considering the alternative alleles, significant SNPs caused RL to be two days shorter, R1 two days later and R8 two days earlier, respectively. We identified association peaks acting independently over R1 and R8, suggesting that trait-specific minor effect loci are also involved in controlling R1 and R8. From the 111 genes highly associated with the three peaks detected in this study, we selected six candidate genes as the most likely cause of R1, R8, and RL variation. High correspondence was observed between a modifying variant SNP at position 04:39294836 in GmFulb and an association peak on Chr 4. Further studies using map-based cloning and fine mapping are necessary to elucidate the role of the candidates we identified for soybean maturity and adaptation to different latitudes and to be effectively used in the marker-assisted breeding of cultivars with optimal yield-related traits.

Molecular breeding for improvement of photothermal adaptability in soybean

Soybean (Glycine max (L.) Merr.) is a typical short-day and temperate crop that is sensitive to photoperiod and temperature. Responses of soybean to photothermal conditions determine plant growth and development, which affect its architecture, yield formation, and capacity for geographic adaptation. Flowering time, maturity, and other traits associated with photothermal adaptability are controlled by multiple major-effect and minor-effect genes and genotype-by-environment interactions. Genetic studies have identified at least 11 loci (E1-E4, E6-E11, and J) that participate in photoperiodic regulation of flowering time and maturity in soybean. Molecular cloning and characterization of major-effect flowering genes have clarified the photoperiod-dependent flowering pathway, in which the photoreceptor gene phytochrome A, circadian evening complex (EC) components, central flowering repressor E1, and FLOWERING LOCUS T family genes play key roles in regulation of flowering time, maturity, and adaptability to photothermal conditions. Here, we provide an overview of recent progress in genetic and molecular analysis of traits associated with photothermal adaptability, summarizing advances in molecular breeding practices and tools for improving these traits. Furthermore, we discuss methods for breeding soybean varieties with better adaptability to specific ecological regions, with emphasis on a novel strategy, the Potalaization model, which allows breeding of widely adapted soybean varieties through the use of multiple molecular tools in existing elite widely adapted varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01406-z.

Genetic relationships and genome selection signatures between soybean cultivars from Brazil and United States after decades of breeding

Soybean is one of the most important crops worldwide. Brazil and the United States (US) are the world's two biggest producers of this legume. The increase of publicly available DNA sequencing data as well as high-density genotyping data of multiple soybean germplasms has made it possible to understand the genetic relationships and identify genomics regions that underwent selection pressure during soy domestication and breeding. In this study, we analyzed the genetic relationships between Brazilian (N = 235) and US soybean cultivars (N = 675) released in different decades and screened for genomic signatures between Brazilian and US cultivars. The population structure analysis demonstrated that the Brazilian germplasm has a narrower genetic base than the US germplasm. The US cultivars were grouped according to maturity groups, while Brazilian cultivars were separated according to decade of release. We found 73 SNPs that differentiate Brazilian and US soybean germplasm. Maturity-associated SNPs showed high allelic frequency differences between Brazilian and US accessions. Other important loci were identified separating cultivars released before and after 1996 in Brazil. Our data showed important genomic regions under selection during decades of soybean breeding in Brazil and the US that should be targeted to adapt lines from different origins in these countries.

Heterosis and Differential DNA Methylation in Soybean Hybrids and Their Parental Lines

Heterosis is an important biological phenomenon and is widely applied to increase agricultural productivity. However, the underlying molecular mechanisms of heterosis are still unclear. Here we constructed three combinations of reciprocal hybrids of soybean, and subsequently used MethylRAD-seq to detect CCGG and CCWGG (W = A or T) methylation in the whole genome of these hybrids and their parents at the middle development period of contemporary seed. We were able to prove that changes in DNA methylation patterns occurred in immature hybrid seeds and the parental variation was to some degree responsible for differential expression between the reciprocal hybrids. Non-additive differential methylation sites (DMSs) were also identified in large numbers in hybrids. Interestingly, most of these DMSs were hyper-methylated and were more concentrated in gene regions than the natural distribution of the methylated sites. Further analysis of the non-additive DMSs located in gene regions exhibited their participation in various biological processes, especially those related to transcriptional regulation and hormonal function. These results revealed DNA methylation reprogramming pattern in the hybrid soybean, which is associated with phenotypic variation and heterosis initiation.

Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes

BACKGROUND

With advances in next-generation sequencing technologies, an unprecedented amount of soybean accessions has been sequenced by many individual studies and made available as raw sequencing reads for post-genomic research.

RESULTS

To develop a consolidated and user-friendly genomic resource for post-genomic research, we consolidated the raw resequencing data of 1465 soybean genomes available in the public and 91 highly diverse wild soybean genomes newly sequenced. These altogether provided a collection of 1556 sequenced genomes of 1501 diverse accessions (1.5 K). The collection comprises of wild, landraces and elite cultivars of soybean that were grown in East Asia or major soybean cultivating areas around the world. Our extensive sequence analysis discovered 32 million single nucleotide polymorphisms (32mSNPs) and revealed a SNP density of 30 SNPs/kb and 12 non-synonymous SNPs/gene reflecting a high structural and functional genomic diversity of the new collection. Each SNP was annotated with 30 categories of structural and/or functional information. We further identified paired accessions between the 1.5 K and 20,087 (20 K) accessions in US collection as genomic "equivalent" accessions sharing the highest genomic identity for minimizing the barriers in soybean germplasm exchange between countries. We also exemplified the utility of 32mSNPs in enhancing post-genomics research through in-silico genotyping, high-resolution GWAS, discovering and/or characterizing genes and alleles/mutations, identifying germplasms containing beneficial alleles that are potentially experiencing artificial selection.

CONCLUSION

The comprehensive analysis of publicly available large-scale genome sequencing data of diverse cultivated accessions and the newly in-house sequenced wild accessions greatly increased the soybean genome-wide variation resolution. This could facilitate a variety of genetic and molecular-level analyses in soybean. The 32mSNPs and 1.5 K accessions with their comprehensive annotation have been made available at the SoyBase and Ag Data Commons. The dataset could further serve as a versatile and expandable core resource for exploring the exponentially increasing genome sequencing data for a variety of post-genomic research.

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.

The Synchronized Efforts to Decipher the Molecular Basis for Soybean Maturity Loci E1, E2, and E3 That Regulate Flowering and Maturity

The general concept of photoperiodism, i.e., the photoperiodic induction of flowering, was established by Garner and Allard (1920). The genetic factor controlling flowering time, maturity, or photoperiodic responses was observed in soybean soon after the discovery of the photoperiodism. E1, E2, and E3 were named in 1971 and, thereafter, genetically characterized. At the centennial celebration of the discovery of photoperiodism in soybean, we recount our endeavors to successfully decipher the molecular bases for the major maturity loci E1, E2, and E3 in soybean. Through systematic efforts, we successfully cloned the E3 gene in 2009, the E2 gene in 2011, and the E1 gene in 2012. Recently, successful identification of several circadian-related genes such as PRR3a, LUX, and J has enriched the known major E1-FTs pathway. Further research progresses on the identification of new flowering and maturity-related genes as well as coordinated regulation between flowering genes will enable us to understand profoundly flowering gene network and determinants of latitudinal adaptation in soybean.

Selection of GmSWEET39 for oil and protein improvement in soybean

Soybean [Glycine max (L.) Merr.] was domesticated from wild soybean (G. soja Sieb. and Zucc.) and has been further improved as a dual-use seed crop to provide highly valuable oil and protein for food, feed, and industrial applications. However, the underlying genetic and molecular basis remains less understood. Having combined high-confidence bi-parental linkage mapping with high-resolution association analysis based on 631 whole sequenced genomes, we mapped major soybean protein and oil QTLs on chromosome15 to a sugar transporter gene (GmSWEET39). A two-nucleotide CC deletion truncating C-terminus of GmSWEET39 was strongly associated with high seed oil and low seed protein, suggesting its pleiotropic effect on protein and oil content. GmSWEET39 was predominantly expressed in parenchyma and integument of the seed coat, and likely regulates oil and protein accumulation by affecting sugar delivery from maternal seed coat to the filial embryo. We demonstrated that GmSWEET39 has a dual function for both oil and protein improvement and undergoes two different paths of artificial selection. A CC deletion (CC-) haplotype H1 has been intensively selected during domestication and extensively used in soybean improvement worldwide. H1 is fixed in North American soybean cultivars. The protein-favored (CC+) haplotype H3 still undergoes ongoing selection, reflecting its sustainable role for soybean protein improvement. The comprehensive knowledge on the molecular basis underlying the major QTL and GmSWEET39 haplotypes associated with soybean improvement would be valuable to design new strategies for soybean seed quality improvement using molecular breeding and biotechnological approaches.

Impacts of genomic research on soybean improvement in East Asia

It has been commonly accepted that soybean domestication originated in East Asia. Although East Asia has the historical merit in soybean production, the USA has become the top soybean producer in the world since 1950s. Following that, Brazil and Argentina have been the major soybean producers since 1970s and 1990s, respectively. China has once been the exporter of soybean to Japan before 1990s, yet she became a net soybean importer as Japan and the Republic of Korea do. Furthermore, the soybean yield per unit area in East Asia has stagnated during the past decade. To improve soybean production and enhance food security in these East Asian countries, much investment has been made, especially in the breeding of better performing soybean germplasms. As a result, China, Japan, and the Republic of Korea have become three important centers for soybean genomic research. With new technologies, the rate and precision of the identification of important genomic loci associated with desired traits from germplasm collections or mutants have increased significantly. Genome editing on soybean is also becoming more established. The year 2019 marked a new era for crop genome editing in the commercialization of the first genome-edited plant product, which is a high-oleic-acid soybean oil. In this review, we have summarized the latest developments in soybean breeding technologies and the remarkable progress in soybean breeding-related research in China, Japan, and the Republic of Korea.

Nested association mapping of important agronomic traits in three interspecific soybean populations

KEY MESSAGE

Glycine soja germplasm can be used to successfully introduce new alleles with the potential to add valuable new genetic diversity to the current elite soybean gene pool. Given the demonstrated narrow genetic base of the US soybean production, it is essential to identify beneficial alleles from exotic germplasm, such as wild soybean, to enhance genetic gain for favorable traits. Nested association mapping (NAM) is an approach to population development that permits the comparison of allelic effects of the same QTL in multiple parents. Seed yield, plant maturity, plant height and plant lodging were evaluated in a NAM panel consisting of 392 recombinant inbred lines derived from three biparental interspecific soybean populations in eight environments during 2016 and 2017. Nested association mapping, combined with linkage mapping, identified three major QTL for plant maturity in chromosomes 6, 11 and 12 associated with alleles from wild soybean resulting in significant increases in days to maturity. A significant QTL for plant height was identified on chromosome 13 with the allele increasing plant height derived from wild soybean. A significant grain yield QTL was detected on chromosome 17, and the allele from Glycine soja had a positive effect of 166 kg ha-1; RIL's with the wild soybean allele yielded on average 6% more than the lines carrying the Glycine max allele. These findings demonstrate the usefulness and potential of alleles from wild soybean germplasm to enhance important agronomic traits in a soybean breeding program.

Characterization of Two Growth Period QTLs Reveals Modification of PRR3 Genes During Soybean Domestication

Soybean yield is largely dependent on growth period. We characterized two growth period quantitative trait loci, Gp11 and Gp12, from a recombinant inbred population generated from a cross of wild (W05) and cultivated (C08) soybean. Lines carrying Gp11C08 and Gp12C08 tend to have a shorter growth period and higher expression of GmFT2a and GmFT5a. Furthermore, multiple interval mapping suggests that Gp11 and Gp12 may be genetically interacting with the E2 locus. This is consistent with the observation that GmFT2a and GmFT5a are activated by Gp11C08 and Gp12C08 at ZT4 in the recessive e2 but not the dominant E2 background. Gp11 and Gp12 are duplicated genomic regions each containing a copy of the soybean ortholog of PSEUDO RESPONSE REGULATOR 3 (GmPRR3A and GmPRR3B). GmPRR3A and GmPRR3B from C08 carry mutations that delete the CCT domain in the encoded proteins. These mutations were selected during soybean improvement and they alter the subcellular localization of GmPRR3A and GmPRR3B. Furthermore, GmPRR3A and GmPRR3B can interact with TOPLESS-related transcription factors, suggesting that they function in a transcription repressor complex. This study addresses previously unexplored components of the genetic network that probably controls the growth period of soybean and puts these loci into context with the well-characterized growth period-regulating E loci.

Accelerating Soybean Breeding in a CO2-Supplemented Growth Chamber

Soybean (Glycine max) is the most important dicot crop worldwide, and is increasingly used as a model legume due to the wide availability of genomic soybean resources; however, the slow generation times of soybean plants are currently a major hindrance to research. Here, we demonstrate a method for accelerating soybean breeding in compact growth chambers, which greatly shortens the generation time of the plants and accelerates breeding and research projects. Our breeding method utilizes commonly used fluorescent lamps (220 µmol m-2 s-1 at the canopy level), a 14 h light (30°C)/10 h dark (25°C) cycle and carbon dioxide (CO2) supplementation at >400 p.p.m. Using this approach, the generation time of the best-characterized elite Japanese soybean cultivar, Enrei, was shortened from 102-132 d reported in the field to just 70 d, thereby allowing up to 5 generations per year instead of the 1-2 generations currently possible in the field and/or greenhouse. The method also facilitates the highly efficient and controlled crossing of soybean plants. Our method uses CO2 supplementation to promote the growth and yield of plants, appropriate light and temperature conditions to reduce the days to flowering, and the reaping and sowing of immature seeds to shorten the reproductive period greatly. Thus, the appropriate parameters enable acceleration of soybean breeding in the compact growth chambers commonly used for laboratory research. The parameters used in our method could therefore be optimized for other species, cultivars, accessions and experimental designs to facilitate rapid breeding in a wide range of crops.

Allelic Variation and Distribution of the Major Maturity Genes in Different Soybean Collections

Soybean time of flowering and maturity are genetically controlled by E genes. Different allelic combinations of these genes determine soybean adaptation to a specific latitude. The paper describes the first attempt to assess adaptation of soybean genotypes developed and realized at Institute of Field and Vegetable Crops, Novi Sad, Serbia [Novi Sad (NS) varieties and breeding lines] based on E gene variation, as well as to comparatively assess E gene variation in North-American (NA), Chinese, and European genotypes, as most of the studies published so far deal with North-American and Chinese cultivars and breeding material. Allelic variation and distribution of the major maturity genes (E1, E2, E3, and E4) has been determined in 445 genotypes from soybean collections of NA ancestral lines, Chinese germplasm, and European varieties, as well as NS varieties and breeding lines. The study showed that allelic combinations of E1-E4 genes significantly determined the adaptation of varieties to different geographical regions, although they have different impacts on maturity. In general, each collection had one major E genotype haplogroup, comprising over 50% of the lines. The exceptions were European varieties that had two predominant haplogroups and NA ancestral lines distributed almost evenly among several haplogroups. As e1-as/e2/E3/E4 was the most common genotype in NS population, present in the best-performing genotypes in terms of yield, this specific allele combination was proposed as the optimal combination for the environments of Central-Eastern Europe.