
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237583
71357
10.1038/s41598-024-71357-8
Article
Genome-wide association study provides new insight into the underlying mechanism of drought tolerance during seed germination stage in soybean
Aleem Muqadas 123
Razzaq Muhammad Khuram 4
Aleem Maida 5
Yan Wenliang 1
Sharif Iram 6
Siddiqui Manzer H. 7
Aleem Saba 8
Iftikhar Muhammad Sarmad 29
Karikari Benjamin 10
Ali Zulfiqar 2
Begum Naheeda 1
Zhao Tuanjie tjzhao@njau.edu.cn

1
1 https://ror.org/05td3s095 grid.27871.3b 0000 0000 9750 7019 Key Laboratory of Biology and Genetics Improvement of Soybean, Ministry of Agriculture/Zhongshan Biological Breeding Laboratory (ZSBBL)National Innovation Platform for Soybean Breeding and Industry-Education Integration/State Key Laboratory of Crop Genetics & Germplasm Enhancement and UtilizationCollege of Agriculture, Nanjing Agricultural University, Nanjing, 210095 China
2 https://ror.org/054d77k59 grid.413016.1 0000 0004 0607 1563 Department of Plant Breeding and Genetics, University of Agriculture, Faisalabad, Pakistan
3 https://ror.org/054d77k59 grid.413016.1 0000 0004 0607 1563 Center for Advanced Studies in Agriculture and Food Security (CAS-AFS), University of Agriculture Faisalabad, Faisalabad, 38040 Pakistan
4 https://ror.org/00yh88643 grid.444934.a 0000 0004 0608 9907 Faculty of Agriculture and Veterinary Sciences, Superior University, Lahore, Pakistan
5 https://ror.org/054d77k59 grid.413016.1 0000 0004 0607 1563 Department of Botany, University of Agriculture, Faisalabad, Pakistan
6 Cotton Research Station, Faisalabad, Pakistan
7 https://ror.org/02f81g417 grid.56302.32 0000 0004 1773 5396 Department of Botany and Microbiology, College of Science, King Saud University, 11451 Riyadh, Saudi Arabia
8 Barani Agricultural Research Station, Fatehjang, Pakistan
9 https://ror.org/00rqy9422 grid.1003.2 0000 0000 9320 7537 School of Agriculture and Food Sciences, The University of Queensland, St Lucia, QLD Australia
10 https://ror.org/052nhnq73 grid.442305.4 0000 0004 0441 5393 Department of Crop Science, Faculty of Agriculture, Food and Consumer Sciences, University for Development Studies, PO Box TL 1882, Tamale, Ghana
5 9 2024
5 9 2024
2024
14 2076515 3 2024
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Drought is one of the major environmental issues that reduce crop yield. Seed germination is a crucial stage of plant development in all crop plants, including soybean. In soybean breeding, information about genetic mechanism of drought tolerance has great importance. However, at germination stage, there is relatively little knowledge on the genetic basis of soybean drought resistance. The objective of this work was to find the quantitative trait nucleotides (QTNs) linked to drought tolerance related three traits using a genome-wide association study (GWAS), viz., germination rate (GR), root length (RL), and whole seedling length (WSL), using germplasm population of 240 soybean PIs with 34,817 SNPs genotype data having MAF > 0.05. It was observed that heritability (H2) for GR, WSL, and RL across both environments (2020, and 2019) were high in the range of 0.76–0.99, showing that genetic factors play a vital role in drought tolerance as compared to environmental factors. A number of 23 and 27 QTNs were found to be linked to three traits using MLM and mrMLM, respectively. Three significant QTNs, qGR8-1, qWSL13-1, and qRL-8, were identified using both MLM and mrMLM methods among these QTNs. QTN8, located on chromosome 8 was consistently linked to two traits (GR and RL). The area (± 100 Kb) associated with this QTN was screened for drought tolerance based on gene annotation. Fifteen candidate genes were found by this screening. Based on the expression data, four candidate genes i.e. Glyma08g156800, Glyma08g160000, Glyma08g162700, and Glyma13g249600 were found to be linked to drought tolerance regulation in soybean. Hence, the current study provides evidence to understand the genetic constitution of drought tolerance during the germination stage and identified QTNs or genes could be utilized in molecular breeding to enhance the yield under drought stress.

Keywords

Soybean
Genome wide association study
SNPs
QTNs
Drought
Germination
Subject terms

Biotechnology
Genetics
Plant sciences
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Soybean (Glycine max (L.) Merr.) is among the world's five most important crops, providing food, feed, protein, and edible oil. However, it is subjected to several abiotic stresses like drought, salinity, and heat stress1. Among these stresses, drought stress is going to be a major threat to crop plants. In the current scenario, increased temperature, extreme variability of seasonal rainfall, and prolonged dry season due to climate change have contributed to the water scarcity challenge2. Water scarcity is becoming a major constraint for crop production throughout the world3. In soybean, drought stress mostly occurs during the whole growing season which leads to considerable yield loss and quality deterioration mainly in arid and semi-arid zones. Some studies have reported that drought stress could cause 24–50% yield losses in soybean4,5 The research studies have also shown that during juvenile and reproductive stages soybean is sensitive to drought stress6,7 In most of the soybean growing regions, irrigation is not a possible option due to lack of water availability and the high cost of irrigation. Therefore, the most economical way of sustainable soybean cultivation is to develop soybean genotypes having the ability to germinate and survived under drought stress condition8.

The most apparent symptoms of drought stress in soybean are reduction in germination percentage, cell elongation and expansion, hypocotyl length, shoot and root fresh weight, increased root length, leaf flipping, leaf clamping and initiation of early flowering5. Drought stress also disturbs carbon assimilation and plant phenology9. The germination and seedling stages are the most sensitive to drought stress in soybean and other crops like rice10, maize11 and wheat12. To combat with the challenge of drought stress, plants have adopted various morphological, physiological and molecular mechanisms including the robust root system, reduced leaf area and stomatal conductance, maintenance of turgor pressure, early maturation, activation of drought responsive genes and transcriptional factors e.g. DREBs, WRKY, ZIP, ZFP, NAC, MYB and ERFs13–16. It is predicted that in future, plants will face drought stress more frequently due to high temperature extremes17. Therefore, efforts are necessary to explore the genetic base and genetic mechanism regulating the drought tolerance in soybean to combat climate change.

Drought tolerance is a polygenic trait influenced by genotypic and environmental factors18. Due to complex genetic architecture of drought tolerance, less genetic diversity, segregation distortion, linkage drag, and distortion and suppression of recombination conventional breeding fails to meet the requirement for stabilizing yield potential in plants under variable conditions19–22. Advances in genetic and genomic approaches have assisted in the dissection of genetic control of various complex traits which have facilitated mapping of numerous quantitative trait loci (QTLs) in different crops, including soybean23,24. However, QTL mapping provides inadequate genetic diversity as reported in Oryza sativa25. To far, numerous genetic maps have been created to investigate QTLs that affect drought tolerance in soybeans. But in practical breeding, confirmation and verifications of QTLs are recommended prior to their application in the breeding program26.

Advancement in next-generation sequencing (NGS) has facilitated exploiting marker assisted selection (MAS) to enhance the genetic gains in an economical way. Genome wide association studies (GWAS) by utilizing genome wide single nucleotide polymorphism (SNPs) have the potential to detect variants that could be used directly for crop plant improvement23,27. GWAS could overwhelm the constraints of traditional QTL mapping and can be utilized to investigate the minor gene effects controlling complex traits in crop plants28. GWAS has been successfully used in breeding of various crops like cotton, wheat, barely, rice and corn13,23,29,30.

Association mapping accuracy relies on marker density and mapping populations31. Among all the approaches utilized for GWAS, the mixed linear model (MLM) is applied widely for association mapping32. Furthermore, advanced MLM- based methods have also been established like compressed MLM (cMLM), which saves time for computation and enhances the statistical performance. Bonferroni correction is indeed for such tests since these approaches rely on single-locus genome-wide scanning32. Currently, a more refined method is purposed known as random-SNP impact MLM (mrMLM), which combines the advantages of both multiple locus model and single locus genome scanning. The mrMLM revealed significantly better accuracy and power in the detection of QTN and affected the estimation24.

Several GWAS have been conducted and reported in soybean33,34. Still, literature related to GWAS for drought tolerance at germination is limited35. Germination is a critical step contributing toward final soybean production. Cultivars with a better germination percentage during drought stress produced higher yields36. The information is scarce about molecular genetics controls drought tolerance during the germination stage. As a result, it's critical to identify the QTLs linked to drought tolerance in soybeans during the germination stage.

The present research includes marker-trait association depending upon phenotypic and genotypic information from 240 worldwide germplasm accessions with wide diversity. The 240 accessions were screened for drought tolerance at germination stage with popular method for phenotyping. We used both single and multi-locus models viz., MLM and mrMLM to detect the significant SNPs/QTNs, as well as for the prediction of candidate genes linked to drought tolerance during germination in soybean. Physiological mechanisms associated with drought tolerance at the germination stage may be akin to those that help plants to survive under drought stress during seedling and reproductive stages. The objective of present study is to unravel the genetic constitution of drought tolerance in 240 soybean accessions at the germination stage using GWAS.

