Multiple-Race Stem Rust Resistance Loci Identified in Durum Wheat Using Genome-Wide Association Mapping

Mapping and Candidate Gene Analysis of an All-Stage Stem Rust Resistance Gene in Durum Wheat Landrace PI 94701

Puccinia graminis f. sp. tritici (Pgt), the causal agent of wheat stem rust, poses a significant threat to global wheat production. Genetic resistance offers a cost-effective and sustainable solution. The durum wheat landrace PI 94701 was previously hypothesized to carry two stem rust resistance (Sr) genes, but their chromosomal locations were unknown. In this study, we mapped and characterized an all-stage Sr gene in PI 94701, temporarily designated as SrPI94701. In seedling tests, SrPI94701 was effective against all six Pgt races tested. Using a large segregating population, we mapped SrPI94701 on chromosome arm 5BL within a 0.17-cM region flanked by markers pku69124 and pku69228, corresponding to 1.04 and 2.15 Mb genomic regions in the Svevo and Chinese Spring reference genomes. Within the candidate region, eight genes exhibited differential expression between the Pgt-inoculated resistant and susceptible plants. Among them, two nucleotide-binding leucine-rich repeat (NLR) genes, TraesCS5B03G1334700 and TraesCS5B03G1335100, showed high polymorphism between the parental lines and were upregulated in Pgt-inoculated resistant plants. However, the flanking and completely linked markers developed in this study could not accurately predict the presence of SrPI94701 in a survey of 104 wheat accessions. SrPI94701 is a promising resource for enhancing stem rust resistance in wheat breeding programs.

Genotype × environment interaction: trade-offs between the agronomic performance and stability of durum (Triticum turgidum) wheat to stem-rust resistance in Kenya

Stem rust significantly threatens durum wheat production, often resulting in substantial yield losses. To better understand resistance mechanisms and the stability of durum lines in stem rust-prone environments, this study evaluated 49 durum genotypes over three seasons at the Kenya Agricultural and Livestock Research Organization in Njoro. Utilizing 7 × 7 alpha lattice design, we assessed adult-plant resistance, monitored disease progression through final disease score (FDS) and area under the disease progress curve (AUDPC), and evaluated agronomic performance. Statistical analyses revealed significant seasonal and genotypic effects on FDS, AUDPC, spike length, and grain yield (p≤0.01; p≤0.001), with important genotype-by-season interactions (p≤0.05; p≤0.001). Broad-sense heritability for AUDPC was high at 0.91 and moderate at 0.35 for kernels per spike, underscoring the genetic basis of these traits. Notably, we observed negative correlations between disease parameters and agronomic traits, suggesting potential trade-offs. GGE biplot analysis singled out the first season (main season of 2019) as crucial for evaluating stem rust resistance and identified several durum lines, such as G45 and G48, as consistently resistant across all conditions. Furthermore, this analysis highlighted G45, G48, G176 and G189 as the highest yielding and most stable lines. The discovery of these resistant and high-performing genotypes is critical for enhancing durum breeding programs, helping to mitigate the impact of stem rust and improve yield stability.

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.

A new wild emmer wheat panel allows to map new loci associated with resistance to stem rust at seedling stage

Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is a major wheat disease worldwide. A collection of 283 wild emmer wheat [Triticum turgidum L. subsp. dicoccoides (Körn. ex Asch. & Graebn.) Thell] accessions, representative of the entire Fertile Crescent region where wild emmer naturally occurs, was assembled, genotyped, and characterized for population structure, genetic diversity, and rate of linkage disequilibrium (LD) decay. Then, the collection was employed for mapping Pgt resistance genes, as a proof of concept of the effectiveness of genome-wide association studies in wild emmer. The collection was evaluated in controlled conditions for reaction to six common Pgt pathotypes (TPMKC, TTTTF, JRCQC, TRTTF, TTKSK/Ug99, and TKTTF). Most resistant accessions originated from the Southern Levant wild emmer lineage, with some showing a resistance reaction toward three to six tested races. Association analysis was conducted considering a 12K polymorphic single-nucleotide polymorphisms dataset, kinship relatedness between accessions, and population structure. Eleven significant marker-trait associations (MTA) were identified across the genome, which explained from 17% to up to 49% of phenotypic variance with an average 1.5 additive effect (based on the 1-9 scoring scale). The identified loci were either effective against single or multiple races. Some MTAs colocalized with known Pgt resistance genes, while others represent novel resistance loci useful for durum and bread wheat prebreeding. Candidate genes with an annotated function related to plant response to pathogens were identified at the regions linked to the resistance and defined according to the estimated small LD (about 126 kb), as typical of wild species.

