Genome-Wide Association Studies to Identify Loci and Candidate Genes Controlling Kernel Weight and Length in a Historical United States Wheat Population

Reviewing the essential roles of remote phenotyping, GWAS and explainable AI in practical marker-assisted selection for drought-tolerant winter wheat breeding

Marker-assisted selection (MAS) plays a crucial role in crop breeding improving the speed and precision of conventional breeding programmes by quickly and reliably identifying and selecting plants with desired traits. However, the efficacy of MAS depends on several prerequisites, with precise phenotyping being a key aspect of any plant breeding programme. Recent advancements in high-throughput remote phenotyping, facilitated by unmanned aerial vehicles coupled to machine learning, offer a non-destructive and efficient alternative to traditional, time-consuming, and labour-intensive methods. Furthermore, MAS relies on knowledge of marker-trait associations, commonly obtained through genome-wide association studies (GWAS), to understand complex traits such as drought tolerance, including yield components and phenology. However, GWAS has limitations that artificial intelligence (AI) has been shown to partially overcome. Additionally, AI and its explainable variants, which ensure transparency and interpretability, are increasingly being used as recognised problem-solving tools throughout the breeding process. Given these rapid technological advancements, this review provides an overview of state-of-the-art methods and processes underlying each MAS, from phenotyping, genotyping and association analyses to the integration of explainable AI along the entire workflow. In this context, we specifically address the challenges and importance of breeding winter wheat for greater drought tolerance with stable yields, as regional droughts during critical developmental stages pose a threat to winter wheat production. Finally, we explore the transition from scientific progress to practical implementation and discuss ways to bridge the gap between cutting-edge developments and breeders, expediting MAS-based winter wheat breeding for drought tolerance.

GWAS in a Collection of Bulgarian Old and Modern Bread Wheat Accessions Uncovers Novel Genomic Loci for Grain Protein Content and Thousand Kernel Weight

Genetic enhancement of grain production and quality is a priority in wheat breeding projects. In this study, we assessed two key agronomic traits-grain protein content (GPC) and thousand kernel weight (TKW)-across 179 Bulgarian contemporary and historic varieties and landraces across three growing seasons. Significant phenotypic variation existed for both traits among genotypes and seasons, and no discernible difference was evident between the old and modern accessions. To understand the genetic basis of the traits, we conducted a genome-wide association study with MLM using phenotypic data from the crop seasons, best linear unbiased estimators, and genotypic data from the 25K Infinium iSelect array. As a result, we detected 16 quantitative trait nucleotides (QTNs) associated with GPC and 15 associated with TKW, all of which passed the false discovery rate threshold. Seven loci favorably influenced GPC, resulting in an increase of 1.4% to 8.1%, while four loci had a positive impact on TKW with increases ranging from 1.9% to 8.4%. While some loci confirmed previously published associations, four QTNs linked to GPC on chromosomes 2A, 7A, and 7B, as well as two QTNs related to TKW on chromosomes 1B and 6A, may represent novel associations. Annotations for proteins involved in the senescence-associated nutrient remobilization and in the following buildup of resources required for seed germination have been found for selected putative candidate genes. These include genes coding for storage proteins, cysteine proteases, cellulose-synthase, alpha-amylase, transcriptional regulators, and F-box and RWP-RK family proteins. Our findings highlight promising genomic regions for targeted breeding programs aimed at improving grain yield and protein content.

Identification and functional analysis of a chromosome 2D fragment harboring TaFPF1 gene with the potential for yield improvement using a late heading wheat mutant

KEY MESSAGE

A chromosome fragment influencing wheat heading and grain size was identified using mapping of m406 mutant. The study of TaFPF1 in this fragment provides more insights into wheat yield improvement. In recent years, wheat production has faced formidable challenges driven by rapid population growth and climate change, emphasizing the importance of improving specific agronomic traits such as heading date, spike length, and grain size. To identify potential genes for improving these traits, we screened a wheat EMS mutant library and identified a mutant, designated m406, which exhibited a significantly delayed heading date compared to the wild-type. Intriguingly, the mutant also displayed significantly longer spike and larger grain size. Genetic analysis revealed that a single recessive gene was responsible for the delayed heading. Surprisingly, a large 46.58 Mb deletion at the terminal region of chromosome arm 2DS in the mutant was identified through fine mapping and fluorescence in situ hybridization. Thus, the phenotypes of the mutant m406 are controlled by a group of linked genes. This deletion encompassed 917 annotated high-confidence genes, including the previously studied wheat genes Ppd1 and TaDA1, which could affect heading date and grain size. Multiple genes in this region probably contribute to the phenotypes of m406. We further investigated the function of TaFPF1 using gene editing. TaFPF1 knockout mutants showed delayed heading and increased grain size. Moreover, we identified the direct upstream gene of TaFPF1 and investigated its relationship with other important flowering genes. Our study not only identified more genes affecting heading and grain development within this deleted region but also highlighted the potential of combining these genes for improvement of wheat traits.

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 Analysis for Hybrid Breeding in Wheat

