Association mapping of agronomic traits on chromosome 2A of wheat

Association mapping is a method to test the association between molecular markers and quantitative trait loci (QTL) based on linkage disequilibrium (LD). In this study, the collection of 108 wheat germplasm accessions form China were evaluated for their plant heights, spike length, spikelets per spike, grains per spike, thousand kernel weight and spikelets density in 3 years at three locations. And they were genotyped with 85 SSR markers and 40 EST-SSR markers. The population structure was inferred on the basis of unlinked 48 SSR markers and 40 EST-SSR markers. The extent of LD on chromosome 2A was 2.3 cM. Association of 37 SSR loci on chromosomes 2A with six agronomic traits was analysed with a mixed linear model. A total of 14 SSR loci were significantly associated with agronomic traits. Some of the associated markers were located in the QTL region detected in previous linkage mapping analysis. Our results demonstrated that association mapping can enhance QTL information and achieves higher resolution with short LD extent.     

Polymorphisms in cinnamoyl CoA reductase (CCR) are associated with variation in microfibril angle in Eucalyptus spp

Linkage disequilibrium (LD) mapping using natural populations results in higher resolution of marker-trait associations compared to family-based quantitative trait locus (QTL) studies. Depending on the extent of LD, it is possible to identify alleles within candidate genes associated with a trait. Analysis of a natural mutant in Arabidopsis has shown that mutations in cinnamoyl CoA reductase (CCR), a key lignin gene, affect physical properties of the secondary cell wall such as stiffness and strength. Using this gene, we tested whether LD mapping could identify alleles associated with microfibril angle (MFA), a wood quality trait affecting stiffness and strength of wood. We identified 25 common single-nucleotide polymorphism (SNP) markers in the CCR gene in Eucalyptus nitens. Using single-marker and haplotype analyses in 290 trees from a E. nitens natural population, two haplotypes significantly associated with MFA were found. These results were confirmed in two full-sib families of E. nitens and Eucalyptus globulus. In an effort to understand the functional significance of the SNP markers, we sequenced the cDNA clones and identified an alternatively spliced variant from the significant haplotype region. This study demonstrates that LD mapping can be used to identify alleles associated with wood quality traits in natural populations of trees.     

Combined GWAS and QTL mapping revealed candidate genes and SNP network controlling recovery and tolerance traits associated with drought tolerance in seedling winter wheat

To date, very little research on drought tolerance has been conducted at the seedling stage in winter wheat. In this study, two types of traits, namely tolerance and recovery traits, associated with drought tolerance were scored in biparental mapping population (BPP) and association mapping population (A-set). The results of this study revealed no or weak significant correlation between the two types of traits. Based on GWAS and QTL mapping analyses, all QTLs associated with recovery traits were completely different from those associated with tolerance traits except one QTL in each population that was found to be associated with one tolerance trait and one recovery trait. The analysis of SNP and gene networks confirmed the results of combined GWAS and QTL mapping. One SNP marker located on the 2B chromosome (S2B_26494801) was found to be associated with recovery traits in both populations. The results of this study provided new information on understanding and improving drought tolerance in winter wheat.     

Molecular and comparative mapping of genes governing spike compactness from wild emmer wheat

The development and morphology of the wheat spike is important because the spike is where reproduction occurs and it holds the grains until harvest. Therefore, genes that influence spike morphology are of interest from both theoretical and practical stand points. When substituted for the native chromosome 2A in the tetraploid Langdon (LDN) durum wheat background, the Triticum turgidum ssp. dicoccoides chromosome 2A from accession IsraelA confers a short, compact spike with fewer spikelets per spike compared to LDN. Molecular mapping and quantitative trait loci (QTL) analysis of these traits in a homozygous recombinant population derived from LDN × the chromosome 2A substitution line (LDNIsA-2A) indicated that the number of spikelets per spike and spike length were controlled by linked, but different, loci on the long arm of 2A. A QTL explaining most of the variation for spike compactness coincided with the QTL for spike length. Comparative mapping indicated that the QTL for number of spikelets per spike overlapped with a previously mapped QTL for Fusarium head blight susceptibility. The genes governing spike length and compactness were not orthologous to either sog or C, genes known to confer compact spikes in diploid and hexaploid wheat, respectively. Mapping and sequence analysis indicated that the gene governing spike length and compactness derived from wild emmer could be an ortholog of the barley Cly1/Zeo gene, which research indicates is an AP2-like gene pleiotropically affecting cleistogamy, flowering time, and rachis internode length. This work provides researchers with knowledge of new genetic loci and associated markers that may be useful for manipulating spike morphology in durum wheat.     

