Wild Soybean Accessions Research Articles

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

Key messageIn 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.

GmHs1-1 and GmqHS1 Simultaneously Contribute to the Domestication of Soybean Hard-Seededness.

Seed physical dormancy (hard-seededness) is an interesting ecological phenomenon and important agronomic trait. The loss of seed coat impermeability/hard-seededness is a key target trait during the domestication of leguminous crops which allows seeds to germinate rapidly and uniformly. In this study, we examined the mutation of quantitative trait locus (QTL) genes, GmHs1-1 and GmqHS1, in 18 wild soybean (G. soja) and 23 cultivated soybean (G. max) accessions. The sequencing results indicate that a G-to-T substitution in GmqHS1 and a C-to-T substitution in GmHs1-1 occurred in all 23 cultivated soybean accessions but not in any of the 18 wild soybean accessions. The mutations in the two genes led to increased seed coat permeability in cultivated soybean. Therefore, we provide evidence that two genes, GmHs1-1 and GmqHS1, simultaneously contribute to the domestication of hard-seededness in soybeans. This finding is of great significance for genetic analysis and improved utilization of the soybean hard-seededness trait.

Genetic Control of Tolerance to Drought Stress in Wild Soybean (Glycine soja) at the Vegetative and the Germination Stages.

Drought stress, which is becoming more prevalent due to climate change, is a significant abiotic factor that adversely impacts crop production and yield stability. Cultivated soybean (Glycine max), a versatile crop for humans and animals, exhibits sensitivity to drought, resulting in reduced growth and development under drought conditions. However, few genetic studies have assessed wild soybean's (Glycine soja) response to drought stress. In this work, we conducted a genome-wide association study (GWAS) and analysis of wild soybean accessions to identify loci responsible for drought tolerance at the vegetative (n = 187) and the germination stages (n = 135) using the available resequencing data. The GWAS analysis of the leaf wilting score (LWS) identified eight single-nucleotide polymorphisms (SNPs) on chromosomes 10, 11, and 19. Of these, wild soybeans with both SNPs on chromosomes 10 (adenine) and 11 (thymine) produced lower LWS, indicating that these SNPs have an important role in the genetic effect on LWS for drought tolerance at the vegetative stage. At the germination stage, nine SNPs associated with five phenotypic measurements were identified on chromosomes 6, 9, 10, 13, 16, and 17, and the genomic regions identified at the germination stage were different from those identified for the LWS, supporting our previous finding that there may not be a robust correlation between the genes influencing phenotypes at the germination and vegetative stages. This research will benefit marker-assisted breeding programs aimed at enhancing drought tolerance in soybeans.

Research advances on the hard seededness trait of soybean and the underlying regulatory mechanisms.

Soybean is one of the world's most economically significant crops and is widely utilized as an essential source of vegetable protein and edible oil. Cultivated soybean is domesticated from its annual counterpart, wild soybean, which is considered valuable germplasm for soybean breeding. However, wild soybean accessions generally produce seeds with impermeable coats, a trait known as hard seededness (HS), which is beneficial for long-term seed survival but is undesirable for the uniform water absorption and germination of seeds, thus limiting the utilization of wild soybeans in breeding. In addition, moderate HS can isolate the embryo from the surrounding environment and is thus beneficial for long-term seed storage and germplasm preservation. The HS trait is primarily associated with the structure and chemical composition of the seed coat. Moreover, its development is also influenced by various environmental conditions, such as water and temperature. Genetic analysis has revealed that HS of soybean is a complex quantitative trait controlled by multiple genes or minor quantitative trait loci (QTL), with many QTLs and several causal genes currently identified. Investigating the physiological and molecular mechanisms underlying this trait is crucial for soybean breeding, production, and food processing. For this article, the literature was reviewed and condensed to create a well-rounded picture of the current understanding of internal and external factors, QTLs, causal genes, and the regulatory mechanisms related to the HS of soybean, with the aim of providing reference for future research and utilization of this trait.

Identification of the domestication gene GmCYP82C4 underlying the major quantitative trait locus for the seed weight in soybean.