Disclosure of markers that are significantly associated with plant traits can help develop new varieties with desirable properties. This study determined the genome-wide associations based on DArTseq markers for six agronomic traits assessed in eight environments for wheat. Moreover, the association study for heterosis and analysis of the effects of markers grouped by linkage disequilibrium were performed based on mean values over all experiments. All results were validated using data from post-registration trials. GWAS revealed 1273 single nucleotide polymorphisms with biologically significant effects. Most polymorphisms were predicted to be modifiers of protein translation, with only two having a more pronounced effect. Markers significantly associated with the considered set of features were clustered within chromosomes based on linkage disequilibrium in 327 LD blocks. A GWAS for heterosis revealed 1261 markers with significant effects.

Genome-wide association mapping of bread wheat genotypes for sustainable food security and yield potential under limited water conditions

Determining the genetic basis of yield and water deficient tolerance in wheat is vital for wheat breeding programs. Herein, a genome-wide association study (GWAS) was performed for water deficient and yield-related attributes on wheat genotypes with high-density Illumina 90K Infinium SNP array. Major yield and drought-related attributes were phenotyped on a panel of Pakistani and foreign accessions grown in non-stressed and water deficient stressed environments during two crop cycles. Among all accessions, highly significant variations were shown in studied environments for examined characters. Water deficient conditions, reduced the wheat yield and had strong and positive correlation among relative water content and grain yield per plant. Population structure analyses based on 90,000 SNP data, classify the accessions into 4 sub-populations. Marker-trait association analyses (MTA) revealed that 134 significant SNPs were linked with yield and drought tolerance attributes. Pleotropic loci RAC875_s117925_244 and RAC875_c16333_340 located on chromosome 5A and 2A respectively, were significantly linked with relative water contents (RWC), cell membrane thermo-stability (CMT), grain per spike (GPS), spikelet per spike (SPS) and grain yield per plant (GYP). The markers Ra_c58279_684, BobWhite_c23828_341 and IAAV3414 located on chromosomes 2A, 6B and 7B respectively, showed pleotropic effects for RWC, GPS and GYP under both environments. The current experiment not only validated several MTAs reported in other studies but also discovered novel MTAs which significant under drought-stressed conditions. A total of 171 candidate genes were recognized that could be cloned and functionally characterized for the respective associated traits. For RWC and CMT, total 11 and 3 associated SNPs were mapped on coding DNA sequence (CDS) of the identified candidate genes. Isolation and characterization of the candidate genes herein mapped SNPs will be useful in discovering novel genes underpinning drought tolerance in bread wheat to fulfill the wheat demand and sustainable food security under limited water conditions.

Overexpression of TaSTT3b-2B improves resistance to sharp eyespot and increases grain weight in wheat

STAUROSPORINE AND TEMPERATURE SENSITIVE3 (STT3) is a catalytic subunit of oligosaccharyltransferase, which is important for asparagine-linked glycosylation. Sharp eyespot, caused by the necrotrophic fungal pathogen Rhizoctonia cerealis, is a devastating disease of bread wheat. However, the molecular mechanisms underlying wheat defense against R. cerealis are still largely unclear. In this study, we identified TaSTT3a and TaSTT3b, two STT3 subunit genes from wheat and reported their functional roles in wheat defense against R. cerealis and increasing grain weight. The transcript abundance of TaSTT3b-2B was associated with the degree of wheat resistance to R. cerealis and induced by both R. cerealis and exogenous jasmonic acid (JA). Overexpression of TaSTT3b-2B significantly enhanced resistance to R. cerealis, grain weight, and JA content in transgenic wheat subjected to R. cerealis stress, while silencing of TaSTT3b-2B compromised resistance of wheat to R. cerealis. Transcriptomic analysis showed that TaSTT3b-2B affected the expression of a series of defense-related genes and JA biosynthesis-related genes, as well as genes coding starch synthase and sucrose synthase. Application of exogenous JA elevated expression levels of the abovementioned defense- and grain weight-related genes, and rescuing the resistance of TaSTT3b-2B-silenced wheat to R. cerealis, while pretreatment with sodium diethyldithiocarbamate, an inhibitor of JA synthesis, attenuated the TaSTT3b-2B-mediated resistance to R. cerealis, suggesting that TaSTT3b-2B played critical roles in regulating R. cerealis resistance and grain weight via JA biosynthesis. Altogether, this study reveals new functional roles of TaSTT3b-2B in regulating plant innate immunity and grain weight, and illustrates its potential application value for wheat molecular breeding.

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.