Results

Phenotypic characterization of the 240 germplasm

The drought tolerance of 240 soybean accessions was assessed on the base of three germination linked traits i.e. RL, GR and WSL. The values of mean, standard deviation, range, ANOVA (F-value) and estimated h2 for all traits of 240 lines in two different years 2018 and 2019 are shown in Table 2-1. During 2018, the average value of GR, WSL and RL were 0.30, 0.32 and 0.43 and their range varied from 0–1.25, 0.11–1.75 to 0.0028–0.96 respectively. However, the mean comparison revealed that, during 2019, the average values of GR and WSL were higher than the mean values in 2018 whereas RL growth pattern was found contradictory to it (higher in 2018). These findings depicted that, in soybean, drought tolerance is influenced by environment. The results of ANOVA indicated highly significant difference between genotypes (G) and environments (E). Furthermore, significant genotype× environment (G×E) interaction was detected for all three studied traits (Table 1). Furthermore, h2 for GR, WSL and RL across both environments were high in the range of 0.76–0.99, showing that genetic factors play a vital role in drought tolerance as compared to environmental factors.Table 1 Descriptive statistics, heritability (H2) and F values from ANOVA of three drought tolerance indices in 2 years (2018 and 2019).

Trait	Env	Mean ± SD	Range	H2	F-values from ANOVA	
Genotype	Env	G × E	
GR	2018	0.30 ± 0.19	0–1.25	0.99	360.3**	351.69**	69.99**	
2019	0.32 ± 0.25	0.02–2.7	
WSL	2018	0.32 ± 0.15	0.11–1.75	0.99	77.22**	702.15**	79.54**	
2019	0.36 ± 0.14	0.01–1.05	
RL	2018	0.43 ± 0.19	0.003–0.96	0.76	136.14**	945.039**	72.76**	
2019	0.36 ± 0.29	0.004–1.02	
Env environment, SD standard deviation, 2018 one environment, 2019s environment, G × E genotype by environment interaction.

**Significance at P < 0.01, H2 broad sense heritability.

SNPs distribution, population structure and LD estimation

In this study, out of a total 42,509 SNPs, 17.98% SNPs having MAF ≤ 0.05 were extracted for analyses, while the leftover 34,817 SNPs with MAF > 0.05 were used to investigate the genetic variability in 240 soybean genotypes and other subsequent analyses. These targeted SNPs were randomly located on 20 chromosomes of G. max (Fig. 1A). On an average, each chromosome (Chr) carried 1762 SNPs. Maximum number of SNPs (2922) were mapped on Chr.18 while minimum number of SNPs (1223) were found on Chr. 12. Furthermore, the SNPs distribution density on 20 chromosomes is presented in Fig. 1B that varied from 25.87 SNPs/Mb on Chr.01 to 51.01 SNPs/Mb on Chr.13.Fig. 1 Distribution of 34,817 SNPs across 240 soybean PIs. (A) Number of SNPs per chromosome (B) SNP density distribution (SNPs/Mb).

Besides this, three different methods, principal component analysis (PCA), phylogenetic analysis, and bayesian model-based methods were employed to build the population structure of 240 soybean accessions. STRUCTURE software was used to compute Delta K (∆K) value (Fig. 2A; K = 1–10, Fig. 2B), and the result revealed that 240 soybean PIs were grouped into three subpopulations i.e. K = 3. Moreover, PCA analysis and phylogenetic tree showed consistent results which agree with the findings of population structure analysis (Fig. 2C, D).Fig. 2 An examination of the population structure of 240 soybean lines. (A) ΔK plot calculated for K = 1–10 (B) structure software was used to determine the population structure. Three subpopulations are represented by three colors (red, green, and blue) (C) the tested 240 lines’ maximum likelihood neighbor-joining tree (D) principal component analysis plot.

Furthermore, the r2 values were calculated for those SNPs pairs that resided within 10 Mb region and the LD decay rate was established based on the regression of negative natural logarithm. In G. max, LD decay was comparatively higher37 than other crop plants like A. thaliana38, O. sativa39, Z. mays40. In addition, our results showed that r2 values decreased gradually as the distance increases and At 1.76 Mb, the LD decay distance was measured, and r2 dropped to half of its highest value (Fig. 3).Fig. 3 Linkage disequilibrium (LD) decay across the genome in a population of 240 soybean lines. The squared correlation coefficient was used to get the mean LD decay rates.

Association mapping for traits via MLM

The MLM model was utilised in GWAS research to identify important QTNs related to GR, WSL, and RL. Results highlighted that a total of 4, 9 and 10 QTNs were found to be associated with GR, WSL, and RL, respectively, at a threshold value of − log10P = 4.5 (Table 2 and Fig. 4, Fig. 5, Fig. 6 respectively). Four QTNs associated with GR were mapped on two chromosomes only: Chr.03 and Chr.08. Among these QTNs, two QTN qGR3-1 and qGR8-1 were distributed on Chr.3 and 8 respectively and covered 8.55 to 10.59% phenotypic variation (PV) and they were identified consistently in all three environments (i.e. 2018, 2019 and combined environment) (Table 2, Fig. 4). These QTNs showed significant association with SNP marker Gm03-4741942, Gm08-12368145 and exhibited largest phenotypic variation effect (10.59%) in combined environment. Besides this, qGR3-2 and qGR8-2 were detected in two environments only and they covered 8.99–9.14% of the PV (Table 2).Table 2 Quantitative trait nucleotides (QTNs) associated with germination rate (GR), whole seedling length (WSL) and root length (RL) via MLM.

Trait	QTN	Env. a	Chrb	Marker	Position (bp)	P value	 − log10 Pc	Allele	R2 (%)	
GR	qGR3-1	2018	3	Gm03-4741942	4741942	1.09E-05	4.96	C/T	8.55	
qGR8-1	2018	8	Gm08-12368145	12368145	1.09E-05	4.96	C/T	8.55	
qGR3-1	2019	3	Gm03-4741942	4741942	1.35E-06	5.87	C/T	10.22	
qGR3-2	2019	3	Gm03-4801039	4801039	6.04E-06	5.22	A/G	8.99	
qGR8-1	2019	8	Gm08-12368145	12368145	1.35E-06	5.87	C/T	10.22	
qGR8-2	2019	8	Gm08-12513720	12513720	6.04E-06	5.22	A/G	8.99	
qGR3-1	Combined	3	Gm03-4741942	4741942	9.62E-07	6.02	C/T	10.60	
qGR3-2	Combined	3	Gm03-4801039	4801039	5.82E-06	5.24	A/G	9.14	
qGR8-1	Combined	8	Gm08-12368145	12368145	9.62E-07	6.02	C/T	10.60	
qGR8-2	Combined	8	Gm08-12513720	12513720	5.82E-06	5.24	A/G	9.14	
WSL	qWSL8	2018	8	Gm08-18192704	18192704	1.28E-05	4.89	A/G	8.56	
qWSL13-1	2018	13	Gm13-35394159	35394159	1.92E-05	4.72	A/G	8.22	
qWSL13-2	2018	13	Gm13-35545045	35545045	9.51E-06	5.02	G/T	8.84	
qWSL13-3	2018	13	Gm13-35552015	35552015	1.01E-05	5.00	C/T	8.78	
qWSL16	2018	16	Gm16-3107434	3107434	1.28E-05	4.89	A/G	8.56	
qWSL8	2019	8	Gm08-18192704	18192704	1.33E-05	4.88	A/G	8.52	
qWSL9	2019	9	Gm09-35829245	35829245	1.33E-05	4.88	A/G	8.52	
qWSL12	2019	12	Gm12-2710068	2710068	7.48E-06	5.13	A/G	9.05	
qWSL13-1	2019	13	Gm13-35394159	35394159	2.41E-05	4.62	A/G	8.02	
qWSL18-1	2019	18	Gm18-224527	224527	2.08E-05	4.68	C/T	8.19	
qWSL18-2	2019	18	Gm18-51452390	51452390	7.85E-06	5.11	A/C	9.00	
qWSL8	Combined	8	Gm08-18192704	18192704	1.19E-05	4.93	A/G	8.63	
qWSL9	Combined	9	Gm09-35829245	35829245	1.19E-05	4.93	A/G	8.63	
qWSL12	Combined	12	Gm12-2710068	2710068	7.48E-06	4.92	A/G	9.05	
qWSL13-1	Combined	13	Gm13-35394159	35394159	1.99E-05	4.70	A/G	8.19	
qWSL13-2	Combined	13	Gm13-35545045	35545045	7.56E-06	5.01	G/T	9.04	
qWSL13-3	Combined	13	Gm13-35552015	35552015	7.99E-06	4.98	C/T	8.98	
qWSL16	Combined	16	Gm16-3107434	3107434	1.28E-05	4.92	A/G	8.56	
qWSL18-1	Combined	18	Gm18-224527	224527	2.08E-05	4.62	C/T	8.19	
qWSL18-2	Combined	18	Gm18-51452390	51452390	7.85E-06	4.92	A/C	9.00	
RL	qRL-6	2018	6	Gm-06-12538994	12538994	1.56E-05	4.81	C/T	8.35	
qRL-8	2018	8	Gm08-12368145	12368145	1.56E-05	4.81	C/T	8.35	
qRL15-4	2018	15	Gm-15-4497855	14841750	2.29E-05	4.64	C/T	9.84	
qRL16-1	2018	16	Gm-16-17726124	17726124	1.95E-05	4.71	C/T	8.15	
qRL-8	2019	8	Gm08-12368145	12368145	8.49E-06	5.07	C/T	8.90	
qRL-14	2019	14	Gm-14-28361034	28361034	1.34E-05	4.87	A/G	8.49	
qRL15-2	2019	15	Gm-15-4485355	4485355	2.22E-05	4.65	A/G	8.08	
qRL15-3	2019	15	Gm-15-4494390	4494390	2.16E-05	4.67	G/T	8.07	
qRL15-4	2019	15	Gm-15-4497855	4497855	2.02E-05	4.69	A/G	8.12	
qRL15-5	2019	15	Gm-15-4509223	4509223	2.02E-05	4.69	C/T	8.12	
qRL15-6	2019	15	Gm-15-4522374	4522374	2.02E-05	4.69	A/C	8.12	
qRL16-2	2019	16	Gm-16-23995972	23995972	2.18E-05	4.66	A/G	8.06	
qRL-6	Combined	6	Gm-06-12538994	12538994	8.55E-06	4.70	C/T	8.88	
qRL-8	Combined	8	Gm08-12368145	12368145	8.55E-06	5.07	C/T	8.88	
qRL15-1	Combined	15	Gm-15-4462597	4462597	2.69E-05	4.57	G/T	8.02	
qRL15-2	Combined	15	Gm-15-4485355	4485355	2.17E-05	4.66	A/G	8.11	
qRL15-3	Combined	15	Gm-15-4494390	4494390	2.20E-05	4.66	G/T	8.05	
qRL15-4	Combined	15	Gm-15-4497855	4497855	1.99E-05	4.70	A/G	8.13	
qRL15-5	Combined	15	Gm-15-4509223	4509223	1.99E-05	4.70	C/T	8.13	
qRL15-6	Combined	15	Gm-15-4522374	4522374	1.99E-05	4.70	A/C	8.13	
qRL16-1	Combined	16	Gm-16-17726124	17726124	2.05E-05	4.69	C/T	8.10	
qRL16-2	Combined	16	Gm-16-23995972	23995972	2.34E-05	4.63	A/G	7.99	
aDifferent environments, 2018 and 2019. Combined combined-environments 2018 and 2019.

