Genome-wide association study identifies novel loci and candidate genes for rust resistance in wheat (Triticum aestivum L.).

BACKGROUND: Wheat rusts are important biotic stresses, development of rust resistant cultivars through molecular approaches is both economical and sustainable. Extensive phenotyping of large mapping populations under diverse production conditions and high-density genotyping would be the ideal strategy to identify major genomic regions for rust resistance in wheat. The genome-wide association study (GWAS) population of 280 genotypes was genotyped using a 35 K Axiom single nucleotide polymorphism (SNP) array and phenotyped at eight, 10, and, 10 environments, respectively for stem/black rust (SR), stripe/yellow rust (YR), and leaf/brown rust (LR). RESULTS: Forty-one Bonferroni corrected marker-trait associations (MTAs) were identified, including 17 for SR and 24 for YR. Ten stable MTAs and their best combinations were also identified. For YR, AX-94990952 on 1A + AX-95203560 on 4A + AX-94723806 on 3D + AX-95172478 on 1A showed the best combination with an average co-efficient of infection (ACI) score of 1.36. Similarly, for SR, AX-94883961 on 7B + AX-94843704 on 1B and AX-94883961 on 7B + AX-94580041 on 3D + AX-94843704 on 1B showed the best combination with an ACI score of around 9.0. The genotype PBW827 have the best MTA combinations for both YR and SR resistance. In silico study identifies key prospective candidate genes that are located within MTA regions. Further, the expression analysis revealed that 18 transcripts were upregulated to the tune of more than 1.5 folds including 19.36 folds (TraesCS3D02G519600) and 7.23 folds (TraesCS2D02G038900) under stress conditions compared to the control conditions. Furthermore, highly expressed genes in silico under stress conditions were analyzed to find out the potential links to the rust phenotype, and all four genes were found to be associated with the rust phenotype. CONCLUSION: The identified novel MTAs, particularly stable and highly expressed MTAs are valuable for further validation and subsequent application in wheat rust resistance breeding. The genotypes with favorable MTA combinations can be used as prospective donors to develop elite cultivars with YR and SR resistance.

Genome-wide association studies for phenological and agronomic traits in mungbean (Vigna radiata L. Wilczek).

Mungbean (Vigna radiata L. Wilczek) is one of the important warm-season food legumes, contributing substantially to nutritional security and environmental sustainability. The genetic complexity of yield-associated agronomic traits in mungbean is not well understood. To dissect the genetic basis of phenological and agronomic traits, we evaluated 153 diverse mungbean genotypes for two phenological (days to heading and days to maturity) and eight agronomic traits (leaf nitrogen status using SPAD, plant height, number of primary branches, pod length, number of pods per plant, seeds per pod, 100-seed weight, and yield per plant) under two environmental conditions. A wide array of phenotypic variability was apparent among the studied genotypes for all the studied traits. The broad sense of heritability of traits ranged from 0.31 to 0.95 and 0.21 to 0.94 at the Delhi and Ludhiana locations, respectively. A total of 55,634 genome-wide single nucleotide polymorphisms (SNPs) were obtained by the genotyping-by-sequencing method, of which 15,926 SNPs were retained for genome-wide association studies (GWAS). GWAS with Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK) model identified 50 SNPs significantly associated with phenological and agronomic traits. In total, 12 SNPs were found to be significantly associated with phenological traits across environments, explaining 7%-18.5% of phenotypic variability, and 38 SNPs were significantly associated with agronomic traits, explaining 4.7%-27.6% of the phenotypic variability. The maximum number of SNPs (15) were located on chromosome 1, followed by seven SNPs each on chromosomes 2 and 8. The BLAST search identified 19 putative candidate genes that were involved in light signaling, nitrogen responses, phosphorus (P) transport and remobilization, photosynthesis, respiration, metabolic pathways, and regulating growth and development. Digital expression analysis of 19 genes revealed significantly higher expression of 12 genes, viz. VRADI01G08170, VRADI11G09170, VRADI02G00450, VRADI01G00700, VRADI07G14240, VRADI03G06030, VRADI02G14230, VRADI08G01540, VRADI09G02590, VRADI08G00110, VRADI02G14240, and VRADI02G00430 in the roots, cotyledons, seeds, leaves, shoot apical meristems, and flowers. The identified SNPs and putative candidate genes provide valuable genetic information for fostering genomic studies and marker-assisted breeding programs that improve yield and agronomic traits in mungbean.