Mapping QTL for spike fertility and related traits in two doubled haploid wheat (Triticum aestivum L.) populations

BACKGROUND

In breeding programs, the selection of cultivars with the highest yield potential consisted in the selection of the yield per se, which resulted in cultivars with higher grains per spike (GN) and occasionally increased grain weight (GW) (main numerical components of the yield). In this study, quantitative trait loci (QTL) for GW, GN and spike fertility traits related to GN determination were mapped using two doubled haploid (DH) populations (Baguette Premium 11 × BioINTA 2002 and Baguette 19 × BioINTA 2002).

RESULTS

In total 305 QTL were identified for 14 traits, out of which 12 QTL were identified in more than three environments and explained more than 10% of the phenotypic variation in at least one environment. Eight hotspot regions were detected on chromosomes 1A, 2B, 3A, 5A, 5B, 7A and 7B in which at least two major and stable QTL sheared confidence intervals. QTL on two of these regions (R5A.1 and R5A.2) have previously been described, but the other six regions are novel.

CONCLUSIONS

Based on the pleiotropic analysis within a robust physiological model we conclude that two hotspot genomic regions (R5A.1 and R5A.2) together with the QGW.perg-6B are of high relevance to be used in marker assisted selection in order to improve the spike yield potential. All the QTL identified for the spike related traits are the first step to search for their candidate genes, which will allow their better manipulation in the future.

Exploitation of Tolerance of Wheat Kernel Weight and Shape-Related Traits from Aegilops tauschii under Heat and Combined Heat-Drought Stresses

Kernel weight and shape-related traits are inherited stably and increase wheat yield. Narrow genetic diversity limits the progress of wheat breeding. Here, we evaluated kernel weight and shape-related traits and applied genome-wide association analysis to a panel of wheat multiple synthetic derivative (MSD) lines. The MSD lines harbored genomic fragments from Aegilops tauschii. These materials were grown under optimum conditions in Japan, as well as under heat and combined heat-drought conditions in Sudan. We aimed to explore useful QTLs for kernel weight and shape-related traits under stress conditions. These can be useful for enhancing yield under stress conditions. MSD lines possessed remarkable genetic variation for all traits under all conditions, and some lines showed better performance than the background parent Norin 61. We identified 82 marker trait associations (MTAs) under the three conditions; most of them originated from the D genome. All of the favorable alleles originated from Ae. tauschii. For the first time, we identified markers on chromosome 5D associated with a candidate gene encoding a RING-type E3 ubiquitin-protein ligase and expected to have a role in regulating wheat seed size. Our study provides important knowledge for the improvement of wheat yield under optimum and stress conditions. The results emphasize the importance of Ae. tauschii as a gene reservoir for wheat breeding.

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.

Identification of New Leaf Rust Resistance Loci in Wheat and Wild Relatives by Array-Based SNP Genotyping and Association Genetics

Leaf rust caused by Puccinia triticina is the most widespread rust disease of wheat. As pathogen populations are constantly evolving, identification of novel sources of resistance is necessary to maintain disease resistance and stay ahead of this plant-pathogen evolutionary arms race. The wild genepool of wheat is a rich source of genetic diversity, accounting for 44% of the Lr genes identified. Here we performed a genome-wide association study (GWAS) on a diverse germplasm of 385 accessions, including 27 different Triticum and Aegilops species. Genetic characterization using the wheat 90 K array and subsequent filtering identified a set of 20,501 single nucleotide polymorphic (SNP) markers. Of those, 9,570 were validated using exome capture and mapped onto the Chinese Spring reference sequence v1.0. Phylogenetic analyses illustrated four major clades, clearly separating the wild species from the T. aestivum and T. turgidum species. GWAS was conducted using eight statistical models for infection types against six leaf rust isolates and leaf rust severity rated in field trials for 3-4 years at 2-3 locations in Canada. Functional annotation of genes containing significant quantitative trait nucleotides (QTNs) identified 96 disease-related loci associated with leaf rust resistance. A total of 21 QTNs were in haplotype blocks or within flanking markers of at least 16 known Lr genes. The remaining significant QTNs were considered loci that putatively harbor new Lr resistance genes. Isolation of these candidate genes will contribute to the elucidation of their role in leaf rust resistance and promote their usefulness in marker-assisted selection and introgression.