Phenotypic Variability, Heritability and Associations of Agronomic and Quality Traits in Cultivated Ethiopian Durum Wheat ( Triticum turgidum L. ssp. Durum , Desf.)

: Quality is an important aspect of durum wheat in the processing sector. Thus, recognizing the variability of quality and agronomic traits and their association is fundamental in designing plant breeding programs. This study aimed to assess the variability, heritability, genetic advance, and correlation of some agronomic and quality traits among 420 Ethiopian durum wheat genotypes and to identify the promising genotypes with distinct processing quality attributes to produce superior quality pasta. The ﬁeld experiment was conducted at two locations (Sinana and Chefe Donsa) using an alpha lattice design with two replications. Analysis of variance, chi-square test, and Shannon–Weaver diversity index revealed the existence of highly signiﬁcant ( p < 0.001) variation among genotypes for all studied traits. The broad-sense heritability values were ranging from 46.2% (days to maturity) to 81% (thousand kernel weight) with the genetic advance as a percent of the mean ranging from 1.1% (days to maturity) to 21.2% (grain yield). The phenotypic correlation coefﬁcients for all possible pairs of quantitative traits showed a signiﬁcant ( p < 0.05) association among most paired traits. The gluten content (GC) and grain protein content (GPC) were negatively correlated with grain yield and yield-related traits and positively associated with phenological traits, while yield and phenological traits correlated negatively. The frequency distributions of amber-colored and vitreous kernels, which are preferable characters of durum wheat in processing, were highly dominant in Ethiopian durum wheat genotypes. The identiﬁed top 5% genotypes, which have amber color and vitreous kernel with high GC and GPC content as well as sufﬁcient grain yield, could be directly used by the processing sector and/or as donors of alleles in durum wheat breeding programs.


Introduction
Durum wheat (Triticum turgidum L. ssp. Durum, Desf.) is one of the most commonly cultivated cereal crop globally with a production of 33.8 million tons in 2019/20, which is 5% of the global wheat production [1]. Ethiopia is considered a center of durum wheat diversity [2][3][4][5] where the crop has been cultivated in small-to-medium-scale farming systems in diverse agro-ecological and climatic conditions [6,7]. Although durum wheat has a lower grain yield potential than bread wheat [8], it continues to be the crop of choice for smallholder farmers. Major factors contributing to its preference include its wide adaptation to diverse climatic conditions [9,10], a higher price than bread wheat, millet,

Data Collection
Data were collected for the following phenotypic traits: days to heading (DH), days to maturity (DM), number of spikelets per spike (SPS), spike length (SPL), plant height (PLH), thousand kernel weight (TKW), grain yield (GY), lodging (LDG), spike density (SPD), glume hairiness (GLH), glume color (GLC), kernel color (KCL), vitreousness (VTR), grain moisture content (GMC), gluten content (GC) and grain protein content (GPC). The SPS, PLH, LDG, SPD, GLH, GLC, KCL, and VTR traits data were collected using the descriptor list for wheat [33] whereas DH, DM, and SPL traits data were col- The field experiment was conducted at two sites in Ethiopia, Sinana (07 • 07 N and 40 • 10 E, 2400 m.a.s.l) and Chefe-Donsa (08 • 44 N and 39 • 09 E, 2450 m.a.s.l), during the 2019/2020 main crop growing season using an alpha lattice design with two replications. The soil types at both Sinana and Chefe Donsa research centers were vertisols with minor light soils. The genotypes were randomly assigned to plots within each block. The area of each plot was 0.4 m 2 (2 rows with 1 m length and 0.4 m width) and the spacing between entry, plants, blocks, and replications was 0.4 m, 0.1 m, 1 m, and 2 m, respectively. Fertilizer was applied at a rate of 100 Kg/ha of DAP and 50 Kg/ha of urea at the planting and jointing stages. All other management practices recommended for wheat were also applied to all experimental plots.

Data Collection
Data were collected for the following phenotypic traits: days to heading (DH), days to maturity (DM), number of spikelets per spike (SPS), spike length (SPL), plant height (PLH), thousand kernel weight (TKW), grain yield (GY), lodging (LDG), spike density (SPD), glume hairiness (GLH), glume color (GLC), kernel color (KCL), vitreousness (VTR), grain moisture content (GMC), gluten content (GC) and grain protein content (GPC). The SPS, PLH, LDG, SPD, GLH, GLC, KCL, and VTR traits data were collected using the descriptor list for wheat [33] whereas DH, DM, and SPL traits data were collected following the description provided in Hailu et al. [34]. TKW (the weight of 1000 kernels in grams) was calculated by multiplying the weight of 250 kernels by four, and GY was estimated based on the weight (in grams) of the total seed harvest per plot. GMC, GC, and GPC were

