Abstract

Bread wheat is the world’s leading cereal grain, and more than one-third of the world’s population uses it as a staple food. The bread wheat production in Ethiopia is low compsssared to the national average yield, mainly due to the lack of high-yielding genotypes. This study was conducted during the 2019-2020 growing season to assess genetic variability and estimate the association of traits among bread wheat genotypes. The experiment consists of 49 bread wheat genotypes and is laid out in 7 × 7 simple lattice designs. The results showed significant differences ( < 0.01) among genotypes for most of the studied traits. Moderate genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV) values were estimated for yield plant −1, thousand seed weight, and biomass yield. High heritability coupled with a high GAM was observed for thousand seed weight and yield plant –1. The grain yield showed a highly significant ( < 0.01) correlation with many yield-related traits at the phenotypic and genotypic levels. The biomass yield and the harvest index exerted the highest positive direct effect on grain yield at the genotypic level. The highest intercluster distance was observed between clusters I and IV (D2 = 31.86), followed by clusters II and IV (D2 = 29.21), and clusters II and III (D2 = 28.24), which indicated the chance of selecting a member of these clusters for hybridization. This experiment’s result indicates sufficient genetic variability among the tested genotypes, which provides ample scope for selecting superior and desired genotypes. Best-performed genotypes should be included in the future breeding program for further yield improvement. In conclusion, attention should be given to traits with moderate to high heritability and GAM, exerting a positive direct effect on the grain yield. However, the experiment should be repeated over locations and seasons to draw a definite conclusion.

1. Introduction

Bread wheat (Triticum aestivum L.) is the world’s leading cereal grain where more than one-third of the world’s population uses it as a staple food [1]. It is a self-pollinating annual plant cultivated worldwide and domesticated 10,000 years ago. World wheat production is based almost entirely on two modern species: common or hexaploid bread wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and durum or tetraploid wheat (T. turgidum subsp. durum, 2n = 4x = 28, AABB).

It is one of the key staple crops for global food security, providing more than 35% of the cereal calorie intake in the developing world, 74% in the developed world, and 41% globally from direct consumption [2]. In Ethiopia, wheat is used for making bread, porridge (genfo), local beer (tela), roasted grain (kolo), and boiled grain (nifro) [3]. In Ethiopia, the crop is grown at an altitude ranging from 1500 to 3000 m above sea level, between 6 and 16°N latitude, and 35 and 42°E longitude. However, the most suitable agro-ecological zones are between 1900 and 2700 m.a.s.l. [4]. Soil types used for wheat production vary from well-drained, fertile soils to waterlogged heavy vertisols [5]. The major wheat-producing areas in Ethiopia are located in Arsi, Bale, Shewa, Ilubabor, Western Hararghe, Sidamo, Tigray, Northern Gondar, and Gojam zones [4].

The development of an improved variety capable of producing a better yield under various agro-climatic conditions depends upon the amount of genotypic variability present in a population for the traits [6]. According to Rahman et al. [7], the first step in developing varieties is assessing the genetic variability of available genotypes for the characters of interest. Grain yield is a complex quantitative trait influenced by several yield-contributing traits [8]. Hence, the selection for the desirable genotypes should be based not only on the yield alone and but also on the other yield components. Direct selection for the yield is often misleading in wheat because the wheat yield is polygenically controlled. Therefore, effective utilization of genetic stock in crop improvement and knowledge about genotypic and phenotypic correlations that indicate the degree to which various characters are associated with grain yield are necessary. Correlation and path coefficient analysis could be an important tool to bring information about the appropriate cause-and-effect relationship between the yield and some yield components [9]. Several researchers indicated a positive correlation between the wheat grain yield and yield component traits, such as biomass yield and harvest index, plant height, and thousand seed weight [1012].

Even though the productivity of wheat has increased in the last few years in Ethiopia, still its productivity in Buno Bedele as well as in Ethiopia is very low (2.85 tha−1 and 3.046 tha−1, respectively), as compared to the world’s average (3.53 tha−1) [13]. This is due to many factors, including the unavailability of improved high-yielding and adapted varieties, poor management practices, and biotic and abiotic factors. Therefore, developments of high-yielding varieties resistant to different diseases and adaptable to environments with abiotic stress could be crucial for further grain yield improvement in bread wheat. In this regard, this study was conducted to answer the following questions: (1) Is there sufficient genetic variability among the tested genotypes for further grain yield improvement? (2) What is the degree of association between the grain yield and the other yield components? (3) Which genotypes are the best performed? Several studies have been conducted on bread wheat genotypes in Ethiopia’s different regions and reported genetic variability [1417]. However, this is not enough to bring effective yield improvement because the productivity of this crop was still under potential, especially in the study area, and there was limited information related to the nature and degree of genetic variability on bread wheat genotype, particularly in the study area, generally in the country. It is, therefore, of great importance to assess genetic variability and determine the association among yield and yield-related traits in bread wheat to provide basic information for further improvement of bread wheat yield. Therefore, the overall objective of the study was to assess genetic variability and determine the association among traits in bread wheat genotypes.