Association mapping of drought tolerance-related traits in barley to complement a traditional biparental QTL mapping study

Key message

Association mapping of drought-related traits in barley was used to increase the density of existing QTL maps without recreating mapping populations.

Abstract

We used 109 spring barley genotypes exhibiting high or low drought tolerance to elucidate the associations between diversity array technology sequencing (DArTseq) and single nucleotide polymorphism (SNP) markers and various physiological parameters related to plant responses to drought conditions. The study was performed in controlled conditions (growth chambers), drought tolerance was phenotyped in the four-leaf seedlings. We identified 58 associations including 34 new markers (i.e., 16 DArTseq and 18 SNP markers). The results for three markers were consistent with the data obtained in an earlier traditional biparental QTL mapping study. The regions neighboring markers on linkage group 2H contained the highest number of significant marker–trait associations. Five markers related to the photosynthetic activity of photosystem II were detected on chromosome 4H. The lowest number of associations were observed for the sequences neighboring DArT markers on linkage group 6H. A chromosome 3H region related to water use efficiency and net photosynthesis rate in both biparental QTL, and association study, may be particularly valuable, as these parameters correspond to the ability of plants to remain highly productive under water deficit stress. Our findings confirm that association mapping can increase the density of existing QTL maps without recreating mapping populations.

     

Association mapping of salt tolerance in barley (Hordeum vulgare L.)

A spring barley collection of 192 genotypes from a wide geographical range was used to identify quantitative trait loci (QTLs) for salt tolerance traits by means of an association mapping approach using a thousand SNP marker set. Linkage disequilibrium (LD) decay was found with marker distances spanning 2–8 cM depending on the methods used to account for population structure and genetic relatedness between genotypes. The association panel showed large variation for traits that were highly heritable under salt stress, including biomass production, chlorophyll content, plant height, tiller number, leaf senescence and shoot Na⁺, shoot Cl^? and shoot, root Na⁺/K⁺ contents. The significant correlations between these traits and salt tolerance (defined as the biomass produced under salt stress relative to the biomass produced under control conditions) indicate that these traits contribute to (components of) salt tolerance. Association mapping was performed using several methods to account for population structure and minimize false-positive associations. This resulted in the identification of a number of genomic regions that strongly influenced salt tolerance and ion homeostasis, with a major QTL controlling salt tolerance on chromosome 6H, and a strong QTL for ion contents on chromosome 4H.     

Extent and genome-wide distribution of linkage disequilibrium in commercial maize germplasm

Association mapping is based on linkage disequilibrium (LD) resulting from historical recombinations and helps understanding the genetic basis of complex traits. Many factors affect LD and, therefore, it must be determined empirically in the germplasm under investigation to examine the prospects of successful genome-wide association mapping. The objectives of our study were to (1) examine the extent of LD with simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers in 1,537 commercial maize inbred lines belonging to four heterotic pools, (2) compare the LD patterns determined by these two marker types, (3) evaluate the number of SNP markers needed to perform genome-wide association analyses, and (4) investigate temporal trends of LD. Mean values of the squared correlation coefficient ( $ \bar{R} $ ) were almost identical for unlinked, linked, and adjacent SSR marker pairs. In contrast, $ \bar{R} $ values were lowest for the unlinked SNP loci and highest for the SNPs within amplicons. LD decay varied across the different heterotic pools and the individual chromosomes. The SSR markers employed in the present study are not adequate for association analysis, because of insufficient marker density for the germplasm evaluated. Based on the decay of LD in the various heterotic pools, we would need between 4,000 and 65,000 SNP markers to detect with a reasonable power associations with rather large quantitative trait loci (QTL). A much higher marker density is required to identify QTL with smaller effects. However, not only the total number of markers but also their distribution among and along the chromosomes are primordial for undertaking powerful association analyses.     

Linkage Disequilibrium and Genome-Wide Association Mapping in Tetraploid Wheat (Triticum turgidum L.)