A major quantitative trait locus (QTL) for the hundred-seed weight (HSW) was identified and confirmed in the two distinct soybean populations, and the target gene GmCYP82C4 underlying this locus was identified that significantly associated with soybean seed weight, and it was selected during the soybean domestication and improvement process. Soybean is a major oil crop for human beings and the seed weight is a crucial goal of soybean breeding. However, only a limited number of target genes underlying the quantitative trait loci (QTLs) controlling seed weight in soybean are known so far. In the present study, six loci associated with hundred-seed weight (HSW) were detected in the first population of 573 soybean breeding lines by genome-wide association study (GWAS), and 64 gene models were predicted in these candidate QTL regions. The QTL qHSW_1 exhibits continuous association signals on chromosome four and was also validated by region association study (RAS) in the second soybean population (409 accessions) with wild, landrace, and cultivar soybean accessions. There were seven genes in qHSW_1 candidate region by linkage disequilibrium (LD) block analysis, and only Glyma.04G035500 (GmCYP82C4) showed specifically higher expression in flowers, pods, and seeds, indicating its crucial role in the soybean seed development. Significant differences in HSW trait were detected when the association panels are genotyped by single-nucleotide polymorphisms (SNPs) in putative GmCYP82C4 promoter region. Eight haplotypes were generated by six SNPs in GmCYP82C4 in the second soybean population, and two superior haplotypes (Hap2 and Hap4) of GmCYP82C4 were detected with average HSW of 18.27g and 18.38g, respectively. The genetic diversity of GmCYP82C4 was analyzed in the second soybean population, and GmCYP82C4 was most likely selected during the soybean domestication and improvement process, leading to the highest proportion of Hap2 of GmCYP82C4 both in landrace and cultivar subpopulations. The QTLs and GmCYP82C4 identified in this study provide novel genetic resources for soybean seed weight trait, and the GmCYP82C4 could be used for soybean molecular breeding to develop desirable seed weight in the future.

Assessment of Drought Responses of Wild Soybean Accessions at Different Growth Stages

Drought is a significant abiotic stress that limits crop production. Soybeans [Glycine max (L.) Merr.] are regarded as drought-sensitive. In the present study, the drought responses of wild soybean accessions were evaluated at different growth stages. Based on the leaf wilting index of 411 accessions at the vegetative stage, seven highly tolerant (HT) and 24 tolerant (T) accessions were identified, although most wild soybeans were classified as moderate (M), sensitive (S), and highly sensitive (HS) genotypes. In addition, with selected wild soybeans, stomatal density decreased in HT accessions but increased in HS accessions under drought conditions at the vegetative stage. However, for cultivated soybeans, the stomatal density of the drought-tolerant and drought-sensitive were not significantly different between the two conditions. The expression levels of drought-related transcriptional factors indicated that the HT genotype showed a higher expression level of drought-related genes than that of the HS genotype at the vegetative stage. At the reproductive stages, 12 wild soybeans randomly selected from HT, T, S, and HS based on assessment at the vegetative stage showed consistent drought responses with seed yield, root development, and water status. However, the extent of the detrimental effect of drought on the germination rates and root length of 165 wild soybeans at the germination stage varied depending on the genotype, indicating that there may not be a robust correlation between phenotypic measurements at the germination stage and drought-related assessments at the two growth stages. The information from this study can provide useful breeding materials for the development of drought-tolerant cultivars from wild soybeans.

A novel miR160a-GmARF16-GmMYC2 module determines soybean salt tolerance and adaptation.

Salt stress is a major challenge that has a negative impact on soybean growth and productivity. Therefore, it is important to understand the regulatory mechanism of salt response to ensure soybean yield under such conditions. In this study, we identified and characterized a miR160a-GmARF16-GmMYC2 module and its regulation during the salt-stress response in soybean. miR160a promotes salt tolerance by cleaving GmARF16 transcripts, members of the Auxin Response Factor (ARF) family, which negatively regulates salt tolerance. In turn, GmARF16 activates GmMYC2, encoding a bHLH transcription factor that reduces salinity tolerance by down-regulating proline biosynthesis. Genomic analysis among wild and cultivated soybean accessions identified four distinct GmARF16 haplotypes. Among them, the GmARF16H3 haplotype is preferentially enriched in localities with relatively saline soils, suggesting GmARF16H3 was artificially selected to improve salt tolerance. Our findings therefore provide insights into the molecular mechanisms underlying salt response in soybean and provide valuable genetic targets for the molecular breeding of salt tolerance.