2. Material and Methods

2.1. Description of the Study Site

The field experiment was conducted during the 2019–2020 growing season at Gito kebele, Gechi district, Buno Bedele Zone, Southwestern Ethiopia. Geographically, the experimental site is located at 08° 12′0''N latitude and 36° 18′0''E longitude, at an altitude of 2132 m.a.s.l. The experimental site receives an average annual rainfall and temperature of 1850 mm and 21oC, respectively. The soil at the experimental site had a pH of 4.3, 2.56 CEC(cmolc/kg), and deep, well-drained Nitosol [18]. The field was treated with 4.22 tha−1 limestone, and its pH was raised to 6.5, which is ideal for wheat growth. Major crops grown in Buno Bedele and the study areas are maize, wheat, sorghum, barley, and tef.

2.2. Experimental Materials

The experimental materials comprised 47 advanced bread wheat genotypes and two released bread wheat varieties, namely, Wane and Lemu. The advanced bread wheat genotypes received for this research were originally introduced from CIMMYT by Kulumsa Agricultural Research Center (KARC) at nursery stages for variety development. These materials were received from KARC for this experiment at a national variety trial (NVT) developmental stage. The advanced genotypes are homozygous lines. The two genotypes used as checks were Lemu (ETBW 6861) and Wane (ETBW 6130), released by Kulumsa Agricultural Research Center (KARC) in 2016) (Appendix Table 1).

2.3. Experimental Design and Field Management

The experiment was laid out in a simple lattice design (7 × 7) with two replications. The total area of the experiment was 488 m2 (23.8 m × 20.5 m), and the plot size was 3 m2 (2.5 m × 1.2 m). The distance between blocks, replications, and plots was 0.5, 1.0, and 0.5 m, respectively. Sowing was carried out on July 2019/2020 by manual drilling in six rows, which are 20 cm apart and 2.5 meters in length. The seed rate was 150 kg ha−1 or (45g plot−1). Blend fertilizer (nitrogen, phosphorus, sulphur, and boron (NPSB)) with a rate of 100 kg ha−1 (NPSB: 18.9% N, 37.7% P2O5, 6.95% S, and 0.1% B) was applied just before sowing, and 100 kg ha−1 of UREA was applied in the split form at the time of planting and knee stage of the crop before heading.

2.4. Data Collected

Data were collected based on the average of ten randomly selected plants and plot basis. Ten representative plants were selected from four central rows for data collection, leaving the two rows as border effect. Data were collected from the following parameters: days to 50% heading (DH), days to maturity (DM), awn length (AL) (cm), grain filling period (GFP), total tillers plant−1 (TPP), fertile tillers plant−1 (FPP), plant height (PH) (cm), spike length (SL) (cm), number of kernel spike−1 (NKPS), grains spikelet−1 (GPs), spikelets spike−1 (SPS), yield plant−1 (YPP) (g), thousand seeds weight (TSW), biomass yield (BMY) (tha−1), grain yield (GY) (tha−1), and harvest index (HI).

2.5. Statistical Analysis
2.5.1. Analysis of Variance (ANOVA)

The data collected for each quantitative trait were subjected to the analysis of variance (ANOVA). The analysis of variance (ANOVA) was performed using Proc GLM procedures of SAS version 9.3 [19].

The mathematical model for a simple lattice design iswhere Yil(j) is the observation of the ith genotype grown in the kth block of the replication j, μ is the grand mean, Gi is the effect of the ith genotype, Rj is the effect of the jth replicate, Bkj is the effect of the kth block in the jth replicate and, and ik(j) is the intrablock residual.

2.5.2. Estimation of Variance Components

The phenotypic and genotypic coefficients of variation are estimated according to the methods suggested Ref. [20]:

2.5.3. Heritability in a Broad Sense

Broad-sense heritability (H2) was estimated as the percentage of the ratio of the genotypic variance () to the phenotypic variance (σ2p) and estimated on a genotype mean basis as described by Allard (1999) as . The heritability was categorized as low (0–30%), moderate (30–60%), and high (>60%) [21].

2.5.4. Genetic Advance under Selection

Genetic advance was calculated by the assuming selection of superior 5% of the genotypes estimated following the methods illustrated in [22]:

2.5.5. Genetic Advance as % of the Mean (GAM)

GAM was estimated according to Hanson et al. [23]. , where  is the population mean of a character.

GAM was categorized as low (0–10%), moderate (10–20%), and high (>20%) [22].

2.6. Correlation Analysis