Improving power of genome-wide association studies via transforming ordinal phenotypes into continuous phenotypes

Introduction

Ordinal traits are important complex traits in crops, while genome-wide association study (GWAS) is a widely-used method in their gene mining. Presently, GWAS of continuous quantitative traits (C-GWAS) and single-locus association analysis method of ordinal traits are the main methods used for ordinal traits. However, the detection power of these two methods is low.

Methods

To address this issue, we proposed a new method, named MTOTC, in which hierarchical data of ordinal traits are transformed into continuous phenotypic data (CPData).

Results

Then, FASTmrMLM, one C-GWAS method, was used to conduct GWAS for CPData. The results from the simulation studies showed that, MTOTC+FASTmrMLM for ordinal traits was better than the classical methods when there were four and fewer hierarchical levels. In addition, when MTOTC was combined with FASTmrEMMA, mrMLM, ISIS EM-BLASSO, pLARmEB, and pKWmEB, relatively high power and low false positive rate in QTN detection were observed as well. Subsequently, MTOTC was applied to analyze the hierarchical data of soybean salt-alkali tolerance. It was revealed that more significant QTNs were detected when MTOTC was combined with any of the above six C-GWAs.

Discussion

Accordingly, the new method increases the choices of the GWAS methods for ordinal traits and helps to mine the genes for ordinal traits in resource populations.

Genetic diversity of Ethiopian durum wheat (T. turgidum subsp. durum) landraces under water stressed and non stressed conditions

Ethiopia, being a major center of origin and diversity for durum wheat, possesses a highly variable genetic pool with diverse agroecological adaptations. Wheat landraces are an important source of genetic variation for breeding programs. This study was conducted to study the genotypic diversity of Ethiopian durum wheat genetic resources under two contrasting environments namely drought-stressed and non-stressed. It was carried out on 100 landraces and 4 local checks using an augmented design. Data were collected on 13 traits comprising yield and yield components, phenology, and canopy condition. The analysis of variance revealed significant differences between landraces for different traits with different sources of variation. Several landraces were found to outyield the checks at both environmental conditions. Intermediate to high estimates of the phenotypic coefficient of variation (PCV), genotypic coefficient of variation (GCV), heritability in a broad sense (h2b), and genetic advance in percent of the mean (GAPM) were observed for all the studied traits except for days to flowering at normal, thousands seed weight at stress, and days to maturity, leaf chlorophyll concentration measurement, and canopy temperature measurement at both conditions. The estimation of variability parameters showed that genotypic variation was higher than environmental variation for most traits. The number of tillers, spike length, kernel per spike, and grain yield indicated higher values for h2b and GAPM (74.42% and 20.86; 83.2% and 28.24; 70.79% and 28.0; and 89.54% and 74.71) at normal and (97.87% and 98.22; 71.27% and 28.51; 75.52% and 43.9; and 90.04% and 103.68) at the stressed condition, respectively. Spikelets per spike, kernel per spike, and thousands seed weight were positively correlated with grain yield. Grain yield exhibited a weak negative correlation with days to heading and days to maturity. Principal components analysis revealed that six traits were the major loadings on the first two principal components that describe 37.9% and 41.0% of the total morphological variance at normal and stressed conditions, respectively. Cluster analysis grouped the landraces into six clusters, with each cluster showing variation in performance for different traits under normal and stressed conditions. The intracluster distance was maximum in cluster I (D2 = 7.68) and (D2 = 8.19) at normal and stressed conditions respectively and the intercluster distance was found to be maximum between clusters I and IV (D2 = 11.02) and clusters I and II (D2 = 10.33) at normal and stressed conditions respectively. The presence of significant genetic variability among the evaluated durum wheat landraces suggests an opportunity for improvement of grain yield through the hybridization of genotypes from different clusters and subsequent selection. Genotypes with superior agronomic traits that outperform the best checks are identified as potential parents for yield improvement programs for moisture stress.

