Association of Agronomic Traits with SNP Markers in Durum Wheat (Triticum turgidum L. durum (Desf.))

Genome-Wide Association Study of Leaf Rust and Stem Rust Seedling and Adult Resistances in Tetraploid Wheat Accessions Harvested in Kazakhstan

Leaf rust (LR) and stem rust (SR) are diseases increasingly impacting wheat production worldwide. Fungal pathogens producing rust diseases in wheat may cause yield losses of up to 50−60%. One of the most effective methods for preventing such losses is the development of resistant cultivars with high yield potential. This goal can be achieved through complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. The objective of this study was to identify sources of tetraploid wheat resistance to LR and SR races, both at the seedling growth stage in the greenhouse and at the adult plant stage in field experiments, under the conditions of the North Kazakhstan region. A panel consisting of 193 tetraploid wheat accessions was used in a genome-wide association study (GWAS) for the identification of quantitative trait loci (QTLs) associated with LR and SR resistance, using 16,425 polymorphic single-nucleotide polymorphism (SNP) markers in the seedling and adult stages of plant development. The investigated panel consisted of seven tetraploid subspecies (Triticum turgidum ssp. durum, ssp. turanicum, ssp. turgidum, ssp. polonicum, ssp. carthlicum, ssp. dicoccum, and ssp. dicoccoides). The GWAS, based on the phenotypic evaluation of the tetraploid collection’s reaction to the two rust species at the seedling (in the greenhouse) and adult (in the field) stages, revealed 38 QTLs (p < 0.001), comprising 17 for LR resistance and 21 for SR resistance. Ten QTLs were associated with the reaction to LR at the seedling stage, while six QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. Eleven QTLs were associated with SR response at the seedling stage, while nine QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. A comparison of these results with previous LR and SR studies indicated that 11 of the 38 QTLs are presumably novel loci. The QTLs identified in this work can potentially be used for marker-assisted selection of tetraploid and hexaploid wheat for the breeding of new LR- and SR-resistant cultivars.

Regions of Chromosome 2A of Bread Wheat (Triticum aestivum L.) Associated with Variation in Physiological and Agronomical Traits under Contrasting Water Regimes

Understanding the genetic architecture of drought tolerance is of great importance for overcoming the negative impact of drought on wheat yield. Earlier, we discovered the critical role of chromosome 2A for the drought-tolerant status of wheat spring cultivar Saratovskaya 29. A set of 92 single-chromosome recombinant double haploid (SCRDH) lines were obtained in the genetic background of Saratovskaya 29. The lines carry fragments of chromosome 2A from the drought-sensitive cultivar Yanetzkis Probat. The SCRDH lines were used to identify regions on chromosome 2A associated with the manifestation of physiological and agronomical traits under distinct water supply, and to identify candidate genes that may be associated with adaptive gene networks in wheat. Genotyping was done with Illumina Infinium 15k wheat array using 590 SNP markers with 146 markers being polymorphic. In four identified regions of chromosome 2A, 53 out of 58 QTLs associated with physiological and agronomic traits under contrasting water supply were mapped. Thirty-nine candidate genes were identified, of which 18 were transcription factors. The region 73.8-78.1 cM included the largest number of QTLs and candidate genes. The variation in SNPs associated with agronomical and physiological traits revealed among the SCRDH lines may provide useful information for drought related marker-assisted breeding.

Linkage disequilibrium patterns, population structure and diversity analysis in a worldwide durum wheat collection including Argentinian genotypes

Background

Durum wheat (Triticum turgidum L. ssp. durum Desf. Husn) is the main staple crop used to make pasta products worldwide. Under the current climate change scenarios, genetic variability within a crop plays a crucial role in the successful release of new varieties with high yields and wide crop adaptation. In this study we evaluated a durum wheat collection consisting of 197 genotypes that mainly comprised a historical set of Argentinian germplasm but also included worldwide accessions.

Results

