Genome-wide association mapping for adult resistance to powdery mildew in common wheat

Genetic mapping of loci affecting seedling and adult-plant resistance to powdery mildew derived from two CIMMYT wheat lines

MAIN CONCLUSION

PM3 and PM8 alleles carried by two CIMMYT wheat lines confer powdery mildew resistance in seedlings and/or adult plants. A stage-specific epistatic interaction was observed between PM3 and PM8. Powdery mildew is an important foliar disease of wheat. Major genes for resistance, which have been widely used in wheat breeding programs, are typically effective against only limited numbers of virulence genes of the pathogen. The main aim of this study was to map resistance loci in wheat lines 7HRWSN58 and ZWW09-149 from the International Maize and Wheat Improvement Center (CIMMYT). Doubled haploid populations (Magenta/7HRWSN58 and Emu Rock/ZWW09-149) were developed and grown in controlled environment experiments and inoculated with a composite of Blumeria graminis f.sp. tritici isolates that had been collected at various locations in Western Australia. Plants were assessed for powdery mildew symptoms (percentage leaf area diseased) on seedlings and adult plants. Populations were subjected to genotyping-by-sequencing and assayed for known SNPs in the resistance gene PM3. Linkage maps were constructed, and markers were anchored to the wheat reference genome sequence. In both populations, there were asymptomatic lines that exhibited no symptoms. Among symptomatic lines, disease severity varied widely. In the Magenta/7HRWSN58 population, most of the observed variation was attributed to the PM3 region of chromosome 1A, with the allele from 7HRWSN58 conferring resistance in seedlings and adult plants. In the Emu Rock/ZWW09-149 population, two interacting quantitative trait loci were mapped: one at PM3 and the other on chromosome 1B. The Emu Rock/ZWW09-149 population was confirmed to segregate for a 1BL·1RS translocation that carries the PM8 powdery mildew resistance gene from rye. Consistent with previous reports that PM8-derived resistance can be suppressed by PM3 alleles, the observed interaction between the quantitative trait loci on chromosomes 1A and 1B indicated that the PM3 allele carried by ZWW09-149 suppresses PM8-derived resistance from ZWW09-149, but only at the seedling stage. In adult plants, the PM8 region conferred resistance regardless of the PM3 genotype. The resistance sources and molecular markers that were investigated here could be useful in wheat breeding.

Genome-Wide Association Study Reveals Novel Powdery Mildew Resistance Loci in Bread Wheat

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-2021, 2021-2022, and 2022-2023), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K Infinium iSelect assay to obtain genome-wide single-nucleotide polymorphism (SNP) marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 marker-trait associations (MTAs) located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p ≤ 0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-2021 and 2021-2022. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interaction-related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Distinct Clades of Protein Phosphatase 2A Regulatory B'/B56 Subunits Engage in Different Physiological Processes