Statistical Analysis
Analysis of variance and variance components: prior to statistical analysis, the collected quantitative data were checked for normal distribution using the Shapiro-Wilk test through the shapiro.test function of the Stats package [36] in R software. The effects of genotype, replication, block, location, and genotype by location interaction were tested via analysis of variance (ANOVA) using the aov function of the Stats package [36] based on the alpha lattice design linear additive model.
where P ijk is the phenotypic value of ith genotype under jth replication and kth block, µ is the grand mean, τ i is the effect of ith genotypes, π j is the effect of jth replication, b k is the effect of k th block and ξ ijk is residual error. The phenotypic, genotypic, and genotype by location interaction variances were computed as described by Singh and Chaudhary [37]: where MS g is the mean square of genotype, MS g×l is the mean square due to genotype by location interaction, MS e is the error mean square, r is the number of replications and l is the number of locations. The genotypic and phenotypic coefficient of variation was calculated as described in Johnson et al. [38]: Genotypic coefficient of variation = GCV = 100 × σ g /µ (5) Phenotypic coefficient of variation = PCV = 100 × σ p/ µ where σ g and σ p are genotypic and phenotypic standard deviations, respectively, and µ is the grand mean. The GCV and PCV values were categorized as high (>20%), medium (10-20%) and low (<10%) as classified by Burton and DeVane [39]. Heritability and genetic advance: broad-sense heritability (H b ) was calculated using the formula presented in Allard [40]: where σ 2 g and σ 2 p are genotypic and phenotypic variance, respectively. The values of broadsense heritability were classified as low (H b < 25%), moderate (25% < H b < 50%), and high (H b > 50%).
The expected genetic advance and genetic advance as a percent of the mean were estimated as described in Johnson et al. [38]: GA as percent of mean (GAM) = (GA/µ) × 100 (9) where σ g and σ p are genotypic and phenotypic standard deviations, respectively, k is a constant value at a selection intensity of 5% (k = 2.06); H b is broad-sense heritability, GA is genetic advance and µ is the grand mean. The GAM is considered low, moderate, and high when its values are below 10%, 10-20%, and above 20%, respectively, as described in Falconer and Mackay [41].
Correlation between traits: The Pearson correlation coefficient (r) between quantitative phenotypic traits was computed using the chart.Correlation function of the Performance-Analytics package [42] in R.
Frequency distribution, chi-square test and Shannon diversity index of qualitative data: the qualitative data were transformed into frequency data and subjected to the chi-square test using prob.table and chisq.test function of the Stats package [36]. The standardized Shannon-Weaver diversity index ("H ") was calculated from the frequency data of qualitative traits as described in Eticha et al. [43]: where n is the number of phenotypic classes of the trait and pi is the proportion of ith phenotypic class. H was estimated for all qualitative traits and all genotypes. The Shannon diversity index was regarded as high (H ≥ 0.60), intermediate (0.40 ≤ H < 0.60) and low (H < 0.40) as described in Eticha et al. [43].

Phenotypic Variation among the Genotypes
A wide range of variation was observed for most of the quantitative phenotypic traits ( Table 1). The mean GPC and GC values of the genotypes were 13.4% and 25.7%, respectively ( Table 1). The optimal GC for pasta and related products is 28% and above, according to Sissons [15]. Interestingly, 97 landraces of the 386 landraces (25.1%) had a GC value 28% and above, with genotypes 364, 150,191,287,354,263,282,285,294,67,197, and 276 having GC values above 32%. On the other hand, all cultivars had a gluten content of below 28%. Likewise, 285 genotypes had a grain protein content (GPC) of 13%, which is the lower limit for pasta production [44], with genotypes 208,190,386,191,364,282,373,354,285,247,224,110,263,287,205,25,272, and 115 having a GPC of above 15%. In total, 208 genotypes (49.5%) produced grain yields exceeding the mean value of 336.4 g of seeds per plot (gm/pt), of which 48 genotypes (11.4%) produced grain yields exceeding 400 g per plot (Table 1). One hundred and seventy-four genotypes (41.4%) had a TKW above the mean value of 42.4 g, of which 34 genotypes had a TKW above 50 g (Table 1). DH = days to heading; DM = days to maturity; SPS = the number of spikelets per spike; SPL = spike length; PLH = plant height; TKW = thousand kernel weight; GY = grain yield; GC = gluten content and GPC = grain protein content.