Unravelling the Genetic Basis of Moisture Deficit Stress Tolerance in Wheat for Seedling Vigour-Related Traits and Root Traits Using Genome-Wide Association Study.

A key abiotic stress that negatively affects seed germination, plant development, and crop yield is moisture deficit stress. Achieving higher vigour and uniform germination under stress conditions is essential for crop establishment and productivity and to enhance the yield. Hence, revealing wheat's capacity to withstand moisture deficit stress during seed germination and early growth stages is fundamental in improving its overall performance. However, the genetic regulation of moisture deficit stress tolerance during the seed germination phase remains largely unexplored. In this study, a total of 193 wheat genotypes were subjected to simulated moisture deficit stress using PEG-6000 (-0.4 MPa) during the seed germination stage. The induced moisture deficit stress significantly reduced various seedling-vigour-related traits. The genetic regions linked to these traits were found using a genome-wide association study (GWAS). The analysis identified 235 MTAs with a significance -log10(p) value of >4. After applying the Bonferroni correction, the study identified 47 unique single nucleotide polymorphisms (SNPs) that are linked to candidate genes important for the trait of interest. The current study emphasises the effectiveness of genome-wide association studies (GWAS) in identifying promising candidate genes, improving wheat seedling vigour and root traits, and offering essential information for the development of wheat cultivars tolerant to moisture deficit stress.

Marker-assisted selection for transfer of QTLs to a promising line for drought tolerance in wheat (Triticum aestivum L.).

Wheat crop is subjected to various biotic and abiotic stresses, which affect crop productivity and yield. Among various abiotic stresses, drought stress is a major problem considering the current global climate change scenario. A high-yielding wheat variety, HD3086, has been released for commercial cultivation under timely sown irrigated conditions for the North Western Plain Zone (NWPZ) and North Eastern Plain Zone NEPZ of India. Presently, HD3086 is one of the highest breeder seed indented wheat varieties and has a stable yield over the years. However, under moisture deficit conditions, its potential yield cannot be achieved. The present study was undertaken to transfer drought-tolerant QTLs in the background of the variety HD3086 using marker-assisted backcross breeding. QTLs governing Biomass (BIO), Canopy Temperature (CT), Thousand Kernel Weight (TKW), Normalized Difference Vegetation Index (NDVI), and Yield (YLD) were transferred to improve performance under moisture deficit conditions. In BC1F1, BC2F1, and BC2F2 generations, the foreground selection was carried out to identify the plants with positive QTLs conferring drought tolerance and linked to traits NDVI, CT, TKW, and yield. The positive homozygous lines for targeted QTLs were advanced from BC2F2 to BC2F4 via the pedigree-based phenotypic selection method. Background analysis was carried out in BC2F5 and obtained 78-91% recovery of the recurrent parent genome in the improved lines. Furthermore, the advanced lines were evaluated for 2 years under drought stress to assess improvement in MABB-derived lines. Increased GWPS, TKW, and NDVI and reduced CT was observed in improved lines. Seven improved lines were identified with significantly higher yields in comparison to HD3086 under stress conditions.

Genome-Wide Association Study for Grain Protein, Thousand Kernel Weight, and Normalized Difference Vegetation Index in Bread Wheat (Triticum aestivum L.).

Genomic regions governing grain protein content (GPC), 1000 kernel weight (TKW), and normalized difference vegetation index (NDVI) were studied in a set of 280 bread wheat genotypes. The genome-wide association (GWAS) panel was genotyped using a 35K Axiom array and phenotyped in three environments. A total of 26 marker-trait associations (MTAs) were detected on 18 chromosomes covering the A, B, and D subgenomes of bread wheat. The GPC showed the maximum MTAs (16), followed by NDVI (6), and TKW (4). A maximum of 10 MTAs was located on the B subgenome, whereas, 8 MTAs each were mapped on the A and D subgenomes. In silico analysis suggest that the SNPs were located on important putative candidate genes such as NAC domain superfamily, zinc finger RING-H2-type, aspartic peptidase domain, folylpolyglutamate synthase, serine/threonine-protein kinase LRK10, pentatricopeptide repeat, protein kinase-like domain superfamily, cytochrome P450, and expansin. These candidate genes were found to have different roles including regulation of stress tolerance, nutrient remobilization, protein accumulation, nitrogen utilization, photosynthesis, grain filling, mitochondrial function, and kernel development. The effects of newly identified MTAs will be validated in different genetic backgrounds for further utilization in marker-aided breeding.