Protein phosphatase 2A (PP2A) is a strongly conserved and major protein phosphatase in all eukaryotes. The canonical PP2A complex consists of a catalytic (C), scaffolding (A), and regulatory (B) subunit. Plants have three groups of evolutionary distinct B subunits: B55, B' (B56), and B''. Here, the Arabidopsis B' group is reviewed and compared with other eukaryotes. Members of the B'α/B'β clade are especially important for chromatid cohesion, and dephosphorylation of transcription factors that mediate brassinosteroid (BR) signaling in the nucleus. Other B' subunits interact with proteins at the cell membrane to dampen BR signaling or harness immune responses. The transition from vegetative to reproductive phase is influenced differentially by distinct B' subunits; B'α and B'β being of little importance, whereas others (B'γ, B'ζ, B'η, B'θ, B'κ) promote transition to flowering. Interestingly, the latter B' subunits have three motifs in a conserved manner, i.e., two docking sites for protein phosphatase 1 (PP1), and a POLO consensus phosphorylation site between these motifs. This supports the view that a conserved PP1-PP2A dephosphorelay is important in a variety of signaling contexts throughout eukaryotes. A profound understanding of these regulators may help in designing future crops and understand environmental issues.

Novel genomic regions on chromosome 5B controlling wheat powdery mildew seedling resistance under Egyptian conditions

Wheat powdery mildew (PM) causes significant yield losses worldwide. None of the Egyptian wheat cultivars was detected to be highly resistant to such a severe disease. Therefore, a diverse spring wheat panel was evaluated for PM seedling resistance using different Bgt conidiospores collected from Egyptian fields in two growing seasons. The evaluation was done in two separate experiments. Highly significant differences were found between the two experiments suggesting the presence of different isolates populations. Highly significant differences were found among the tested genotypes confirming the ability to improve PM resistance using the recent panel. Genome-wide association study (GWAS) was done for each experiment separately and a total of 71 significant markers located within 36 gene models were identified. The majority of these markers are located on chromosome 5B. Haplotype block analysis identified seven blocks containing the significant markers on chromosome 5B. Five gene models were identified on the short arm of the chromosome. Gene enrichment analysis identified five and seven pathways based on the biological process and molecular functions respectively for the detected gene models. All these pathways are associated with disease resistance in wheat. The genomic regions on 5B seem to be novel regions that are associated with PM resistance under Egyptian conditions. Selection of superior genotypes was done and Grecian genotypes seem to be a good source for improving PM resistance under Egyptian conditions.

Genome-wide association study of yield-related traits in common wheat (Triticum aestivum L.) under normal and drought treatment conditions

The primary goal of modern wheat breeding is to develop new high-yielding and widely adaptable varieties. We analyzed four yield-related agronomic traits in 502 wheat accessions under normal conditions (NC) and drought treatment (DT) conditions over three years. The genome-wide association analysis identified 51 yield-related and nine drought-resistance-related QTL, including 13 for the thousand-grain weight (TGW), 30 for grain length (GL), three for grain width (GW), five for spike length (SL) and nine for stress tolerance index (STI) QTL in wheat. These QTL, containing 72 single nucleotide polymorphisms (SNPs), explained 2.23 - 7.35% of the phenotypic variation across multiple environments. Eight stable SNPs on chromosomes 2A, 2D, 3B, 4A, 5B, 5D, and 7D were associated with phenotypic stability under NC and DT conditions. Two of these stable SNPs had association with TGW and STI. Several novel QTL for TGW, GL and SL were identified on different chromosomes. Three linked SNPs were transformed into kompetitive allele-specific PCR (KASP) markers. These results will facilitate the discovery of promising SNPs for yield-related traits and/or drought stress tolerance and will accelerate the development of new wheat varieties with desirable alleles.

Genome-wide association study reveals quantitative trait loci for waterlogging-triggered adventitious roots and aerenchyma formation in common wheat

Waterlogging severely affects wheat growth and development. Limited availability of oxygen in the root zone negatively affects the metabolism of plants. The formation of adventitious roots (ARs) and root cortical aerenchyma (RCA) are the most important adaptive trait contributing to plants' ability to survive in waterlogged soil conditions. This study used a genome-wide association study (GWAS) approach with 90K single nucleotide polymorphisms (SNPs) in a panel of 329 wheat genotypes, to reveal quantitative trait loci (QTL) conferring ARs and RCA. The wheat genotypes exposed to waterlogging were evaluated for ARs and RCA in both field and glasshouse over two consecutive years. Six and five significant marker-trait associations (MTAs) were identified for ARs and RCA formation under waterlogging, respectively. The most significant MTA for AR and RCA was found on chromosome 4B. Two wheat cultivars with contrasting waterlogging tolerance (tolerant: H-242, sensitive: H-195) were chosen to compare the development and regulation of aerenchyma in waterlogged conditions using staining methods. Results showed that under waterlogging conditions, H₂O₂ signal generated before aerenchyma formation in both sensitive and tolerant varieties with the tolerant variety accumulating more H₂O₂ and in a quicker manner compared to the sensitive one. Several genotypes which performed consistently well under different conditions can be used in breeding programs to develop waterlogging-tolerant wheat varieties.

Fine mapping of QPm.caas-3BS, a stable QTL for adult-plant resistance to powdery mildew in wheat (Triticum aestivum L.)

A stable QTL QPm.caas-3BS for adult-plant resistance to powdery mildew was mapped in an interval of 431 kb, and candidate genes were predicted based on gene sequences and expression profiles. Powdery mildew is a devastating foliar disease occurring in most wheat-growing areas. Characterization and fine mapping of genes for powdery mildew resistance can benefit marker-assisted breeding. We previously identified a stable quantitative trait locus (QTL) QPm.caas-3BS for adult-plant resistance to powdery mildew in a recombinant inbred line population of Zhou8425B/Chinese Spring by phenotyping across four environments. Using 11 heterozygous recombinants and high-density molecular markers, QPm.caas-3BS was delimited in a physical interval of approximately 3.91 Mb. Based on re-sequenced data and expression profiles, three genes TraesCS3B02G014800, TraesCS3B02G016800 and TraesCS3B02G019900 were associated with the powdery mildew resistance locus. Three gene-specific kompetitive allele-specific PCR (KASP) markers were developed from these genes and validated in the Zhou8425B derivatives and Zhou8425B/Chinese Spring population in which the resistance gene was mapped to a 0.3 cM interval flanked by KASP14800 and snp_50465, corresponding to a 431 kb region at the distal end of chromosome 3BS. Within the interval, TraesCS3B02G014800 was the most likely candidate gene for QPm.caas-3BS, but TraesCS3B02G016300 and TraesCS3B02G016400 were less likely candidates based on gene annotations and sequence variation between the parents. These results not only offer high-throughput KASP markers for improvement of powdery mildew resistance but also pave the way to map-based cloning of the resistance gene.

Genome-wide association study reveals a genomic region on 5AL for salinity tolerance in wheat

Soil salinity is a major threat to crop productivity and quality worldwide. In order to reduce the negative effects of salinity stress, it is important to understand the genetic basis of salinity tolerance. Identifying new salinity tolerance QTL or genes is crucial for breeders to pyramid different tolerance mechanisms to improve crop adaptability to salinity. Being one of the major cereal crops, wheat is known as a salt-sensitive glycophyte and subject to substantial yield losses when grown in the presence of salt. In this study, both pot and tank experiments were conducted to investigate the genotypic variation present in 328 wheat varieties in their salinity tolerance at the vegetative stage. A Genome-Wide Association Studies (GWAS) were carried out to identify QTL conferring salinity tolerance through a mixed linear model. Six, five and eight significant marker-trait associations (MTAs) were identified from pot experiments, tank experiments and average damage scores, respectively. These markers are located on the wheat chromosomes 1B, 2B, 2D, 3A, 4B, and 5A. These tolerance alleles were additive in their effects and, when combined, increased tolerance to salinity. Candidate genes identified in these QTL regions encoded a diverse class of proteins involved in salinity tolerance in plants. A Na⁺/H⁺ exchanger and a potassium transporter on chromosome 5A (IWB30519) will be of a potential value for improvement of salt tolerance of wheat cultivars using marker assisted selection programs. Some useful genotypes, which showed consistent tolerance in different trials, can also be effectively used in breeding programs.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Genotyping Platforms for Genome-Wide Association Studies: Options and Practical Considerations

Genome-wide association studies (GWAS) in crops requires genotyping platforms that are capable of producing accurate high density genotyping data on hundreds of plants in a cost-effective manner. Currently there are multiple commercial platforms available that are being effectively used across crops. These platforms include genotyping arrays such as the Illumina Infinium arrays and the Applied Biosystems Axiom Arrays along with a variety of resequencing methods. These methods are being used to genotype tens of thousands of markers up to millions of markers on GWAS panels. They are being used on crops with simple genomes to crops with very complex, large, polyploid genomes. Depending on the crop and the goal of the GWAS, there are several options and practical considerations to take into account when selecting a genotyping technology to ensure that the right coverage, accuracy, and cost for the study is achieved.

Mining for New Sources of Resistance to Powdery Mildew in Genetic Resources of Winter Wheat

Genetic pathogen control is an economical and sustainable alternative to the use of chemicals. In order to breed resistant varieties, information about potentially unused genetic resistance mechanisms is of high value. We phenotyped 8,316 genotypes of the winter wheat collection of the German Federal ex situ gene bank for Agricultural and Horticultural Crops, Germany, for resistance to powdery mildew (PM), Blumeria graminis f. sp. tritici, one of the most important biotrophic pathogens in wheat. To achieve this, we used a semi-automatic phenotyping facility to perform high-throughput detached leaf assays. This data set, combined with genotyping-by-sequencing (GBS) marker data, was used to perform a genome-wide association study (GWAS). Alleles of significantly associated markers were compared with SNP profiles of 171 widely grown wheat varieties in Germany to identify currently unexploited resistance conferring genes. We also used the Chinese Spring reference genome annotation and various domain prediction algorithms to perform a domain enrichment analysis and produced a list of candidate genes for further investigation. We identified 51 significantly associated regions. In most of these, the susceptible allele was fixed in the tested commonly grown wheat varieties. Eleven of these were located on chromosomes for which no resistance conferring genes have been previously reported. In addition to enrichment of leucine-rich repeats (LRR), we saw enrichment of several domain types so far not reported as relevant to PM resistance, thus, indicating potentially novel candidate genes for the disease resistance research and prebreeding in wheat.

Mapping Powdery Mildew (Blumeria graminis f. sp. tritici) Resistance in Wild and Cultivated Tetraploid Wheats

Wheat is the most widely grown crop and represents the staple food for one third of the world’s population. Wheat is attacked by a large variety of pathogens and the use of resistant cultivars is an effective and environmentally safe strategy for controlling diseases and eliminating the use of fungicides. In this study, a collection of wild and cultivated tetraploid wheats (Triticum turgidum) were evaluated for seedling resistance (SR) and adult plant resistance (APR) to powdery mildew (Blumeria graminis) and genotyped with a 90K single nucleotide polymorphism (SNP) array to identify new sources of resistance genes. The genome-wide association mapping detected 18 quantitative trait loci (QTL) for APR and 8 QTL for SR, four of which were identical or at least closely linked to four QTL for APR. Thirteen candidate genes, containing nucleotide binding sites and leucine-rich repeats, were localized in the confidence intervals of the QTL-tagging SNPs. The marker IWB6155, associated to QPm.mgb-1AS, was located within the gene TRITD1Av1G004560 coding for a disease resistance protein. While most of the identified QTL were described previously, five QTL for APR (QPm.mgb-1AS, QPm.mgb-2BS, QPm.mgb-3BL.1, QPm.mgb-4BL, QPm.mgb-7BS.1) and three QTL for SR (QPm.mgb-3BL.3, QPm.mgb-5AL.2, QPm.mgb-7BS.2) were mapped on chromosome regions where no resistance gene was reported before. The novel QTL/genes can contribute to enriching the resistance sources available to breeders.

Appraising the Genetic Architecture of Kernel Traits in Hexaploid Wheat Using GWAS

Kernel morphology is one of the major yield traits of wheat, the genetic architecture of which is always important in crop breeding. In this study, we performed a genome-wide association study (GWAS) to appraise the genetic architecture of the kernel traits of 319 wheat accessions using 22,905 single nucleotide polymorphism (SNP) markers from a wheat 90K SNP array. As a result, 111 and 104 significant SNPs for Kernel traits were detected using four multi-locus GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, and pLARmEB) and three single-locus models (FarmCPU, MLM, and MLMM), respectively. Among the 111 SNPs detected by the multi-locus models, 24 SNPs were simultaneously detected across multiple models, including seven for kernel length, six for kernel width, six for kernels per spike, and five for thousand kernel weight. Interestingly, the five most stable SNPs (RAC875_29540_391, Kukri_07961_503, tplb0034e07_1581, BS00074341_51, and BobWhite_049_3064) were simultaneously detected by at least three multi-locus models. Integrating these newly developed multi-locus GWAS models to unravel the genetic architecture of kernel traits, the mrMLM approach detected the maximum number of SNPs. Furthermore, a total of 41 putative candidate genes were predicted to likely be involved in the genetic architecture underlining kernel traits. These findings can facilitate a better understanding of the complex genetic mechanisms of kernel traits and may lead to the genetic improvement of grain yield in wheat.