Correlation between variables was estimated using the R software (R version 4.1.0, 2021) package “Variability” version 0.1.0 and the method described by Miller et al. [24], and Kashiani and Saleh [25] from the corresponding variance and covariance components as follows:

2.6.1. Path Coefficient Analysis

Using R software “Variability” package version 0.1.0, path coefficient analysis was performed as suggested in Ref. [26] with the phenotypic and genotypic correlation coefficients as :

Path coefficient was rated as very high (>1.0), high (0.30–0.99), moderate (0.2–0.29), low (0.1–0.19), and negligible (0.00–0.09) [27].

2.7. Cluster Analysis

The cluster analysis was performed by using R software. Cut points, where local peaks of the pseudo F statistic joined with small values of the pseudo t2 statistic followed by a larger pseudo t2 for the next cluster fusion, were examined to decide the number of clusters. Then, the 49 genotypes were clustered with the average linkage method.

2.8. Genetic Divergence

Genetic divergence between clusters was determined using Mahalanobis D2 statistics [28]. The D2 analysis was performed based on the R software program’s mean values of all traits. The distance between any two groups will be estimated in matrix notation from the following relationship:

2.9. Principal Component Analysis

The principal component analysis was performed using R software (R version 4.1.0, 2021, package “FactoMineR” version 2.4) and computed using the formula suggested by Saporta and Niang [29]:where f2i is the squared factor scores of each observation and λl is the eigenvalue of the lth component.

3. Results and Discussion

3.1. Analysis of Variance

The present study showed that there were highly significant differences ( ≤ 0.01) among genotypes for days to 50% heading, days to 90% maturity, grain filling period, plant height, spike length, number of spikelets spike−1, number of grains spike−1, total tillers plant−1, fertile tillers plant−1, yield plant−1, thousand seed weight, grain yield, and biomass yield. There were significant differences ( ≤ 0.05) among genotypes for awn length and harvest index (Table 2). A complete ANOVA table is presented in Appendix Table 3. The highly significant variation among genotypes for the majority of quantitative traits indicated the presence of sufficient variability, which provides ample scope for selecting superior and desired genotypes by plant breeders for further improvement.

Similarly, Refs. [3034] reported highly significant differences ( < 0.01) for days to 50% heading, days to 90% maturity, number of productive tillers plant−1, number of spikelets spike−1, plant height, spike length, grain yield plant−1, thousand seed weight, biomass yield, and number of kernel spike−1 in their previous studies on bread wheat genotypes.

3.2. Mean and Range
3.2.1. Days to Heading, Grain Filling Period, and Maturity

Days to 50% heading ranged from 61 (BW 174466) to 76 (ETBW 9066) and days to 90% maturity from 98(BW 174466) to 119(ETBW 9066), with a mean of 68 and 106, respectively (Table 4). The maximum days to 90% maturity was recorded for genotype ETBW 9066 (119 days), followed by BW 174467, ETBW 9594 (118 days), ETBW 9635 (116 days), ETBW 9627 (113 days), ETBW 9185 (113 days), and BW 174465 (112 days), while the minimum days to 90% maturity was recorded for BW 174466 (98 days) followed by ETBW 9653(99 days) and ETBW 9653(99 days). The grain filling period ranged from 34 days (BW 174462) to 46 days (ETBW 9594). The long grain filling period was recorded for ETBW 9594 (46 days) followed by ETBW 9185, ETBW 9627, and BW 174468, which had (45 days), whereas the shortest grain filling period was recorded for BW 174462 (34 days) followed by BW 174464 (35 days) and ETBW 9304 (35 days) (Table 5). This indicates that the tested genotypes were in the early to medium maturing category. Therefore, genotypes characterized by early heading and early maturity could be promising, especially when the breeding objective is related to developing varieties resistant to drought stress. In agreement with this result, variation in days to 50% heading and days to 90% maturity has been reported by several authors [31, 35, 36].

3.2.2. Plant Height, Spike Length

Plant height ranged from 65.28 cm (BW 174465) to 87.275 cm (ETBW 9619) with a mean of 74.63 cm. The highest plant height was recorded for ETBW 9619 (87.23 cm) followed by ETBW 9626 (81.24 cm), ETBW 9615 (80.58 cm), ETBW 9594 (80.49 cm), and ETBW 9185 (80.21 cm), whereas the lowest plant height was recorded for BW 174465 (65.28 cm) followed by ETBW 9086 (66.19 cm), ETBW 6130 (67.98 cm), BW 174468 (69.55 cm), and ETBW 9128 (69.75 cm). Spike length ranged from 6.65 cm (ETBW 6130) to 9.36 (ETBW 9594) with a mean value of 8.29 cm. The highest spike length was recorded for ETBW 9594 (9.36 cm) followed by ETBW 9627(9.345 cm), ETBW 9605 (9.27 cm), ETBW 9623 (9.135 cm), and ETBW 9089 (8.9 cm), whereas the lowest spike length recorded for ETBW 6130 (6.65 cm) followed by BW 174116 (6.675 cm), ETBW 9617 (7.21 cm), and BW 174460 (7.415 cm) (Table 5). In agreement with this result, several authors reported similar ranges of variation in bread wheat for spike length and plant height [16, 32].