Mapping of the QTLs governing grain micronutrients and thousand kernel weight in wheat (Triticum aestivum L.) using high density SNP markers.

Biofortification is gaining importance globally to improve human nutrition through enhancing the micronutrient content, such as vitamin A, iron, and zinc, in staple food crops. The present study aims to identify the chromosomal regions governing the grain iron concentration (GFeC), grain zinc concentration (GZnC), and thousand kernel weight (TKW) using recombinant inbred lines (RILs) in wheat, developed from a cross between HD3086 and HI1500. The experiment was conducted in four different production conditions at Delhi viz., control, drought, heat, and combined heat and drought stress and at Indore under drought stress. Grain iron and zinc content increased under heat and combined stress conditions, while thousand kernel weight decreased. Medium to high heritability with a moderate correlation between grain iron and zinc was observed. Out of 4,106 polymorphic markers between the parents, 3,407 SNP markers were used for linkage map construction which spanned over a length of 14791.18 cm. QTL analysis identified a total of 32 chromosomal regions governing the traits under study, which includes 9, 11, and 12 QTLs for GFeC, GZnC, and TKW, respectively. A QTL hotspot was identified on chromosome 4B which is associated with grain iron, grain zinc, and thousand kernel weight explaining the phenotypic variance of 29.28, 10.98, and 17.53%, respectively. Similarly, common loci were identified on chromosomes 4B and 4D for grain iron, zinc, and thousand kernel weight. In silico analysis of these chromosomal regions identified putative candidate genes that code for proteins such as Inositol 1,3,4-trisphosphate 5/6-kinase, P-loop containing nucleoside triphosphate hydrolase, Pleckstrin homology (PH) domains, Serine-threonine/tyrosine-protein kinase and F-box-like domain superfamily proteins which play role in many important biochemical or physiological process. The identified markers linked to QTLs can be used in MAS once successfully validated.

Molecular Mapping of Biofortification Traits in Bread Wheat (Triticum aestivum L.) Using a High-Density SNP Based Linkage Map.

A set of 188 recombinant inbred lines (RILs) derived from a cross between a high-yielding Indian bread wheat cultivar HD2932 and a synthetic hexaploid wheat (SHW) Synthetic 46 derived from tetraploid Triticum turgidum (AA, BB 2n = 28) and diploid Triticum tauschii (DD, 2n = 14) was used to identify novel genomic regions associated in the expression of grain iron concentration (GFeC), grain zinc concentration (GZnC), grain protein content (GPC) and thousand kernel weight (TKW). The RIL population was genotyped using SNPs from 35K Axiom® Wheat Breeder's Array and 34 SSRs and phenotyped in two environments. A total of nine QTLs including five for GPC (QGpc.iari_1B, QGpc.iari_4A, QGpc.iari_4B, QGpc.iari_5D, and QGpc.iari_6B), two for GFeC (QGfec.iari_5B and QGfec.iari_6B), and one each for GZnC (QGznc.iari_7A) and TKW (QTkw.iari_4B) were identified. A total of two stable and co-localized QTLs (QGpc.iari_4B and QTkw.iari_4B) were identified on the 4B chromosome between the flanking region of Xgwm149-AX-94559916. In silico analysis revealed that the key putative candidate genes such as P-loop containing nucleoside triphosphatehydrolase, Nodulin-like protein, NAC domain, Purine permease, Zinc-binding ribosomal protein, Cytochrome P450, Protein phosphatase 2A, Zinc finger CCCH-type, and Kinesin motor domain were located within the identified QTL regions and these putative genes are involved in the regulation of iron homeostasis, zinc transportation, Fe, Zn, and protein remobilization to the developing grain, regulation of grain size and shape, and increased nitrogen use efficiency. The identified novel QTLs, particularly stable and co-localized QTLs are useful for subsequent use in marker-assisted selection (MAS).

QTL mapping: insights into genomic regions governing component traits of yield under combined heat and drought stress in wheat.