We assessed the genetic diversity, population structure and linkage disequilibrium (LD) patterns in this collection using a 35 K SNP array. The level of polymorphism was considered, taking account of the frequent and rare allelic variants. A total of 1547 polymorphic SNPs was located within annotated genes. Genetic diversity in the germplasm collection increased slightly from 1915 to 2010. However, a reduction in genetic diversity using SNPs with rare allelic variants was observed after 1979. However, larger numbers of rare private alleles were observed in the 2000–2009 period, indicating that a high reservoir of rare alleles is still present among the recent germplasm in a very low frequency. The percentage of pairwise loci in LD in the durum genome was low (13.4%) in our collection. Overall LD and the high (r² > 0.7) or complete (r² = 1) LD presented different patterns in the chromosomes. The LD increased over three main breeding periods (1915–1979, 1980–1999 and 2000–2020).

Conclusions

Our results suggest that breeding and selection have impacted differently on the A and B genomes, particularly on chromosome 6A and 2A. The collection was structured in five sub-populations and modern Argentinian accessions (cluster Q4) which were clearly differentiated. Our study contributes to the understanding of the complexity of Argentinian durum wheat germplasm and to derive future breeding strategies enhancing the use of genetic diversity in a more efficient and targeted way.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07519-z.

Molecular Mapping and Genomics of Grain Yield in Durum Wheat: A Review

Durum wheat is the most relevant cereal for the whole of Mediterranean agriculture, due to its intrinsic adaptation to dryland and semi-arid environments and to its strong historical cultivation tradition. It is not only relevant for the primary production sector, but also for the food industry chains associated with it. In Mediterranean environments, wheat is mostly grown under rainfed conditions and the crop is frequently exposed to environmental stresses, with high temperatures and water scarcity especially during the grain filling period. For these reasons, and due to recurrent disease epidemics, Mediterranean wheat productivity often remains under potential levels. Many studies, using both linkage analysis (LA) and a genome-wide association study (GWAS), have identified the genomic regions controlling the grain yield and the associated markers that can be used for marker-assisted selection (MAS) programs. Here, we have summarized all the current studies identifying quantitative trait loci (QTLs) and/or candidate genes involved in the main traits linked to grain yield: kernel weight, number of kernels per spike and number of spikes per unit area.

SNP-based association study of kernel architecture in a worldwide collection of durum wheat germplasm

Durum wheat, genetic resource with favorable alleles is considered as natural gene pool for wheat breeding. Kernel size and weight are important factors affecting grain yield in crops. Here, association analysis was performed to dissect the genetic constitution of kernel-related traits in 150 lines collected from 46 countries and regions using a set of EST-derived and genome-wide SNP markers with five consecutive years of data. Total 109 significant associations for eight kernel-related traits were detected under a mix linear model, generating 54 unique SNP markers distributed on 13 of 14 chromosomes. Of which, 19 marker-trait associations were identified in two or more environments, including one stable and pleiotropic SNP BE500291_5_A_37 on chromosome 5A correlated with six kernel traits. Although most of our SNP loci were overlapped with the previously known kernel weight QTLs, several novel loci for kernel traits in durum were reported. Correlation analysis implied that the moderate climatic variables during growth and development of durum are needed for the large grain size and high grain weight. Combined with our previous studies, we found that chromosome 5A might play an important role in durum growth and development.

Genome-wide association and epistatic interactions of flowering time in soybean cultivar

Genome-wide association studies (GWAS) have enabled the discovery of candidate markers that play significant roles in various complex traits in plants. Recently, with increased interest in the search for candidate markers, studies on epistatic interactions between single nucleotide polymorphism (SNP) markers have also increased, thus enabling the identification of more candidate markers along with GWAS on single-variant-additive-effect. Here, we focused on the identification of candidate markers associated with flowering time in soybean (Glycine max). A large population of 2,662 cultivated soybean accessions was genotyped using the 180k Axiom® SoyaSNP array, and the genomic architecture of these accessions was investigated to confirm the population structure. Then, GWAS was conducted to evaluate the association between SNP markers and flowering time. A total of 93 significant SNP markers were detected within 59 significant genes, including E1 and E3, which are the main determinants of flowering time. Based on the GWAS results, multilocus epistatic interactions were examined between the significant and non-significant SNP markers. Two significant and 16 non-significant SNP markers were discovered as candidate markers affecting flowering time via interactions with each other. These 18 candidate SNP markers mapped to 18 candidate genes including E1 and E3, and the 18 candidate genes were involved in six major flowering pathways. Although further biological validation is needed, our results provide additional information on the existing flowering time markers and present another option to marker-assisted breeding programs for regulating flowering time of soybean.