Analysis of Variance (ANOVA)
The ANOVA has revealed significant (p < 0.05) effects of most phenotypic traits on Ethiopian durum wheat genotypes with R 2 values ranging from 0.67 (DM) to 0.89 (DH) ( Table 2). Highly significant (p < 0.001) differences among genotypes were observed for all phenotypic traits. The test locations (environments), as well as the genotype by location interaction (genotype by environment interaction; G × E), also had significant (p < 0.05) effects on all traits ( Table 2). DH = days to heading; DM = days to maturity; SPS = the number of spikelets per spike; SPL = spike length; PLH = plant height; TKW = thousand kernel weight; GY = grain yield; GC = gluten content and GPC = grain protein content. df = degree freedom; GEN = genotype; LOC = location; REP = replication; BLK = block; LSD = the least significant difference and CV = coefficient of variation. "***" = significant at p < 0.001; "*" = significant at p < 0.05; "ns" = non-significant at p < 0.05.

Broad-Sense Heritability and Genetic Advance
The values of broad-sense heritability were high to intermediate for all traits, which ranged from 42.6% for DM to 81% for TKW (Table 1). Within this range, high heritability values (H b > 50%) were observed for TKW (81%), DH (76.5%), SPL (68.5%), GC (55.6%), and PLH (52.9%), whereas moderate heritability values (25% < H b < 50%) were recorded for the remaining traits. GA is a useful selection parameter in addition to heritability for crop improvement. The GAM values for most of the traits were low except for GY (21.2%), SPL (19.1%), and GC (13.4%), which can be regarded as moderate (Table 1).

Frequency Distribution and Diversity of Qualitative Traits
For data analysis, five of the six qualitative traits (SPD, GLC, GLH, KCL, and VTR) were classified into three phenotypic classes, while LDG was classified into four classes (Table 3). In terms of lodging, 47.8% of the genotypes were lodging types with medium to complete lodging. The vast majority (92.3%) of the genotypes had hairless glumes. Interestingly, about 70% and 76.4% of the genotypes had amber kernel color and vitreous kernel, respectively, which are among the most important processing quality traits. The chi-square test revealed highly significant (p < 0.05) variation among the genotypes in terms of all qualitative traits with Shannon diversity index (H′) ranging from 0.27 (GLH) to 0.96 (SPD) ( Table 3). Hence, the diversity of each trait is high (H′ > 0.60) for all qualitative traits except for GLH (H′ = 0.27). Test locations (environments) had a significant effect (p < 0.05) on LDG and GLC but not on the other qualitative traits (Table 3).

Frequency Distribution and Diversity of Qualitative Traits
For data analysis, five of the six qualitative traits (SPD, GLC, GLH, KCL, and VTR) were classified into three phenotypic classes, while LDG was classified into four classes (Table 3). In terms of lodging, 47.8% of the genotypes were lodging types with medium to complete lodging. The vast majority (92.3%) of the genotypes had hairless glumes. Interestingly, about 70% and 76.4% of the genotypes had amber kernel color and vitreous kernel, respectively, which are among the most important processing quality traits. The chi-square test revealed highly significant (p < 0.05) variation among the genotypes in terms of all qualitative traits with Shannon diversity index (H ) ranging from 0.27 (GLH) to 0.96 (SPD) ( Table 3). Hence, the diversity of each trait is high (H > 0.60) for all qualitative traits except for GLH (H = 0.27). Test locations (environments) had a significant effect (p < 0.05) on LDG and GLC but not on the other qualitative traits (Table 3).