Drought and heat frequently co-occur during crop growth leading to devastating yield loss. The knowledge of the genetic loci governing component traits of yield under combined drought and heat stress is essential for enhancing the climate resilience. The present study employed a mapping population of 180 recombinant inbred lines (RILs) derived from a cross between GW322 and KAUZ to identify quantitative trait loci (QTLs) governing the component traits of yield under heat and combined stress conditions. Phenotypic evaluation was conducted across two consecutive crop seasons (2021-2022 and 2022-2023) under late sown irrigation (LSIR) and late sown restricted irrigation (LSRI) conditions at the Indian Council of Agricultural Research Institute-Indian Agricultural Research Institute (ICAR-IARI), New Delhi. Various physiological and agronomic traits of importance were measured. Genotyping was carried out with 35K SNP Axiom breeder's genotyping array. The linkage map spanned a length of 6769.45 cM, ranging from 2.28 cM/marker in 1A to 14.21 cM/marker in 5D. A total of 35 QTLs were identified across 14 chromosomes with 6B containing the highest (seven) number of QTLs. Out of 35 QTLs, 16 were major QTLs explaining the phenotypic variance greater than 10%. The study identified eight stable QTLs along with two hotspots on chromosomes 6B and 5B. Five QTLs associated with traits thousand-grain weight (TGW), normalized difference vegetation index (NDVI), and plant height (PH) were successfully validated. Candidate genes encoding antioxidant enzymes, transcription factors, and growth-related proteins were identified in the QTL regions. In silico expression analysis highlighted higher expression of transcripts TraesCS2D02G021000.1, TraesCS2D02G031000, TraesCS6A02G247900, and TraesCS6B02G421700 under stress conditions. These findings contribute to a deeper understanding of the genetic architecture underlying combined heat and drought tolerance in wheat, providing valuable insights for wheat improvement strategies under changing climatic conditions.

Marker assisted improvement for leaf rust and moisture deficit stress tolerance in wheat variety HD3086.

There is a significant yield reduction in the wheat crop as a result of different biotic and abiotic stresses, and changing climate, among them moisture deficit stress and leaf rust are the major ones affecting wheat worldwide. HD3086 is a high-yielding wheat variety that has been released for commercial cultivation under timely sown irrigated conditions in the Indo-Gangetic plains of India. Variety HD3086 provides a good, stable yield, and it is the choice of millions of farmers in India. It becomes susceptible to the most prevalent pathotypes 77-5 and 77-9 of Puccinia triticina (causing leaf rust) in the production environment and its potential yield cannot be realized under moisture deficit stress. The present study demonstrates the use of a marker-assisted back cross breeding approach to the successful transfer of leaf rust resistance gene Lr24 and QTLs linked to moisture deficit stress tolerance in the background of HD3086. The genotype HI1500 was used as a donor parent that possesses leaf rust-resistant gene Lr24, which confers resistance against the major pathotypes found in the production environment. It possesses inbuilt tolerance under abiotic stresses with superior quality traits. Foreground selection for gene Lr24 and moisture deficit stress tolerance QTLs linked to Canopy temperature (CT), Normal Differential Vegetation Index (NDVI) and Thousand Kernel Weight (TKW) in different generations of the backcrossing and selection. In BC2F2, foreground selection was carried out to identify homozygous lines based on the linked markers and were advanced following pedigree based phenotypic selection. The selected lines were evaluated against P. triticina pathotypes 77-5 and 77-9 under controlled conditions. Recurrent parent recovery of the selected lines ranged from 78-94%. The identified lines were evaluated for their tolerance to moisture stress under field conditions and their resistance to rust under artificial epiphytotic conditions for two years. In BC2F5 generation, eight positive lines for marker alleles were selected which showed resistance to leaf rust and recorded an improvement in component traits of moisture deficit stress tolerance such as CT, NDVI, TKW and yield compared to the recurrent parent HD3086. The derived line is named HD3471 and is nominated for national trials for testing and further release for commercial cultivation.

Identification of genomic regions of wheat associated with grain Fe and Zn content under drought and heat stress using genome-wide association study.

Wheat is the staple food crop of global importance for its grain nutrient quality. Grain iron and zinc content of the wheat grain is an important quantitatively inherited trait that is influenced by the environmental factors such as drought and heat stress. Phenotypic evaluation of 295 advanced breeding lines from the wheat stress breeding program of IARI was carried out under timely sown irrigated (IR), restricted irrigated, and late-sown conditions at New Delhi during the cropping season of 2020-21, and grain iron (GFeC) and zinc (GZnC) contents were estimated from both control and treatments. A statistically significant increase in GFeC and GZnC was observed under stress conditions compared to that of the control. Genotyping was carried out with the SNPs from the 35K Axiom Breeder's array, and marker-trait association was identified by GWAS analysis. Of the 23 MTAs identified, seven were linked with GFeC and sixteen were linked with GZnC. In silico analysis revealed a few important transcripts involved in various plant metabolism, growth, and development activities such as auxin response factor, root UVB sensitive proteins, potassium transporter, glycosyl transferase, COBRA, and F-box-like domain. The identified MTAs can be used for molecular breeding after validation and also for rapid development of micronutrient-rich varieties of wheat to mitigate hidden hunger.