Development of molecular markers for invasive alien plants in Korea: a case study of a toxic weed, Cenchrus longispinus L., based on next generation sequencing data

BACKGROUND

Genomic data play an important role in plant research because of its implications in studying genomic evolution, phylogeny, and developing molecular markers. Although the information of invasive alien plants was collected, the genomic data of those species have not been intensively studied.

METHODS

We employ the next generation sequencing and PCR methods to explore the genomic data as well as to develop and test the molecular markers.

RESULTS

In this study, we characterize the chloroplast genomes (cpDNA) of Cenchrus longispinus and C. echinatus, of which the lengths are 137,144 and 137,131 bp, respectively. These two newly sequenced genomes include 78 protein-coding genes, 30 tRNA, and four rRNA. There are 56 simple single repeats and 17 forward repeats in the chloroplast genome of C. longispinus. Most of the repeats locate in non-coding regions. However, repeats can be found in infA, ndhD, ndhH, ndhK, psbC, rpl22, rpoC2, rps14, trnA-UGC, trnC-GCA, trnF-GAA, trnQ-UUG, trnS-UGA, trnS-GCU, and ycf15. The phylogenomic analysis revealed the monophyly of Cenchrus but not Panicum species in tribe Paniceae. The single nucleotide polymorphism sites in atpB, matK, and ndhD were successfully used for developing molecular markers to distinguish C. longispinus and related taxa. The simple PCR protocol for using the newly developed molecular markers was also provided.

Genetic analysis of threshability and other spike traits in the evolution of cultivated emmer to fully domesticated durum wheat

Genetic mutations in genes governing wheat threshability were critical for domestication. Knowing when these genes mutated during wheat evolution will provide more insight into the domestication process and lead to further exploitation of primitive alleles for wheat improvement. We evaluated a population of recombinant inbred lines derived from a cross between the durum variety Rusty and the cultivated emmer accession PI 193883 for threshability, rachis fragility, and other spike-related traits. Quantitative trait loci (QTL) associated with spike length, spikelets per spike, and spike compactness were primarily associated with known genes such as the pleiotropic domestication gene Q. Interestingly, rachis fragility was not associated with the Q locus, suggesting that this trait, usually a pleiotropic effect of the q allele, can be influenced by the genetic background. Threshability QTL were identified on chromosome arms 2AS, 2BS, and 5AL corresponding to the tenacious glume genes Tg^2A and Tg^2B as well as the Q gene, respectively, further demonstrating that cultivated emmer harbors the primitive non-free-threshing alleles at all three loci. Genetic analysis indicated that the effects of the three genes are mostly additive, with Q having the most profound effects on threshability, and that free-threshing alleles are necessary at all three loci to attain a completely free-threshing phenotype. These findings provide further insight into the timeline and possible pathways of wheat domestication and evolution that led to the formation of modern day domesticated wheats.

Haplotype Loci Under Selection in Canadian Durum Wheat Germplasm Over 60 Years of Breeding: Association With Grain Yield, Quality Traits, Protein Loss, and Plant Height

Durum wheat was introduced in the southern prairies of western Canada in the late nineteenth century. Breeding efforts have mainly focused on improving quality traits to meet the pasta industry demands. For this study, 192 durum wheat lines were genotyped using the Illumina 90K Infinium iSelect assay, and resulted in a total of 14,324 polymorphic SNPs. Genetic diversity changed over time, declining during the first 20 years of breeding in Canada, then increased in the late 1980s and early 1990s. We scanned the genome for signatures of selection, using the total variance Fst-based outlier detection method (Lositan), the hierarchical island model (Arlequin) and the Bayesian genome scan method (BayeScan). A total of 407 outliers were identified and clustered into 84 LD-based haplotype loci, spanning all 14 chromosomes of the durum wheat genome. The association analysis detected 54 haplotype loci, of which 39% contained markers with a complete reversal of allelic state. This tendency to fixation of favorable alleles corroborates the success of the Canadian durum wheat breeding programs over time. Twenty-one haplotype loci were associated with multiple traits. In particular, hap_4B_1 explained 20.6, 17.9 and 16.6% of the phenotypic variance of pigment loss, pasta b∗ and dough extensibility, respectively. The locus hap_2B_9 explained 15.9 and 17.8% of the variation of protein content and protein loss, respectively. All these pleiotropic haplotype loci offer breeders the unique opportunity for further improving multiple traits, facilitating marker-assisted selection in durum wheat, and could help in identifying genes as functional annotations of the wheat genome become available.

