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.

Recurrent stem rust epidemics have resulted in large-scale production losses (Singh et al., 2015), posing serious and imminent threat to food and nutrition security in Kenya and other durum wheat production regions.While durum adoption in Kenya lags behind Ethiopia, food insecurity urges diversification with durum as key cereal because of increasing consumption of durum products.The emergence of race TTKSK (Ug99) in Africa has rendered the widely used stem rust resistance gene Sr31 (Zhang et al., 2017) and multiple additional Sr genes ineffective (Pretorius et al., 2000).It was anticipated that this race would continue to spread to Iran, South Africa, Yemen, and East Africa (Singh et al., 2015;Patpour et al., 2020).
The global threat posed by Ug99 stem rust strain and its variants, which have compromised previously effective resistance genes such as Sr31, highlights the critical need for ongoing genetic innovation (Pretorius et al., 2000;Singh et al., 2015).Over 13 virulent variants of Ug99 have been identified, exposing up to 95% of global wheat cultivars to stem rust (Newcomb et al., 2016;Mergessa et al., 2020;Bhavani et al., 2022).The rapid evolution of these pathogens emphasizes the importance of a proactive approach in breeding programs that incorporates a broader base of resistance genes from both tetraploid and hexaploid wheat varieties (McIntosh et al., 1995;Park, 2015).
Novel resistance genes are being incorporated into durum wheat breeding programs to address global imperative, which goes beyond regional concerns.Since Puccinia graminis f. sp.tritici (Pgt) does not consider national borders, stem rust poses persistent threat to food security worldwide.Reducing the reliance on fungicides globally is essential for lowering production costs and boosting the economic resilience of farming communities across the globe.This can be achieved through strategic integration of these resistance genes into breeding programs (Singh et al., 2011;Park, 2008).We focused on comprehensive phenotypic analysis of durum wheat genotypes in order to discover and define these novel sources of resistance.This method offers useful foundation for selection in breeding programs by enabling the direct observation of resistance under various environmental conditions (Chao et al., 2017;Bhavani et al., 2022).
Our focus was to identify durum wheat cultivars that demonstrate durable resistance to both current and emerging Pgt races by fusing conventional breeding techniques with cutting-edge genomic technologies.In light of the danger posed by emerging infections, such initiatives are crucial for maintaining wheat output and enhancing yield stability (Chaves et al., 2013;Singh et al., 2016).This approach is in line with sustainable farming techniques, which are essential for ensuring future food security, and it also lessens the effects of deadly diseases like stem rust (Ellis et al., 2014;Nirmala et al., 2017).
Collaboratively, this effort is being carried out collaboratively with international organizations such as CIMMYT, the Kenya Agriculture and Livestock Research Organization, and the Ethiopia Institute of Agricultural Research.These collaborations have resulted in over 500 new durum wheat varieties adapted to wide range of environmental conditions, demonstrating the power of collaborative research and development (Macharia and Ngina, 2017;Bhavani et al., 2019).Furthermore, our understanding of genotype-environment interactions has revolutionized with the application of complex statistical models such as the Additive Main Effects and Multiplicative Interaction (AMMI) and the Genotype Main Effect (Yan and Tinker, 2006) and Genotype by Environment Interaction (GGE) biplots (Gauch, 2006(Gauch, , 2013)).These models have improved our capacity to choose genotypes that reliably yield high and stable yields under various climatic conditions (Mohammadi and Amri, 2013;Heidari et al., 2016;Kendal et al., 2019).
In addition to providing important new insights into the genetic underpinnings of stem rust resistance, this research aimed to support the sustainability of durum wheat production worldwide, promising food security in the face of changing biotic threats.Through series of comprehensive field trials conducted in Njoro, Kenya, and analysis of CIMMYT's large durum germplasm pool, we hoped to discover new genetic resources that would strengthen the regional and global's durum wheat supply in the face of stem rust's persistent challenge.