The Performance of Selected Genotype
First, genotypes having vitreous grains with amber kernel color were selected (Table S2), and then the top 5% of the total genotypes analyzed were selected based on their GC and GPC content by also taking TKW and grain yield into consideration. Because of the highly negative correlation between GC and GPC with grain yield and grain yield-related traits, genotypes with better TKW and GY were prioritized. The difference between the grand mean and the mean of the selected genotypes was significant (p < 0.05) for GC and GPC, but not for the other quantitative traits ( Table 4). The GC mean value for the top 5% genotypes (30.5%) was greater by 4.8% than the grand mean of all genotypes (25.7%) ( Table 4), thus meeting the highest grain quality grade as per the description of Sissons [15]. Interestingly, the mean GPC for the top 5% genotypes was 14.4%, exceeding the minimum threshold (13%) required for high processing quality [44]. It should be noted that none of the released cultivars are part of these top 5% genotypes. 1.0 6.9 2.43 ** DH = days to heading; DM = days to maturity; SPS = the number of spikelets per spike; SPL = spike length; PLH = plant height; TKW = thousand kernel weight; GY = grain yield; GC = gluten content and GPC = grain protein content. t-value (0.05) (df = 41) = 1.68. "***" = significant at p < 0.001; "**" = significant at p < 0.01; "ns" = non-significant at p < 0.05.

