Introduction
Wheat is a commodity with a high market value that generates income for farmers in Ethiopia, the largest wheat producer in sub-Saharan Africa. Two types of wheat are predominantly grown in Ethiopia: tetraploid wheat, indigenous to the country, and hexaploid wheat, a recent introduction to the country (Getachew and Worede, Reference Getachew and Worede1991; Alamerew et al., Reference Alamerew, Chebotar, Huang and Roder2004). The majority of tetraploid wheat species grown by farmers are mixtures of landraces varying in botanical and morphological features. For millennia, farmers have preferred to grow mixtures of tetraploid wheat landraces to add variety to their diet and to reduce the risks of losses due to new disease or pest outbreaks or due to unusual environmental conditions (Pecetti and Damania, Reference Pecetti and Damania1996) because of useful alleles they possess as compared to hexaploid wheat. Ethiopia is a centre of diversity for the cultivated tetraploid wheat (2n = 4x = 28) (Vavilov, Reference Vavilov1929; Abate, Reference Abate2018; Brasesco et al., Reference Brasesco, Asgedom, Sommacal and Casari2019). Despite this potential, Ethiopia remains a net importer of wheat due to the huge gap between production and consumption (Abate, Reference Abate2018; Brasesco et al., Reference Brasesco, Asgedom, Sommacal and Casari2019) emanating from very low national yield (Haile et al., Reference Haile, Hammer, Ayele, Nachit and Röder2013a, Reference Haile, Hammer, Badebo, Singh and Roder2013b) and increased demand for wheat (Zegeye et al., Reference Zegeye, Alamirew and Tolossa2020). The mean productivity of wheat in Ethiopia is 3 tons per hectare (t/ha) (CSA 2021) which is below the attainable yield for the crop which can be up to 5 t/ha (Alemu et al., Reference Alemu, Zegeye, Kassa, Asnake, Solomon and Asefa2019; Zegeye et al., Reference Zegeye, Alamirew and Tolossa2020; Nigus et al., Reference Nigus, Shimelis, Mathew and Abady2022). Hence, low productivity continues to be the major challenge facing wheat production in Ethiopia (Alemu et al., Reference Alemu, Zegeye, Kassa, Asnake, Solomon and Asefa2019). Accordingly, continuous enhancement of wheat yield is crucially needed. This necessitates the development of wheat varieties that are high yielding with the required quality and are stable under biotic and abiotic stresses. This requires a continuous supply of new germplasm as a source of desirable genes and gene variants (Asmamaw et al., Reference Asmamaw, Keneni and Tesfaye2019) for which tetraploid wheat landraces are a valuable source.
Tetraploid wheat has been under cultivation in Ethiopia for thousands of years and has acquired a diverse set of characteristics and enormous genetic variability. However, the diversity present in tetraploid wheat through domestication has not been fully evaluated (Negisho et al., Reference Negisho, Shibru, Id, Ordon and Id2021). Landraces are dynamic population(s) of cultivated species that have historic origin with distinct identity and locally adapted in association with traditional farming systems (Ceccarelli, Reference Ceccarelli2016). Landraces possess wide genetic diversity and are underutilised primary sources of desirable genes for economically important traits (Teklu and Hammer, Reference Teklu and Hammer2008; Haile et al., Reference Haile, Hammer, Badebo, Singh and Roder2013b; Muleta et al., Reference Muleta, Rouse, Rynearson, Chen, Buta and Pumphrey2017). Tetraploid wheat holds diverse alleles for disease resistance that can reverse the existing genetic diversity erosion in established, elite cultivars Muleta et al. (Reference Muleta, Rouse, Rynearson, Chen, Buta and Pumphrey2017). The development of new varieties is critical to enhancing the yield of wheat, and use of landrace populations is a viable strategy to improve yield and yield stability as well as resistance to biotic and abiotic stresses (Abbasabad et al., Reference Abbasabad, Mohammadi, Moghaddam and Kamali2016).
Crop improvement substantially depends on the extent of genetic variability existing within the species and their crop wild relatives. The genetic variation that exists among plant populations is a basic requirement for their efficient improvement and also serves as an evidence to prove whether the population of such plants can withstand unpredictable changes in the environment (Nandwani, Reference Nandwani2019). Crop breeding programmes depend on the availability of large germplasm collections, which are invaluable source of parental strains for hybridization and subsequent development of improved varieties (Asins and Carbonell, Reference Asins and Carbonell1989). More than 7000 tetraploid wheat landraces were maintained at the Biodiversity Institute of Ethiopia (https://ebi.gov.et/). However, only limited portions of the collections were characterized (Negisho et al., Reference Negisho, Shibru, Id, Ordon and Id2021) using morphological markers (Getachew and Worede, Reference Getachew and Worede1991; Belay et al., Reference Belay, Tesemma, Becker and Merker1993; Bechere et al., Reference Bechere, Belay, Mitiku and Merker1996; Belay et al., Reference Belay, Bechere, Mitiku, Merker and Tsegaye1997; Kebebew et al., Reference Kebebew, Tsehaye and McNeilly2001; Alamerew et al., Reference Alamerew, Chebotar, Huang and Roder2004; Eticha et al., Reference Eticha, Bekele, Belay and Börner2005; Faris et al., Reference Faris, Arnulf, Getachew and Eva2006; Teklu and Hammer, Reference Teklu and Hammer2008; Tsegaye et al., Reference Tsegaye, Dessalegn, Dessalegn and Share2012; Mengistu et al., Reference Mengistu, Kiros and Pè2015; Asmamaw et al., Reference Asmamaw, Keneni and Tesfaye2020). Therefore, more information is needed about phenotypic and genetic variation present in Ethiopian tetraploid wheat landraces through morphological characterization to reveal the potential input of the landraces for breeding (Teklu and Hammer, Reference Teklu and Hammer2008).
Genetic variability among tetraploid wheat genotypes can be estimated based on quantitative traits (Azene et al., Reference Azene, Menzir and Dejene2020). Furthermore, knowledge of the naturally occurring diversity in tetraploid wheat landraces helps to identify diverse groups of genotypes to be incorporated in the breeding programme (Azene et al., Reference Azene, Menzir and Dejene2020). To inform breeding, it is important to estimate heritability and genetic advance (Pandey and Tiwari, Reference Pandey and Tiwari1983). Heritability denotes the proportion of phenotypic variance that is due to genetic reasons (Singh, Reference Singh1990). Genetic advance provides a prior quantitative estimate of the magnitude of the progress that can be achieved through selection (Panse and Sukhatme, Reference Panse and Sukhatme1957). In this study, we assessed the genetic variation, heritability, genetic advance and correlation of phenotypic traits in 196 tetraploid wheat germplasm to inform future breeding programmes.
Materials and methods
Planting materials
A total of 196 tetraploid wheat (Triticum turgidum spp.) genotypes, representing 174 landraces and 22 varieties, collected from different parts of the country were used in this study (Supplementary Table S1). The landraces accessions used were originated from the different Ethiopian wheat producing regions: Shewa, Jima, Bale, Tigray, Wello, Gonder, Gojam, Agaw Awi (Fig. 1). The released varieties and 40 of the landraces were obtained from Debrezeit Agricultural Research Centre (www.eiar.gov.et). The remaining 134 landraces were obtained from Sinana Agricultural Research Centre (http://www.iqqo.org). The landraces were collected by the Biodiversity Institute of Ethiopia (IBE).
Figure 1. Geographical map of Ethiopia, indicating areas of the collection of the tetraploid wheat landraces and field trial sites of the research.
Methods
Each genotype was grown in plots of two rows of 1 metre long and 20 cm inter-row spacing with two replicates per accession in a simple lattice design at Debrezeit and Sinana during 2020. Debrezeit Agricultural Research Centre is located at an altitude of approximately 1900 m above sea level, with latitude of 80 44′ N and longitude of 380 85′ E. Sinana, located at an altitude of 2400 metre above sea level, has a range of mean annual rainfall of 563–1018 mm with minimum and maximum temperature of 7.9 °C and 24.3 °C, respectively. All agronomic recommendations were used as recommended: 100 kg urea and 150 kg DAP per hectare and three times hand weeding was applied. Ten plants were randomly selected and tagged for phenotypic data collection. The data were collected for 11 quantitative morphological traits using the descriptors for wheat (IBPGR, 1985). Days to booting, days to maturity, seed yield and biomass yield were recorded on a plot basis.
Plant height
Height of plant at maturity, measured in cm from ground to top of spike.
Spike length
The length of spike from the base of spike to the tip of spike measured in cm.
Days to booting
Counted as days from sowing to 50% of plants in booting stage.
Days to maturity
Counted as days from sowing to 50% of plants at physiological maturity.
Number of spikelet per spike
The average number of spikelet per spike from five typical spikes randomly selected from a growing accession.
Number of seeds per spikelet
The average number of seeds from a spikelet – obtained from the central portion of the five randomly selected typical spikes.
Number of effective tiller per plant
The number of tillers bearing spike from the five randomly selected plants.
Biomass yield per plot
Dry weight of the above ground wheat per plots taken at harvest.
Grain yield per plot
The grain weight of all plants grown per plot taken from each genotype, moisture content adjusted to 12.5%.
Thousand kernel weight
Thousand grains were counted from each genotypes harvested from a plot and their weight in gram was recorded.
Harvest index
The ratio of dried grain weight per plot divided by above ground biomass at 12.5% moisture content.
Phenotypic data analysis: Analysis of variance (ANOVA) for each location was carried out using the PBIB test in R by considering genotypes and block as fixed and random factors, respectively. In each location, the observed phenotypic response of the ith genotype in the jth replication and lth sub-block was computed using the following model:
where y ijl = the observed phenotype, μ = the grand mean, g i = fixed effect of the ith genotype, y j = effect of the jth replication, bl(j) = random effect of the lth block nested within the jth replication and εijl = random error term.
The combined ANOVA across the locations was executed by considering genotype as a fixed effect and the block, and location as random effects according to the following model:
where Y ijkl = observed response of genotype l and replication j of block k of location i; μ = grand mean; g l = fixed effect of genotype l; r ijk = effect of replication j in location i; l i = random effect of location i that is ~NID(0, δ2e); b ijkl + random effect of block k nested within replication j in location i that is ~NID(0, δ2b); (gl)il = random effect of the interaction between genotype l and location i that is ~NID(0, δ2gl) + εijkl = random residual effect that is ~NID(0, δ2ε).
Homogeneity of the error mean square (MS) was tested from individual ANOVA at Sinana and Debrezeit was checked following the F-max technique of Hartley (Reference Hartley1950) described as: maximum F statistics (F max) = Larger error mean square (MSE)/smaller error mean square (MSE). The error variance is declared as homogenous if the larger MSE is not three times greater than the smaller MSE (Gomez and Gomez, Reference Gomez and Gomez1984). After deciding that the error variances were homogeneous for specific triat, the combined ANOVA was performed. The statistical significance between genotypes is decided based on P-value that corresponds to the F statistics. If P-value is less than the specified alpha (α) level, the null hypothesis is rejected and the difference between the genotypes concluded significant; however, if the P-value is not less than the specified alpha (α) level, the null hypothesis is accepted and we conclude that the genotypes showed statistically no significant difference (https://statisticalpoint.com/anova-f-value-p-value).
The variability of each quantitative trait was estimated by simple statistical measures such as mean, range, phenotypic and genotypic variances and coefficient of variation. Phenotypic, genotypic, environmental and genotype by environment interaction coefficient of variation, broad sense heritability (H 2), and genetic advance were the parameters assessed. The phenotypic and genotypic variation and coefficient of variations were calculated using the formula suggested by (Singh and Chaudhary, Reference Singh and Chaudhary1985) and (Allard, Reference Allard1960).given below:
Genotypic variance (δ2g)
where MSg stands for the mean square of genotype, MSgl represents the mean square due to genotype by environment interaction, l is the number of locations and r stands for number of replications.
Environmental variance (δ2e) = MSe where MSe = combined error mean square.
Genotype by environment interaction variance (δ2gl)
where MSgl = mean square due to genotype by environment interaction and MSe = combined error mean square.
Phenotypic variance (δ2p)
Estimates of coefficient of variation will be obtained as follows.
Phenotypic coefficient of variation (PCV)
where PCV = phenotypic coefficient of variation, δ2p = phenotypic variance and μ = population mean for the trait considered.
Genotypic coefficient of variation (GCV)
where GCV = genotypic coefficient of variation, δ2 g = genotypic variance and μ = population mean for the trait considered.
Environmental coefficient of variations (ECV)
Coefficient of variation due to genotype by environment interaction was computed by the formula
where, δ2gl = genotypic by environment interaction variance and μ = population mean for the trait considered.
Broad sense heritability (H2) and genetic advance
Heritability in broad sense, for the two locations, was estimated based on the formula given by (Allard, Reference Allard1960). H 2 = (δ2g/δ2p) × 100, where δ2p = δ2 g + (δ2gl/l) + (δ2e/rl) where δ2e = error variance, l = the number of locations and r = the number of replications. Expected genetic advance under selection was calculated with the formula of (Allard, Reference Allard1960), at 5% selection intensity, as: GA = (K) (δp) (H 2), where GA = expected genetic advance, K represents a selection differential that varies based on the selection intensity and is equal to 2.056 if one chooses 5% of the genotypes, δp stands for phenotypic standard deviation and H 2 represents broad sense heritability. Genetic advance as percentage of the mean will be calculated as GA (% of mean) = () x100%, where GA = genetic advance and μ population mean for the trait considered.
Correlation coefficient analysis
Pearson correlation analysis of quantitative traits was performed for quantitative traits using R software (version 4.1.1) (R Development Core Team, 2018).
Results
Combined ANOVA
A total of 196 tetraploid wheat (T. turgidum spp.), including 174 landraces and 22 varieties, were assessed for the genetic variation, genetic advance, heritability and correlation of their eleven phenotypic traits at Debrezeit and Sinana Agricultural Research centres during 2020. After checking for homogeneity of the error MS from individual ANOVA at Sinana and Debrezeit during 2020 following the formula of (Gomez and Gomez, Reference Gomez and Gomez1984) described as: F- larger error MS/smaller error MS, the error variances were homogeneous for all the traits studied. Table 1 illustrates the variance results from a pooled analysis of eleven phenotypic traits for 196 genotypes of tetraploid wheat genotypes at Sinana and Debrezeit during 2020.
Table 1. Mean square of combined ANOVA of quantitative traits
NB: SdpSpklet, number of seeds per spikelet, number of spikelet per spike; SpL, spike length (cm); TPP, number of effective tillers per plant; PH, Plant height (cm); TKW, weight of thousand kernels per plant (gm). SY, seed yield (T/ha); BY, biomass yield (T/ha); HI, harvest index; DB, days to booting; DM, days to maturity; gen, genotype; Loc, location; Rep:block, block within replication; SE, standard error; CV, coefficient of variation, and Lsd, least significant difference.
For all traits other than the number of effective tillers per plant, the mean squares resulting from genotypes and genotype by location interaction differed significantly among genotypes (P < 0.001). All traits showed highly significant variation (P < 0.001) across locations.
Patterns of quantitative traits variation
The mean value of the phenotypic traits of the accessions along with their pedigree is given Supplementary Table S1. The respective mean and range values for the number of seeds per spikelet (3 and 2–5), the number of spikelets per spike (18 and 15–21), the number of effective tillers per plant (6 and 3–9), spike length (8.20 cm and 5.57–10.38 cm), plant height (99.86 cm and 81.45–117.86 cm), thousand kernel weight (34.24 g and 23.6–44.9 g), seed yield (3.25 t/ha and 1.33–5.98 t/ha), biomass yield (11.23 t/ha and 3.75–22.5 t/ha), days to booting (73 and 153 61–81), days to maturity (123 and 117–130) and harvest index (33.20% and 13.55–72.21%) were given in Table 2. The coefficient of variation was high for grain yield (31.34), biomass yield (31.03), harvest index (33.95) and the number of effective tillers per plant (23.91) (Table 2).
Table 2. The mean, minimum, maximum and range values of quantitative traits at the entire genotypes level
Remarks: SDpsp, the number of seeds per spikelet; SpPSp, the number of spikelet per spike; SPL, spike length (cm); TPP, the number of effective tillers per plant; PH, Plant height (cm); TKW, weight of thousand kernels per plant (gm). SY, seed yield (t/ha); BY, biomass yield (t/ha); HI, harvest index; DB, days to booting; DM, days to maturity; Min, minimum, and Max, maximum.
Phenotypic and genotypic coefficient of variation
The values of phenotypic coefficient of variation (PCV = the variation due to genotype and environment), genotypic coefficient of variation (GCV = the variation due genotype only), genotype by environment interaction coefficient of variation (GECV = the variation due to the interaction of genotype and environment), broad sense heritability (H 2 = how much a variation in a trait is due to genetic factors) and genetic advance (explains the degree of gain obtained in a trait under a particular selection pressure) are given in Table 3. PCV and GCV below 10%, 10–20% and above 20% were respectively regarded as low, intermediate and high (Burton and DeVane, Reference Burton and DeVane1953). The values of PCV and GCV were low for the number of seeds per spikelet (6.98, 5.06), plant height (8.31, 6.65), moisture content (2.66, 2.11), days to booting (5.25, 3.43) and days to maturity; were intermediate for the number of spikelets per spike (17, 15.64), spike length (14.74, 12.01) and thousand kernel weight (15.1, 11.58) and were high for seed yield (37.18, 24.03), biomass yield (36, 24.18) and harvest index (41.11, 29.82). High PCV (20.79) and intermediate GCV (14.93) values were obtained for the number of effective tillers per plant.
Table 3. Variability, heritability and genetic advance
NB: PH, plant height; DB, days to booting; DM, days to maturity; TKW, thousand kernel weight; TPP, tillers per plant; SYTPH, seed yield (t/ha); HI, harvest index; BY, biomass yield (t/ha); spl, spike length (cm); SDpSp, seed per spikelet; Spkltpsp, spikelet per spike; GV, genetic variance; EV, environmental variance; GxEV, genotype by environment interaction variance; PV, phenotypic variance; GCV, genotypic coefficient of variation; PCV, phenotypic coefficient of variation; ECV, environmental coefficient of variation; GECV, genotype by environment interaction coefficient of variation,H2, broad sense heritability; GA, genetic advance and GAM, genetic advance as a percentage of mean.
The difference between PCV and GCV was 0.85, 1.36, 1.66, 1.82, 1.92, 2.73, 3.52, 5.86, 11.28, 11.82.13.15 for days to maturity, number of spikelet per spike, plant height, days to booting, number of seed per spikelet, spike length, thousand kernel weight, number of effective tillers per plant, harvest index, biomass yield and seed yield respectively. The observed environmental coefficient of variation were high for seed yield (31.38), biomass yield (31.04), harvest index (33.95) and the number of effective tillers per plant (23.9); were intermediate for the number of spikelet per spike (10.54), spike length (13.47) and thousand kernel weight (13.46); and were low for the number of seed per spikelet (7.52), plant height (7.04), moisture content (2.49), days to booting (4.24) and days to maturity (2).
The value of genotype by environment interaction variation was high for seed yield (33.42), biomass yield (30.69) and harvest index (32.01); was intermediate for the number of effective tillers per plant (11.57) and was low for the number of spikelet per spike (5.77), the number of seed per spikelet (4.27), spike length (7.32), plant height (4.99), thousand kernel weight (9.85), moisture content (1.63), days to booting (4.24) and days to maturity (2.97).
Broad sense heritability (H2) and genetic advance
H2 values <40%, 40–80%, and > 80% were categorized as low, medium and high, respectively (Mesele et al., Reference Mesele, Wassu and Tadesse2015). Estimates of heritability (H2) ranged from 41.78% to 84.62% for seed yield and the number of spikelets per spike, respectively (Table 2). High value of broad sense heritability was observed for the number of spikelet per spike (84.62) and medium value of heritability was recorded for spike length (66.44), plant height (64.07), days to booting (42.72), days to maturity (58.65), the number of seeds per spikelet (52.53), the number of effective tillers per plant (51.61), thousand kernel weight (58.85), biomass yield (45.11), seed yield (41.78) and harvest index (52.64).
The values of genetic advance as a percentage of mean >20%, 10–20% and <10% were categorized as high, intermediate and low, respectively (Johnson et al., Reference Johnson, Robinson and Comstock1955). Genetic advance as the percentage of mean was low for the number of seed per spikelet (7.54%), days to booting (4.61%), and days to maturity (4.36%); intermediate for plant height (10.95) and thousand kernel weight (18.27), and high for spike length (20.13), the number of spikelets per spike (29.57), the number of effective tillers per plant (22.06), seed yield (31.94), biomass yield (33.39) and harvest index (44.49). The observed genetic advance and broad sense heritability (H2) were high for the number of spikelet per spike.
Correlation analysis
The result of Pearson correlation coefficient was given in Table 4. Grain yield showed a highly significant (P < 0.001) negative correlation with days to booting (−0.36***) and days to maturity (−0.31***) and a highly significant (P < 0.001) positive correlation with thousand kernel weight (0.43***), biomass yield (0.31***) and harvest index (0.49***). Grain yield, on the other hand, showed a significant positive (0.01) positive correlation with plant height (0.18*), spike length (0.15*) and the number of spikelet per spike (0.15*) and a positive correlation with the number of seeds per spikelet (0.14) and the number of effective tillers per plant (0.05).
Table 4. Pearson correlation of quantitative traits of 174 tetraploid wheat land races and 22 improved cultivars of Ethiopia
PH, plant height; DB, days to booting; DM, days to maturity; TKW, thousand kernel weight; TPP, number of effective tillers per plant; SY, Seed yield (t/ha); HI, harvest index; BY, biomass yield (t/ha); SPL, spike length (cm); SDpSpiklet, seed per spikelet; spkletPSsp, spikelet per spike.
Discussion
Genetic variation, genetic advance, heritability and correlation analysis was carried out for 196 genotypes evaluated at Sinana and Debrezeit during 2020 based on eleven phenotypic traits. The genotypes, locations and genotypes by location interaction showed significant variation for the majority of the traits evaluated (Table 1). There is significant (P < 0.001) variation among the genotypes for all traits other than the number of effective tillers per plant indicating the presence of genetic variation among the genotypes which in turn suggests that selection of lines can be effective in improving both yield and quality traits (Azene et al., Reference Azene, Menzir and Dejene2020). In line with this study, Azene et al. (Reference Azene, Menzir and Dejene2020) reported significant variation among genotypes of durum wheat in Ethiopia. The significant (P < 0.001) variation acorss locations and genotype by location interaction suggest that the significant phenotypic variation among the tetraploid wheat genotypes is influenced by the environmental factors such as weather and farming practices, such as soil characteristics, field management or weather, affect how genes are expressed, which may help to explain the situation (Yao et al., Reference Yao, Liu, Liu and Li2008; Persaud et al., Reference Persaud, Persaud, Gobind, Khan, Subramanian and Corredor2022). Further investigations will be required to ascertain stability of traits over several years to assess their suitability for crossing with other desirable traits in a breeding programme.
PCV and GCV below 10%, 10–20%, and above 20% were regarded as low, intermediate and high (Burton and DeVane, Reference Burton and DeVane1953). PCV is a measure of variation due to genetic and environmental factors and GCV is a measure of the relative variability of a trait due to genetic differences among individuals. The difference between PCV and GCV was high for seed yield, biomass yield and harvest index and low for days to maturity, number of spikelets per spike, plant height, days to booting, number of seed per spikelet, spike length, thousand kernel weight and number of effective tillers per plant implying that seed yield, biomass yield and harvest index were influenced by the environment whereas the remaining traits were mainly due to genetic factors. Arega et al. (Reference Arega, Hussein and Singh2010) reported similar result on days to maturity, plant height and spike length; however, their result disagrees with the present result on the number of effective tillers per plant, the number of spikelets per spike, biomass yield, thousand kernel weight, harvest index and grain yield. Additionally, Abebe and Desta (Reference Abebe and Desta2017) reported similar result on days to maturity, number of effective tillers per plant, biomass yield and harvest index; however, their work disagrees with the present study on spike length, grain yield and plant height. Moreover, Azene et al. (Reference Azene, Menzir and Dejene2020) reported similar result on PCV and GCV values of days to maturity, spikelet per spike, thousand kernel weight and spike length; however, their result disagree with the present study on the GCV and PCV values of other traits. Furthermore, Meles et al. (Reference Meles, Mohammed and Tsehaye2017) reported similar results for plant height and days to maturity and thousand kernel weight; however, their result disagrees with the present result on effective tillers per plant, the number of spikelet per spike, spike length, the number of effective tillers per plant, thousand kernel weight, harvest index, grain yield and biomass yield.
The number of spikelet per spike, the number of effective tillers per plant, seed yield and harvest index gave high values of genetic advance and heritability. High heritability accompanied with high genetic advance is an indication of additive gene effects (Johnson et al., Reference Johnson, Robinson and Comstock1955) and hence, high genetic gain from selection of the number of spikelets per spike would be expected. The work by Arega et al. (Reference Arega, Hussein and Singh2011) on durum wheat Ethiopia agrees with the present study on GAM values of plant height (17.4), thousand kernel weight, and days to maturity while it disagrees with this study for the other traits. Moreover, Mesele et al. (Reference Mesele, Wassu and Tadesse2015) reported similar result on biomass yield, seed yield, spike length, thousand kernel weight and plant height on bread wheat of Ethiopia; however, their report disagrees with the present study on the number of effective tillers per plant, harvest index, seed per spikelet and days to maturity. Furthermore, Azene et al. (Reference Azene, Menzir and Dejene2020) also reported similar to the present study on days to maturity, harvest index and thousand kernel weight; however, their report disagrees with present study on plant height, spikelet per spike, spike length, biomass yield and grain yield.
Seed yield showed positive association with thousand kernel weight, biomass yield, harvest index, plant height, spike length, the number of spikelet per spike, seeds per spikelet and effective tillers per plant. This implies that there might be common gene(s) that control seed yield and these traits, which indicates that improving either one or more of these traits could result in high seed yield (Arega et al., Reference Arega, Hussein and Singh2010). According to Kearsey and Pooni (Reference Kearsey and Pooni1996), the positive association of these traits with seed yield might be due to either the presence of strong coupling of genes or pleiotropic genes controlling the traits in the same direction. In line with present study, Arega et al. (Reference Arega, Hussein and Singh2010) reported that seed yield had significant association with biomass yield, plant height, thousand kernel weight and harvest index at both phenotypic and genotypic level. This result disagrees with the work of (Baye et al., Reference Baye, Berihun, Bantayehu and Derebe2020). Similarly, Azene et al. (Reference Azene, Menzir and Dejene2020) reported that there was a highly significant positive correlation of seed yield with thousand kernel weight, biomass yield and harvest index and a positive correlation with spike length. Grain yield showed a highly significant (P < 0.001) negative correlation with days to booting (−0.36***) and days to maturity (−0.31***). In line with this study, Ayer et al. (Reference Ayer, Sharma, Ojha, Paudel and Dhakal2017) reported non-significant negative correlation of days to booting and days to maturity with seed yield, highly significant positive correlation of seed yield with spike length, highly significant negative correlation with 1000 grain weight and harvest index.
The present study result provided preliminary indications that tetraploid wheat of Ethiopia hold huge genetic variation, which could be used as potential input in the breeding programme. Most of the traits studied showed positive association with seed yield implying that improving either one or more of these traits could result in high seed yield. Further work is needed to evaluate these and additional traits under different environmental conditions to assess their suitability for contributing to future breeding programmes.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S1479262123001089
Acknowledgements
We acknowledge Oromia Agricultural Research Institute for financial support. We also extend our appreciation to Sinana Agricultural Research centre and Debrezeit Agricultural Research centre for providing the germplasm and research area. We also extend our sincere appreciation to Tesfaye Tadesse and Demisew Nigusse for their support in data collection.
Author's contributions
MF designed and implemented the experiment, analysed the data and wrote the manuscript. All authors read and approved the final manuscript.
Funding statement
This study was funded by the Oromia Agricultural Research Institute, one of the regional research institutes in Ethiopia.
Competing interest
None.
Ethical standards
Not applicable.
Consent for publication
Not applicable.
Availability of data and materials
All data produced during this study were included this article.
Open access