bChromosome.

c − log10 transformed P-value.

Fig. 4 Manhattan plots and quantile–quantile (QQ) plots were used in the population-based genome-wide association studies for (A) GR-2018. (B) GR-2019. (C) GR-Combined.

Fig. 5 Manhattan plots and quantile–quantile (QQ) plots were used in the population-based genome-wide association studies for (A) WSL-2018. (B) WSL-2019. (C) WSL-Combined.

Fig. 6 Manhattan plots and quantile–quantile (QQ) plots were used in the population-based genome-wide association studies for (A) RL-2018. (B) RL-2019. (C) RL-Combined.

For WSL, a total of 9 QTNs were identified on six different chromosomes i.e. Chr.08, 09, 12, 13, 16 and 18 (Table 2, Fig. 5). Among these, qWSL-8 and qWSL13-1 were found in all the three environments and categorized as stable QTNs linked with WSL. These two QTNs accounted for 8.01%–8.63% phenotypic variations, respectively. Two QTNs linked to WSL viz., qWSL13-2 and qWSL13-3 were identified in two environments and explicated 8.84–9.77% of the PV respectively. However, qWSL-16–1 and qWSL-9–1 were detected during 2018 and 2019, respectively (Table 2).

We used the MLM approach to identify 10 QTNs in RL, and all these QTNs were located on 5 distinct chromosomes viz., Chr.06, 08, 15, 14 and 16 (Table 2). Out of these QTNs, only qRL-8 and qRL15-4 were continually detected in all environments (2018, 2019 and Combined), elucidating 8.12–9.83% of the PV and these QTNs were linked SNP marker Gm08-12368145 on Chr.08 and Gm15-4497855 on Chr.15, respectively. Other QTNs, qRL16-1 and qRL-14 were detected in 2018 and 2019 environments, explaining 8.14% and 8.49% of the phenotypic variance, respectively. The three remaining QTNs associated with RL were only found in one environment and accounted for 8.08 percent–8.12 percent of the PV (Table 2 and Fig. 6).

Association mapping for traits via mrMLM

Multi-locus GWAS analysis was performed by a mrMLM approach to ensure the QTNs dependability detected by the MLM model and to identify further QTNs related with GR, WSL and RL. As shown in Table 2-3, a total of 8, 8 and 11 QTNs were linked with GR, WSL and RL respectively (Table 3). For GR, one QTN qGR8-1 linked with the Gm08-12368145 marker was constantly found in all environments, and this QTN showed broad PV effects ranging from 8.40 to 9.95%. The remaining seven GR QTNs were only found in one environment, accounting for 4.04 to 9.95% of the PV.Table 3 Quantitative trait nucleotides (QTNs) associated with germination rate (GR), whole seedling length (WSL) and root length (RL) via mrMLM.

Trait ID	QTN	Env.a	Chr.b	SNP	Position (bp)	QTN effect	LOD score	 − log10Pc	R2 (%)	MAF	Genotype	
GR	qGR4	2018	4	Gm04-8805064	8805064	0.084	6.329	7.173	7.903	0.168	A	
qGR8-1	2018	8	Gm08-12368145	12368145	0.073	7.536	8.415	9.9567	0.383	A	
qGR10	2018	10	Gm10-2912266	2912266	 − 0.052	3.314	4.029	4.3644	0.269	G	
qGR19-1	2018	19	Gm19-45138884	45138884	 − 0.055	3.096	3.797	4.2803	0.208	C	
qGR6	2019	6	Gm06-11913935	11913935	 − 0.075	4.288	5.054	8.0296	0.351	T	
qGR8-1	2019	8	Gm08-12368145	12368145	0.071	5.418	6.230	8.4029	0.383	A	
qGR9	2019	9	Gm09-35829245	35829245	 − 0.075	5.289	6.097	6.0674	0.156	G	
qGR16	2019	16	Gm16-33529724	33529724	0.077	5.452	6.266	7.6347	0.185	A	
qGR2-2	Combined	2	Gm02-2861258	2861258	0.066	3.907	4.654	7.1115	0.311	T	
qGR8-1	Combined	8	Gm08-12368145	12368145	0.062	6.741	7.597	8.9977	0.383	A	
qGR11	Combined	11	Gm11-9062549	9062549	 − 0.059	3.532	4.259	4.0415	0.198	G	
WSL	qWSL13-1	2018	13	Gm13-35394159	35394159	 − 0.069	4.592	5.372	8.068	0.109	A	
qWSL14	2018	14	Gm14-47884956	47884956	0.041	4.553	5.331	6.2969	0.324	A	
qWSL16	2018	16	Gm16-18589836	18589836	0.043	5.166	5.969	4.9986	0.191	G	
qWSL18-2	2018	18	Gm18-51452390	51452390	 − 0.025	4.073	4.828	2.6056	0.477	A	
qWSL20-1	2018	20	Gm20-36122830	36122830	 − 0.040	4.709	5.493	5.7337	0.324	A	
qWSL12	2019	12	Gm12-2710068	2710068	 − 0.063	3.655	4.388	6.4491	0.097	G	
qWSL13-1	2019	13	Gm13-35394159	35394159	 − 0.062	4.367	5.136	7.1721	0.109	A	
qWSL18	2019	18	Gm18-224527	224527	0.050	4.625	5.405	6.0322	0.162	T	
qWSL13-1	Combined	13	Gm13-35394159	35394159	0.033	4.288	5.054	8.0296	0.351	A	
qWSL18-1	Combined	18	Gm18-224527	224527	0.052	6.882	7.743	6.5292	0.162	T	
qWSL20-2	Combined	20	Gm20-34381893	34381893	 − 0.053	5.721	6.545	8.5594	0.214	A	
RL	qRL6	2018	6	Gm06-12538994	12538994	 − 0.072	4.716	5.500	6.0961	0.130	C	
qRL-8	2018	8	Gm08-12368145	12368145	0.061	6.624	7.477	8.4509	0.315	G	
qRL15-1	2018	15	Gm15-14841750	14841750	 − 0.089	5.418	6.230	13.4029	0.132	T	
qRL16-2	2018	16	Gm16-17726124	17726124	0.057	4.591	5.370	5.9049	0.223	C	
qRL19	2018	19	Gm19-31536655	31536655	 − 0.062	3.246	3.956	5.0791	0.128	T	
qRL7-1	2019	7	Gm07-6015318	6015318	 − 0.102	7.625	8.508	18.6021	0.278	C	
qRL7-2	2019	7	Gm07-6604493	6604493	0.062	6.741	7.597	7.9977	0.383	A	
qRL-8	2019	8	Gm08-12368145	12368145	0.055	4.033	4.786	5.6986	0.316	G	
qRL14	2019	14	Gm14-28361034	28361034	0.058	6.385	7.230	7.451	0.483	G	
qRL20	2019	20	Gm20-40765691	40765691	0.061	5.564	6.382	6.0109	0.237	G	
qRL7-1	Combined	7	Gm07-6015318	6015318	 − 0.077	6.501	7.350	11.33	0.278	C	
qRL7-2	Combined	7	Gm07-6604493	6604493	0.075	7.555	8.435	12.5244	0.383	A	
qRL-8	Combined	8	Gm08-12368145	12368145	0.070	7.782	8.668	9.8464	0.316	G	
qRL13	Combined	13	Gm13-29961821	29961821	 − 0.047	3.077	3.778	3.7954	0.254	T	
qRL16-1	Combined	16	Gm16-5586574	5586574	0.046	3.334	4.050	4.6205	0.376	G	
MAF minor allele frequency.