Experimental site
The experiment was conducted at Kenya Agricultural and Livestock Research Organization (KALRO), Njoro, for three seasons, i.e off-season 2019 (Season 1), main-season 2019 (Season 2) and off-season 2020 (Season 3).The site is located at 35°55′ 60″ E, 0°19′ 60″ S with an elevation of 2185 m above sea level in lower highland agro-ecological zone III (LH3) with predominant welldrained mollic andosols soils.The specific location was at KALRO wheat rust phenotyping platform, which is situated at 0.20513°S latitude; 35.56801°E longitude and 2171m elevation, as shown in Figure 1.The mean annual precipitation is 1000mm with minimum and maximum temperature of 9°C and 22°C respectively (KALRO Meteorological Station No. 903502) (Shewry, 2009;Shewry and Hey, 2015).The weather conditions experienced during the growing period are as shown in Figure 2.

Plant materials
A set of 49 durum genotypes derived from CIMMYT collection lines along with stem rust susceptible checks (Cacuke and PBW343) were evaluated as represented in Table 1.They were selected based on stem rust severity, morphology (height, spike and uniformity in particular), grain yield and biomass production.

Experimental procedure
This study was conducted in a field where canola (Brassica napus) was previously grown.A suitable seedbed was prepared using disc plough and harrow operations.Each line was sown in 2-row plots of 0.2 m × 0.75 m.The plots were arranged in 7 × 7 alpha-lattice design with 0.2 m spacing between rows and 0.5 m spacing between blocks.This allowed control of variability within blocks and replications, enhancing the reliability of comparisons among genotypes across environments/seasons. Three replicates were used.At the time of sowing, di-ammonium phosphate (DAP) fertilizer was applied at a rate of 125 kg ha-1 to provide 22.5 kg N and 25.1 kg P. Spreader rows containing susceptible wheat genotypes and Sr24-carrying lines were sown between replicates, creating quad rows around the experimental unit.Sr24-carrying lines included served as a reference for comparative analysis of genotype susceptibility, and buffering each unit to standardize disease exposure, thus ensuring vigorous assessment of resistance under field conditions.To initiate disease infection, the spreader rows were artificially inoculated with mixture of stem rust urediniospores of different races (TTKSK, TTKST, TTKTT, TTKTK, and TTTTF).Urediniospores harvested from the disease nursery at KALRO-Njoro were suspended in distilled water containing drop of Tween 20 at a concentration of 1×10 5 spores per ml and were inoculated into the spreader rows using a syringe.The inoculation process was repeated multiple times.Additional, urediniospores were sprayed on the spreader rows to enhance disease development.Overhead irrigation was employed in off-season nurseries to maintain a humid environment conducive to disease development, operating for four hours each during the morning and evening.Preemergence herbicide Stomp 455C (pendimethalin) was applied after planting, while post-emergence herbicide Buctril MC (bromoxynil octanoate) and MCPA ethyl-hexyl ester were applied at specific growth stages.Calcium ammonium nitrate (CAN) was top-dressed twice at tillering and stem elongation growth stage to supply nitrogen.Systemic insecticide Thunder OD 145 containing imidacloprid and beta-cyfluthrin was applied at the tillering and ear emergence stages.Foliar diseases were not controlled in this trial.Sprinkler and drip irrigation systems were installed to provide additional moisture during dry periods.

Data collection
The severity of stem rust disease was assessed using the modified Cobb's scale, which involved determining the percentage of stem area covered by rust pustules (Peterson et al., 1948).Infection response was scored based on pustule size, chlorosis, and necrosis on the stem (Roelfs, 1992).Data collection began when approximately 50% of the test genotypes had headed and susceptible checks displayed 50% disease severity.Stem rust severity was recorded at 7-8-day intervals throughout the study.Response classes included "0" for no visible infection, "R" for resistance, "MR" for moderate resistance, "MS" for moderate susceptibility, and "S" for susceptibility (Peterson et al., 1948).In cases where there was an overlapping infection response on a single genotype, the most frequent response was recorded first, followed by the less frequent one.Stem rust severity was scored three times per season at 7-day intervals.The area under disease progress curve (AUDPC) was calculated to quantify the disease data, and the final severity score   was also used in the analysis (Gregory and Finney, 1977).
Where y i =disease observation (severity) at the i th genotype, t i =time (days) at the i th disease observation, n=total number of observations.
Phenological traits, such as plant height, spike length, seeds per spike, 1000-kernel weight, and grain yield, were measured for all durum genotypes.Heading was determined when 50% of the spike emerged from the boot.At maturity, plant height was measured from the soil level to the tip of the spikes, excluding awns, in a random sample of 3 plants per plot.Spike length was measured from base to tip, excluding awns, in a random sample of three spikes per plot.The number of seeds per spike was calculated using three random spikes per plot.Grain yield was estimated by harvesting the plots at the base and recorded in kilograms per plot area (Kg m −2 ), then converted to t ha −1 as: where; GY=grain yield and PA=plot area.Using a Contador seed counter (brand Pfeuffer, Serial number: 14176107), 1,000 kernels were counted from threshed grains and weighed to estimate 1000-Kernel weight.