A QTL allele from wild soybean enhances protein content without reducing the oil content

Abstract Soybean is one of the chief crops producing protein and oil for human consumption. Wild soybean, the ancestor of cultivated soybean, possesses high seed protein content; therefore, it is a valuable genetic resource that could enhance protein content in the cultivated varieties. To identify the genes responsible for increasing protein content in wild soybean, a population comprising 113 BC4F6 chromosome segment substitution lines (CSSL) was developed from a cross between soybean cultivar ‘Jackson’ and wild soybean accession JWS156-1. The CSSL population was cultivated in the field conditions for 3 years (2018, 2019 and 2020), and the seeds harvested from each line were analysed for protein and oil contents by InfraTec Nova instrument. Quantitative trait locus (QTL) analysis with 243 SSR markers identified 12 QTLs associated with seed protein, oil and protein + oil contents. Among these QTLs, qPro8 and qPro19, two major and stable QTLs for protein content, were detected on chromosomes 8 and 19, respectively. No QTL for oil content was detected in the vicinity of qPro19, indicating that qPro19 did not influence the seed oil content. The effect of qPro19 was validated using near-isogenic lines (NILs) of qPro19. By introducing the qPro19 allele from wild soybean into another soybean variety, ‘Tachiyutaka’, a BC4 line, T-678, that showed enhanced seed protein content, without reducing the seed oil content. This study implied that the qPro19 allele from wild soybean could be a potential genetic resource for breeding programmes aimed to improve soybean seed quality.

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.

QTL Mapping of Soybean (Glycine max) Vine Growth Habit Trait.

The vine growth habit (VGH) is a notable property of wild soybean plants that also holds a high degree of importance in domestication as it can preclude using these wild cultivars for breeding and improving domesticated soybeans. Here, a bulked segregant analysis (BSA) approach was employed to study the genetic etiology of the VGH in soybean plants by integrating linkage mapping and population sequencing approaches. To develop a recombinant inbred line (RIL) population, the cultivated Zhongdou41 (ZD41) soybean cultivar was bred with ZYD02787, a wild soybean accession. The VGH status of each line in the resultant population was assessed, ultimately leading to the identification of six and nine QTLs from the BSA sequencing of the F4 population and F6-F8 population sequence mapping, respectively. One QTL shared across these analyzed generations was detected on chromosome 19. Three other QTLs detected by BSA-seq were validated and localized to the 90.93 kb, 2.9 Mb, and 602.08 kb regions of chromosomes 6 and 13, harboring 14, 53, and 4 genes, respectively. Three consistent VGH-related QTLs located on chromosomes 2 and 19 were detected in a minimum of three environments, while an additional six loci on chromosomes 2, 10, 13, and 18 were detected in at least two environments via ICIM mapping. Of all the detected loci, five had been reported previously whereas seven represent novel QTLs. Together, these data offer new insights into the genetic basis of the VGH in soybean plants, providing a rational basis to inform the use of wild accessions in future breeding efforts.

Single‐grain‐based widely targeted metabolomics profiling of sixty‐four accessions of Japanese wild soybean (Glycin soja Sieb. Et Zucc.)

SummarySoybean (Glycine max L. Merr.) is a major crop worldwide, providing essential dietary proteins and edible oils for agriculture, food, animal feed and biofuel industries. However, its limited genetic variation makes it vulnerable to harsh environmental conditions, and its chemical diversity was strongly reduced after long domestication. Wild soybean (Glycine soja Sieb. et Zucc.), the ancestor of the G. max, possesses a high level of genetic diversity and thus present a precious reservoir for improving cultivated soybean. Characterising metabolomes of different G. soja accessions is crucial in exploring metabolite accumulation patterns for developing metabolic breeding strategies. This study successfully revealed the metabolomes of sixty‐four seed lines of G. soja collected throughout Japan by a widely targeted metabolomics approach, therefore, providing a good source of information for further studies to induce the production of bioactive compounds and the defence mechanism to stress.

The AP2/ERF transcription factor TOE4b regulates photoperiodic flowering and grain yield per plant in soybean.

