Genotyping-by-sequencing (GBS) is a high-throughput genotyping approach that is starting to be used in several crop species, including bread wheat. Anchoring GBS tags on chromosomes is an important step towards utilizing them for wheat genetic improvement. Here we use genetic linkage mapping to construct a consensus map containing 28644 GBS markers. Three RIL populations, PBW343 × Kingbird, PBW343 × Kenya Swara and PBW343 × Muu, which share a common parent, were used to minimize the impact of potential structural genomic variation on consensus-map quality. The consensus map comprised 3757 unique positions, and the average marker distance was 0.88 cM, obtained by calculating the average distance between two adjacent unique positions. Significant variation of segregation distortion was observed across the three populations. The consensus map was validated by comparing positions of known rust resistance genes, and comparing them to wheat reference genome sequences recently published by the International Wheat Genome Sequencing Consortium, Rye and Ae. tauschii genomes. Three well-characterized rust resistance genes (Sr58/Lr46/Yr29, Sr2/Yr30/Lr27, and Sr57/Lr34/Yr18) and 15 published QTLs for wheat rusts were validated with high resolution. Fifty-two per cent of GBS tags on the consensus map were successfully aligned through BLAST to the right chromosomes on the wheat reference genome sequence. The consensus map should provide a useful basis for analyzing genome-wide variation of complex traits. The identified genes can then be explored as genetic markers to be used in genomic applications in wheat breeding.

Tan spot, caused byPyrenophora tritici-repentis, is a serious foliar disease of wheat in Kazakhstan with reported yield losses as high as 50% during epidemic years. Here, we report the evaluation of a collection of 191 hexaploid spring and winter wheat lines for tan spot resistance and its underlying genetic architecture using genome-wide association study (GWAS). Our wheat collection comprised candidate varieties from Kazakhstan, Russia, and CIMMYT. It was genotyped using the DArTseq technology and phenotyped for resistance to tan spot at seedling and adult plant stages in Kazakhstan. DArTseq SNPs revealed high genetic diversity (average polymorphic information content = 0.33) in the panel and genome-wide linkage disequilibrium decay at 22 Mb (thresholdr²= 0.1). Principal component analysis revealed a clear separation of Eurasian germplasm from CIMMYT and IWWIP lines. GWAS identified 34 marker-trait associations (MTA) for resistance to tan spot and the amount of phenotypic variation explained by these MTA ranged from 4% to 13.7%. Our results suggest the existence of novel valuable resistant alleles on chromosomes 3BS, and 5DL and 6AL for resistance to Race 1 and Race 5, respectively, in addition to known genestsn1andtsc2.On chromosome 6AL, a genomic region spanning 3 Mb was identified conferring resistance to both Race 1 and Race 5. Epistatic interaction of associated loci was revealed on chromosomes 1B, 5B, 7B, 5A, and 6A contributing to additional variation of 3.2–11.7%. Twenty-five lines with the best allele combinations of SNPs associated with resistance to both races have been identified as candidates for future variety release and breeding. The results of the present study will be further validated in other independent genetic backgrounds to be able to use markers in breeding.