Statistical analysis
Statistical analyses were performed to evaluate the significance of seasonal and genotypic effects on disease severity and agronomic traits.The genotype by environment (GxE) interactions were analyzed using GEA-R (Genotype x Environment Analysis with R), which allows for the visualization and interpretation of GxE interactions (Pacheco-Gil et al., 2015).Broad-sense heritability estimates were calculated to determine the genetic influence on these traits.Negative correlations between disease parameters and agronomic traits were identified to explore potential trade-offs between resistance and agronomic performance.

Data analysis
The disease and agronomic data underwent two-step combined analysis of variance.First, individual seasons were analyzed, and then combined analysis was performed using 7 × 7 alpha-lattice design.The general linear model procedure (GLM) of SAS version 9.4 (SAS Institute, N.C, 2009) was used for the analysis.
where Y ijkl =observation of the experimental units; µ =overall mean; S i =effects due to i th season; R j(i) =effects due to the j th replicate nested in the i th season; B k(ji) =effects due to the k th block nested in the j th replicate and the i th season; G l = effects due to l th genotype in the i th season; GS li = effects due to the interaction of the j th genotype and i th season and e ijklm = residual.
Tukey's honestly significant difference (Tukey, 1949) was used to compare means of the agronomic and disease traits.The mixed model procedure of SAS was employed to analyze both the disease and agronomic data, allowing for correction of experimental effects such as seasons, replicates, and blocks variability.The mixed model estimated the random effects of experiment, replicates, blocks, and seasons variables, which were then removed from the trait values of each genotype.This process resulted in Best Linear Unbiased Predictors (BLUPs) or adjusted means (LSMeans) for each genotype.The mean BLUPs from all seasons were used for subsequent analyses, as they provide estimates of genotype effects that closely align with true genotypic effects.Broad-sense heritability was estimated based on genotype mean values.
where, s 2 g =genotypic variance, r=number of replicates, s=number of seasons, s 2 gs =genotype × season interaction variance and s 2 e =residual variance.
The GGE biplot, utilizing principal components scores from singular value partitioning (SVP), was employed to depict the impact of genotype and genotype by environment interaction on stem rust severity using GEA-R Software (Angela et al., 2015).The average environment coordination (AEA) from the biplots assisted in identifying the top-performing genotypes and selecting stable and high-performing genotypes.Polygons were utilized to determine the best-performing genotypes in each season.The GGE model is given by; where, Y ij =the mean performance of i t h genotype in j th environment, m= the grand mean, E j = main effect of j th environment, k=multiplicative term for genotype by environment interaction, b ik =genotypic score, z jk =environmental score, e ij =error term.
Correlations among the traits evaluated were computed from the BLUPs to establish the association among all the agronomic and disease traits.The correlation coefficients were further used to determine direct and indirect effects of each disease trait to agronomic traits.The principal component analysis (PCA) was conducted using GENSTAT (Version 16), correlation matrix to determine the relationship between durum traits collected in this study.The PCA allows explanation of the effect of multiple traits that are acting in almost similar ways.