Photoperiod-mediated flowering determines the phenological adaptability of crops including soybean (Glycine max). A genome-wide association study (GWAS) identified a new flowering time locus, Time of flowering 13 (Tof13), which defined a gene encoding an AP2/ERF transcription factor. This new transcription factor, which we named TOE4b, is localized in the nucleus. TOE4b has been selected for soybean latitude adaptability. The existing natural variant TOE4bH4 was rare in wild soybean accessions but occurred more frequently in landraces and cultivars. Notably, TOE4bH4 improved high-latitude adaptation of soybean to some extent. The gene-edited TOE4b knockout mutant exhibited earlier flowering, conversely, TOE4b overexpression delayed flowering time. TOE4b is directly bound to the promoters and gene bodies of the key flowering integration factor genes FT2a and FT5a to inhibit their transcription. Importantly, TOE4b overexpression lines in field trials not only showed late flowering but also altered plant architecture, including shorter internode length, more internodes, more branches and pod number per plant, and finally boosted grain yield per plant by 60% in Guangzhou and 87% in Shijiazhuang. Our findings therefore identified TOE4b as a pleiotropic gene to increase yield potential per plant in soybean, and these results provide a promising option for breeding a soybean variety with an idealized plant architecture that promotes high yields.

Differential SW16.1 allelic effects and genetic backgrounds contributed to increased seed weight after soybean domestication.

Although seed weight has increased following domestication from wild soybean (Glycine soja) to cultivated soybean (Glycine max), the genetic basis underlying this change is unclear. Using mapping populations derived from chromosome segment substitution lines of wild soybean, we identified SW16.1 as the causative gene underlying a major quantitative trait locus controlling seed weight. SW16.1 encodes a nucleus-localized LIM domain-containing protein. Importantly, the GsSW16.1 allele from wild soybean accession N24852 had a negative effect on seed weight, whereas the GmSW16.1 allele from cultivar NN1138-2 had a positive effect. Gene expression network analysis, reverse-transcription quantitative polymerase chain reaction, and promoter-luciferase reporter transient expression assays suggested that SW16.1 regulates the transcription of MT4, a positive regulator of seed weight. The natural variations in SW16.1 and other known seed weight genes were analyzed in soybean germplasm. The SW16.1 polymorphism was associated with seed weight in 247 soybean accessions, showing much higher frequency of positive-effect alleles in cultivated soybean than in wild soybean. Interestingly, gene allele matrix analysis of the known seed weight genes revealed that G. max has lost 38.5% of the G. soja alleles and that most of the lost alleles had negative effects on seed weight. Our results suggest that eliminating negative alleles from G. soja led to a higher frequency of positive alleles and changed genetic backgrounds in G. max, which contributed to larger seeds in cultivated soybean after domestication from wild soybean. Our findings provide new insights regarding soybean domestication and should assist current soybean breeding programs.

A Genome-Wide Association Study of Protein, Oil, and Amino Acid Content in Wild Soybean (Glycine soja).

Soybean (Glycine max L.) is a globally important source of plant proteins, oils, and amino acids for both humans and livestock. Wild soybean (Glycine soja Sieb. and Zucc.), the ancestor of cultivated soybean, could be a useful genetic source for increasing these components in soybean crops. In this study, 96,432 single-nucleotide polymorphisms (SNPs) across 203 wild soybean accessions from the 180K Axiom® Soya SNP array were investigated using an association analysis. Protein and oil content exhibited a highly significant negative correlation, while the 17 amino acids exhibited a highly significant positive correlation with each other. A genome-wide association study (GWAS) was conducted on the protein, oil, and amino acid content using the 203 wild soybean accessions. A total of 44 significant SNPs were associated with protein, oil, and amino acid content. Glyma.11g015500 and Glyma.20g050300, which contained SNPs detected from the GWAS, were selected as novel candidate genes for the protein and oil content, respectively. In addition, Glyma.01g053200 and Glyma.03g239700 were selected as novel candidate genes for nine of the amino acids (Ala, Asp, Glu, Gly, Leu, Lys, Pro, Ser, and Thr). The identification of the SNP markers related to protein, oil, and amino acid content reported in the present study is expected to help improve the quality of selective breeding programs for soybeans.