Association mapping is a powerful tool for the identification of quantitative trait loci through the exploitation of the differential decay of linkage disequilibrium (LD) between marker loci and genes of interest in natural and domesticated populations. Using a sample of 230 tetraploid wheat lines (Triticum turgidum ssp), which included naked and hulled accessions, we analysed the pattern of LD considering 26 simple sequence repeats and 970 mostly mapped diversity array technology loci. In addition, to validate the potential for association mapping in durum wheat, we evaluated the same genotypes for plant height, heading date, protein content, and thousand-kernel weight. Molecular and phenotypic data were used to: (i) investigate the genetic and phenotypic diversity; (ii) study the dynamics of LD across the durum wheat genome, by investigating the patterns of LD decay; and (iii) test the potential of our panel to identify marker–trait associations through the analysis of four quantitative traits of major agronomic importance. Moreover, we compared and validated the association mapping results with outlier detection analysis based on population divergence. Overall, in tetraploid wheat, the pattern of LD is extremely population dependent and is related to the domestication and breeding history of durum wheat. Comparing our data with several other studies in wheat, we confirm the position of many major genes and quantitative trait loci for the traits considered. Finally, the analysis of the selection signature represents a very useful complement to validate marker–trait associations.     

Kernel weight per spike: what contributes to it at the individual QTL level?

Spike length (SL), spikelet number (SPN) per spike, kernel number per spike (KNPS), and thousand-kernel weight (TKW) have strong genetic associations with kernel weight per spike (KWPS) in wheat. To investigate their genetic relationships at the individual quantitative trait locus (QTL) level, both unconditional and conditional QTL mapping for KWPS with respect to SL, SPN, KNPS, and TKW were conducted. Two related F_8:9 recombinant inbred line populations, comprising 485 and 229 lines, respectively, were used. The trait phenotypic performances of each population were evaluated in four different environments. Unconditional QTL mapping analysis identified 22 putative additive QTL for KWPS, five of which were stable QTL, and only QKwps-WJ-1B.2 showed significant additive-by-environment interaction effects. In comparison with unconditional QTL mapping analysis, conditional QTL mapping analysis indicated that, at the QTL level, KNPS and TKW contributed more to KWPS than did SL and SPN. Any unconditional QTL for KWPS detected in this study was associated with at least one of its four related traits. The present study will provide assistance in the understanding of the genetic relationships between KWPS and its related traits.     

Genome-wide linkage mapping of yield-related traits in three Chinese bread wheat populations using high-density SNP markers

Key message

We identified 21 new and stable QTL, and 11 QTL clusters for yield-related traits in three bread wheat populations using the wheat 90 K SNP assay.

Abstract

Identification of quantitative trait loci (QTL) for yield-related traits and closely linked molecular markers is important in order to identify gene/QTL for marker-assisted selection (MAS) in wheat breeding. The objectives of the present study were to identify QTL for yield-related traits and dissect the relationships among different traits in three wheat recombinant inbred line (RIL) populations derived from crosses Doumai?×?Shi 4185 (D?×?S), Gaocheng 8901?×?Zhoumai 16 (G?×?Z) and Linmai 2?×?Zhong 892 (L?×?Z). Using the available high-density linkage maps previously constructed with the wheat 90 K iSelect single nucleotide polymorphism (SNP) array, 65, 46 and 53 QTL for 12 traits were identified in the three RIL populations, respectively. Among them, 34, 23 and 27 were likely to be new QTL. Eighteen common QTL were detected across two or three populations. Eleven QTL clusters harboring multiple QTL were detected in different populations, and the interval 15.5–32.3 cM around the Rht-B1 locus on chromosome 4BS harboring 20 QTL is an important region determining grain yield (GY). Thousand-kernel weight (TKW) is significantly affected by kernel width and plant height (PH), whereas flag leaf width can be used to select lines with large kernel number per spike. Eleven candidate genes were identified, including eight cloned genes for kernel, heading date (HD) and PH-related traits as well as predicted genes for TKW, spike length and HD. The closest SNP markers of stable QTL or QTL clusters can be used for MAS in wheat breeding using kompetitive allele-specific PCR or semi-thermal asymmetric reverse PCR assays for improvement of GY.

     

Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars

Association mapping is a method for detection of gene effects based on linkage disequilibrium (LD) that complements QTL analysis in the development of tools for molecular plant breeding. In this study, association mapping was performed on a selected sample of 95 cultivars of soft winter wheat. Population structure was estimated on the basis of 36 unlinked simple-sequence repeat (SSR) markers. The extent of LD was estimated on chromosomes 2D and part of 5A, relative to the LD observed among unlinked markers. Consistent LD on chromosome 2D was <1 cM, whereas in the centromeric region of 5A, LD extended for approximately 5 cM. Association of 62 SSR loci on chromosomes 2D, 5A, and 5B with kernel morphology and milling quality was analyzed through a mixed-effects model, where subpopulation was considered as a random factor and the marker tested was considered as a fixed factor. Permutations were used to adjust the threshold of significance for multiple testing within chromosomes. In agreement with previous QTL analysis, significant markers for kernel size were detected on the three chromosomes tested, and alleles potentially useful for selection were identified. Our results demonstrated that association mapping could complement and enhance previous QTL information for marker-assisted selection.     