Combined analysis of variance
The combined analysis of variance across seasons demonstrated significant effects of seasons on stem rust severity (p ≤ 0.001) and area under disease progress curve (AUDPC) (p ≤ 0.01), as well as on various agronomic traits (Table 2).Genotypes also exhibited significant effects (p ≤ 0.001) on both disease parameters and agronomic traits (Table 2).Furthermore, significant interaction (p ≤ 0.001) between genotype and season was observed for stem rust severity, AUDPC, and spike length, indicating genotypic performance across different seasons (Table 2).
The three test seasons significantly influenced the performance of genotypes in terms of agronomic traits and the occurrence and progression of the disease.Season 1 exhibited higher stem rust severity compared to seasons 2 and 3 (Table 3), with disease progression also more rapid in Season 1 than in other seasons.The higher disease pressure in Season 1, likely due to more favorable environmental conditions for pathogen development such as higher humidity, resulted in stressed plants reallocating resources towards defense rather than growth, thereby limiting spike length (SPL).
Season 3 resulted in the longest spikes, followed by season 2 and season 1. Plant growth was similar in seasons 1 and 3, but the growth rate was limited in season 1 (Table 3), suggesting that the environmental stresses affecting spike development did not similarly impact overall plant growth rate.Season 3 favored kernels per spike (KPS) and yield, which were comparable to those in season 1 and season 2. Additionally, the thousand kernel weight (TKW) increased from season 2 to season 3.Although disease incidence was lower in Season 2, kernel weight remained reduced, likely affected by suboptimal abiotic conditions such as transient drought or heat stress during the crucial grain-filling period, which adversely impacted kernel development.

Broad-sense heritability (H 2 ) estimates
Broad sense heritability varied across traits, ranging from 0.90 for AUDPC to 0.35 for KPS, indicating the extent to which genetic factors influenced disease and agronomic traits over environmental factors (Table 4).Traits such as final disease score, AUDPC, spike length, plant height, kernel weight, and yield showed high heritability estimates (Table 4).However, KPS exhibited the lowest heritability estimate of 0.35, suggesting stronger influence of environmental effects on this trait.The environmental variance for KPS was nine times greater than the genotypic variance, while the residual variance was six times greater than the genotypic variance (Table 4).

Correlation analysis among agronomic and disease traits
At the phenotypic level, both stem rust severity and AUDPC exhibited significant negative correlations with all agronomic traits (Table 5).The significant negative correlations suggest an association between stem rust severity, AUDPC, and agronomic traits in the overall phenotype.Nevertheless, it is worth mentioning that the phenotypic correlation between spike length (rp=-0.23),plant height (rp=-0.26),and AUDPC was not statistically significant (Table 5).This implies that spike length and plant height may not have a substantial direct impact on the severity or progression of stem rust.

Effects of genotypes, environments, and genotype × environment interactions of stem rust severity and yield
The GGE biplots for stem rust severity (Figure 3) and yield (Figure 4) provided comprehensive analysis of durum wheat genotypes' performance across different seasons.The disease biplot explained 89.80% of the variation, with PC1 accounting for 76.54% and PC2 for 13.26%, while the yield biplot explained 87.80% of the variation, with PC1 at 69.95% and PC2 at 17.85%.Genotype 45 stood out in both biplots, showing strong resistance to stem rust, particularly in Season 2, and high yield performance, though with variable stability.Genotype 48 also demonstrated strong performance, especially in Season 1, with better stability in both disease resistance and yield (Figures 3A, 4A).Genotypes 108, 1171, and 1184 showed significant interactions with environmental factors, indicating distinct responses that could be leveraged for targeted breeding.Figures 3B and 4B highlight the similarity between Seasons 1 and 3, suggesting that genotypes were performing well in Season 1 were likely to do similarly in Season 3, while Season 2 differed significantly, indicating that genotypes might respond differently due to unique environmental factors like temperature and humidity.Figures 3C and 4C identified top performers in specific environments: Genotype 45 excelled in Season 2, and Genotype 48 excelled in Season 1.Other genotypes, such as 108 and 1168, also showed good performance under certain conditions, highlighting