Genetic Diversity of Korean Wild Soybean Core Collections and Genome-Wide Association Study for Days to Flowering.

The utilization of wild soybean germplasms in breeding programs increases genetic diversity, and they contain the rare alleles of traits of interest. Understanding the genetic diversity of wild germplasms is essential for determining effective strategies that can improve the economic traits of soybeans. Undesirable traits make it challenging to cultivate wild soybeans. This study aimed to construct a core subset of 1467 wild soybean accessions of the total population and analyze their genetic diversity to understand their genetic variations. Genome-wild association studies were conducted to detect the genetic loci underlying the time to flowering for a core subset collection, and they revealed the allelic variation in E genes for predicting maturity using the available resequencing data of wild soybean. Based on principal component and cluster analyses, 408 wild soybean accessions in the core collection covered the total population and were explained by 3 clusters representing the collection regions, namely, Korea, China, and Japan. Most of the wild soybean collections in this study had the E1e2E3 genotype according to association mapping and a resequencing analysis. Korean wild soybean core collections can provide helpful genetic resources to identify new flowering and maturity genes near the E gene loci and genetic materials for developing new cultivars, facilitating the introgression of genes of interest from wild soybean.

Genome-Wide Association Study for Agronomic Traits in Wild Soybean (Glycine soja)

The agronomic traits of soybean are important because they are directly or indirectly related to its yield. Cultivated soybean (Glycine max (L.) Merr) has lost genetic diversity during domestication and selective breeding. However, wild soybean (G. soja) represents a useful breeding material because it has a diverse gene pool. In this study, a total of 96,432 single-nucleotide polymorphisms (SNPs) across 203 wild soybean accessions from the 180K Axiom® Soya SNP array were employed in the association analysis. Wild soybean accessions were divided into four clusters based on their genetic distance using ADMIXTURE, principal component analysis, and neighbor-joining clusters. The linkage disequilibrium decayed rapidly in wild soybean. A genome-wide association study was conducted for days to flowering (DtF), days to maturity (DtM), the number of pods (NoP), and the 100-seed weight (100SW), which are major agronomic traits for wild soybean accessions. A total of 22 significant SNPs were found to be associated with DtF, DtM, and the 100SW. Based on the detected SNP markers, Glyma.12g210400, a gene related to DtF, Glyma.17g115300, a gene related to DtM, and Glyma.14g140200, a gene related to the 100SW, were selected as candidate genes. The SNP markers related to agronomic traits identified in this study are expected to help improve the quality of soybean cultivars through selective breeding.

The genetic basis of high-latitude adaptation in wild soybean

In many plants, flowering time is influenced by daylength as an adaptive response. In soybean (Glycine max) cultivars, however, photoperiodic flowering reduces crop yield and quality in high-latitude regions. Understanding the genetic basis of wild soybean (Glycine soja) adaptation to high latitudes could aid breeding of improved cultivars. Here, we identify the Tof4 (Time of flowering 4) locus, which encodes by an E1-like protein, E1La, that represses flowering and enhances adaptation to high latitudes in wild soybean. Moreover, we found that Tof4 physically associates with the promoters of two important FLOWERING LOCUS T (FT2a and FT5a) and with Tof5 to inhibit their transcription under long photoperiods. The effect of Tof4 on flowering and maturity is mediated by FT2a and FT5a proteins. Intriguingly, Tof4 and the key flowering repressor E1 independently but additively regulate flowering time, maturity, and grain yield in soybean. We determined that weak alleles of Tof4 have undergone natural selection, facilitating adaptation to high latitudes in wild soybean. Notably, over 71.5% of wild soybean accessions harbor the mutated alleles of Tof4 or a previously reported gain-of-function allele Tof5H2, suggesting that these two loci are the genetic basis of wild soybean adaptation to high latitudes. Almost no cultivated soybean carries the mutated tof4 allele. Introgression of the tof4-1 and Tof5H2 alleles into modern soybean or editing E1 family genes thus represents promising avenues to obtain early-maturity soybean, thereby improving productivity in high latitudes.

Genome-wide association study uncovers major genetic loci associated with flowering time in response to active accumulated temperature in wild soybean population