Advances in integrated genomic selection for rapid genetic gain in crop improvement: a review.

MAIN CONCLUSION: Genomic selection and its importance in crop breeding. Integration of GS with new breeding tools and developing SOP for GS to achieve maximum genetic gain with low cost and time. The success of conventional breeding approaches is not sufficient to meet the demand of a growing population for nutritious food and other plant-based products. Whereas, marker assisted selection (MAS) is not efficient in capturing all the favorable alleles responsible for economic traits in the process of crop improvement. Genomic selection (GS) developed in livestock breeding and then adapted to plant breeding promised to overcome the drawbacks of MAS and significantly improve complicated traits controlled by gene/QTL with small effects. Large-scale deployment of GS in important crops, as well as simulation studies in a variety of contexts, addressed G × E interaction effects and non-additive effects, as well as lowering breeding costs and time. The current study provides a complete overview of genomic selection, its process, and importance in modern plant breeding, along with insights into its application. GS has been implemented in the improvement of complex traits including tolerance to biotic and abiotic stresses. Furthermore, this review hypothesises that using GS in conjunction with other crop improvement platforms accelerates the breeding process to increase genetic gain. The objective of this review is to highlight the development of an appropriate GS model, the global open source network for GS, and trans-disciplinary approaches for effective accelerated crop improvement. The current study focused on the application of data science, including machine learning and deep learning tools, to enhance the accuracy of prediction models. Present study emphasizes on developing plant breeding strategies centered on GS combined with routine conventional breeding principles by developing GS-SOP to achieve enhanced genetic gain.

Genome-wide association mapping for component traits of drought and heat tolerance in wheat.

Identification of marker trait association is a prerequisite for marker-assisted breeding. To find markers linked with traits under heat and drought stress in bread wheat (Triticum aestivum L.), we performed a genome-wide association study (GWAS). GWAS mapping panel used in this study consists of advanced breeding lines from the IARI stress breeding programme produced by pairwise and complex crosses. Phenotyping was done at multi locations namely New Delhi, Karnal, Indore, Jharkhand and Pune with augmented-RCBD design under different moisture and heat stress regimes, namely timely sown irrigated (IR), timely sown restricted irrigated (RI) and late sown (LS) conditions. Yield and its component traits, viz., Days to Heading (DH), Days to Maturity (DM), Normalized Difference Vegetation Index (NDVI), Chlorophyll Content (SPAD), Canopy temperature (CT), Plant Height (PH), Thousand grain weight (TGW), Grain weight per spike (GWPS), Plot Yield (PLTY) and Biomass (BMS) were phenotyped. Analysis of variance and descriptive statistics revealed significant differences among the studied traits. Genotyping was done using the 35k SNP Wheat Breeder's Genotyping Array. Population structure and diversity analysis using filtered 10,546 markers revealed two subpopulations with sufficient diversity. A large whole genome LD block size of 7.15 MB was obtained at half LD decay value. Genome-wide association search identified 57 unique markers associated with various traits across the locations. Twenty-three markers were identified to be stable, among them nine pleiotropic markers were also identified. In silico search of the identified markers against the IWGSC ref genome revealed the presence of a majority of the SNPs at or near the gene coding region. These SNPs can be used for marker-assisted transfer of genes/QTLs after validation to develop climate-resilient cultivars.

Genome-wide association study for grain yield and component traits in bread wheat (Triticum aestivum L.).