Principal component analysis of durum wheat traits
The Principal Component Analysis (PCA) biplot, was used to capture the performance similarities among genotypes based on multiple trait analyses.This biplot explained 72.24% of the total variance, with Principal Component 1 (PC1) accounting for 49.96% and Principal Component 2 (PC2) accounting for 22.28% (Figure 5).Genotypes with high yield were strongly differentiated along PC1, indicating that yield was a major determinant of overall performance.Traits associated with PC2, such as lower Area Under Disease Progress Curve (AUDPC) values, highlighted disease resistance as key for genotype differentiation.Genotypes exhibiting high yield and plant height were identified as key candidates for high-yield breeding programs, while those with lower AUDPC values were vital for breeding disease-resistant varieties.Most genotypes clustered around the origin, showing balanced performance and stability across multiple traits, suggesting reliability under diverse environments.

Discussion
This study found significant seasonal variations in response of durum wheat genotypes to stem rust infections, confirming the substantial phenotypic diversity among the genotypes tested.These variations highlighted the role of environmental factors in expression of stem rust severity, as measured by FDS and AUDPC.The study's ability to capture these variations is emphasized by our rigorous approach and extensive evaluation conducted over several seasons.Similarly, findings of Shaw et al. (2008) and Mideksa et al. (2018), noted that environmental  conditions, such as rainfall and temperature, play role in progression of stem rust, with temporal factors having significant impact on disease dynamics.However, our focus on specific geographical locations might limit the generalizability of these results to other regions experiencing different climatic conditions.Furthermore, our findings that seasonal conditions contribute to variations in disease severity emphasized the complex influence of environmental factors, as noted by Nopsa and Pfender's (2014) observation that environmental conditions can sometimes overshadow genetic resistance.Other insights also from Sgrò and Hoffmann (2004), highlighted the significant impact of environmental variations on phenotypic trait expression and genetic correlations.Season 1 exhibited significantly greater disease severity compared to Seasons 2 and 3, suggesting that the environmental conditions during this period exacerbated disease symptoms.This also aligns with Velaśquez et al.'s (2018) findings, which demonstrated how weather conditions intricately affect plant-pathogen interactions.A significant negative correlation was found between area under the disease progression curve (AUDPC), final disease severity (FDS), and yield components.These findings highlight the importance of effective resistance genes in maintaining high yields even in the face of disease, demonstrating the impact of genetic resistance in agricultural settings.Other studies, Ogutu and Charimbu (2020); Msundi et al. (2021) and Wanyoike et al., 2022 also discovered that both AUDPC and FDS had negative impact on grain yield, emphasizing the importance of genetic resistance in mitigating disease-induced yield reductions.
The race-specific gene Sr13 has been essential in providing resistance against Ug99 race of stem rust and its variants in durum wheat (Singh et al., 2015).Despite the effectiveness of Sr13 (Zhang et al., 2017), the emergence of non-Ug99 races such as TTKSK, JRCQC, and TTRTF, which have developed virulence against alleles like Sr13b and Sr9e (Olivera et al., 2012), underscores the dynamic nature of pathogen evolution and the challenges it poses (Olivera et al., 2019).These races have developed virulence against alleles like Sr13b and Sr9e, which differ from Sr13 primarily in their genetic makeup and the resulting functional responses to the pathogen.For instance, Sr13b might confer resistance through slightly altered protein configuration that affects its interaction with the pathogen, leading to differential susceptibility when compared to Sr13.The presence of Sr13 and related resistance alleles in the genetic backgrounds of durum wheat genotypes we assessed is crucial, as it continuously necessitates updates in breeding strategies to effectively counter these evolving threats.Mergessa et al. (2020) drew attention to the necessity of constant innovation in breeding strategies, highlighting the implications of our findings for continued adaptability in the face of changing disease threats.The Genotype main effect and Genotype × Environment interaction biplots depicts the performance of 49 durum wheat genotypes in terms of stem rust severity across three seasons,with subfigures showing the GGE biplot for main effects and interactions (A), the radial relationship among environments (B), polygon view highlighting top performers by environment (C), and plot of mean performance versus stability (D).
Our findings revealed that approximately 40.82% of durum wheat genotypes, exhibiting up to 30% stem rust severity, were categorized as resistant to moderately resistant.These partially resistant genotypes underscore the effectiveness of Sr genes, under field conditions.However, to fully utilize Sr gene-mediated resistance in breeding programs, additional research into the genetic backgrounds and environmental interactions is likely necessary given the heterogeneity in resistance levels between genotypes.This resistance is crucial not only for maintaining yield but also for reducing the selection pressure that drives the evolution of virulence in pathogen populations.However, the variability in resistance levels across genotypes suggests further exploration into the genetic backgrounds and environmental interactions is necessary to fully leverage Sr gene-mediated resistance in breeding programs.Studies by Olivera et al. (2021) also demonstrated similar levels of resistance in USDA durum wheat accessions, noting that partial resistance can effectively slow the adaptation of pathogens by reducing selection pressure for virulent mutations.Research by Burdon and Laine (2019) on the evolutionary dynamics of plant-pathogen interactions suggested that partial resistance, involving spectrum of resistance genes, may provide more sustainable approach to managing disease than relying on complete resistance from single genes.These models propose that such strategies not only mitigate the risk of resistance breakdown but also contribute to maintaining crop yield under various stress conditions.
Our research showed genetic influences on key traits such as plant height and final disease score, with broad-sense heritability estimates indicating strong genetic basis.However, the low heritability of kernels per spike (KPS) suggested significant environmental impact, complicating selection based solely on this trait.Furthermore, we found negative correlations between disease severity (AUDPC) and key agronomic traits like KPS, kernel weight, and grain yield.These correlations, as also described by Chiko et al. (2022), imply that lower disease severity and increased biomass are linked to higher yields and grain weight.Our findings agree with those of Willocquet et al. (2018Willocquet et al. ( , 2022)), who found that increased disease severity not only reduces plant height, but also redirects resources from reproductive structures to pathogen spore growth during critical growth phases.This redirection has significant impact on grain quality and yield, highlighting the importance of selecting for traits that improve resource allocation to reproductive structures in order to reduce stem rust impacts and maintain yield under disease pressure.Such physiological impacts suggest that breeding for traits that enhance resource allocation to reproductive structures could mitigate the effects of stem rust, thus preserving yield under disease pressure.
Durum wheat production in Kenya is relatively low compared to other countries, such as Ethiopia (Tidiane Sall et al., 2019).