Flowering time and active accumulated temperature (AAT) are two key factors that limit the expanded production especially for soybean across different regions. Wild soybean provides an important germplasm for functional genomics study in cultivar soybean. However, the studies on genetic basis underlying flowering time in response to AAT especially in wild soybean were rarely reported. In this study, we used 294 wild soybean accessions derived from major soybean production region characterized by different AAT in Northeast of China. Based on genome-wide association study (GWAS), we identified 96 SNPs corresponded to 342 candidate genes that significantly associated with flowering time recorded in two-year experiments. Gene Ontology enrichment analysis suggests that the pathways of photosynthesis light reaction and actin filament binding were significantly enriched. We found three lead SNPs with -log10(p-value) > 32 across the two-year experiments, i.e., Chr02:9490318, Chr04:8545910 and Chr09:49553555. Linkage disequilibrium block analysis shows 28 candidate genes within the genomic region centered on the lead SNPs. Among them, expression levels of three genes (aspartic peptidase 1, serine/threonine-protein kinase and protein SCAR2-like) were significantly differed between two subgroups possessing contrasting flowering time distributed at chromosome 2, 4 and 9, respectively. There are 6, 7 and 3 haplotypes classified on the coding regions of the three genes, respectively. Collectively, accessions with late flowering time phenotype are typically derived from AAT zone 1, which is associated with the haplotypic distribution and expression levels of the three genes. This study provides an insight into a potential mechanism responsible for flowering time in response to AAT in wild soybean, which could promote the understanding of genetic basis for other major crops.

Whole-genome resequencing reveals signature of local adaptation and divergence in wild soybean.

Global climate change has threatened world crop production and food security. Decoding the adaptive genetic basis of wild relatives provides an invaluable genomic resource for climate-smart crop breedinG.Here, we performed whole-genome sequencing of 185 diverse wild soybean (Glycine soja) accessions collected from three major agro-ecological zones in China to parse the genomic basis of local adaptation in wild soybean. The population genomic diversity pattern exhibited clear agro-ecological zone-based population structure, and multiple environmental factors were observed to contribute to the genetic divergence. Demographic analysis shows that wild soybeans from the three ecological zones diverged about 1 × 105years ago, and then the effective population sizes have undergone different degrees of expansions. Genome-environment association identified multiple genes involved in the local adaptation, such as flowering time and temperature-related genes. A locus containing two adjacent MADS-box transcription factors on chromosome 19 was identified for multiple environmental factors, and it experienced positive selection that enables the adaptation to high-latitude environment. This study provides insights into the genetic mechanism of ecological adaptation in wild soybean that may facilitate climate-resilient soybean breeding.

Genome-Wide Association Studies Reveal Novel Loci for Herbivore Resistance in Wild Soybean (Glycine soja).

The production of soybean [Glycine max (L.) Merr.] is seriously threatened by various leaf-feeding insects, and wild soybean [Glycine soja Sieb. & Zucc.] has a greater resistance capacity and genetic diversity. In this study, a natural population consisting of 121 wild soybean accessions was used for detecting insect resistance genes. The larval weight (LW) of the common cutworm (CCW), the resistance level (RL) and the index of damaged leaf (IDL) were evaluated as resistance indicators to herbivores. An association synonymous SNP AX-94083016 located in the coding region of the respiratory burst oxidase gene GsRbohA1 was identified by genome-wide association study (GWAS) analyses. The overexpression of GsRbohA1 in soybean hairy roots enhanced resistance to CCW. One SNP in the promoter region cosegregated with AX-94083016 contributing to soybean resistance to CCW by altering GsRbohA1 gene expression and reactive oxygen species (ROS) accumulation. Two major haplotypes, GsRbohA1A and GsRbohA1G, were identified based on the SNP. The resistant haplotype GsRbohA1A predominates in wild soybeans, although it has been gradually lost in landraces and cultivars. The nucleotide diversity around GsRbohA1 is much lower in landraces and cultivars than in its ancestors. In conclusion, a new resistant haplotype, GsRbohA1A, was identified in wild soybean, which will be a valuable gene resource for soybean insect resistance breeding through introducing into improvement lines, and it offers a strategy for exploring resistance gene resources from its wild relatives.