Associations of canopy leaf traits with SNP markers in durum wheat (Triticum turgidum L. durum (Desf.))

The canopy leaves including the top three, i.e., the flag, the 2nd and 3rd from the top, are important for photosynthesis and grain yield of wheat. Molecular markers associated with traits of these leaves should be helpful for the high-yielding breeding. In this study, 1366 single nucleotide polymorphisms (SNP) markers covering the whole genome of durum wheat were used to genotype 150 cultivars collected from 46 countries and regions in the world. Leaf length, leaf width and chlorophyll content of the top three leaves were measured, respectively, in three consecutive years. Association analyses were performed on the leaf traits and SNP markers. A total of 120 SNP marker associations were detected on 13 of the 14 chromosomes. Among these markers, 83 were associated with the canopy leaf traits, 10 with 1000-grain weight, and 29 with kernel number per spike. This study is helpful for better understanding the potential and genetic basis of functional leaves, and facilitates pyramiding of the favorable alleles using marker assisted selection for ideal plant-type and high photosynthesis efficiency in durum wheat breeding.

Wheat genetic resources in the post-genomics era: promise and challenges

Background

Wheat genetic resources have been used for genetic improvement since 1876, when Stephen Wilson (Transactions and Proceedings of the Botanical Society of Edinburgh 12: 286) consciously made the first wide hybrid involving wheat and rye in Scotland. Wide crossing continued with sporadic attempts in the first half of 19th century and became a sophisticated scientific discipline during the last few decades with considerable impact in farmers' fields. However, a large diversity of untapped genetic resources could contribute in meeting future wheat production challenges.

Perspectives and Conclusion

Recently the complete reference genome of hexaploid (Chinese Spring) and tetraploid (Triticum turgidum ssp. dicoccoides) wheat became publicly available coupled with on-going international efforts on wheat pan-genome sequencing. We anticipate that an objective appraisal is required in the post-genomics era to prioritize genetic resources for use in the improvement of wheat production if the goal of doubling yield by 2050 is to be met. Advances in genomics have resulted in the development of high-throughput genotyping arrays, improved and efficient methods of gene discovery, genomics-assisted selection and gene editing using endonucleases. Likewise, ongoing advances in rapid generation turnover, improved phenotyping, envirotyping and analytical methods will significantly accelerate exploitation of exotic genes and increase the rate of genetic gain in breeding. We argue that the integration of these advances will significantly improve the precision and targeted identification of potentially useful variation in the wild relatives of wheat, providing new opportunities to contribute to yield and quality improvement, tolerance to abiotic stresses, resistance to emerging biotic stresses and resilience to weather extremes.

Environmental and Genetic Variation for Yield-Related Traits of Durum Wheat as Affected by Development

Phenology has a profound effect on adaptation and productivity of crops. The impact of phenology on tillering and fertility traits of durum wheat (Triticum turgidum L. subsp. durum Desf.) was evaluated with the aim of specifying which group of flowering genes (Vrn, Ppd, or eps) was involved in their control. A recombinant inbred line population was grown under four contrasting conditions of vernalization and daylength. Phenotyping was carried out according to robust phenological models dissecting both phenological and yield related traits. Whole-genome profiling was performed using the DArT-Seq technology. The genetic variability for tillering was mainly related to the genetic variability for vernalization sensitivity, as shown by the many quantitative trait loci (QTLs) identified in non-vernalized plants associated to both tillering and phenological traits. No effects of photoperiod sensitivity on spikelet number were detected in short-day-grown plants, apparently because of limited genetic variability in photoperiod sensitivity of the population. Earliness per se was involved in control of spikelet number via final leaf number, with these traits genetically correlated and sharing some QTLs. Chaff weight and number of kernels per g chaff were negatively associated and related to anthesis date under most conditions. QTL mapping uncovered novel loci involved in phenological control of tillering and fertility traits, and confirmed the presence of several well-established loci. Phenotyping of both phenology and kernel number according to a robust physiological model amplified the possibility of identifying genetic factors underlying their variations. Also, isolating known flowering gene cues by manipulation of environmental conditions provided the opportunity for each group of genes to be expressed without confounding effects of the others. This information helps to predict the consequences of either genetic manipulation of flowering genes and changes in environmental conditions on the potential yield of durum wheat.