A B D C FIGURE 4
The Genotype main effect and Genotype × Environment interaction biplots illustrates the stability of 49 durum wheat genotypes' grain yield across three seasons,featuring GGE biplot for relationship and discrimination by seasons (A), radial plot of the relationship among environments (B), polygon view identifying top-performing genotypes in each environment (C), and plot assessing mean yield versus stability (D).
Therefore, choosing high-yielding, stable genotypes for use in breeding programs is essential to enhance durum wheat productivity.A comprehensive understanding of genotype-byenvironment interactions (GEI) is crucial for setting breeding goals, selecting optimal testing locations, and making variety recommendations (Yan and Hunt, 2001).This understanding is also vital for assessing the adaptability and stability of different genotypes.Our research showed that the major source of variability came from G × E interactions, which suggested significant environmental influence on genotype performance.This implies that genotypes developed for specific environments might not perform consistently well in other settings, highlighting the need to broaden the scope of evaluation trials across various ecological conditions.Additionally, it's important to cultivate identified genotypes in their ideal environments to maximize the benefits of positive GEI (Khan et al., 2021).Thus, identifying suitable environments for breeding and choosing strategic sites for genotype testing are key priorities.
In our study, we used the GGE biplot method to evaluate the interaction of three environments with genotype performance, focusing on stem rust severity and grain yield.The analysis pinpointed genotypes that were stable and specifically adapted to different environmental conditions, which is vital for effective breeding strategies.The results underscore a crucial link between yield stability and disease resistance in the evaluated genotypes.For example, Genotype 45 displayed strong resistance to stem rust in the second season and also showed high yield performance, suggesting that breeding for disease resistance does not necessarily detract from yield potential.By selecting for both traits, robust varieties suitable for various conditions can be developed.Genotype 48 was notable for its stability in both disease resistance and yield, reinforcing the potential for breeding high-performing, disease-resistant wheat.The significant GxE interactions observed highlight the importance of testing genotypes in multiple settings.These interactions reveal the unique performance of specific genotypes, like 1171 and 1184, under varying environmental conditions, providing essential insights for breeders in selecting adaptable and resilient varieties.Furthermore, the observed genotype × environment interactions emphasize the need for breeding programs to consider both resistance and yield stability.This reflects the complex interplay of genetics and environmental influences, further complicated by shifting climatic conditions that can exacerbate disease severity and affect water availability.Hence, developing genotypes that are not only resistant to pathogens but also resilient to abiotic stresses such as drought and heat is critical (Ramegowda and Senthil-Kumar, 2015;Pandey et al., 2017).This approach is essential, particularly in areas prone to extreme weather events, highlighting the need for multi-environment trials to capture the full performance potential of genotypes, as supported by other findings from Tekdal and Kendal (2018) and Madahana et al. (2021).
Principal component analysis (PCA) on the other hand, is a multivariate technique that helps to understand the interrelationships among multiple traits by transforming them Principal Component Analysis show performance similarities among the evaluated 49 durum genotype traits through multiple trait analysis.
into a set of new, uncorrelated variables called principal components (Mishra et al., 2017;Greenacre et al., 2022).This technique is particularly valuable in plant breeding as it enables the identification of key traits contributing to genetic diversity and performance across different environments.PCA provides visual representation of the patterns of similarity among genotypes, allowing breeders to pinpoint those with desirable characteristics, such as high yield and disease resistance, while maintaining stability under various environmental conditions.Here, PCA offered insights into the performance similarities among durum wheat genotypes based on multiple trait analyses.The PCA biplot explained 72.24% of the total variance, with Principal Component 1 (PC1) accounting for 49.96% and Principal Component 2 (PC2) for 22.28%.This analysis differentiated genotypes based on key traits, emphasizing the importance of yield and disease resistance.For instance, genotypes with high yield were strongly associated with PC1, highlighting yield as a major determinant of overall performance.On the other hand, traits related to PC2, such as lower Area Under Disease Progress Curve (AUDPC) values, highlighted the importance of disease resistance in genotype differentiation.The clustering of genotypes around the origin indicated balanced performance and stability across multiple traits, suggesting their reliability under diverse environmental conditions.This aligns with the findings of Alemu et al. (2020), who emphasized the importance of genetic variability and trait associations in Ethiopian durum wheat landraces.