Genomic regions governing days to heading (DH), grain filling duration (GFD), grain number per spike (GNPS), grain weight per spike (GWPS), plant height (PH), and grain yield (GY) were investigated in a set of 280 diverse bread wheat genotypes. The genome-wide association studies (GWAS) panel was genotyped using a 35K Axiom Array and phenotyped in five environments. The GWAS analysis showed a total of 27 Bonferroni-corrected marker-trait associations (MTAs) on 15 chromosomes representing all three wheat subgenomes. The GFD showed the highest MTAs (8), followed by GWPS (7), GY (4), GNPS (3), PH (3), and DH (2). Furthermore, 20 MTAs were identified with more than 10% phenotypic variation. A total of five stable MTAs (AX-95024590, AX-94425015, AX-95210025 AX-94539354, and AX-94978133) were identified in more than one environment and associated with the expression of DH, GFD, GNPS, and GY. Similarly, two novel pleiotropic genomic regions with associated MTAs i.e. AX-94978133 (4D) and AX-94539354 (6A) harboring co-localized QTLs governing two or more traits were also identified. In silico analysis revealed that the SNPs were located on important putative candidate genes such as F-box-like domain superfamily, Lateral organ boundaries, LOB, Thioredoxin-like superfamily Glutathione S-transferase, RNA-binding domain superfamily, UDP-glycosyltransferase family, Serine/threonine-protein kinase, Expansin, Patatin, Exocyst complex component Exo70, DUF1618 domain, Protein kinase domain involved in the regulation of grain size, grain number, growth and development, grain filling duration, and abiotic stress tolerance. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection (MAS).

Genetic dissection of grain iron and zinc, and thousand kernel weight in wheat (Triticum aestivum L.) using genome-wide association study.

Genetic biofortification is recognized as a cost-effective and sustainable strategy to reduce micronutrient malnutrition. Genomic regions governing grain iron concentration (GFeC), grain zinc concentration (GZnC), and thousand kernel weight (TKW) were investigated in a set of 280 diverse bread wheat genotypes. The genome-wide association (GWAS) panel was genotyped using 35 K Axiom Array and phenotyped in five environments. The GWAS analysis showed a total of 17 Bonferroni-corrected marker-trait associations (MTAs) in nine chromosomes representing all the three wheat subgenomes. The TKW showed the highest MTAs (7), followed by GZnC (5) and GFeC (5). Furthermore, 14 MTAs were identified with more than 10% phenotypic variation. One stable MTA i.e. AX-95025823 was identified for TKW in both E4 and E5 environments along with pooled data, which is located at 68.9 Mb on 6A chromosome. In silico analysis revealed that the SNPs were located on important putative candidate genes such as Multi antimicrobial extrusion protein, F-box domain, Late embryogenesis abundant protein, LEA-18, Leucine-rich repeat domain superfamily, and C3H4 type zinc finger protein, involved in iron translocation, iron and zinc homeostasis, and grain size modifications. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection. The identified SNPs will be valuable in the rapid development of biofortified wheat varieties to ameliorate the malnutrition problems.

Genetic dissection of marker trait associations for grain micro-nutrients and thousand grain weight under heat and drought stress conditions in wheat.

Introduction: Wheat is grown and consumed worldwide, making it an important staple food crop for both its calorific and nutritional content. In places where wheat is used as a staple food, suboptimal micronutrient content levels, especially of grain iron (Fe) and zinc (Zn), can lead to malnutrition. Grain nutrient content is influenced by abiotic stresses, such as drought and heat stress. The best method for addressing micronutrient deficiencies is the biofortification of food crops. The prerequisites for marker-assisted varietal development are the identification of the genomic region responsible for high grain iron and zinc contents and an understanding of their genetics. Methods: A total of 193 diverse wheat genotypes were evaluated under drought and heat stress conditions across the years at the Indian Agricultural Research Institute (IARI), New Delhi, under timely sown irrigated (IR), restricted irrigated (RI) and late sown (LS) conditions. Grain iron content (GFeC) and grain zinc content (GZnC) were estimated from both the control and treatment groups. Genotyping of all the lines under study was carried out with the single nucleotide polymorphisms (SNPs) from Breeder's 35K Axiom Array. Result and Discussion: Three subgroups were observed in the association panel based on both principal component analysis (PCA) and dendrogram analysis. A large whole-genome linkage disequilibrium (LD) block size of 3.49 Mb was observed. A genome-wide association study identified 16 unique stringent marker trait associations for GFeC, GZnC, and 1000-grain weight (TGW). In silico analysis demonstrated the presence of 28 potential candidate genes in the flanking region of 16 linked SNPs, such as synaptotagmin-like mitochondrial-lipid-binding domain, HAUS augmin-like complex, di-copper center-containing domain, protein kinase, chaperonin Cpn60, zinc finger, NUDIX hydrolase, etc. Expression levels of these genes in vegetative tissues and grain were also found. Utilization of identified markers in marker-assisted breeding may lead to the rapid development of biofortified wheat genotypes to combat malnutrition.