Variations of Quantitative Traits among Genotypes
A thorough understanding of genetic variation within a crop gene pool, as well as genetic relationships between germplasm of different origins, will facilitate the development of new cultivars with desirable traits. Our understanding of the heritability and correlations of their traits is also imperative. These can be gained through phenotypic studies based on agro-morphological and quality traits. In the present study, nine quantitative and six qualitative traits of durum wheat were analyzed, revealing a wide range of phenotypic variation among the genotypes in terms of most of the traits studied. This suggests the presence of high genetic variation within the Ethiopian durum wheat gene pool. Previous studies on durum wheat also reported a wide phenotypic variation for GPC, TKW, and GY [30,[45][46][47]. Among the 420 genotypes studied, 285 of them meet the minimum threshold value of 13% GPC content necessary for making high-quality pasta and related products [44]. Interestingly, 68 of the 285 genotypes performed well at both experimental sites and were therefore stable across environments, suggesting their abilities to adapt to various agro-ecological conditions.
The grain gluten content of durum wheat is also crucial in the processing for making pasta and related products because gluten confers elasticity and viscosity to the dough [17]. Ninety-seven genotypes had grain gluten content of above 28% with the highest values recorded in genotypes 364,150,191,287,354,263,282,285,294,67,197, and 276, which were also characterized by having high grain protein content, delayed days to maturity, high plant height, long spike, high thousand kernel weight, and intermediate grain yield. All cultivars evaluated in this study had less than 28% grain gluten. This is mainly because this trait has not been seriously taken into consideration during the breeding process, which focused mainly on grain yield improvements. Hence, crossbreeding that involves genotypes with high gluten content, on the one hand, and high-yielding genotypes, on the other hand, is crucial in order to develop high-yielding cultivars with desirable levels of grain gluten content. Although the genotype by location interaction (G×E) was significant for GC, 24 landraces scored GC above 28% at both locations, suggesting a wide adaptation of these genotypes. Therefore, they should be considered in breeding programs for durum wheat as sources of desirable agro-morphological and quality traits.
The highly significant (p < 0.01) phenotypic variation among genotypes for most of the studied traits indicates the presence of an adequate amount of genetic variation among characterized genotypes. In line with the results of the present study, Mengistu and Pe [27] found significant differences among 49 Ethiopian durum wheat genotypes (comprising 45 landraces and 4 cultivars) in DM, PLH, TKW, and GY. Similarly, Wolde et al. [47] and Mengistu et al. [48] reported a significant variation among 68 and 289 Ethiopian durum wheat genotypes, respectively, in DM, PLH, TKW, and GY. Studies on durum wheat genetic resources of non-Ethiopian origins have also demonstrated significant variation among genotypes in GC, GPC, GY, TKW, and PLH [46,49,50]. These indicate that durum wheat has the potential for better adaptation to different agro-ecological conditions.
The highly significant (p < 0.01) G×E interaction for all quantitative traits obtained in this study and similar results of previous reports for GPC, TKW, and GY [50,51] indicate that genotypes differ in their adaptability to different agro-ecological conditions, with some genotypes having excellent adaptation to specific environments.
For all traits, except GY, relatively small differences were observed between genotypic and phenotypic coefficients of variation. This indicates that genotypes have a more profound influence than environmental factors on the phenotypic expressions of these traits in durum wheat. Similarly, relatively smaller differences between GCV and PCV values were reported by Mohammed et al. [52] for TKW and GPC, and Wolde et al. [47] and Alemu et al. [53] for DH, DM, PLH, and TKW. These suggest a lower effect of environmental factors on phenological and quality traits. On the other hand, the relatively large differences between GCV and PCV values for GY obtained in the present study and that of Alemu et al. [53] indicate that grain yield is highly influenced by environmental factors.

Selection Efficiency of Quantitative Traits in Durum Wheat
The determination of the heritability and genetic advance of traits of interest in crops allows the selection of suitable parents for crossbreeding [24], prediction of the performance of progenies in successive generations, and estimation of the degree of genetic gains obtained from the selection [25]. Assuming 5% selection intensity, high heritability combined with high to moderate genetic advance indicates the presence of additive gene effects [54,55] while high heritability coupled with low genetic advance suggests the presence of non-additive gene effects for a given trait [55,56]. The results of the present study indicate the presence of strong additive gene effects for GC, TKW, and GY, which suggests the importance of undertaking effective phenotype-based selection at the early generation of breeding during cultivar development. The additive gene effect was also previously reported for TKW [46,47,51,52]. The influence of the additive gene effect on GY identified in this study was also reported by Wolde et al. [47], and hence early generation selection could be effective. In contrast, Alemu et al. [53], as well as Mengistu and Pe [27] reported a non-additive gene effect on GY, which could be attributed to the small numbers of genotypes (64 and 49 durum wheat landraces, respectively) they used.
In the present study, plant height was negatively correlated with grain yield, which may suggest a trade-off between the two traits in durum wheat. However, PLH was positively correlated with both grain gluten and protein contents (Figure 2). It may be difficult, therefore, to select genotypes using PLH as a selection criterion for crossbreeding with the aim of producing high-yielding cultivars with desirable levels of grain gluten and protein content. Hence, introducing dwarfing genes through modern breeding approaches could improve the grain yield of genotypes along with desirable levels of grain gluten and protein content [57].
The significant influence of the environment on the expression of protein traits lowered the heritability of both GC (55.6%) and GPC (44.8%) ( Table 1). Based on the moderate heritability and low genetic advance obtained for GPC in this study, it appears that there is a non-additive gene effect, which is in accordance with previous studies [52,53], so it should be noted that selection for GPC may be complicated both in the early and advanced generations. The high heritability with a moderate GAM value of GC indicates the presence of an additive gene effect on the trait. This suggests the possibility of undertaking effective genotype selection through phenotypic expression at early generation and enables the development of effective breeding strategies for improving GC in durum wheat. Previously, there were limited numbers of researches conducted to determine the heritability value of GC in durum wheat. For instance, Brankovic et al. [45] and Mohammed et al. [52] reported high and moderate, respectively, heritability and moderate GA for wet gluten content.

Associations of Quantitative Traits in Durum Wheat
An understanding of the associations between different traits is a crucial step in crop breeding for both improving target traits simultaneously and indirectly selecting other desirable traits. The significant negative correlation between GPC and GY obtained in this study is in agreement with that of Alemu et al. [53], Sourour et al. [49], and Bilgin et al. [58].
However, Grahmann et al. [59] found a significant positive correlation, which deviates from the commonly reported negative correlation between the two traits. The significant negative correlation between GPC and TKW obtained in the present study has also been reported previously [59][60][61]. In contrast, Mohammed et al. [52] found a highly significant positive correlation (p < 0.01) between these traits. Sourour et al. [49] reported a significant positive correlation (p < 0.05) of GC with PLH, SPL, SPS, and GPC, while negatively correlated with GY. Similarly, Bilgin et al. [58] found a significant positive correlation (p < 0.05) of GC with GPC and a significant negative correlation (p < 0.05) with TKW and GY. These results suggest that grain yield tends to correlate negatively with processing quality traits in durum wheat landraces. Hence, breeding for the improvement of processing quality traits may lead to a decrease in grain yield. Hence, exploring the diverse durum wheat gene pool to identify genotypes with high GPC and GC along with moderate to high TKW and GY that can be used in a breeding program is a crucial step for avoiding the trade-off between the two groups of traits.

Distribution and Variability of Qualitative Traits among Genotypes
The proportion of phenotypic classes of the studied qualitative traits was variable among Ethiopian durum wheat genotypes ( Table 3). The existence of technologically important phenotypic characters in high frequency in Ethiopian durum wheat genotypes allows them to be valuable for the processing sector. Among these traits, amber kernel color and vitreous kernel classes are dominant in the studied genotypes. Similarly, Hailu et al. [34] reported the dominance of amber kernels among tetraploid wheat. However, Mengistu et al. [30] described the rarity of amber kernels in Ethiopian durum wheat landraces. The distinct genotypes with a purple kernel color observed in this study were also reported by Bechere et al. [28], which accounted for about 16% of their collections. In addition, Zeven [62] identified that genotypes with purple kernels are endemic to Ethiopia, and hence classified as Triticum turgidum sub. aethopicum [43].
The GLH traits showed lower diversity than other qualitative traits with the dominance of the hairless subclass. Similar findings were reported in previous works [28,43,63], which indicate the majority of Ethiopian landraces have hairless glumes. In contrast, a high level of polymorphism of GLH in Ethiopian durum wheat genotypes was reported by Bekele [64]. Negassa [9] suggested that the hairy glume phenotypic class of durum wheat may be associated with resistance to powdery mildew. The observed higher frequencies of lax and intermediate spike subclasses than that of dense spike subclass for SPD (Table 3) were contrary to the results obtained by Bechere et al. [28] and Hailu et al. [34] who reported the dominance of a dense spike in the Ethiopian durum wheat collections. The highest diversity index of SPD creates opportunities for breeders to develop suitable varieties for various environmental conditions. For instance, the lax spike is associated with resistance to spike-related diseases [43,63]. In agreement with this study, Jain et al. [65] reported the rarity of black glume color in Ethiopian durum wheat genotypes which is associated with gluten quality [9].
The observed significant genotypic differences, based on the chi-square test, and high Shannon diversity index for all qualitative traits, except for GLH, suggest the existence of high variation within Ethiopian durum wheat genotypes for the traits targeted in this study ( Table 3). The significant impact of test sites on LDG and GLC traits identify that the growing environment had significant effects on the performance of durum wheat genotypes for these traits, while location had no significant effect on the remaining qualitative traits, implying that these traits are more governed by genetic than environmental factors (Table 3).

Promising Landraces for Processing Quality and Future Breeding
In this study, all genotypes that were ranked in the top 5% were landraces (Table S2), and none of the cultivars made it to the top 5%, mainly due to their low GC and GPC. The results may imply a trade-off between grain yield and some quality traits, as the cultivars were bred primarily for grain yield and related traits. The results also suggest the superiority of durum wheat landraces over cultivars with regard to processing quality traits. The mean value of the grand population and the selected genotypes showed significant differences for GC and GPC, as they are the criterion for the selection of genotypes (Table 4). However, grain yield and related traits, such as TKW, which correlated negatively with the prioritized quantitative quality traits (GC and GPC), appeared to have non-significant differences in the mean value of the grand population and the selected top 5% genotypes (Tables 4 and S2). These landraces with optimal GC and GPC along with sufficient grain yield could be used in wheat breeding programs to improve both groups of traits simultaneously.

Conclusions
This study revealed the existence of substantial genetic variation among Ethiopian durum wheat genotypes for agro-morphological and quality traits. Some genotypes performed better than others for several phenotypic traits across the experimental sites. Hence, further studies on top performers to validate their stability over multiple environments are required. The moderate to high heritability with moderate genetic advance for TKW, SPL, GY, and GC suggests that the selection of genotypes for these traits in early generations could be effective. Ethiopian durum wheat may add value to the processing sector due to its significant proportions of distinct phenotypic characters such as amber color and vitreous kernels. Protein-related traits showed a significant positive correlation (p < 0.05) with phenological traits and a significant negative correlation (p < 0.05) with yield-related traits. The results may suggest breeding for the improvement of protein-related traits might negatively affect GY in durum wheat. Therefore, the top 5% of genotypes that have vitreous kernels with amber color, high GC and GPC along with optimal TKW and GY are highly valuable resources that should be used as donors of desirable alleles in durum wheat breeding programs. Given that only two locations were used in this study, traits that showed significance (p < 0.05) for genotypes by environmental interaction should be studied further across environments to determine the consistency of their performance.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/agronomy12071714/s1, Table S1: Description of Ethiopian durum wheat genotypes used in the phenotypic diversity study. Table S2: Genotypes ranked in the top 5% based on their values in gluten content, grain protein content, thousand-kernel weight and grain yield over combined locations.