Wheat Stem Rust Detection and Race Characterization in Tunisia

Climate changes over the past 25 years have led to conducive conditions for invasive and transboundary fungal disease occurrence, including the re-emergence of wheat stem rust disease, caused by Puccinia graminis f.sp. tritici (Pgt) in East Africa, Europe, and the Mediterranean basin. Since 2018, sporadic infections have been observed in Tunisia. In this study, we investigated Pgt occurrence at major Tunisian wheat growing areas. Pgt monitoring, assessment, and sampling from planted trap nurseries at five different locations over two years (2021 and 2022) revealed the predominance of three races, namely TTRTF (Clade III-B), TKKTF (Clade IV-F), and TKTTF (Clade IV-B). Clade III-B was the most prevalent in 2021 as it was detected at all locations, while in 2022 Pgt was only reported at Beja and Jendouba, with the prevalence of Clade IV-B. The low levels of disease incidence during these two years and Pgt population diversity suggest that this fungus most likely originated from exotic incursions and that climate factors could have caused disease establishment in Tunisia. Further evaluation under the artificial disease pressure of Tunisian wheat varieties and weather-based modeling for early disease detection in the Mediterranean area could be helpful in monitoring and predicting wheat stem rust emergence and epidemics.

Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein

Rice blast is a worldwide fungal disease that seriously affects the yield and quality of rice. Identification of resistance genes against rice blast disease is one of the effective ways to control this disease. However, panicle blast resistance genes, which are useful in the fields, have rarely been studied due to the difficulty in phenotypic identification and the environmental influences. Here, panicle blast resistance-3 (Pb3) was identified by a genome-wide association study (GWAS) based on the panicle blast resistance phenotypes of 230 Rice Diversity Panel I (RDP-I) accessions with 700,000 single-nucleotide polymorphism (SNP) markers. A total of 16 panicle blast resistance loci (PBRLs) within three years including one repeated locus PBRL3 located in chromosome 11 were identified. In addition, 7 genes in PBRL3 were identified as candidate genes by haplotype analysis, which showed significant differences between resistant and susceptible varieties. Among them, one nucleotide-binding domain and Leucine-rich Repeat (NLR) gene Pb3 was highly conserved in multiple resistant rice cultivars, and its expression was significantly induced after rice blast inoculation. Evolutionary analysis showed that Pb3 was a typical disease resistance gene containing coiled-coil, NB-ARC, and LRR domains. T-DNA insertion mutants and CRISPR lines of Pb3 showed significantly reduced panicle blast resistance. These results indicate that Pb3 is a panicle blast resistance gene and GWAS is a rapid method for identifying panicle blast resistance in rice.

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.

Meta-QTLs for multiple disease resistance involving three rusts in common wheat (Triticum aestivum L.)