3.2.3. Number of Spikelets Spike−1, Kernel Spike−1, and Grains Spikelets−1

The number of spikelets spike−1 ranged from 12.6 (ETBW 9647) to 16.95 (BW 174465) with a mean of 14.57. The highest number of spikelets spike−1 was recorded from BW 174465 (26.95) followed by ETBW 9613 (16.8), ETBW 9615 (16.7), ETBW 9185 (15.95), and ETBW 9313 (15.85), whereas the lowest number of spikelets spike−1 was recorded from ETBW 9647 (12.6) followed by ETBW 9623 (12.75), BW 174116 (12.9), and ETBW 9128 (13.1). Kernel spike−1 ranged from 30.95 (BW 174116) to 50.45 (ETBW 9185) with a mean of 42.69. The highest number of kernel spike−1 was recorded for genotypes ETBW 9185(50), BW 174468 (50), and ETBW 9626(49), whereas the lowest number of kernel spike−1 was recorded for BW 174116 (31) followed by ETB W9647 (34) and ETBW 9623(35). The number of grains spikelet−1 ranged from 2.5 (ETBW 9128) to 3.75 (ETBW 9613) with a mean of 3.257. The highest number of grains spikelets−1 was recorded for ETBW 9613 (3.75) followed by ETBW 9617 (3.675), BW 172070 (3.6), and BW 174469 (3.55), whereas the lowest number of grains spikelets−1 recorded for ETBW 9128 (2.5) followed by BW 174463 (2.65), ETBW 9651, and ETBW 9089 (2.9) (Table 5).

3.2.4. Grain Yield, Yield Plant−1, Biomass Yield, and Harvest Index

Grain yield ranged from 2.435 tha−1 (ETBW 9102) to 5.0915 tha−1(ETBW 9613) with a mean value of 3.33 tha−1. In line with this result, Refs. [16, 31] reported a wide variation in the grain yield per hectare, which ranged from 2.115 tha−1 (Menze) to 5.955 tha−1 (Alidoro) and 2.5883 tha−1(ETBW 8511) to 4.6833 tha−1 (ETBW 7871) in bread wheat. The highest grain yield was recorded for genotypes ETBW 9613 (5.0915 tha−1) followed by ETBW 9605 (4.8857 tha−1), ETBW 9626 (4.8846 tha−1), ETBW 9594 (4.8018 tha−1), ETBW 9185(4.6444 tha−1), BW 172070 (4.1782 tha−1), ETBW 9619(4.1243 tha−1), ETBW 9304 (3.9992 tha−1), ETBW 9636 (3.8383 tha−1), and BW 174116 (3.7198 tha−1). The above genotypes were high yielders than the other genotypes including checks (Wane (3.2644 tha−1) and Lemu (3.511 tha−1)). The lowest grain yield was recorded for ETBW 9102 (2.435 tha−1) followed by ETBW 9635 (2.6213 tha−1), BW 174467 (2.622 tha−1), and ETBW 9624 (2.6567 tha−1) (Table 5). Therefore, the above genotypes performed better than all other genotypes including check varieties could be selected for further evaluation to develop high-yielding varieties. Generally, the difference in the mean performance of genotypes for yield and yield-related traits indicates the difference in genetic potential of the genotypes, which confirms the presence of a wide range of genetic variability among genotypes.

Yield plant−1 ranged from 1.02 g (BW 174116) to 1.885 g (ETBW 9185) with a mean of 1.496 g. The highest yield plant−1 was recorded for ETBW 9185 (1.885 g) followed by ETBW 9615 (1.865 g), ETBW 9613 (1.83 g), ETBW 9651 (1.825 g), and ETBW 9594 (1.8), whereas the lowest yield plant−1 was recorded for BW 174116 (1.02 g) followed by ETBW 9647 (1.03 g) and ETBW 9653 (1.11 g). Biomass yield showed a wide range of variation from 6.138 tha−1 (ETBW 9647) to 12.407 tha−1 (ETBW 9185) with a mean of 8.819 tha−1. The highest biomass yield was observed for ETBW 8185 (12.407 tha1), ETBW 9613 (11.563 tha−1), BW 172070 (11.563 tha−1), ETBW 9605 (11.238), ETBW 9601 (11.063 tha−1), and ETBW 9618 (10.188 tha−1), whereas the lowest biomass yield was observed for ETBW 9647(6.138 tha−1) followed by ETBW 9624 (6.39 t ha-14), BW 174470 (6.662 tha−1), BW 174470 (6.728 tha−1), and BW 174464 (6.99 tha−1). The harvest index ranged from 0.28 (ETBW 9635) to 0.48 (ETBW 9626) with mean of 0.377. The highest harvest index was recorded for ETBW 9624 (0.48) followed by BW 174116 (0.475), ETBW 9595 (0.45), ETBW 9613 (0.44), and ETBW 9605 (0.4385), while the lowest harvest index was recorded for ETBW 9635 (0.28), ETBW 9618 (0.2935), ETBW 9601 (0.295), ETBW 9233 (0.298), and ETBW 9623 (0.305) (Table 5).