Yellow rust (YR) or stripe rust, caused by Puccinia striformis f. sp tritici Eriks (Pst), is a major challenge to resistance breeding in wheat. A genome wide association study (GWAS) was performed using 22,415 single nucleotide polymorphism (SNP) markers and 591 haplotypes to identify genomic regions associated with resistance to YR in a subset panel of 419 pre-breeding lines (PBLs) developed at International Center for Maize and Wheat Improvement (CIMMYT). The 419 PBLs were derived from an initial set of 984 PBLs generated by a three-way crossing scheme (exotic/elite1//elite2) among 25 best elites and 244 exotics (synthetics, landraces) from CIMMYT’s germplasm bank. For the study, 419 PBLs were characterized with 22,415 high-quality DArTseq-SNPs and phenotyped for severity of YR disease at five locations in Mexico. A population structure was evident in the panel with three distinct subpopulations, and a genome-wide linkage disequilibrium (LD) decay of 2.5 cM was obtained. Across all five locations, 14 SNPs and 7 haplotype blocks were significantly (P < 0.001) associated with the disease severity explaining 6.0 to 14.1% and 7.9 to 19.9% of variation, respectively. Based on average LD decay of 2.5 cM, identified 14 SNP–trait associations were delimited to seven quantitative trait loci in total. Seven SNPs were part of the two haplotype blocks on chromosome 2A identified in haplotypes-based GWAS. In silico analysis of the identified SNPs showed hits with interesting candidate genes, which are related to pathogenic process or known to regulate induction of genes related to pathogenesis such as those coding for glunolactone oxidase, quinate O-hydroxycinnamoyl transferase, or two-component histidine kinase. The two-component histidine kinase, for example, acts as a sensor in the perception of phytohormones ethylene and cytokinin. Ethylene plays a very important role in regulation of multiple metabolic processes of plants, including induction of defense mechanisms mediated by jasmonate. The SNPs linked to the promising genes identified in the study can be used for marker-assisted selection.

Terminal drought stress poses a big challenge to sustain wheat grain production in rain-fed environments. This study aimed to utilize the genetically diverse pre-breeding lines for identification of genomic regions associated with agro-physiological traits at terminal stage drought stress in wheat. A total of 339 pre-breeding lines panel derived from three-way crosses of ‘exotics × elite × elite’ lines were evaluated in field conditions at Obregon, Mexico for two years under well irrigated as well as drought stress environments. Drought stress was imposed at flowering by skipping the irrigations at pre and post anthesis stage. Results revealed that drought significantly reduced grain yield (Y), spike length (SL), number of grains spikes⁻¹ (NGS) and thousand kernel weight (TKW), while kernel abortion (KA) was increased. Population structure analysis in this panel uncovered three sub-populations. Genome wide linkage disequilibrium (LD) decay was observed at 2.5 centimorgan (cM). The haplotypes-based genome wide association study (GWAS) identified significant associations of Y, SL, and TKW on three chromosomes; 4A (HB10.7), 2D (HB6.10) and 3B (HB8.12), respectively. Likewise, associations on chromosomes 6B (HB17.1) and 3A (HB7.11) were found for NGS while on chromosome 3A (HB7.12) for KA. The genomic analysis information generated in the study can be efficiently utilized to improve Y and/or related parameters under terminal stage drought stress through marker-assisted breeding.

The two recombinant inbred line (RIL) populations developed by crossing Almaly × Avocet S (206 RILs) and Almaly × Anza (162 RILs) were used to detect the novel genomic regions associated with adult plant resistance (APR) and seedling or all-stage resistance (ASR) to yellow rust (YR) and leaf rust (LR). The quantitative trait loci (QTLs) were detected through multi-year phenotypic evaluations (2018–2020) and using high-throughput DArTseq genotyping technology. RILs exhibited significant genetic variation with p < 0.001, and the coefficient of variation ranged from 9.79% to 47.99% for both LR and YR in all Environments and stages of evaluations. The heritability is quite high and ranged between 0.47 and 0.98. We identified nine stable QTLs for YR APR on chromosomes 1B, 2A, 2B, 3D, and 4D and four stable QTLs for LR APR on chromosomes 2B, 3B, 4A, and 5A. Furthermore, in silico analysis revealed that the key putative candidate genes such as cytochrome P450, protein kinase-like domain superfamily, zinc-binding ribosomal protein, SANT/Myb domain, WRKY transcription factor, nucleotide sugar transporter, and NAC domain superfamily were in the QTL regions and probably involved in the regulation of host response toward pathogen infection. The stable QTLs identified in this study are useful for developing rust-resistant varieties through marker-assisted selection (MAS).

Crown rot (CR), caused by various Fusarium species, is a major disease in many cereal-growing regions worldwide. Fusarium culmorum is one of the most important species, which can cause significant yield losses in wheat. A set of 126 advanced International Maize and Wheat Improvement Center (CIMMYT) spring bread wheat lines were phenotyped against CR for field crown, greenhouse crown and stem, and growth room crown resistance scores. Of these, 107 lines were genotyped using Diversity Array Technology (DArT) markers to identify quantitative trait loci linked to CR resistance by genome-wide association study. Results of the population structure analysis grouped the accessions into three sub-groups. Genome wide linkage disequilibrium was large and declined on average within 20 cM (centi-Morgan) in the panel. General linear model (GLM), mixed linear model (MLM), and naïve models were tested for each CR score and the best model was selected based on quarantine-quarantine plots. Three marker-trait associations (MTAs) were identified linked to CR resistance; two of these on chromosome 3B were associated with field crown scores, each explaining 11.4% of the phenotypic variation and the third MTA on chromosome 2D was associated with greenhouse stem score and explained 11.6% of the phenotypic variation. Together, these newly identified loci provide opportunity for wheat breeders to exploit in enhancing CR resistance via marker-assisted selection or deployment in genomic selection in wheat breeding programs.