In wheat, multiple disease resistance meta-QTLs (MDR-MQTLs) and underlying candidate genes for the three rusts were identified which may prove useful for development of resistant cultivars. Rust diseases in wheat are a major threat to global food security. Therefore, development of multiple disease-resistant cultivars (resistant to all three rusts) is a major goal in all wheat breeding programs worldwide. In the present study, meta-QTLs and candidate genes for multiple disease resistance (MDR) involving all three rusts were identified using 152 individual QTL mapping studies for resistance to leaf rust (LR), stem rust (SR), and yellow rust (YR). From these 152 studies, a total of 1,146 QTLs for resistance to three rusts were retrieved, which included 368 QTLs for LR, 291 QTLs for SR, and 487 QTLs for YR. Of these 1,146 QTLs, only 718 QTLs could be projected onto the consensus map saturated with 2, 34,619 markers. Meta-analysis of the projected QTLs resulted in the identification of 86 MQTLs, which included 71 MDR-MQTLs. Ten of these MDR-MQTLs were referred to as the 'Breeders' MQTLs'. Seventy-eight of the 86 MQTLs could also be anchored to the physical map of the wheat genome, and 54 MQTLs were validated by marker-trait associations identified during earlier genome-wide association studies. Twenty MQTLs (including 17 MDR-MQTLs) identified in the present study were co-localized with 44 known R genes. In silico expression analysis allowed identification of several differentially expressed candidate genes (DECGs) encoding proteins carrying different domains including the following: NBS-LRR, WRKY domains, F-box domains, sugar transporters, transferases, etc. The introgression of these MDR loci into high-yielding cultivars should prove useful for developing high yielding cultivars with resistance to all the three rusts.

Mapping and Characterization of a Wheat Stem Rust Resistance Gene in Durum Wheat "Kronos"

Wheat stem (or black) rust is one of the most devastating fungal diseases, threatening global wheat production. Identification, mapping, and deployment of effective resistance genes are critical to addressing this challenge. In this study, we mapped and characterized one stem rust resistance (Sr) gene from the tetraploid durum wheat variety Kronos (temporary designation SrKN). This gene was mapped on the long arm of chromosome 2B and confers resistance to multiple virulent Pgt races, such as TRTTF and BCCBC. Using a large mapping population (3,366 gametes), we mapped SrKN within a 0.29 cM region flanked by the sequenced-based markers pku4856F2R2 and pku4917F3R3, which corresponds to 5.6- and 7.2-Mb regions in the Svevo and Chinese Spring reference genomes, respectively. Both regions include a cluster of nucleotide binding leucine-repeat (NLR) genes that likely includes the candidate gene. An allelism test failed to detect recombination between SrKN and the previously mapped Sr9e gene. This result, together with the similar seedling resistance responses and resistance profiles, suggested that SrKN and Sr9e may represent the same gene. We introgressed SrKN into common wheat and developed completely linked markers to accelerate its deployment in the wheat breeding programs. SrKN can be a valuable component of transgenic cassettes or gene pyramids that includes multiple resistance genes to control this devastating disease.

Genome-wide association mapping of seedling and adult plant response to stem rust in a durum wheat panel

Many of the major stem rust resistance genes deployed in commercial wheat (Triticum spp.) cultivars and breeding lines become ineffective over time because of the continuous emergence of virulent races. A genome-wide association study (GWAS) was conducted using 26,439 single nucleotide polymorphism (SNP) markers and 280 durum wheat [Triticum turgidum L. subsp. Durum (Desf.) Husnot] lines from CIMMYT to identify genomic regions associated with seedling resistance to races TTKSK, TKTTF, JRCQC, and TTRTF and field resistance to TKTTF and JRCQC. The phenotypic data analysis across environments revealed 61-91 and 59-77% of phenotypic variation was explained by the genotypic component for seedling and adult plant response of lines, respectively. For seedling resistance, mixed linear model (MLM) identified eight novel and nine previously reported quantitative trait loci (QTL) while a fixed and random model circulating probability unification (FarmCPU) detected 12 novel and eight previously reported QTL. For field resistance, MLM identified 12 novel and seven previously reported loci while FarmCPU identified seven novel and nine previously reported loci. The regions of Sr7a, Sr8155B1, Sr11, alleles of Sr13, Sr17, Sr22/Sr25, and Sr49 were identified. Novel loci on chromosomes 3B, 4A, 6A, 6B, 7A, and 7B could be used as sources of resistance to the races virulent on durum wheat. Two large-effect markers on chromosome 6A could potentially be used to differentiate resistant haplotypes of Sr13 (R1 and R3). Allelism tests for Sr13, breaking the deleterious effect associated with Sr22/Sr25 and retaining the resistance allele at the Sr49 locus, are needed to protect future varieties from emerging races.