3.3. Estimation of Components of Variation

In the present study, the genotypic coefficient of variation (GCV) ranged from 4.28 for plant height to 15.48 for the grain yield. In contrast, the phenotypic coefficient of variation (PCV) values ranged from 4.91 for days to 90% maturity to 22.35 for the grain yield (Table 4). The GCV and PCV values were categorized as low (<10%), moderate (10 to 20%), and high (>20%) [37, 38].The estimates of GCV and PCV values were moderate for the number of total tillers plant−1 (11.07 and 15.6), the number of fertile tillers plant−1 (13.77 and 18.8), yield plant−1 (13.72 and 17.4), thousand seed weight (14.8 and 16.67), and biomass yield (12.32 and 18.93) and moderate GCV and high PCV (15.48 and 22.35) for grain yield (Table 4), respectively. This indicates that the genotype could be reflected by the phenotype and the effectiveness of selection based on the phenotypic performance of these characters. On the other hand, it indicates the existence of substantial variability, ensuring ample scope for their improvement through selection. In agreement with this result, Refs. [31, 34] reported moderate GCV and PCV for the grain yield and above-ground biomass yield.

The lowest GCV and PCV were estimated for days to 50% heading (4.55 and 5.07), day 90% maturity (4.49 and 4.91), grain filling period (6.84 and 8.45), plant height (4.28 and 6.1), spike length (5.91 and 7.87), awn length(5.01 and 9.52), and number spikelets spike−1 (6.7 and 8.04), respectively. This indicates the environmental factors had more influence on the expression of these characters than the genetic factors, suggesting the limited scope for improvement of these traits through direct selection for better-performing genotypes.

3.3.1. Estimation of Broad-Sense Heritability and Genetic Advance