Conclusion
This research conducted comprehensive analysis of genotypeenvironment interactions impacting durum wheat in Kenya, shedding light on the trade-offs between stem rust resistance and agronomic performance.Our findings identified specific genotypes such as G45 and G48, that consistently demonstrate resistance to stem rust across various environmental conditions, marking them as ideal candidates for breeding programs aimed at enhancing crop resilience and productivity.The high resistance to stem rust ensures lower disease incidence and severity, which translates to healthier plants and higher yields.Moreover, their stable performance across diverse environments implied robust adaptability, crucial for resilience against climatic variability.Additionally, the application of genotype and genotype by environment (GGE) biplot analysis had been instrumental in identifying these genotypes that not only perform well agronomically but also withstand the pressure of stem rust.To enhance food security and promote sustainable agriculture in Kenya and potentially in similar regions globally, it is essential to focus on breeding programs that balance yield stability and disease resistance.Leveraging the genetic diversity within durum wheat will enable breeders to develop varieties tailored to Kenya's diverse agro-ecological zones, ultimately ensuring robust and resilient crop performance across varying conditions.

TABLE 1
Pedigree crosses of selected durum germplasm used in the experiment.

TABLE 3
Effect of seasons on stem rust severity, Area under disease progress curve, spike length, plant height, kernels per spike, 1000-kernel weight and yield of durum genotypes.
Means followed by the same letters are not significantly different at p ≤ 0.05.AUDPC, Area under disease progress curve; SPKL, spike length; PH, plant height; KPS, kernels per spike; TKW, 1000kernel weight.

TABLE 4
Heritability, Variance components, and means of stem rust severity, area under disease progress curve and yield and yield components of durum genotypes evaluated over three seasons.

TABLE 5
Phenotypic (bold)and genotypic correlation of final disease score, area under disease progress curve and agronomic traits of durum traits bold=phenotypic correlation.