Comparative LD mapping using single SNPs and haplotypes identifies QTL for plant height and biomass as secondary traits of drought tolerance in maize

Drought often delays developmental events so that plant height and above-ground biomass are reduced, resulting in yield loss due to inadequate photosynthate. In this study, plant height and biomass measured by the Normalized Difference Vegetation Index (NDVI) were used as criteria for drought tolerance. A total of 305 lines representing temperate, tropical and subtropical maize germplasm were genotyped using two single nucleotide polymorphism (SNP) chips each containing 1536 markers, from which 2052 informative SNPs and 386 haplotypes each constructed with two or more SNPs were used for linkage disequilibrium (LD) or association mapping. Single SNP- and haplotype-based LD mapping identified two significant SNPs and three haplotype loci [a total of four quantitative trait loci (QTL)] for plant height under well-watered and water-stressed conditions. For biomass, 32 SNPs and 12 haplotype loci (30 QTL) were identified using NDVIs measured at seven stages under the two water regimes. Some significant SNP and haplotype loci for NDVI were shared by different stages. Comparing significant loci identified by single SNP- and haplotype-based LD mapping, we found that six out of the 14 chromosomal regions defined by haplotype loci each included at least one significant SNP for the same trait. Significant SNP haplotype loci explained much higher phenotypic variation than individual SNPs. Moreover, we found that two significant SNPs (two QTL) and one haplotype locus were shared by plant height and NDVI. The results indicate the power of comparative LD mapping using single SNPs and SNP haplotypes with QTL shared by plant height and biomass as secondary traits for drought tolerance in maize.     

Genome-wide association mapping of yield and yield components of spring wheat under contrasting moisture regimes

Key message

A stable QTL that may be used in marker-assisted selection in wheat breeding programs was detected for yield, yield components and drought tolerance-related traits in spring wheat association mapping panel.

Abstract

Genome-wide association mapping has become a widespread method of quantitative trait locus (QTL) identification for many crop plants including wheat (Triticum aestivum L.). Its benefit over traditional bi-parental mapping approaches depends on the extent of linkage disequilibrium in the mapping population. The objectives of this study were to determine linkage disequilibrium decay rate and population structure in a spring wheat association mapping panel (n = 285–294) and to identify markers associated with yield and yield components, morphological, phenological, and drought tolerance-related traits. The study was conducted under fully irrigated and rain-fed conditions at Greeley, CO, USA and Melkassa, Ethiopia in 2010 and 2011 (five total environments). Genotypic data were generated using diversity array technology markers. Linkage disequilibrium decay rate extended over a longer genetic distance for the D genome (6.8 cM) than for the A and B genomes (1.7 and 2.0 cM, respectively). Seven subpopulations were identified with population structure analysis. A stable QTL was detected for grain yield on chromosome 2DS both under irrigated and rain-fed conditions. A multi-trait region significant for yield and yield components was found on chromosome 5B. Grain yield QTL on chromosome 1BS co-localized with harvest index QTL. Vegetation indices shared QTL with harvest index on chromosome 1AL and 5A. After validation in relevant genetic backgrounds and environments, QTL detected in this study for yield, yield components and drought tolerance-related traits may be used in marker-assisted selection in wheat breeding programs.     

Mapping QTL for grain yield, yield components, and spike features in a doubled haploid population of bread wheat

A doubled haploid (DH) population derived from a cross between the Japanese cultivar 'Fukuho-kumogi' and the Israeli wheat line 'Oligoculm' was used to map genome regions involved in the expression of grain yield, yield components, and spike features in wheat (Triticum aestivum L). A total of 371 markers (RAPD, SSR, RFLP, AFLP, and two morphological traits) were used to construct the linkage map that covered 4190 cM of wheat genome including 28 linkage groups. The results of composite interval mapping for all studied traits showed that some of the quantitative trait loci (QTL) were stable over experiments conducted in 2004 and 2005. The major QTL located in the Hair-Xpsp2999 interval on chromosome 1A controlled the expression of grains/spike (R(2) = 12.9% in 2004 and 22.4% in 2005), grain weight/spike (R(2) = 21.4% in 2004 and 15.8% in 2005), and spike number (R(2) = 15.6% in 2004 and 5.4% in 2005). The QTL for grain yield located on chromosomes 6A, 6B, and 6D totally accounted for 27.2% and 31.7% of total variation in this trait in 2004 and 2005, respectively. Alleles inherited from 'Oligoculm' increased the length of spikes and had decreasing effects on spike number. According to the data obtained in 2005, locus Xgwm261 was associated with a highly significant spike length QTL (R(2) = 42.33%) and also the major QTL for spikelet compactness (R(2) = 26.1%).     