Multi-Locus Genome-Wide Association Studies for 14 Main Agronomic Traits in Barley

The agronomic traits, including morphological and yield component traits, are important in barley breeding programs. In order to reveal the genetic foundation of agronomic traits of interest, in this study 122 doubled haploid lines from a cross between cultivars "Huaai 11" (six-rowed and dwarf) and "Huadamai 6" (two-rowed) were genotyped by 9680 SNPs and phenotyped 14 agronomic traits in 3 years, and the two datasets were used to conduct multi-locus genome-wide association studies. As a result, 913 quantitative trait nucleotides (QTNs) were identified by five multi-locus GWAS methods to be associated with the above 14 traits and their best linear unbiased predictions. Among these QTNs and their adjacent genes, 39 QTNs (or QTN clusters) were repeatedly detected in various environments and methods, and 10 candidate genes were identified from gene annotation. Nineteen QTNs and two genes (sdw1/denso and Vrs1) were previously reported, and eight candidate genes need to be further validated. The Vrs1 gene, controlling the number of rows in the spike, was found to be associated with spikelet number of main spike, spikelet number per plant, grain number per plant, grain number per spike, and 1,000 grain weight in multiple environments and by multi-locus GWAS methods. Therefore, the above results evidenced the feasibility and reliability of genome-wide association studies in doubled haploid population, and the QTNs and their candidate genes detected in this study are useful for marker-assisted selection breeding, gene cloning, and functional identification in barley.

Dissecting the old Mediterranean durum wheat genetic architecture for phenology, biomass and yield formation by association mapping and QTL meta-analysis

Association mapping was used to identify genome regions affecting yield formation, crop phenology and crop biomass in a collection of 172 durum wheat landraces representative of the genetic diversity of ancient local durum varieties from the Mediterranean Basin. The collection was genotyped with 1,149 DArT markers and phenotyped in Spanish northern and southern locations during three years. A total of 245 significant marker trait associations (MTAs) (P<0.01) were detected. Some of these associations confirmed previously identified quantitative trait loci (QTL) and/or candidate genes, and others are reported for the first time here. Eighty-six MTAs corresponded with yield and yield component traits, 70 to phenology and 89 to biomass production. Twelve genomic regions harbouring stable MTAs (significant in three or more environments) were identified, while five and two regions showed specific MTAs for northern and southern environments, respectively. Sixty per cent of MTAs were located on the B genome and 29% on the A genome. The marker wPt-9859 was detected in 12 MTAs, associated with six traits in four environments and the mean across years. To refine QTL positions, a meta-analysis was performed. A total of 477 unique QTLs were projected onto a durum wheat consensus map and were condensed to 71 meta-QTLs and left 13 QTLs as singletons. Sixty-one percent of QTLs explained less than 10% of the phenotypic variance confirming the high genetic complexity of the traits analysed.

Solar Radiation-Associated Adaptive SNP Genetic Differentiation in Wild Emmer Wheat, Triticum dicoccoides

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.

SNP-based association analysis for seedling traits in durum wheat (Triticum turgidum L. durum (Desf.))

