Information on genetic progress achieved over time from a breeding program is absolutely essential to develop effective and efficient breeding strategies. Thirty-seven improved lowland sorghum varieties released between 1976 and 2016 and promising advanced lines were evaluated to estimate the genetic progresses made in 40 years of sorghum breeding in Ethiopia. The study was conducted at 2 environments during 2018 cropping seasons in a randomized complete block design with 3 replications. Records taken on grain yield and yield attributes were subjected to statistical analysis. Combined analysis of variance revealed highly significant differences among the genotypes and the test environments for most of the traits, the G×E interaction effects being significant for grain yield. Regression analysis revealed an increase in estimated average annual rate in grain yield potential of 12.2 kg ha-1 year-1 with annual relative genetic change of 0.60% year-1 over the last 40 years of sorghum improvement. Increasing trends along variety release year were also evident for biomass yield, grain yield production per day, biomass production rate and seed growth rate. Stepwise regression analysis revealed that seed growth rate was the most important character, which greatly contributed to the improvement in grain yield. Grain yield was positively correlated with biomass yield, biomass production rate, grain yield production per day, seed growth rate, and thousand seed weight. It is, therefore, strategically advisable that breeding efforts in the future should give due emphasize traits such as seed growth rate.
Performance of the genotypes
The difference between the 37 sorghum genotypes in the combined analysis was statistically significant for DTE, HWT, GFP, PW, GYHA, HI, GYPDAY and SGR (Table 4). The difference between the genotypes was significant at both locations for these traits except for DTE (Table 3). The average grain yield of all sorghum varieties was 2448.03 kg ha-1, which ranged from 1861.2 kg ha-1 for 76T1#23 to 3190.3 kg ha-1for P9534A/Gambella 1107 (elite genotype in NVT 2016). The recently advanced elite genotype P9534A/Gambella 1107 was the first best yielder (Table 8) among the 37 varieties. The superiority of the higher yielder variety, P9534A/Gambella 1107 represents 1152.7 kg ha-1 or 36.1% increment (Table 5) over the average of the first two older varieties (Gambella 1107 and 76T1#23) released in 1976.
The G x E interaction was highly significant for all traits except for DTF, DTM, FLW, NLPP and GL, indicating the inconsistency of performance of the genotypes over the 2 locations for most of the traits (Table 4). Due to the low G x E Mean square against which they were tested, the difference between the genotypes for DTF, DTM, FLW, NLPP and GL was significant. These five traits also had the lowest variance due to rank change (from 13 to 62 with mean variance of 45.4) (results not presented). For PHT, NTPP, FLL, PE, PL, and NSPP, the G x E interaction was significant, but did not mask the difference between the genotypes; mean square of genotypes was also significant, indicating the difference between the 37 genotypes. The variance of rank changes for these traits was also not very high and ranged between 39 and 86 with a mean of 64.7. Most of the interaction seems to be not that due to lack of correlation but due heterogeneity of variances and the performance over locations for PHT, NTPP, FLL, PE, PL, and NSPP is also relatively consistent. For the majority of the remaining traits, DTE, HWT, GFP, PW, TSW, GYHA, BYHA, HI, GYPDAY, BPR and SGR, where the G x E interaction was highly significant while the difference between genotypes was non-significant. For example, for HWA, GFP, GYHA, HI, and SGR, where the G x E was highly significant. The variance of ranks was high for these traits and ranged between 81 and 127 with mean of 101.5. There was large rank change of the varieties over the two locations for these traits. This means that some varieties performed best at Sheraro only while others performed best only at Miesso for this trait. For example, the 37 genotypes can be classified into four categories according to their adaptation (i) widely adapted genotypes with above-average grain yield at both locations (9, 16, 19, 21, 22, 27, 28,31, 32, 34 and 37), (ii) genotypes adapted to Sheraro (3, 4, 5, 7, 13, 14, 25, 26, 29 and 35), (iii) genotypes adapted to Miesso (8, 11, 12, 20, 23, 24, 33 and 36) and (iv) genotypes not adapted to any of the locations [(1, 2, 6, 10, 15, 17, 18 and 30) (Figure 1).
All the 11 highest yielding and superior genotypes were of 2000 and 2010s release (4th and 5th decades) (6 are advanced genotypes in the pipeline in 2016 and 5 were among those released in the 2000s). If we divide the 37 sorghum genotypes into 5 decades of release, that is, 1970s, 1980s, 1990s, 2000s and 2010s, mean grain yield of the 5 decades was 1872, 1846, 2095, 2190 and 2451 kg ha-1 (Table 5). Genotypes of the recent decade had a yield advantage of 607.1 kg ha-1 or 23.0% (Table 5). It can, therefore, be concluded that grain yield has been improved over the 40 years of the national sorghum improvement for the dry lowlands of Ethiopia. It is also interesting to note that the three varieties released in the second and third decades (1986 and 1998) still give high grain yield at Sheraro (Table 6). The following is the average yield of these four groups when tested at Sheraro, at Miesso and when averaged over the 2 locations
The 11 genotypes with narrow adaptation to Sheraro had a yield advantage of 473.2 Kg/ha (19.9%) over the mean of the 8 genotypes adapted to Miesso, when tested at Sheraro (2848.0 vs. 2374.8 Kg/ha, respectively). The 11 genotypes had a disadvantage of 739.1 Kg/ha (28.6%) as compared to the mean of the 8 genotypes specifically adapted to Miesso, when tested at Miesso (1848.2 vs. 2587.3 Kg/ha), respectively (Table 2). The widely adapted genotypes gave the highest yields at both locations and over both locations; however, their advantage over both the Sheraro-adapted and Miesso-adapted genotypes is the highest when tested over both locations. They were superior to Sheraro-adapted genotypes by 1.2% at Sheraro, to Miesso- adapted genotypes by 4.2% at Miesso and by 17.1 and 10.8% to the two groups when tested over locations. This indicates the importance of testing the varieties across locations and over years to check their stability for use as reliable genetic materials for crop improvement in a specific location.


In this study, the result of significant interaction of variety by location for grain yield is in contrary to the finding of Hailu et al. (2009) who reported no variety x environment interaction for grain yield.
Genetic improvement in sorghum grain yield
Mean grain yields of varieties released in the years such as 1986, 1998, 2000, 2002, 2003, 2005, 2007, 2009, 2011,2012, 2013, 2014 and 2016 exceeded that of the average of the first released older varieties. Moreover, Seredo, Dinkmash by 300.2 kg ha-1(12.84%), Meko by 587.0 kg ha-1(22.36%), Abshir, Gobiye by 217.5 kg ha-1(9.65%), Teshale, Yeju and Birhan by 408.8 kg ha-1(16.71%), Abuare by 338.85 kg ha-1 (14.26%), Hormat by 457.2 kgha-1 (18.33%), Macia, Red swazi and Raya by 194.23 kgha-1 (8.70%) (Table 5). The least and highest increases were 194.3 kg ha-1 (8.71%) and 648.7 kgha-1(24.15%), respectively, over varieties released in 1976 (Table 5). These indicated that there was a gradual increase in grain yield across years of release although this increment was not consistent over the years. For example; varieties 6. Abshir, released in 2000, 10. Birhan released in 2002, 15. Raya released in 2007 and 18. Gedo-1 released in 2007 was among the lowest yielding varieties. Among the advanced genotypes that were in the pipeline in 2016, genotype 30. PU209A/PU304 (ESH-4) was very low yielding (Table 5). This emphasizes the care that should be taken in the process of variety release. The candidate advanced genotypes should be subjected to rigorous multi-environment and multi-season testing before release. When the candidates are tested only in one season over very few sites, the best ones (widely adapted high yielding genotypes) may be missed. For example, among the 11 elite advanced genotypes tested in 2016, only two (27) 2005MI5064 (Argeti) and (30) PU209A/PU304 (ESH-4)) were released, while genotypes 34. IESV23007DL and 36. ETSC 300001, that gave higher yield than varieties 27 and 30 (Figure 1) were missed from release and are testing in national variety trial (NVT). The average rate of increase in yield potential per year of release over the last 40 years period from the slope of linear regression was 12.2**kg ha-1 year-1 (Figure 2). Although there was significant G x E interaction, the direction of yield improvement was positive at both locations; 7.02 kg ha-1 year-1 at Sheraro and 17.3* kg ha-1 year-1 at Miesso.