Identification and validation of a major cadmium accumulation locus and closely associated SNP markers in North Dakota durum wheat cultivars

Durum wheat has the tendency of accumulating more cadmium (Cd), a biotoxic heavy metal, in seeds than other commonly grown cereals, thus posing a serious food safety/public health concern. This could have serious negative impact on the national pasta industry and the international export market of durum wheat. The phenotyping for selecting low Cd lines is expensive and time consuming. The use of markers could be a more sustainable approach for selecting lines with low Cd levels. Here, a RIL population developed from a cross between Grenora (high Cd) × Haurani (low Cd) and two association mapping panels consisting of advanced breeding lines from the North Dakota durum wheat breeding program were used to identify QTL and associated markers for Cd. A major QTL, with Haurani contributing low Cd uptake allele and explaining 54.3 % phenotypic variation, was detected on chromosome 5BL. Association mapping using 2010 collection validated the results of linkage mapping and identified major QTL on 5BL. The 2009 collection, showed the presence of a major QTL on chromosome 2B. The SNP marker associated with major QTL on 5BL was converted to user friendly KASPar assay. The major QTL and associated KASPar marker were further validated using another RIL population developed from a cross of Strongfield (low Cd) and Alkabo (high Cd). The development of suitable marker assay, associated with major Cd uptake QTL, would help the selection for low Cd accumulating lines, minimizing the costly phenotypic evaluation for this important trait.     

QTL mapping of potato chip color and tuber traits within an autotetraploid family

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid crop species, and this creates challenges for traditional line development and molecular breeding. Recent availability of a single-nucleotide polymorphism (SNP) array with 8303 features and software packages for linkage and association mapping in autotetraploid species present new opportunities for the identification of genomic regions that contribute to high-value traits in cultivated potato. A biparental tetraploid potato population was evaluated across three field seasons and storage trials in order to identify quantitative trait loci (QTL) for multiple tuber traits including fried chip color after 5.5–7.2 °C storage. Tetra-allelic dosage information was used to construct a genetic linkage map that covered 1041 cM and contained 2095 SNP markers with a median marker interval of 0.4 cM. A total of 41 QTL were identified for flower color, tuber yield, tuber number per plant, tuber weight, tuber size, and chip color after various storage regimes. Moderate effect QTL for chip color at 3 months were identified that co-localized with candidate genes vacuolar invertase (VInv), invertase inhibitor (INH2), and apoplastic invertase (Inv _ap -b). A separate QTL for chip color after 6 months of storage was identified in the short arm of chromosome 2, and this locus may contribute to variation in senescent sweetening resistance. QTL for tuber weight, length, and width co-localized with a known QTL for plant maturity on chromosome 5. Genome-wide association mapping using a polyploid model detected the tuber size QTL and identified a number of candidate SNPs, but was unable to detect markers significantly associated with chip color.     

Mixed model association scans of multi-environmental trial data reveal major loci controlling yield and yield related traits in Hordeum vulgare in Mediterranean environments

An association panel consisting of 185 accessions representative of the barley germplasm cultivated in the Mediterranean basin was used to localise quantitative trait loci (QTL) controlling grain yield and yield related traits. The germplasm set was genotyped with 1,536 SNP markers and tested for associations with phenotypic data gathered over 2?years for a total of 24?year?×?location combinations under a broad range of environmental conditions. Analysis of multi-environmental trial (MET) data by fitting a mixed model with kinship estimates detected from two to seven QTL for the major components of yield including 1000 kernel weight, grains per spike and spikes per m², as well as heading date, harvest index and plant height. Several of the associations involved SNPs tightly linked to known major genes determining spike morphology in barley (vrs1 and int-c). Similarly, the largest QTL for heading date co-locates with SNPs linked with eam6, a major locus for heading date in barley for autumn sown conditions. Co-localization of several QTL related to yield components traits suggest that major developmental loci may be linked to most of the associations. This study highlights the potential of association genetics to identify genetic variants controlling complex traits.     

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

BackgroundIncreasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.ResultsThe populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.ConclusionsThis study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.     

Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

Genome-Wide Association Studies Identifies Seven Major Regions Responsible for Iron Deficiency Chlorosis in Soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.