In the present study, 150 accessions of worldwide originated durum wheat germplasm (Triticum turgidum spp. durum) were observed for major seedling traits and their growth. The accessions were evaluated for major seedling traits under controlled conditions of hydroponics at the 13th, 20th, 27th and 34th day-after germination. Biomass traits were measured at the 34th day-after germination. Correlation analysis was conducted among the seedling traits and three field traits at maturity, plant height, grain weight and 1000-grain weight observed in four consecutive years. Associations of the measured seedling traits and SNP markers were analyzed based on the mixed linear model (MLM). The results indicated that highly significant genetic variation and robust heritability were found for the seedling and field mature traits. In total, 259 significant associations were detected for all the traits and four growth stages. The phenotypic variation explained (R2) by a single SNP marker is higher than 10% for most (84%) of the significant SNP markers. Forty-six SNP markers associated with multiple traits, indicating non-neglectable pleiotropy in seedling stage. The associated SNP markers could be helpful for genetic analysis of seedling traits, and marker-assisted breeding of new wheat varieties with strong seedling vigor.

Genome-Wide Single-Nucleotide Polymorphisms in CMS and Restorer Lines Discovered by Genotyping Using Sequencing and Association with Marker-Combining Ability for 12 Yield-Related Traits in Oryza sativa L. subsp. Japonica

Heterosis or hybrid vigor is closely related with general combing ability (GCA) of parents and special combining ability (SCA) of combinations. The evaluation of GCA and SCA facilitate selection of parents and combinations in heterosis breeding. In order to improve combining ability (CA) by molecular marker assist selection, it is necessary to identify marker loci associated with the CA. To identify the single nucleotide polymorphisms (SNP) loci associated with CA in the parental genomes of japonica rice, genome-wide discovered SNP loci were tested for association with the CA of 18 parents for 12 yield-related traits. In this study, 81 hybrids were created and evaluated to calculate the CA of 18 parents. The parents were sequenced by genotyping by sequencing (GBS) method for identification of genome-wide SNPs. The analysis of GBS indicated that the successful mapping of 9.86 × 10⁶ short reads in the Nipponbare reference genome consists of 39,001 SNPs in parental genomes at 11,085 chromosomal positions. The discovered SNPs were non-randomly distributed within and among the 12 chromosomes of rice. Overall, 20.4% (8026) of the discovered SNPs were coding types, and 8.6% (3344) and 9.9% (3951) of the SNPs revealed synonymous and non-synonymous changes, which provide valuable knowledge about the underlying performance of the parents. Furthermore, the associations between SNPs and CA indicated that 362 SNP loci were significantly related to the CA of 12 parental traits. The identified SNP loci of CA in our study were distributed genome wide and caused a positive or negative effect on the CA of traits. For the yield-related traits, such as grain thickness, days to heading, panicle length, grain length and 1000-grain weight, a maximum number of positive SNP loci of CA were found in CMS A171 and in the restorers LC64 and LR27. On an individual basis, some of associated loci that resided on chromosomes 2, 5, 7, 9, and 11 recorded maximum positive values for the CA of traits. From our results, we suggest that heterosis in japonica rice would be improved by pyramiding the favorable SNP loci of CA and eliminating the unfavorable loci from parental genomes.

Genome-Wide Association Mapping of Leaf Rust Response in a Durum Wheat Worldwide Germplasm Collection

Leaf rust (caused by Erikss. []) is increasingly impacting durum wheat ( L. var. ) production with the recent appearance of races with virulence to widely grown cultivars in many durum producing areas worldwide. A highly virulent race on durum wheat was recently detected in Kansas. This race may spread to the northern Great Plains, where most of the US durum wheat is produced. The objective of this study was to identify sources of resistance to several races from the United States and Mexico at seedling stage in the greenhouse and at adult stage in field experiments. Genome-wide association study (GWAS) was used to identify single-nucleotide polymorphism (SNP) markers associated with leaf rust response in a worldwide durum wheat collection of 496 accessions. Thirteen accessions were resistant across all experiments. Association mapping revealed 88 significant SNPs associated with leaf rust response. Of these, 33 SNPs were located on chromosomes 2A and 2B, and 55 SNPs were distributed across all other chromosomes except for 1B and 7B. Twenty markers were associated with leaf rust response at seedling stage, while 68 markers were associated with leaf rust response at adult plant stage. The current study identified a total of 14 previously uncharacterized loci associated with leaf rust response in durum wheat. The discovery of these loci through association mapping (AM) is a significant step in identifying useful sources of resistance that can be used to broaden the relatively narrow leaf rust resistance spectrum in durum wheat germplasm.