Overall increase in grain yield over the older varieties was estimated to be 433.9 kgha-1 (17.56%) considering all varieties in the trial, whereas 1152.7 kgha-1 (36.1%) was obtained from variety P9534A/Gambella 1107 (Table 5). Hence, grain yield was found to increase substantially with the release of improved varieties (Figure 2). This agrees with the findings of Karmakar and Bhatnagar (1996) which reported a significant increase in grain yield of new soybean [Glycine max (L.) Merrill] cultivars over the older ones. Likewise, Mihret et al. (2015) reported that a significant increase in grain yield of new sorghum [Sorghum bicolor (L.) Moench] varieties over the early released cultivars. Similarly, Yifru and Hailu (2005) in tef [Eragrostis tef (Zucc.) Trotter], Tibebu (2011) in Chickpea (Cicer arietinum L.)., Kebere et al. (2006) in haricot bean (Phaseolu vulgaris L.) and Wondimu (2010) in food barley (Hordeum vulgare L.) who reported a substantial increment in grain yield of modern cultivars over the older ones.
There was no indication of a yield potential plateau in sorghum over the period studied indicating that the opportunity for breeders to further improve yields exists, and that continued progress towards that end may be expected. The average relative annual gain in grain yield of varieties since 1976 was 0.60% per year, or about 0.20% for the whole period of 40 years (Table 10). Present results indicated that plant breeders have made substantial progress for over the past 40 years in improving the yields of sorghum varieties in Ethiopia although; a yield fluctuation was occurring during the release of some of the varieties, Red Swazi (1991.1kg ha-1) and Raya (1977.1kg ha-1) released in (2007) showed yield reduction while, others showed a yield increment (Table 5).
Genetic improvement in biomass yield, harvest index and plant height
Mean biomass yields of varieties released in the years such as 1986, 1998, 2000, 2002, 2005, 2007, 2009, 2011, 2013, and 2016 exceeded that of the average of the first released older varieties. Moreover, Dinkmash by 431.8 kg ha-1(5.96%), Meko by 1400.1 kg ha-1(17.05%), Abshir, Gobiye by 20.3kg ha-1(0.30%), Teshale, Yeju and Birhan by 742.7kg ha-1(10.20%), Hormat by 1034.9 kgha-1 (13.18%), Macia, Red swazi and Raya by 430.8 kgha-1 (5.95%), Gedo, Melkam, ESH-1 &2 by 497.5 kgha-1 (6.81%), Messay, Chare by 6.58 kgha-1 (8.81%),Melkamash-79 by 669.6 kg ha-1 (8.95%), 2005MI504, 2005MI205 and ETSC300002 by 705.1 kg ha-1 (9.38%) (Table 9). The least and highest increases were 20.3kg ha-1(0.30%) and 1034.9 kgha-1 (13.18%, respectively, over varieties released in 1976 (Table 5). These indicated that there was a gradual increase in biomass yield across years of release although this increment was not consistent over years. The regression of the mean biomass yields of variety on the year of release indicated that there was 13.59** kg ha-1 year-1 average annual rate of increase (Figure 3). The increasing rate was 15.5 and 11.64*** kg ha-1 year-1 at Sheraro and at Miesso, respectively. There was positive significant trend of improvement in biomass over the last 40 years of sorghum improvement. The relative annual biomass yield increment in sorghum varieties was estimated to be 0.20% per year for the last 40 years (Table 10). The present result agrees with the findings of Daniel and Parzies (2011) in the study of genetic improvement of Sesame (Sesamum indicum L.) in Ethiopia. The authors reported higher biomass yield in recently developed varieties than in older ones. Similarly, Hailu et al. (2009) indicated that fodder yield of early-maturing soybean varieties can show a positive trend although not significant, the linear regression of fodder yield of variety means on year of release showed an increasing trend (22.81 kg ha-1 year-1) during a 16 years period. Fano et al. (2016) also indicated that biomass yield in tef was greater in newer varieties and linearly related to variety age which positively and significantly correlated to grain yield. Contrary to these findings, Sinha et al. (1981) reported that breeding had failed to raise the biomass of wheat and the grain yield improvement was solely due to the result of higher harvest index and similarly, Wondimu (2010) reported a non-significant trend in biomass yield in food barley breeding program.
Although the difference between sorghum genotypes was statistically non-significant in the analysis due to the high G x E interaction, this difference was significant at Miesso. At this location, varieties released in the fourth and fifth decades had the highest harvest index (Table 7).
In combined analysis, the highest harvest index was recorded for varieties, 34(IESV23007DL) released in 2016 (42.0%), P9534A/Gambella 1107(41.6%) released in 2016, 21(ESH-2) (40.8%) released in 2009, and 19(Melkam) (38.6%) released in 2009 (Table 8).
Genotypes 1(Gambella1107) and 2(76T1# 23), which were released in 1976 had harvest index of 33.2 and 28.6% (Table 8), respectively; at Miesso higher harvest index is associated with higher grain yield. Harvest index was increased by 0.12% year-1 in the combined analysis (Figure 4). It was increased by 0.25%* at Sheraro, but declined by 0.059% at Miesso (Table 10); this is the consequence of the highly significant G x E interaction. The present finding is similar to the findings of Tafesse et al. (2011) that reported newer sesame varieties had high harvest index; harvest index was improved at an annual rate of 0.97% year-1 over a period of 47 years. Similarly, Wondimu (2010) showed that newer food barley varieties developed in Ethiopia had higher harvest index and the regression slope of the trait over years of release was 0.004. Jin et al. (2010) reported that harvest index increased significantly with year of release, averaging 0.40% per year, rising from 0.31 to 0.38 for soybean cultivars released from 1950 to 2006 in Northeast China. In the same way, yield potential improvement in bread wheat produced marked positive change in harvest index (0.42% year-1) in Ethiopia (Amsal, 1994). In contrary, Fano et al. (2016) in tef, Kebere, et al. (2006) in haricot bean and Tamene (2008) in faba bean who have reported that harvest index was not steadily changed with the year of release of the varieties in the respective crops they investigated. Likewise, Era et al. (2009) also reported that soybean varieties did not show significant differences for harvest index over the period of the genetic improvement. Besides the increment in biomass, there was also a consistent gradual increment in plant height from the older to the newer varieties.


Genetic improvement in growth parameters
PHT was increased by 1.05*** cm year-1 (Figure 5).
The same tendency was observed at Sheraro (0.90 cm year-1) and at Miesso (1.20 cm year-1). NTPP was reduced by 0.003* tillers plant-1 year-1 (Table 10). It was reduced by 0.005* and by 0.001 tillers plant-1 year-1 at Sheraro and at Miesso, respectively. PE was extended (increased) by 0.04* cm year-1. This increment was 0.07* cm and 0.01 cm year-1 at Sheraro and at Miesso, respectively. PW was improved by 0.02* cm year-1 over the two locations. This improvement was 0.003 cm year-1 and 0.03* cm year-1 at Sheraro and at Miesso, respectively. There was a tendency for number of leaves plant-1 to increase over time. It increased by 0.01, 0.02* and 0.02 leaves plant-1 year-1, at Sheraro, Miesso and over locations (Table 10).
PL and FLW had non-significant positive slope; there was a tendency to increase over time, but this increment was statistically non-significant. FLL remained unchanged over time, although there was a tendency of decrease in FLL at Sheraro and in the combined analysis.
A combined analysis averaged over both locations indicated that there was highly significant (p ≤ 0.01) and significant difference among locations, varieties and location x variety interaction effect in plant height (Table 4). These highly significant differences observed among varieties for plant height agrees with different authors (Saleem et al., 2002; Fikru, 2004; Melese, 2005). Among those recently released varieties, minimum plant height was observed in genotypes ‟Seredo (128.0 cm), while, ETSC300002‟ advanced in 2016 (210.5 cm) exhibited maximum plant height. As it was estimated from regression of variety means against year of release, the annual rate of gain, 0.79 cm ha-1 year-1 was different from zero (Table 12). This indicated that yield potential improvement program had markedly affected plant height. Similarly, Yifru and Hailu (2005) reported that plant height was higher for the modern tef varieties than the older ones, even though the relative genetic gain over the past 35 years of breeding, was low (0.4285 cm per year) and was not significantly (p<0.05) different from zero. Similarly, Amsal (1994) indicated that the newest varieties were significantly taller than the older ones but, it did not show relation with year of variety release. On the contrary, Donmez et al. (2001) reported that modern varieties showed significantly decreased plant height and reduced lodging in winter wheat varieties. Similar reports were presented by different researchers in different crops (Wondimu, 2010 in barley; Mihret et al., 2015 in sorghum, Kebere et al., 2006 in haricot bean).
Yield attributes
HWHA increased by 2.47 tones ha-1 year-1. At Sheraro, it decreased by 5.06 tones ha-1 year-1 while, at Miesso it increased by 10.0* tones ha-1 year-1. NSPP increased significantly over the last 40 years by 13.7, 17.4* and 15.5* seeds panicle-1 year-1 at Sheraro, at Miesso and over locations, respectively. However, TSW was reduced by 0.01 g per thousand weight year-1 at both locations and over locations, although this reduction was statistically not significant. Seed length (size) also reduced by 0.001, 0.0001 and 0.001 mm year-1 at Sheraro, at Miesso and over locations (Table 10).
Location mean squares from combined analysis of variance were significant (p ≤ 0.05) for panicle width, number of productive tillers per plant, number of seeds per panicle, grain length and thousand seed weight (Table 4). Likewise, Abebe (1985), Fikru (2004), Melese (2005) and Temesgen (2007), found significant difference in the above yield components traits among tested genotypes in different crops. The mean number of grain (seed) per panicle, panicle width, number of productive tillers per plant and grain length increased significantly over the 40 years period from 2546.50 to 2975.69 (an increase of 14.39%), from 8.67 to 8.84, from 1.03 to 1.04 and from 2.30 to 2.32, respectively. However, over the same period, thousand seed weight decreased from 28.82 to 28.13 cm (by 2.39%).
Generally, older varieties had lower number of grains per panicle and productive tillers than the newer and high yielding varieties. Similar trend was reported by Amasal (1994) who reported 0.438% grains gain in the number of grains per panicle. This difference is reflected in the linear regression coefficient that showed a significant (p ≤ 0.05) increase in number of grains per panicle with annual rate of gain of 429.19 or by 0.61% year-1 as compared to the older variety for the last 40 years in Sorghum varieties improvement program (Table 12). Similarly, Demissew (2010) reported a linear regression of mean which is highly significant increment with a relative genetic gain of 0.61% year-1 grain per panicle. Panicle length showed an increasing trend with years of variety release, which indicated that newer varieties had longer panicle length, higher number of productive tillers per plant and higher number of grains per panicle than the older ones. Linear regression of variety means against year of variety release showed significant (p ≤ 0.05) increment trend in panicle length with relative annual genetic increment of 0.12% (Table 12). Similar reports were published for progress in seed length from soybean breeding in the USA during the period between 1902 and 1977 (Specht and Williams, 1984), in durum wheat improvement (Tafese, 2011) and chickpea breeding in Ethiopia (Tibebu, 2011).

Analysis of variance revealed highly significant differences among genotypes for 1000 seed weight. The values for 1000 seed weight ranged from 28.8 to 32.1 g with a mean value of 30.5 g. Accordingly, Meko (released in 1998) and mean of Gambella 1107 and 76T1#23 (oldest varieties) exhibited maximum and minimum 1000 seed weight of 32.1 and 8.82 g, respectively. In line with this, Tamene (2008) reported annual rate of genetic progress of 8 g thousand seed weight year-1 with relative genetic gain of 1.06 % in fababean. In the same way, Mihret et al. (2015) reported a significant increase in thousand seed weighs with a relative annual rate of gain was 0.94% year-1 over 39 year. However, a highly significant (P ≤ 0.01) decrease in thousand seed weight with a relative annual reduction of 0.96% was obtained in soybean varieties (Demissew, 2010). But Kebere et al. (2006) reported non-significant change in hundred seed weight of haricot bean varieties released between 1972 and 1998.
Results of the present study of the studied genotypes for Panicle width, which ranged from7.33 to10.83cm with mean value of 9.1 and coefficient of variations of 5.71. Highest Panicle width was depicted by genotype ETSC300002 (10.83 cm), while the lowest exhibited by Birhan (7.33 cm). Linear regression coefficient showed increment in Panicle width with annual rate of gain of 0.02*cm year-1 or by 0.004% year-1 relative increase as compared to the older varieties for the last 40 years in sorghum varieties improvement program (Table 12).
Phonological development traits (flowering and maturity)
The analysis of variance of the phonological traits at individual location revealed that there was significant difference (P ≤ 0.05) among varieties in days to flowering. DTE was increased by 0.03* and 0.02* days year-1 at Sheraro and over locations but was reduced by 0.002 days year-1 at Miesso. Days to flowering were increased by 0.06, 0.11 and 0.09 days year-1, at Sheraro, Miesso and over locations. These changes were significant at 10% probability level. Days to maturity was also increased by 0.09, 0.10- and 0.10-days year-1 at Sheraro, Miesso and over locations and these changes were significant at 10% probability level. However, the grain filling period was decreased by 0.11* and 0.04 at Miesso and over locations, but was increased by 0.03 days year-1 at Sheraro (Table 10).
The regression analysis of days to flowering against the year of release indicated a significant annual genetic gain of 0.06 days y-1 at Sheraro and a 0.11 days year-1 gain at Miesso. In addition, days to physiological maturity showed a significant positive trend with the year of variety release in both locations (Table 10). The relative annual genetic gain since 1976 was found to be 0.11% (Sheraro) and 0.16% (Miesso) for days to flowering and 0.09% (Sheraro) and 0.09% (Miesso) for physiological maturity (Table 11).
A combined analysis of variance across the two test locations showed significant (p ≤ 0.01) differences among locations, among varieties and location by variety interaction for days to maturity (Table 4). In line with these results, Tigist (2003) and Ketema (2007) reported the presence of significant difference among genotypes for days to flowering. Most of the recently released varieties were the earliest in flowering and physiological maturity. Variety Raya (released in 2007) is the variety that took longer period (79 to 125) to flower and mature and the highest yielder ‟P9534A/Gambella1107 ‟ is among the varieties that flower and mature early. This shift towards early maturity by decreasing the flowering and maturation time without significantly reducing the grain filling period is important to escape from terminal moisture stress in mid and low altitude areas. In the same way, Wondimu (2010) in barley, Tafese (2011) in sesame, Mihret et al. (2015) in sorghum and Tibebu (2011) in Desi type chickpea in Ethiopia reported a decrease in days to flowering. However, Hailu et al. (2009) observed insignificant yield increment with delayed flowering and maturity in soybean genotypes. Similarly, in a study on haricot bean and durum wheat in Ethiopia, Kebere et al. (2006) and Yifru and Hailu (2005) also reported a non-significant increase in days to maturity.

Biomass production rate, seed growth rates and grain production per day
Biomass production rate, seed growth rate, and grain yield production per day showed significant (p≤0.05) difference among varieties in both locations. At Sheraro genotypes observed to produce the highest biomass production rate, seed growth rate and grain yield production per day of 90.0; 65.1, 27.2 kg ha-1day-1 at Sheraro and 54.8; 54.8, 21.1 kg ha-1day-1 at Miesso, respectively (Table 10). Low seed growth rate might be due to early termination of rain which caused lower biomass yield and grain yield. Most of the older varieties produced a higher biomass production rate, seed growth rate and grain yield production per day than the recent varieties at both locations. The relative annual gain of 0.13% per year for biomass production rate, 0.44% per year for seed growth rate and 1.80% for grain production per day (Table 12) was high, indicating that these characters were effectively and significantly improved as a result of the 40 years period of grain yield potential improvement. This agrees with the investigation of Kebere et al. (2006) on haricot bean and Fano et al. (2016) on tef. This data indicated that, biomass production rate; seed growth rate and grain production per day from the five decades of plant breeding and selection was increased by 11.61, 18.86 and 34.17%. The annual genetic gain as estimated from the regression coefficient was 0.22, 0.20 and 0.09 kg ha-1 day-1 year-1 for biomass production rate, seed growth rate and grain yield production per day respectively (Table 10).
According to Kusmenoglu and Muehlbauer (1998), increased seed yield has been obtained through development of cultivars with shorter vegetative and generative growth periods, and greater rates of crop and seed growth. In agreement with the finding of the present study, Pedersen et al. (1998) and Demissew (2010) reported significant increase in seed growth rate in soybean. Similarly, Kebere et al. (2006), Tamene (2008) and Tibebu (2011) reported significant increase in seed growth rate and biomass production rates of haricot bean, fababean and chickpea varieties released in Ethiopia, respectively.
However, Wondimu (2010) showed that there was a significant increase only in seed growth rate in barley. In contrast, Yifru and Hailu (2005) observed non-significant increases in both seed growth rate and biomass yields of tef genotypes over 35 years of breeding and selection. From this, it can be concluded that substantial improvement was apparent in the rate of biomass production, seed growth rate and grain production per day due to grain yield improvement.
BYPR increased by 0.07, 0.37**, and 0.22 kg ha-1 day-1 at Sheraro, Miesso and over locations, for the last 40 years. GYPDAY also increased by 0.05, 0.14* and 0.09* kg ha-1 day-1 at Sheraro, Miesso and over locations. SGR increased by 0.14, 0.55** and 0.34 kg ha-1 day-1 at Sheraro, Miesso and over locations, respectively (Table 10). The highest improvements were achieved at Miesso.
Correlation analysis
Biomass yield, harvest index and plant height
The correlation coefficients of grain yield, thousand seed weight and biomass yield with all the traits studied are presented in Table 12. The results of correlation analysis indicated that grain yield showed a highly significant (p ≤ 0.01) and positive association with biomass yield (r= 0.57***), harvest index (r=0.86***), grain yield production per day(r=0.94**), biomass production rate (r=0.37***), and seed growth rate(r=0.92***). Moreover, biomass yield showed significant positive correlation with biomass production rate (r=0.84***), grain yield production per day (r=0.61***) and seed growth rate(r=0.55***), but non-significant association with all other traits (Table 14) In agreement with the present study, Singh et al. (1990) on chickpea found that biological yield and harvest index had significant positive association with grain yield and therefore simultaneous selection for these two traits would lead to high seed yield. Conversely, Yifru and Hailu (2005) on tef, (Kebere et al., 2006) on haricot bean, Tamene (2008) on faba bean, Hailu et al. (2009 and Demissew (2010) on soybean found a highly significant positive correlations between grain yield and biomass yield, but no significant correlation between grain yield and harvest index. However, Amsal (1994) on bread wheat and Wondimu (2010) on food barley reported a significant and positive association between harvest index and grain yield and a non-significant association between biomass and grain yield. The association between grain yield and plant height was positive (r=0.24). Different authors (Wondimu, 2010 on food barley and Jin et al., 2010 on soybean) found a significant correlation of grain yield and plant height. In contrary, Yifru and Hailu (2005), Kebere et al. (2006), Tamene (2008), Hailu et al. (2009) observed no relation between grain yield and plant height in tef, haricot bean, faba bean, and soybean respectively. In general, grain yield in the modern varieties appears to be associated more with a higher partitioning efficiency to the grain sink than the production of a higher biomass. This indicated that partitioning efficiency may serve as an index for identifying varieties with higher seed yield.
Yield components
Significant and negative correlation was observed between grain yield and number of tillers per plant, grain yield and panicle length, while the association of grain yields with panicle width, head weight and number of seeds per panicle were positive (Table 14). This indicates that these characters are important traits used as indirect selection criteria in breeding for improving grain yield in sorghum. Similar results were also reported by Saleem et al. (2002) in chickpea found that there was significant and negative association of grain yield with spike length. Likewise, Majumder et al., (2008); Degewione and Alamerew (2013) reported that positive and non-significant correlation of grain yield with number of kernel per spike and 1000 seed weight in chickpea.
Phonological traits
Days to flowering and days to maturity showed a non-significant and positive association (r = 0.08) and (r=0.03) with grain yield respectively (Table 14). This is in agreement with the investigation of Tibebu (2011) who reported a non-significant association of days to flowering and days to maturity with grain yield in chickpea. However, Singh et al. (1990) in chickpea, Hailu et al. (2009) and Demissew (2010) in soybean reported strong positive correlations of grain yield with days to flowering and days to maturity. In contrast, Fikru (2004), reported a negative association between days to flowering and days to maturity with grain yield in wheat. According to Amsal (1994) in wheat and Fano et al. (2016) in tef, days to flowering and days to maturity were association poorly with grain yield. Kebere et al. (2006) in haricot bean reported lack of correlation between grain yield and these phonological traits.
Productivity traits
Biomass production rate (r=0.37**), grain yield production per day (r=0.94**) and seed growth rate (r=0.92**) showed a highly significant (P ≤ 0.01) positive relation with grain yield (Table 14). This clearly showed that improvement in these traits was markedly concurrent with the yield improvement achieved in the past and can further be exploited in future breeding program. Similar results were reported by Kebere et al. (2006) in haricot bean and Demissew (2010) in soybean. Likewise, DeBruin and Pedersen (2009) found positive relation between grain yield and crop growth rate during seed set, seed growth rate, grain yield production per day. Moreover, positive correlation between grain yield with grain yield production per day and biomass production rate (Yifru and Hailu, 2005) on Tef were reported.
Yield components
Non-significant and negative correlation was observed between grain yield and plant height and grain filling period, while the association of grain yields with panicle length and panicle width was positive. This indicates that these characters are important traits used as indirect selection criteria in breeding for improving grain yield in sorghum. Similar results were also reported by Saleem et al. (2002) in chickpea found that there was significant and negative association of grain yield with spike length. Likewise, Majumder et al., (2008); Degewione and Alamerew (2013) reported that positive and non-significant correlation of grain yield with number of kernel per spike and 1000 seed weight in chickpea.
Phonological traits
Days to flowering and days to maturity showed a non-significant and positive association (r = 0.20) and (r=0.13) with grain yield respectively (Table 14). This is in agreement with the investigation of Fikru (2004), reported a negative association between days to flowering and days to maturity with grain yield in wheat, Tibebu (2011) who reported a non-significant association of days to flowering and days to maturity with grain yield in chickpea. However, Singh et al. (1990) in chickpea, Hailu et al. (2009) and Demissew (2010) in soybean reported strong positive correlations of grain yield with days to flowering and days to maturity. According to Amsal (1994) in wheat and Fano et al. (2016) in tef, days to flowering and days to maturity were association poorly with grain yield. Kebere et al. (2006) in haricot bean reported lack of correlation between grain yield and these phonological traits.
Productivity traits
Biomass production rate(r = 0.22NS) revealed a non-significant and negative association with grain yield, while grain yield production per day (r=0.92***) and seed growth rate (r =0.95***) showed a highly significant (P ≤ 0.01) and positive relation with grain yield (Table 14). This clearly showed that improvement in these traits was markedly concurrent with the yield improvement achieved in the past and can further be exploited in future breeding program. Similar results were reported by Kebere et al. (2006) in haricot bean and Demissew (2010) in soybean. Likewise, DeBruin and Pedersen (2009) found positive relation between grain yield and crop growth rate during seed set, seed growth rate, grain yield production per day. Moreover, positive correlation between grain yield with grain yield production per day and biomass production rate (Yifru and Hailu, 2005) on Tef were reported.
Stepwise regression analyses using grain yield as dependant variable and other traits independent variables indicated that, seed growth rate, harvest index, grain production per day and biomass yield production rate are traits which contributed to gain in grain yield. Particularly, 49% of the variation in grain yield was explained by seed growth rate and was the single most important trait that contributed most to the variation in grain yield among others, where as 27, 23 and 20% variation in grain yield were contributed by grain yield production per day, biomass production rate and harvest index respectively (Table 13). This illustrates that the improvement in grain yield was achieved by combination of different factors.
According to Wondimu (2010) results of a stepwise regression analysis of grain yield on selected yield components revealed that harvest index, biomass yield and seed yield per day altogether accounted for 46, 73 and 74% of the variation in grain production respectively. Amsal, (1994) also reported number of grains per meter square alone accounted for most of the variation (>68%) in grain yield while number of gains per meter square, 1000-seed weight, plant height, biomass yield collectively contributed for more than 93% variation in wheat grain yield. About 96% of the variation in fababean grain yield was explained by economic (seed) growth rate, whereas economic growth rate, number of pods per plant, harvest index and biomass together accounted for 99% of the variation in grain yield (Tamene, 2008).