aDifferent environments, 2018 and 2019; Combined combined environment 2018 and 2019.

bChromosome.

c − log10 transformed P-value.

In case of WSL, 8 QTNs were identified, scattered on 6 different chromosomes. Among them, only single QTN qWSL-13 was mapped in all three environments with larger effect, while qWSL-18 was identified in two environments (2019 and combined). However, the rest of 6 QTNs were identified in single environment only. For RL, a total of eleven QTNs were detected, explaining 3.7–18.6% of the PV. Out of these, qRL-8-1 associated with RL in all environments (2018, 2019 and combined). Two QTNs qRL7-1 and qRL7-2 were observed in two environments (2019 and combined). However, the remaining eight QTNs were environment specific QTNs (Table 3). In short, a total of 23 and 27 QTNs linked with all studied traits of drought tolerance were observed via MLM and mrMLM model, respectively (Tables 2 and 3).

Comparison between MLM and mrMLM results

In comparison with the MLM and mrMLM, For MLM, six QTNs, viz., qGR3-1, qGR8-1, qWSL8, qWSL13-1, qRL-8 and qRL15-4 were detected in all three environments. As predicted, for mrMLM, three major QTN qGR8-1, qWSL13-1 and qRL-8 were linked with GR, WSL and RL were detected. These three critical QTNs were detected by both MLM and mrMLM. Nevertheless, the qGR3-1 and qWSL8-1 detected by MLM could not be confirmed by mrMLM. Furthermore, qRL6 and qRL14 were also detected by both MLM and mrMLM in 2018 and 2019, respectively (Table 3). Based on above mentioned results, two major and stable QTN i.e. qGR8-1(Gm08-12368145), qWSL13-1(Gm13-35394159) and qRL8 (Gm08-12368145) were identified in all studied environments via MLM and mrMLM methods. The QTNs, qRL8 and qGR8-1 were renamed as QTN8. QTN qRL15-4 (Gm15-4497855) was identified in all three environments in MLM but not in mrMLM. Thus, these major QTN on Chr.8 and 13 were potentially employed as the hotspot area for fine mapping and analysis of candidate genes linked with drought tolerance in soybean.

Candidate gene prediction