The estimated values of heritability in the present study ranged from 28.54% for awn length to 83.6% for days to maturity (Table 4). Robinson et al. [21] classified heritability values as low (0–30%), moderate ((30–60%), and high (>60%). Based on this, high heritability was observed for days to 90% maturity (83.6%), days to 50% heading (80.6%), thousand seed weight (78.92%), number of spikelets spike−1 (69.48%), grain filling period (65.5%), and yield plant−1 (62.1%) (Table 4). This indicated that the variation observed was mainly under genetic control and was less influenced by the environment, and there will be the possibility of progress from the selection. However, only high broad-sense heritability does not always provide a high prediction of genetic gain to ensure effective selection for improvement. Rather, higher heritability coupled with a higher estimate of GCV and GAM [39].

This indicates that for these traits to be effective for selection to ensure grain yield improvement, they must be supported with higher GCV and GAM estimates. This result was in close agreement with the finding of Ref. [40], which reported high heritability for days to 90% maturity, days to 50% heading, number of spikelets spike−1, number of kernel spike−1, and thousand seed weight in bread wheat. Moderate broad-sense heritability was estimated for spike length (56.37%), number of kernel spike−1 (56.2%), number of fertile tillers plant−1(53.62%), plant height (49.13%), grain yield (47.94%), biomass yield (42.4%), number of grains spikelet−1(38.97%), harvest index (30.27%), and number of total tillers plant−1(50.32%) (Table 4). This indicates that the variation observed was mainly due to the influence of the environment. In line with this result, Refs. [10, 41] reported moderate heritability for grain yield, plant height, spike length, number of kernel spike−1, number of spikelets spike−1, and biomass yield.

In this study, genetic advance as percent of the mean (GAM) ranged from 5.6% (awn length) to 27.1% for (thousand seed weight), which indicates that selecting the top 5% of the base population could result in advance of 5.6% to 27.1% percent over the base population mean (Table 4). In the present study, high heritability coupled with high genetic advance as a percent of the mean was observed for thousand seed weight (78.92% and 27.1%) and yield plant−1 (62.1% and 22.26%), respectively (Table 4). This suggests that environmental factors do not much influence these traits, and substantial improvement for these traits could also be achieved through direct selection. These traits are considered to be governed by additive genes. In conformity with this result, Refs. [30, 32] reported high heritability coupled with high genetic advance as a percentage of the mean for thousand seed weight and yield plant−1. Moderate heritability coupled with high genetic advance as a percent of the mean was observed for the grain yield (47.94%, 22.03), and the number of fertile tiller plant−1 (53.62%, 20.76) indicates that improving these traits through selection would be effective (Table 4). High heritability values associated with low genetic advances as a percentage of mean were observed for days to 50% heading (80.6% and 8.427) and days to 90% maturity (83.6% and 8.455), indicating nonadditive gene expression of these traits (Table 4). Therefore, selection based on these traits may not be effective, but heterosis breeding might be necessary for improving these traits. Similar reports of high heritability coupled with low genetic advances as a mean percentage for days to 50% heading and days to 90% maturity [42].

3.4. Association Analysis
3.4.1. Correlation of Grain Yield with Other Traits at the Genotypic and Phenotypic Level

Grain yield showed a positive and highly significant correlation with plant height, the number of grains spikelet−1, yield plant−1, thousand seed weight, and biomass yield at both phenotypic and genotypic levels. It also exhibited a positive and significant correlation with total tillers and the number of productive tillers at the phenotypic and genotypic levels (Table 6). It showed a highly significant positive correlation with the number of spikelet spike−1, spike length, and harvest index at phenotypic and genotypic levels, respectively (Table 6).

A significant positive correlation of these traits with grain yield suggests that any improvement would result in a substantial increment in grain yield. Likewise, [42] reported that plant height, thousand seed weight, biomass yield, and harvest index showed a significant positive correlation with grain yield at the genotypic and phenotypic levels. Grain yield exhibited a nonsignificant correlation with days to heading, days to maturity, and awn length at both the genotypic and phenotypic levels. In contrast, the grain filling period shows a significant correlation with grain yield at a phenotypic level only, which suggests that an increase in days to 50% heading, days to 90% maturity, awn length, and grain filling period had no significant effect on grain yield increment (Table 6). In agreement with this result, Refs. [42, 43] reported a nonsignificant correlation of grain yield with days to 50% heading and days to 90% maturity.

3.4.2. Correlation between Major Traits at the Genotypic and Phenotypic Level

Days to 50% heading had a significant positive correlation with days to 90% maturity, grain filling period, number of total tillers plant−1, and number of fertile tillers plant−1 at both the genotypic and phenotypic level and a significant positive correlation with the number of grains spike−1, number of spikelet spike−1, and biomass yield at the genotypic level. In agreement with this result, Refs. [4446] reported a positive and highly significant correlation of days to 50% heading with days to 90% maturity, while a negative and significant correlation with harvest index at both genotypic and phenotypic levels. Days to 90% maturity showed positive and highly significant correlation with grain filling period, number of total tillers, number of fertile tillers, number of grains spike−1, and number of spikelet spike−1, and significant positive correlation with biomass yield and yield plant−1 at both genotypic and phenotypic levels. Similar to this result, Ref. [47] reported a highly significant and positive correlation between days to 90% maturity and the number of spikelets spike−1, number of kernel spike−1, and number of grains spikelet−1 at both the genotypic and phenotypic levels.

Plant height showed a positive and highly significant correlation with spike length, yield plant−1, biomass yield, thousand seed weight, and grain yield at both the genotypic and phenotypic levels. This indicates that an increase in plant height leads to increase spike length, yield plant−1, biomass yield, thousand seed weight, and grain yield. The number of grains spikelet−1 showed a positive and highly significant correlation with grain yield at the genotypic and phenotypic levels. The number of grains spikelet−1 showed a positive and highly significant association with biomass yield and a significant positive association with harvest index at the genotypic and phenotypic levels, respectively. The number of kernel spike−1 showed a positive and highly significant association with the number of spikelets spike−1, yield plant−1, thousand seed weight, and biomass yield at both the phenotypic and genotypic levels (Table 6). In agreement with the present findings, Refs. [48, 49] reported that the number of kernel spike−1 showed a significant positive correlation with spike length and thousand seed weight at both the genotypic and phenotypic levels.

3.4.3. Path Coefficient Analysis of Different Traits on Grain Yield at the Genotypic Level

Genotypic path coefficient analysis revealed that biomass yield (rg = 0.88;  ≤ 0.01) exerted the highest positive direct effect on grain yield, followed by harvest index (rg = 0.70;  ≤ 0.01) and thousand seed weight (rg = 0.56;  ≤ 0.01), which showed a good relationship between these traits and grain yield (Table 7). This result indicates that the positive and significant correlation of the biomass yield and the harvest index with grain yield at the genotypic level was due to the direct effect of these characters on the grain yield. Therefore, these traits could be considered the main components of the breeding program, and selection based on these traits may be effective in improving the grain yield of bread wheat. This result is in agreement with the findings of Refs. [36, 42, 43], which reported the highest and most positive direct effect of the biomass yield and the harvest index on the grain yield. The residual genotypic value was low (0.0271), indicating that traits included in the genotypic path analysis explained 97.29% of the variation in the grain yield.

3.4.4. Path Coefficient Analysis of Different Traits on Grain Yield at the Phenotypic Level

Direct effects of one character on the grain yield showed that this character directly affected the grain yield and might be helpful in selection for the improvement of yield. Indirect effects on the grain yield indicated that such traits affected the yield via other component characters. The relationship between yield and yield-related characters through the phenotypic path coefficient reveals that the biomass yield (rp = 0.77;  ≤ 0.01) and the harvest index (rp = 0.72;  ≤ 0.01) exerted the highest positive direct effect on the grain yield at the phenotypic level. In contrast, other traits exert a negligible direct effect on the grain yield (Table 8). These highest direct effects of the biomass yield and the harvest index on the grain yield suggest a true relationship, and selection based on these traits may give better response for the improvement of grain yield. In agreement with the present findings, Refs. [15, 43, 50] reported the biomass yield and the harvest index exerted the highest positive direct effect on the grain yield at the phenotypic level.

Plant height (0.39), spike length (0.28), number of spikelet’s spike−1(0.3), number of kernel spike−1 (0.25), and thousand seed weight (0.48) exerted a moderate to high indirect effect via the biomass yield and negligible direct effect on the grain yield (Table 8). This indicates that the significant positive correlation of these traits with the grain yield at the phenotypic level was due to the indirect effect of these traits via the biomass yield. In agreement with this result, Ref. [43] reported positive indirect effects of plant height and the number of kernel spike−1 on the grain yield. The residual phenotypic value was low (0.0759), which indicated that traits included in phenotypic analysis explained 92.41% of the variation in the grain yield.

3.5. Principal Component Analysis

In the present study, the principal component analysis revealed five principal components with eigenvalues greater than 1 (Table 9). The first five principal components accounted for about 76.7% of the total variation for sixteen studied traits, which showed the presence of genetic variability among the studied bread wheat genotypes.

The first principal component alone explains about 34.41% of total variation due to high positive load from the biomass yield (0.82), yield plant−1 (0.82), thousand seed weight (0.75), number of total tillers (0.68), spike length (0.65), number of fertile tillers plant−1 (0.65), plant height (0.64), number of spikelets spike−1(0.63), number of kernel spike−1 (0.61), days to 90% maturity (0.54), and grain yield (0.53). This indicates their significant importance for this component. On the other hand, other traits are less important to PC1. PC2 accounted for 14.46% of total variation due to high negative load from days to 90% maturity (−0.72), days to 50% heading (−0.63), and grain filling period (−0.53) and positive load from plant height (0.47) and grain yield (0.4). The third PC accounted for about 11.42% of total variation due to a high positive load from the harvest index (0.7) and the grain yield (0.6) and negative load from awn length (−0.47), the number of total tillers plant−1 (−0.42), and number fertile tillers plant−1 (−0.39). The fourth PC accounts for about 8.79% due to a positive load from the number of grains per spikelets (0.50), fertile tillers plant−1(0.39), harvest index(0.37), number of total tiller plant−1(0.36), and negative load from yield plant−1(−0.41). Lastly, PC5 contributed 7.63% for total variability due to load positive from grain filling period (0.47) and awn length (0.38), and negative load from grains spikelets−1(−0.47) (Table 9).

The first three PCs accounted for a cumulative 60.3% of the total variation, which indicates much of the variability that occurred originated from traits included in these PCs. The most important yield and yield attributing traits were present in the first three PCs. Generally, maximum contributions toward the variability of genotypes were due to days to 90% maturity, days to 50% heading, grain filling period, biomass yield, yield plant−1, thousand seed weight, and grain yield. Similar works have also been reported by Dargicho et al. [14], who found that three principal components (PCs) out of 12 contributed 63.2% of the variation. Ref. [51] studied 26 genotypes and reported that five PCs contributed 91% of the variation. The traits far from the PCA biplot origin contributed more to the total variability. Accordingly, the primary traits contributing more to total variability were days to maturity, days to heading, plant height, thousand seed weight, biomass yield, grain yield, and yield plant−1.

And the traits close to the biplot origin that contributed less to total variability were awn length, harvest index, grain filling period, number of grains spike−1, number of grains spikelet−1, and number of fertile tillers plant−1. In addition, the traits nearest to the x-axis contributed more for PC1, whereas those nearest to the y-axis contributed more for PC2 (Figure 1).

3.6. Cluster Analysis

The cluster analysis based on the unweighted pair group method with arithmetic means (UPGMA) clustering method from the Euclidean distance matrix grouped 49 bread wheat genotypes into five major clusters of 2 to 21 genotypes (Figure 2). Cluster II (42.857%) was the largest cluster of 21 genotypes. Cluster I (20.408%) was the second-largest cluster and consisted of 10 genotypes. Cluster IV (18.367%) was the third largest cluster containing nine genotypes, whereas cluster V (4.0816%) was the least and followed by cluster III (14.285%), consisting of 2 and 7 genotypes, respectively (Table 10).

3.6.1. Cluster Means Analysis of the Major Contributed Traits for Variability

In the current study, cluster IV was characterized by having the lowest cluster means for days to 50% heading (64.28) and days to 90% maturity (100.7), which indicates the presence of early heading and early maturing genotypes, which were represented in this cluster (Table 11). Thus, further evaluation of members of this cluster is needed to develop promising early maturing genotypes. Cluster I exhibited the lowest grain filling period (36.67), whereas the highest mean for the grain filling period (43.85) was from cluster I. The highest mean for the number of grains spike−1 (45.88) and the number of spikelets spike−1 (15.15) was recorded from cluster I, whereas the highest mean for the grain yield (3.82) and the biomass yield (9.542) was recorded from clusters V and I, respectively.

This result indicates sufficient scope for genotypic improvement through hybridization between the genotypes from divergent clusters. The genotypes in clusters III and IV might not be selected for further evaluation due to the undesirable character of low yield comparably, while the members of these clusters can be selected for the traits with greater mean values. Generally, cluster I was characterized by having the highest mean values for spike length, number of grains spike−1, number of spikelets spike−1, and biomass yield. In contrast, cluster V was characterized by having the highest mean values for the number of total tillers plant−1, number of fertile tillers plant−1, number of kernel spikelet−1, thousand seed weight, harvest index, and grain yield (Table 11). Therefore, genotypes in these clusters could be used as a parent for crossing to develop high-yielding bread wheat varieties.

3.7. Analysis of the Genetic Divergence
3.7.1. Intra- and Intercluster Distances

The average intra- and intercluster distance D2 values are presented in Table 12. The maximum intracluster distance was recorded in cluster V (12.77), whereas the minimum intracluster distance was recorded in cluster IV (9.2). However, the chi-square test for within (intra) clusters was insignificant for all clusters. This indicates that genotypes within the cluster are less diverged than genotypes of different clusters. On the other hand, the chi-square values for intercluster distances showed that there was a significant acceptable difference between clusters I and IV, clusters III and V, clusters II and III, clusters I and V, and clusters IV and V. The highest intercluster distance was observed between clusters I and IV (D2 = 31.86) and followed by clusters II and IV (D2 = 29.21), clusters II and III (D2 = 28.24), clusters III and V (D2 = 27.92), clusters I and V (26.32), and clusters IV and V (D2 = 26.15) (Table 12). This result indicated the presence of genetically distant genotypes, which could be used in hybridization programs to combine the most desirable traits to develop varieties. Among these, cluster V had a significant distance from cluster I, which was divergent for the grain yield and yield-related traits, indicating crossing members of these clusters could be used for yield improvement. The members of clusters III and IV may not be used due to undesirable lower yields, but they could be used for other traits. Refs. [14, 41, 52, 53] reported that the selection of wheat parents for hybridization should be carried out from two clusters having wider intercluster distances to get maximum variability in segregating generations. The lowest and nonsignificant intercluster distance was observed between clusters II and V (D2 = 18.58 ns), followed by clusters I and III (D2 = 23.93 ns) and clusters III and IV (D2 = 24.10 ns) (Table 12). This indicated that genotypes of these clusters were closely related and crossing of genotypes from those clusters might not give a higher heterotic value.

4. Conclusions

The present study showed highly significant differences ( < 0.01) among genotypes for all the traits except awn length and harvest index, which can be exploited to develop high-yielding varieties. The existence of genetic variability among the studied genotypes for quantitative traits shows the possibility of direct selection of parental genotypes to develop hybrid varieties. Based on the mean performance of genotypes, genotypes ETBW 9613, ETBW 9605, ETBW 9624, ETBW 9594, ETBW 9185, BW 172070, ETBW 9304, BW 174116, ETBW 9066, ETBW 9615, and BW 174466 were performed better than all genotypes including check varieties. Therefore, for grain yield improvement, direct selection of these genotypes can be rewarding. In general, the presence of sufficient variability creates opportunities for further improvement of bread wheat genotypes. The biomass yield followed by the harvest index and thousand seed weight exerted the highest positive direct effect on grain yield indicating that these traits are effective for selection. Therefore, the information generated from this study can be used by breeders to improve their traits of interest. Therefore, it is recommended to evaluate the genotypes for high yielding over location and season to know the stability of the genotypes and to draw a definite conclusion.

Appendix

Planting materials used in the study.

Mean square values from ANOVA.

Mean values of forty-nine bread wheat genotypes for sixteen traits.

Data Availability

The data supporting the findings of this study can be obtained from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this article.