As the stable QTN found on Chr.08 (QTN8) and Chr.13 (qWSL13-1) were constantly detected in all studied traits i.e. GR, WSL and RL by both methods in all environments. For further confirmation, candidate gene prediction analysis has been conducted in the genomic region adjacent to the QTN8-1 and QTN13-1. Based on LD distance, 500 kb up as well as downstream region carrying peak SNPs markers was explored and within the interval of 1.0 Mb, a total of 100 and 139 genes were mapped on Chr.08 and Chr.13, respectively. The haplotype block analysis was carried out with Haploview 4.2 for all the SNPs located within ± 100 kb candidate region (Fig. 7). Based on gene annotation data derived from Soybase (http://www.soybase.org) and the Phytozome (https://phytozome.jgi.doe.gov), 15 genes were identified as possible candidate genes regulating the drought tolerance in soybean (Table2-4). Out of 15 candidate genes, Glyma08g156800, Glyma08g160000, Glyma08g162700 and Glyma13g249600 only showed higher expression in soybean root tissues.Fig. 7 Candidate regions and haplotype blocks of the genome significantly linked to identified major QTNs for GR and RL on Chr.08 and for WSL on Chr.13 in 240 soybean PIs.

Table 4 In two distinct environments, phenotypic correlations (r) between GR, WSL, and RL.

Trait	GR	WSL	RL	
GR	1	0.634	0.797	
WSL	0.634	1	0.582	
RL	0.797	0.582	1	

Discussion

In this study, significant diversity has been observed in 240 soybean PIs regarding the GR, WSL and RL values that are three important drought-linked parameters at germination stage. Results depicted that GWAS analysis using 34,817 SNPs detected QTLs related with drought tolerance (Table 2, Table 3). Furthermore, strong positive correlations of r = 0.79 was observed between GR and RL in the 240 soybean PIs (Table 4). Drought stress impairs germination by limiting water imbibition41. Therefore, for evaluation of soybean genotypes, we screened the population at germination stage by measuring GR, WSL and RL. On the basis of correlation analysis, it was cleared that WSL showed a positive correlation with GR and RL (Table 4).

Two stable QTNs were discovered among the QTNs identified, viz. QTN8, (qGR8-1, qRL-8) and qWSL13-1 were detected in all environments (2018, 2019, and Combined) via MLM, explaining 8.34–10.59% of the PV, whereas these QTNs (QTN8 and qWSL13-1) were also detected by mrMLM. Remarkably, Stable QTN (QTN8) located on Chr.08, was linked with SNP marker Gm08-12368145 for two traits (qGR8-1 and qRL-8), this QTN was discovered in all of the environments examined. Furthermore, mrMLM was used to effectively identify this stable QTN8 on Chr. 8. Another major QTN (qWSL13-1) associated with Gm13-35394159 was detected in all environments (Table 3). Thus, this stable QTN8 detected on Chr.08 in present study is a newly QTN linked with drought-tolerance in soybean.

However, the second stable QTNs viz., qWSL13-1 identified via MLM was detected in all environments, was only linked to a single trait. Therefore, on the base of above mentioned results, it is summarized that in this research all the identified stable QTNs were novel for drought tolerance in soybean and mostly carried larger phenotypic effects. Besides this, QTN8 and qWSL13-1 exerted larger phenotypic effects on drought tolerance. That’s why, we carried out candidate gene prediction within 1.0 Mb region adjacent to the QTN8/Gm08-12368145 and Gm13-35394159. Based on functional annotations, fifteen candidate genes which are directly or indirectly linked with drought tolerance were selected (Table 5).Table 5 Functional annotation of candidate genes.

Genes	Position (bp)	Homolog in Arabidopsis	Functional annotation	
Glyma.08g154900	11,953,876–11,957,663	AT5G13490.1	ADP,ATP Carrier protein 1	
Glyma.08g155200	12,000,644–12,005,368	AT5G13450.1	F-type H + -transporting ATPase subunit O	
Glyma.08g156800	12,115,786–12,130,558	AT2G13680.1	Callose synthase	
Glyma.08g157700	12,213,783–12,222,299	AT2G44350.1	Citrate Synthase	
Glyma.08g158600	12,290,252–12,297,255	AT2G44280.1	Major Facilitator superfamily domain containing protein	
Glyma.08g159800	12,386,775–12,390,398	AT3G04770.2	40S ribosomal protein SA	
Glyma.08g159900	12,394,459–12,399,478	AT1G10670.1	ATP-citrate synthase	
Glyma.08g160000	12,401,898–12,405,089	AT1G79260.1	THAP4 PROT	
Glyma.08g162200	12,635,478–12,640,656	AT1G79270.1	Methylated RNA-binding protein	
Glyma.08g162700	12,736,920–12,739,132	AT1G30870.1	Peroxidase	
Glyma.13g241300	35,130,654–35,134,130	AT2G20760.1	Clathrin light chain protein	
Glyma.13g241500	35,148,816–35,152,426	AT1G16700.1	NADH DEHYDROGENASE IRON-SULFUR PROTEIN 8,	
Glyma.13g243100	35,256,095–35,257,810	AT3G49910.1	60S RIBOSOMAL PROTEIN L26 FAMILY MEMBER	
Glyma.13g244800	35,408,963–35,410,102	AT4G21620.1	Glycine-rich protein	
Glyma.13g249600	35,726,229–35,729,214	AT4G27270.1	1,4-benzoquinone reductase-like	

Up till now, various studies have been carried out in soybean to explore the genetic mechanism that regulates drought tolerance. Several QTLs linked with drought tolerance have been identified across different environments42–44. However, most of the research related to QTLs mapping was carried out by employing linkage mapping analysis on the basis of low-density genetic maps. Drought tolerance is a complex character as controlled by several minor loci with significant environmental influence. To overcome the linkage mapping limitations, LD based GWAS analysis was conducted in this research. Several studies have confirmed that utilization of GWAS analysis could be an alternative and powerful approach for precise detection and defining the genomic regions associated with quantitative traits. Utilization of GWAS is an efficient alternative approach for precise detection and marking the genomic region of loci/genes linked with complex traits16,45.

Linkage disequilibrium refers to non-random association of SNP alleles that are located at different locus. Nonetheless, the amount of genome-wide LD is critical in GWAS research. As it is influenced by several forces viz., population stratification, linkage, genetic drift, and relatedness13,46. For the soybean panel, the LD decay was approximately assessed at 1.76 Mb, which was markedly higher as compared to other plant species i.e. rice (75 kb in indica)39, Arabidopsis thaliana (10 kb)39 and 1.1 kb in sunflower47. In the self-pollinated species like soybean comparatively higher LD value has been recorded as compared to cross-pollinated species48,49. Furthermore, resequencing of 31 soybean genomes, Lam, Xu, Liu, Chen, Yang, Wong, Li, He, Qin and Wang37 stated that both cultivated and wild soybean species showed high LD, while the comparison of both species showed that cultivated soybean demonstrated higher LD level than wild species. Therefore, for effective improvement of the genetic variability of GWAS population genotypes should have diverse genetic background (including wild species, cultivars, and landraces) collected from different geographical regions.

From previous functionally verified genes, transgenic lines with GmNAC8 overexpression show noticeably greater drought tolerance than those with GmNAC8 deletion, indicating that GmNAC8 is a positive regulator of drought tolerance50. GmWRKY54 enhanced drought tolerance by activating genes in the abscisic acid and Ca2+ signaling pathways from functionally confirmed genes. In Arabidopsis, overexpression of WRKY20 from a wild soybean greatly increased drought tolerance51,52. Using a constitutive expression promoter, it was possible to overexpress the LOS5/ABA3 gene in soybeans, encoding the molybdenum cofactor sulfurase, therefore improving drought tolerance and increasing seed output by at least 21% during drought stress53. In plants, U-box (PUB) proteins serve as E3 ligases. Under drought stress conditions, the transgenic lines of Arabidopsis showed reduced plant viability when GmPUB6 was overexpressed54.

Manavalan, Prince, Musket, Chaky, Deshmukh, Vuong, Song, Cregan, Nelson and Shannon43 has reported one QTN associated with drought tolerance linked with marker Satt315 on Chr.08 within a confidence interval of 15 cm. Depending upon the genomic sequence information, it was found that this QTL was linked with Satt315 from SoyBase (http://www.soybase.org), this QTL was present between 6751690 and 6751725 Mb of physical interval on Chr.08. In another study, Abdul Haleem et al. (2011) identified one QTL Q_root_Gm08 (flanked with Satt228 and Satt429 on Gm08) via MLM relevant to the fibrous root system. This QTL was located between 46366169 and 46366092 Mb of physical interval of Chr.08. (QTN8)55 also reported a cluster of 4 QTL associated with seed and root related traits (RFW, RDW, SFW, and SDW) on LG A2 (chromosome 8). This QTL is linked with marker Satt177 within a confidence interval of 15 cM. Charlson et al. (2009) also identified QTL (Sat_319) on MLG A2 contributed slow-wilting. This QTL was located between 4717934 and 4717907 Mb of physical interval of Chr.08. Wang et al., (2012) also identified QTL (I-Satt327 on G2-A2) controlling drought tolerant traits during germination stage. This QTL was located between 20511527 and 20511574 Mb of physical interval of Chr.08 but the physical position of Gm08-12368145 was 12.36 Mb of Chr.08.

Materials and methods

Plant materials and field evaluation

An association mapping population comprised 240 worldwide (Table S2-1) soybean plant introductions (PIs). The collections were obtained from 22 countries across Europe (EU), West Asia and North Africa, East Asia, and the Americas. The list of germplasm accessions and their origins are presented in Supplementary Table 1. Seeds of the 240 germplasm were collected from National Center for Soybean Improvement, Nanjing Agricultural University, Nanjing–China. The seeds obtained were sown in Jiangpu Experimental station of Nanjing Agricultural University in 2018 and 2019. The healthy and uniform seeds were selected for drought tolerance screening in the laboratory.

Germination test

To test the germination percentage, twenty seeds per entry were randomly selected for each of three replications. For evaluation of growth performance and phenotypic diversity related to seedling traits, seeds were kept on wetted filter paper in a 9 cm diameter petri dish. Drought stress was created by the addition of Polyethylene glycol 1000 (PEG-1000, 15% w/v) whereas distilled water (without PEG) was employed in control conditions. The Petri dishes were retained in an incubator at 25 ℃ in the darkness. The experiment was carried out in complete randomized design (CRD). The seeds were categorized as germinated when radicle attained the length of 2 mm. Germination was recorded at 24 h’ intervals for five consecutive days. To evaluate germination traits, we measured germination rate (GR, whole seedling length (WSL) and root length (RL) under drought and control conditions. Germination related traits were scored as per International Seed Testing Association rules (ISTA) (ISTA, 2003):Germination percentage=n/N×100

where n and N represent number of seeds germinated and the total number of seeds sown, respectively.

For measurement of whole seedling length (WSL) and root length (RL), three plants of each accession were taken under control conditions and 15% PEG stress. After five days, RL and WSL (in cm) were measured manually by scaled ruler56. The data for these traits was utilised to create indices, which were calculated as the ratio of GR, WSL, and RL, germination traits, at drought conditions to the control conditions of same trait. The germination rate (GR) was scored by using the following equations:GR=Germination rate under drought conditionGermination rate under normal condition

Data analyses

Using SAS PROC UNIVARIATE tools, descriptive statistics such as mean, standard deviation (SD), and range were calculated for 240 soybean PIs in two environments: 2018 and 2019 respectively. The PROC GLM and PROC CORR programmes of SAS were used to calculate analysis of variance (ANOVA) and Pearson’s correlation for drought traits57. The broad sense heritability of studied traits under drought and normal condition was calculated for combined environment using the below formulaH2=σ2g(σ2g+σ2gen+σ2enr)

For individual environment heritability was calculated by using the following formula:H2=σ2g(σ2g+σ2e)

where σ2g is the genotypic variance, σ2ge is interaction variance between genotypes and environment, σ2e is the error variance, n is the number of environments, r and e represent replications and error respectively58.

Genotyping and quality control

The SNP data for the 240 accessions is from SoySNP50K BeadChip available on SoyBase database (soybase.org/snps/index.php) with 42,449 SNPs59. Briefly, a high-throughput SNP assay, the SoySNP50K iSelect BeadChip containing 42,449 SNPs was utilized for genotyping following procedure reported by60. Chromosomal distributions of SNPs and quality control were estimated following the procedure of Wen et al. (2014). A subset of 34,817 SNPs (MAF > 0.05 and missing rates < 0.25) was selected out of 42,449 SNPs for further analyses.

Principal component, linkage disequilibrium, and population structure analysis

After screening using the PLINK in-depth-pairwise command 7358 SNPs were retrieved for this study61. Using the software STRUCTURE 2.2, a model-based cluster assay was performed to find the genetic structure and define the cluster (gene pools) in the data62. Potential populations (K), the number was from 1 to 10 with four repetitions, and each Q and P-value was estimated. The Delta K method was used to calculate the most likely number of subpopulations63. TASSEL 5.0 software was used to create a neighbor-joining tree (NJT)64. R studio was used to perform principal component analysis (PCA). TASSEL 5.0 was used to estimate Kinship, and further pairwise SNPs linkage disequilibrium (LD) was estimated as a squared correlation coefficient (r2) using RTM-GWASV1.165. The distance of chromosome at which r2 dropped to half of its highest value was used to determine the linkage disequilibrium decay rate.

Marker-trait association and haplotype block analysis (GWAS)

After eliminating SNPs with an MAF < 5%, 34,817 SNPs were remained for the 240 PIs. To reduce false-positives and enhance statistical performance, the principal component (PC) and kinship (K) matrix were calculated. For the SLM (Single locus model) viz., MLM (PCA + K) was carried out in TASSEL5.2 (Bradbury et al., 2007) while for multi-locus model, mrMLM was conducted in R using mrMLMGUI package (version 2.1)66. A threshold level (–log10 (P) value = 4.5; P = 1/34,817) was set to calculate a significant linkage of SNPs with drought tolerance parameters in MLM whereas a critical LOD value of 3 was adopted for mrMLM (Xu et al., 2018; Misra et al., 2018). The SNPs identified by both models and in all environments was regarded as comparatively stable SNPs. The QTNs were named after the nomenclature of McCouch et al. (1997) with minor changes. Furthermore, the abbreviation was started with ‘q’, followed by trait name (GR, germination rate; WSL, whole seedling length; RL, root length) and the chromosome name, succeeded by number of QTN identified on same chromosome. Thus, qGR8-1 designated a locus on chromosome 08 related with GR in all environments. Haplotype block analysis of the comparatively stable SNPs was carried out in Haploview software version 4.2 with default parameters (Barrett et al.37.

Prediction of candidate genes

Possible candidate genes were obtained within major QTN region found on Chr.08 and Chr.13 by either single or multi-locus models in the soybean (William 82) gene models 1.0 (Grant et al., 2010). The functional annotations of reference genes were obtained from SoyBase and phytozome and candidate genes were selected manually based on their expression pattern in roots.

Conclusions

A total of 23 and 27 QTNs linked with GR, WSL and RL were detected using MLM and mrMLM approaches in all three environments (2018, 2019, and Combined). Among these QTNs, six QTNs, viz., qGR3-1, qGR8-1, qWSL8, qWSL13-1, qRL-8 and qRL15-4 were identified in all three environments via MLM and three of them (qGR8-1, qWSL13-1 and qRL-8) were also detected using mrMLM method. On the base of MLM and mrMLM comparative analysis, the most stable QTN8 was linked with RL and GR and it was confirmed through MLM and mrMLM methods in all environments. Eventually, the 1.0 Mb region around QTN8/ Gm08-12368145 and Gm13-35394159 ± 500 kb was utilized for candidate gene prediction. Based on the expression data, Glyma08g156800, Glyma08g16000, Glyma08g162700 and Glyma13g249600 were regarded as most possible candidate genes which regulate drought tolerance in soybean. Although, functional validation is further needed to define its role in drought tolerance. In short, this study broadens our vision to understand the genetic mechanisms regulating drought tolerance in soybean and it will be helpful in future to breed soybean cultivars with increased drought tolerance via gene cloning and marker assisted selection.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71357-8.

Acknowledgements

The authors would like to extend their sincere appreciation to the Researchers Supporting Project number (RSP2024R347), King Saud University, Riyadh, Saudi Arabia.

Author contributions

Author contributions: Conceptualization, T.Z, Muqadas Aleem, MHS, M.K.R; Data curation, Maida Aleem., M.K.R, and W.Y; Formal analysis, Muqadas Aleem, M.S.I, I.S, S.A; Investigation, Maida Aleem., W.Y., Z.A, and B.K.; Methodology, T.Z, BK, and Muqadas Aleem; Project administration, T.Z..; Supervision, T.Z.; Writing—original draft, Muqadas Aleem, NB and Maida A. Writing—review & editing, T.Z., M.A, Z.A, and N.B.

Funding

This research was funded by the National Key Research and Development Program, National Natural Science Foundation of China (Grant Nos. 31571691), the Core Technology Development for Breeding Program of Jiangsu Province (JBGS-2021-014), the Program of Collaborative Innovation Center for Modern Crop Production co-sponsored by Province and Ministry (CIC-MCP), and also by the Researchers Supporting Project number (RSP2024R347), King Saud University, Riyadh, Saudi Arabia.

Data availability

Manuscript and supplementary material has data.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Miransari M Abiotic and biotic stresses in soybean production 2016 Elsevier 53 77
Miransari, M. Abiotic and biotic stresses in soybean production 53–77 (Elsevier, 2016).
2. Payus C Impact of extreme drought climate on water security in North Borneo: Case study of Sabah Water 2020 12 1135 10.3390/w12041135
Payus, C. et al. Impact of extreme drought climate on water security in North Borneo: Case study of Sabah. Water 12, 1135 (2020).10.3390/w12041135
3. Wollenweber B Porter J Schellberg J Lack of interaction between extreme high-temperature events at vegetative and reproductive growth stages in wheat J. Agron. Crop Sci. 2003 189 142 150 10.1046/j.1439-037X.2003.00025.x
Wollenweber, B., Porter, J. & Schellberg, J. Lack of interaction between extreme high-temperature events at vegetative and reproductive growth stages in wheat. J. Agron. Crop Sci. 189, 142–150 (2003).10.1046/j.1439-037X.2003.00025.x
4. Frederick JR Camp CR Bauer PJ Drought-stress effects on branch and mainstem seed yield and yield components of determinate soybean Crop Sci. 2001 41 759 763 10.2135/cropsci2001.413759x
Frederick, J. R., Camp, C. R. & Bauer, P. J. Drought-stress effects on branch and mainstem seed yield and yield components of determinate soybean. Crop Sci. 41, 759–763 (2001).10.2135/cropsci2001.413759x
5. Sadeghipour O Abbasi S Soybean response to drought and seed inoculation World Appl. Sci. J. 2012 17 55 60
Sadeghipour, O. & Abbasi, S. Soybean response to drought and seed inoculation. World Appl. Sci. J. 17, 55–60 (2012).
6. Samarah N Mullen R Anderson I Soluble sugar contents, germination, and vigor of soybean seeds in response to drought stress J. New Seeds 2009 10 63 73 10.1080/15228860902786525
Samarah, N., Mullen, R. & Anderson, I. Soluble sugar contents, germination, and vigor of soybean seeds in response to drought stress. J. New Seeds 10, 63–73 (2009).10.1080/15228860902786525
7. Eck H Mathers A Musick J Plant water stress at various growth stages and growth and yield of soybeans Field Crops Res. 1987 17 1 16 10.1016/0378-4290(87)90077-3
Eck, H., Mathers, A. & Musick, J. Plant water stress at various growth stages and growth and yield of soybeans. Field Crops Res. 17, 1–16 (1987).10.1016/0378-4290(87)90077-3
8. Pennisi E The blue revolution, drop by drop, gene by gene Science 2008 320 171 173 10.1126/science.320.5873.171 18403686
Pennisi, E. The blue revolution, drop by drop, gene by gene. Science 320, 171–173 (2008).18403686 10.1126/science.320.5873.171
9. Zeid I Shedeed Z Response of alfalfa to putrescine treatment under drought stress Biol. Plant. 2006 50 635 10.1007/s10535-006-0099-9
Zeid, I. & Shedeed, Z. Response of alfalfa to putrescine treatment under drought stress. Biol. Plant. 50, 635 (2006).10.1007/s10535-006-0099-9
10. Nayyeripasand L Garoosi GA Ahmadikhah A Genome-wide association study (GWAS) to identify salt-tolerance QTLs carrying novel candidate genes in rice during early vegetative stage Rice 2021 14 1 21 10.1186/s12284-020-00433-0 33409869
Nayyeripasand, L., Garoosi, G. A. & Ahmadikhah, A. Genome-wide association study (GWAS) to identify salt-tolerance QTLs carrying novel candidate genes in rice during early vegetative stage. Rice 14, 1–21 (2021).33409869 10.1186/s12284-020-00433-0
11. Seetharam K Genomic regions associated with heat stress tolerance in tropical maize (Zea mays L.) Sci. Rep. 2021 11 1 11 10.1038/s41598-021-93061-7 33414495
Seetharam, K. et al. Genomic regions associated with heat stress tolerance in tropical maize (Zea mays L.). Sci. Rep. 11, 1–11 (2021).33414495 10.1038/s41598-021-93061-7
12. Ahmed AA Mohamed EA Hussein MY Sallam A Genomic regions associated with leaf wilting traits under drought stress in spring wheat at the seedling stage revealed by GWAS Environ. Exp. Bot. 2021 184 104393 10.1016/j.envexpbot.2021.104393
Ahmed, A. A., Mohamed, E. A., Hussein, M. Y. & Sallam, A. Genomic regions associated with leaf wilting traits under drought stress in spring wheat at the seedling stage revealed by GWAS. Environ. Exp. Bot. 184, 104393 (2021).10.1016/j.envexpbot.2021.104393
13. Zhao K Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa Nat. Commun. 2011 2 1 10 10.1038/ncomms1467
Zhao, K. et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2, 1–10 (2011).10.1038/ncomms1467
14. Kavar T Maras M Kidrič M Šuštar-Vozlič J Meglič V Identification of genes involved in the response of leaves of Phaseolus vulgaris to drought stress Mol. Breed. 2008 21 159 172 10.1007/s11032-007-9116-8
Kavar, T., Maras, M., Kidrič, M., Šuštar-Vozlič, J. & Meglič, V. Identification of genes involved in the response of leaves of Phaseolus vulgaris to drought stress. Mol. Breed. 21, 159–172 (2008).10.1007/s11032-007-9116-8
15. Lacerda JS Importance of zinc for arabica coffee and its effects on the chemical composition of raw grain and beverage quality Crop Sci. 2018 58 1360 1370 10.2135/cropsci2017.06.0373
Lacerda, J. S. et al. Importance of zinc for arabica coffee and its effects on the chemical composition of raw grain and beverage quality. Crop Sci. 58, 1360–1370 (2018).10.2135/cropsci2017.06.0373
16. Lehnert H Serfling A Friedt W Ordon F Genome-wide association studies reveal genomic regions associated with the response of wheat (Triticum aestivum L.) to mycorrhizae under drought stress conditions Front. Plant Sci. 2018 9 1728 10.3389/fpls.2018.01728 30568663
Lehnert, H., Serfling, A., Friedt, W. & Ordon, F. Genome-wide association studies reveal genomic regions associated with the response of wheat (Triticum aestivum L.) to mycorrhizae under drought stress conditions. Front. Plant Sci. 9, 1728 (2018).30568663 10.3389/fpls.2018.01728
17. Coumou D Robinson A Historic and future increase in the global land area affected by monthly heat extremes Environ. Res. Lett. 2013 8 034018 10.1088/1748-9326/8/3/034018
Coumou, D. & Robinson, A. Historic and future increase in the global land area affected by monthly heat extremes. Environ. Res. Lett. 8, 034018 (2013).10.1088/1748-9326/8/3/034018
18. Jiang Y Quantitative trait loci associated with drought tolerance in Brachypodium distachyon Front. Plant Sci. 2017 8 811 10.3389/fpls.2017.00811 28567049
Jiang, Y. et al. Quantitative trait loci associated with drought tolerance in Brachypodium distachyon. Front. Plant Sci. 8, 811 (2017).28567049 10.3389/fpls.2017.00811
19. Tester M Langridge P Breeding technologies to increase crop production in a changing world Science 2010 327 818 822 10.1126/science.1183700 20150489
Tester, M. & Langridge, P. Breeding technologies to increase crop production in a changing world. Science 327, 818–822 (2010).20150489 10.1126/science.1183700
20. Oluoch G QTL mapping for salt tolerance at seedling stage in the interspecific cross of Gossypium tomentosum with Gossypium hirsutum Euphytica 2016 209 223 235 10.1007/s10681-016-1674-6
Oluoch, G. et al. QTL mapping for salt tolerance at seedling stage in the interspecific cross of Gossypium tomentosum with Gossypium hirsutum. Euphytica 209, 223–235 (2016).10.1007/s10681-016-1674-6
21. Paterson AH DeVerna JW Lanini B Tanksley SD Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato Genetics 1990 124 735 742 10.1093/genetics/124.3.735 1968874
Paterson, A. H., DeVerna, J. W., Lanini, B. & Tanksley, S. D. Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato. Genetics 124, 735–742 (1990).1968874 10.1093/genetics/124.3.735
22. Baytar AA Genome-wide association mapping of yield components and drought tolerance-related traits in cotton Mol. Breed. 2018 38 1 16 10.1007/s11032-018-0831-0
Baytar, A. A. et al. Genome-wide association mapping of yield components and drought tolerance-related traits in cotton. Mol. Breed. 38, 1–16 (2018).10.1007/s11032-018-0831-0
23. Abou-Elwafa SF Shehzad T Genetic diversity, GWAS and prediction for drought and terminal heat stress tolerance in bread wheat (Triticum aestivum L.) Genet. Resour. Crop Evol. 2021 68 711 728 10.1007/s10722-020-01018-y
Abou-Elwafa, S. F. & Shehzad, T. Genetic diversity, GWAS and prediction for drought and terminal heat stress tolerance in bread wheat (Triticum aestivum L.). Genet. Resour. Crop Evol. 68, 711–728 (2021).10.1007/s10722-020-01018-y
24. Wang L QTL fine-mapping of soybean (Glycine max L.) leaf type associated traits in two RILs populations BMC Genom. 2019 20 1 15
Wang, L. et al. QTL fine-mapping of soybean (Glycine max L.) leaf type associated traits in two RILs populations. BMC Genom. 20, 1–15 (2019).
25. Cui Y Zhang F Zhou Y The application of multi-locus GWAS for the detection of salt-tolerance loci in rice Front. Plant Sci. 2018 9 1464 10.3389/fpls.2018.01464 30337936
Cui, Y., Zhang, F. & Zhou, Y. The application of multi-locus GWAS for the detection of salt-tolerance loci in rice. Front. Plant Sci. 9, 1464 (2018).30337936 10.3389/fpls.2018.01464
26. Boerma, H. R. & MianMR. In Proceedings of the World Soybean Research Conference VI. 4–7.
27. Yang W Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice Nat. Commun. 2014 5 1 9 10.1038/ncomms6087
Yang, W. et al. Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice. Nat. Commun. 5, 1–9 (2014).10.1038/ncomms6087
28. Korte A Farlow A The advantages and limitations of trait analysis with GWAS: A review Plant Methods 2013 9 1 9 10.1186/1746-4811-9-29 23286457
Korte, A. & Farlow, A. The advantages and limitations of trait analysis with GWAS: A review. Plant Methods 9, 1–9 (2013).23286457 10.1186/1746-4811-9-29
29. Sun H Genome-wide association mapping of stress-tolerance traits in cotton Crop J. 2019 7 77 88 10.1016/j.cj.2018.11.002
Sun, H. et al. Genome-wide association mapping of stress-tolerance traits in cotton. Crop J. 7, 77–88 (2019).10.1016/j.cj.2018.11.002
30. Pham A-T Genome-wide association of barley plant growth under drought stress using a nested association mapping population BMC Plant Biol. 2019 19 134 10.1186/s12870-019-1723-0 30971212
Pham, A.-T. et al. Genome-wide association of barley plant growth under drought stress using a nested association mapping population. BMC Plant Biol. 19, 134 (2019).30971212 10.1186/s12870-019-1723-0
31. Gowda M Relatedness severely impacts accuracy of marker-assisted selection for disease resistance in hybrid wheat Heredity 2014 112 552 561 10.1038/hdy.2013.139 24346498
Gowda, M. et al. Relatedness severely impacts accuracy of marker-assisted selection for disease resistance in hybrid wheat. Heredity 112, 552–561 (2014).24346498 10.1038/hdy.2013.139
32. Yu J A unified mixed-model method for association mapping that accounts for multiple levels of relatedness Nat. Genet. 2006 38 203 208 10.1038/ng1702 16380716
Yu, J. et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 38, 203–208 (2006).16380716 10.1038/ng1702
33. Zeng A Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines Mol. Breed. 2017 37 1 14 10.1007/s11032-017-0634-8 28127252
Zeng, A. et al. Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines. Mol. Breed. 37, 1–14 (2017).28127252 10.1007/s11032-017-0634-8
34. Zatybekov A GWAS of agronomic traits in soybean collection included in breeding pool in Kazakhstan BMC Plant Biol. 2017 17 1 8 10.1186/s12870-017-1125-0 28049439
Zatybekov, A. et al. GWAS of agronomic traits in soybean collection included in breeding pool in Kazakhstan. BMC Plant Biol. 17, 1–8 (2017).28049439 10.1186/s12870-017-1125-0
35. Liu Z Genome-wide association study of soybean seed germination under drought stress Mol. Genet. Genom. 2020 295 661 673 10.1007/s00438-020-01646-0
Liu, Z. et al. Genome-wide association study of soybean seed germination under drought stress. Mol. Genet. Genom. 295, 661–673 (2020).10.1007/s00438-020-01646-0
36. Baalbaki R Zurayk R Bleik M Talhouk S Germination and seedling development of drought tolerant and susceptible wheat under moisture stress Seed Sci. Technol. 1999 27 291 302
Baalbaki, R., Zurayk, R., Bleik, M. & Talhouk, S. Germination and seedling development of drought tolerant and susceptible wheat under moisture stress. Seed Sci. Technol. 27, 291–302 (1999).
37. Schneider R Genome-wide association mapping of partial resistance to Phytophthora sojae in soybean plant introductions from the Republic of Korea BMC Genom. 2016 17 1 14 10.1186/s12864-016-2918-5
Schneider, R. et al. Genome-wide association mapping of partial resistance to Phytophthora sojae in soybean plant introductions from the Republic of Korea. BMC Genom. 17, 1–14 (2016).10.1186/s12864-016-2918-5
38. Lam H-M Resequencing of 31 wild and cultivated soybean genomes identifies patterns of genetic diversity and selection Nat. Genet. 2010 42 1053 1059 10.1038/ng.715 21076406
Lam, H.-M. et al. Resequencing of 31 wild and cultivated soybean genomes identifies patterns of genetic diversity and selection. Nat. Genet. 42, 1053–1059 (2010).21076406 10.1038/ng.715
39. Kim S Recombination and linkage disequilibrium in Arabidopsis thaliana Nat. Genet. 2007 39 1151 1155 10.1038/ng2115 17676040
Kim, S. et al. Recombination and linkage disequilibrium in Arabidopsis thaliana. Nat. Genet. 39, 1151–1155 (2007).17676040 10.1038/ng2115
40. Mather KA The extent of linkage disequilibrium in rice (Oryza sativa L.) Genetics 2007 177 2223 2232 10.1534/genetics.107.079616 17947413
Mather, K. A. et al. The extent of linkage disequilibrium in rice (Oryza sativa L.). Genetics 177, 2223–2232 (2007).17947413 10.1534/genetics.107.079616
41. Wang M Yang W-M Du W-J Construction of a molecular marker linkage map and its use for quantitative trait locus (QTLs) underlying drought tolerance at germination stage in soybean Afr. J. Biotechnol. 2012 11 12830 12838
Wang, M., Yang, W.-M. & Du, W.-J. Construction of a molecular marker linkage map and its use for quantitative trait locus (QTLs) underlying drought tolerance at germination stage in soybean. Afr. J. Biotechnol. 11, 12830–12838 (2012).
42. Carter, T. In World soybean research conference IV; Proceedings. (Asociacion Agrentina de la Soja).
43. Kaya MD Okçu G Atak M Cıkılı Y Kolsarıcı Ö Seed treatments to overcome salt and drought stress during germination in sunflower (Helianthus annuus L.) Eur. J. Agron. 2006 24 291 295 10.1016/j.eja.2005.08.001
Kaya, M. D., Okçu, G., Atak, M., Cıkılı, Y. & Kolsarıcı, Ö. Seed treatments to overcome salt and drought stress during germination in sunflower (Helianthus annuus L.). Eur. J. Agron. 24, 291–295 (2006).10.1016/j.eja.2005.08.001
44. Manavalan LP Identification of novel QTL governing root architectural traits in an interspecific soybean population PLoS One 2015 10 e0120490 10.1371/journal.pone.0120490 25756528
Manavalan, L. P. et al. Identification of novel QTL governing root architectural traits in an interspecific soybean population. PLoS One 10, e0120490 (2015).25756528 10.1371/journal.pone.0120490
45. Abdel-Haleem H Lee G-J Boerma RH Identification of QTL for increased fibrous roots in soybean Theor. Appl. Genet. 2011 122 935 946 10.1007/s00122-010-1500-9 21165732
Abdel-Haleem, H., Lee, G.-J. & Boerma, R. H. Identification of QTL for increased fibrous roots in soybean. Theor. Appl. Genet. 122, 935–946 (2011).21165732 10.1007/s00122-010-1500-9
46. Riedelsheimer C Genome-wide association mapping of leaf metabolic profiles for dissecting complex traits in maize Proc. Natl. Acad. Sci. 2012 109 8872 8877 10.1073/pnas.1120813109 22615396
Riedelsheimer, C. et al. Genome-wide association mapping of leaf metabolic profiles for dissecting complex traits in maize. Proc. Natl. Acad. Sci. 109, 8872–8877 (2012).22615396 10.1073/pnas.1120813109
47. Stich B Linkage disequilibrium in European elite maize germplasm investigated with SSRs Theor. Appl. Genet. 2005 111 723 730 10.1007/s00122-005-2057-x 15997389
Stich, B. et al. Linkage disequilibrium in European elite maize germplasm investigated with SSRs. Theor. Appl. Genet. 111, 723–730 (2005).15997389 10.1007/s00122-005-2057-x
48. Zhu C Gore M Buckler ES Yu J Status and prospects of association mapping in plants Plant Genome 2008 1 1 10.3835/plantgenome2008.02.0089
Zhu, C., Gore, M., Buckler, E. S. & Yu, J. Status and prospects of association mapping in plants. Plant Genome 1, 1 (2008).10.3835/plantgenome2008.02.0089
49. Liu A Burke JM Patterns of nucleotide diversity in wild and cultivated sunflower Genetics 2006 173 321 330 10.1534/genetics.105.051110 16322511
Liu, A. & Burke, J. M. Patterns of nucleotide diversity in wild and cultivated sunflower. Genetics 173, 321–330 (2006).16322511 10.1534/genetics.105.051110
50. Yang C GmNAC8 acts as a positive regulator in soybean drought stress Plant Sci. 2020 293 110442 10.1016/j.plantsci.2020.110442 32081255
Yang, C. et al. GmNAC8 acts as a positive regulator in soybean drought stress. Plant Sci. 293, 110442 (2020).32081255 10.1016/j.plantsci.2020.110442
51. Wei W GmWRKY54 improves drought tolerance through activating genes in abscisic acid and Ca2+ signaling pathways in transgenic soybean Plant J. 2019 100 384 398 10.1111/tpj.14449 31271689
Wei, W. et al. GmWRKY54 improves drought tolerance through activating genes in abscisic acid and Ca2+ signaling pathways in transgenic soybean. Plant J. 100, 384–398 (2019).31271689 10.1111/tpj.14449
52. Luo X Expression of wild soybean WRKY20 in Arabidopsis enhances drought tolerance and regulates ABA signalling J. Exp. Bot. 2013 64 2155 2169 10.1093/jxb/ert073 23606412
Luo, X. et al. Expression of wild soybean WRKY20 in Arabidopsis enhances drought tolerance and regulates ABA signalling. J. Exp. Bot. 64, 2155–2169 (2013).23606412 10.1093/jxb/ert073
53. Li Y Expression of an Arabidopsis molybdenum cofactor sulphurase gene in soybean enhances drought tolerance and increases yield under field conditions Plant Biotechnol. J. 2013 11 747 758 10.1111/pbi.12066 23581509
Li, Y. et al. Expression of an Arabidopsis molybdenum cofactor sulphurase gene in soybean enhances drought tolerance and increases yield under field conditions. Plant Biotechnol. J. 11, 747–758 (2013).23581509 10.1111/pbi.12066
54. Wang N The soybean U-box gene GmPUB6 regulates drought tolerance in Arabidopsis thaliana Plant Physiol. Biochem. 2020 155 284 296 10.1016/j.plaphy.2020.07.016 32795910
Wang, N. et al. The soybean U-box gene GmPUB6 regulates drought tolerance in Arabidopsis thaliana. Plant Physiol. Biochem. 155, 284–296 (2020).32795910 10.1016/j.plaphy.2020.07.016
55. Brensha W Genetic analysis of root and shoot traits in the ‘Essex’ by ‘Forrest’ recombinant inbred line (RIL) population of soybean [Glycine max (L.) Merr.] Plant Genet. Genom. Biotechnol. 2012 1 1 9
Brensha, W. et al. Genetic analysis of root and shoot traits in the ‘Essex’ by ‘Forrest’ recombinant inbred line (RIL) population of soybean [Glycine max (L.) Merr.]. Plant Genet. Genom. Biotechnol. 1, 1–9 (2012).
56. Hetz W Hochholdinger F Schwall M Feix G Isolation and characterization of rtcs, a maize mutant deficient in the formation of nodal roots Plant J. 1996 10 845 857 10.1046/j.1365-313X.1996.10050845.x
Hetz, W., Hochholdinger, F., Schwall, M. & Feix, G. Isolation and characterization of rtcs, a maize mutant deficient in the formation of nodal roots. Plant J. 10, 845–857 (1996).10.1046/j.1365-313X.1996.10050845.x
57. Base, S. 9.3 Procedures Guide: Statistical Procedures (SAS Institute Inc., Cary, NC, 2011).
58. Nyquist WE Baker R Estimation of heritability and prediction of selection response in plant populations Crit. Rev. Plant Sci. 1991 10 235 322 10.1080/07352689109382313
Nyquist, W. E. & Baker, R. Estimation of heritability and prediction of selection response in plant populations. Crit. Rev. Plant Sci. 10, 235–322 (1991).10.1080/07352689109382313
59. Song Q Fingerprinting soybean germplasm and its utility in genomic research Genes Genomes Genet. 2015 5 1999 2006 10.1534/g3.115.019000
Song, Q. et al. Fingerprinting soybean germplasm and its utility in genomic research. Genes Genomes Genet. 5, 1999–2006 (2015).10.1534/g3.115.019000
60. Song Q Development and evaluation of SoySNP50K, a high-density genotyping array for soybean PLoS One 2013 8 e54985 10.1371/journal.pone.0054985 23372807
Song, Q. et al. Development and evaluation of SoySNP50K, a high-density genotyping array for soybean. PLoS One 8, e54985 (2013).23372807 10.1371/journal.pone.0054985
61. Purcell S PLINK: A tool set for whole-genome association and population-based linkage analyses Am. J. Hum. Genet. 2007 81 559 575 10.1086/519795 17701901
Purcell, S. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007).17701901 10.1086/519795
62. Pritchard JK Stephens M Donnelly P Inference of population structure using multilocus genotype data Genetics 2000 155 945 959 10.1093/genetics/155.2.945 10835412
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).10835412 10.1093/genetics/155.2.945
63. Evanno G Regnaut S Goudet J Detecting the number of clusters of individuals using the software structure: A simulation study Mol. Ecol. 2005 14 2611 2620 10.1111/j.1365-294X.2005.02553.x 15969739
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).15969739 10.1111/j.1365-294X.2005.02553.x
64. Bradbury PJ TASSEL: Software for association mapping of complex traits in diverse samples Bioinformatics 2007 23 2633 2635 10.1093/bioinformatics/btm308 17586829
Bradbury, P. J. et al. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635 (2007).17586829 10.1093/bioinformatics/btm308
65. He J An innovative procedure of genome-wide association analysis fits studies on germplasm population and plant breeding Theor. Appl. Genet. 2017 130 2327 2343 10.1007/s00122-017-2962-9 28828506
He, J. et al. An innovative procedure of genome-wide association analysis fits studies on germplasm population and plant breeding. Theor. Appl. Genet. 130, 2327–2343 (2017).28828506 10.1007/s00122-017-2962-9
66. Wang S-B Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology Sci. Rep. 2016 6 1 10 28442746
Wang, S.-B. et al. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. Sci. Rep. 6, 1–10 (2016).28442746
