 Research
 Open access
 Published:
Phenotypic characterization and seed viability test in exsitu conserved Ethiopian cultivated barley (Hordeum vulgare L.) landraces
BMC Plant Biology volume 23, Article number: 613 (2023)
Abstract
Background
Nowadays, most of the Ethiopian barley landraces had been lost from farmer’s field and exclusively found exsitu conserved at the Ethiopian Biodiversity Institute (EBI). Those exsitu conserved are generally believed to be representative of the original population and possess high genetic diversity and important unique genes that are useful for tackling the various biotic and abiotic stresses in the face of the current climate change. Thus, this research was aimed at testing the performance of 150 exsitu conserved landraces that had been collected from Arsi and Bale highlands, Southeastern Ethiopia. The landraces were tested at multiple test locations over two years (2021 and 2022).
Results
All the tested landraces showed a good germination rate regardless of their long storage duration. In addition, performance of all the qualitative traits revealed a varying frequency for each character state. For example, most of the accessions (51.3%) had six kernel row numbers (KRN). All the remaining accessions had two rows (28.7%) and irregular KRN with variable lateral florets (20%). Likewise, some of the quantitative traits considered showed a significant variation among the landraces. However, there observed a significant variation for all the interaction effects in some of the traits considered signifying the importance of considering environment effects while targeting genetic selection and improvement of exsitu conserved germplasms. The phenotypic coefficients of variation (PCV) were considerably high to medium in most of the traits considered including seed yield per hectare (SYPH) but with no associated higher genotypic coefficients of variation (GCV). Moreover, all the traits showed a far greater phenotypic coefficient of variation (PCV) to that of genotypic coefficients of variation (GCV) once again suggesting the pronounced effect of environmental factors to the variation. This was far supported by the significantly higher absolute magnitudes in phenotypic correlation compared to their corresponding genotypic correlation in most of the traits. Low estimates of heritability and genetic advance observed in all the traits considered except seed yield per hectare indicate importance of the trait for selection in Ethiopian barley improvement programs. Clustering patterns of the accessions, in narrow sense, revealed the existence of low divergence among the samples.
Conclusion
Ethiopian barley landraces are promising candidates for further yield improvement and conservation. However, further regular testing and screening should be conducted for the exsitu conserved landraces because of the current erratic climate change. In addition, more robust molecular marker systems could be used to clearly reveal the extents of genetic diversity and to facilitate the breeding and conservation of Ethiopian barley landraces.
Introduction
Ethiopia is one of the world’s richest genetic resource centers for various food crops including cultivated barley [1]. The country is home for several indigenous food crops, which purely constitute landraces that have been maintained over centuries by farmers, mainly through traditional cultivation systems. They are genetic resources that are believed to have considerable breeding value due to their coadapted gene complexes with tolerance or adaptation to diseases and environmental constraints [2]. In addition, they are useful in breeding for marginal conditions [3] as they offer genes responsible for a more stable yield over a wide range of environmental conditions [4, 5].
Cultivated barley (Hordeum vulgare L.) is among the top genetic resources for which the country is known to be the secondary center of diversity, according to Vavilov [1] and claimed to be the center of origin, according to Bekele [6] and Negassa [7]. Until recently, its cultivation is purely traditional and hence, exists largely in landrace form [8]. Moreover, it is among the top neglected cereals regardless of its huge potential for subsidizing household food security [9], and supporting national and international breeding programs targeting improved adaptation potential to enhance resilience to drought, diseases and other biotic and abiotic crises [10]. In this regard, since its start in the 1950s, barley research has gained several successes. So far, more than 36 improved food barley varieties with different unique important characters and breeding objectives have been released [11].
However, in recent days, the cultivation of Ethiopian barley landraces is declining to the extent of total genetic loss because of the preference for and replacement with a limited number of modern, genetically uniform cultivars or other crops suited for high input agriculture [12, 13]. For instance, some earlier morphotypes such as hulled barley, smoothawned types, hullless types, many naked and some rare covered forms are no longer found in Ethiopia [12]. Some are only found ex situ conserved at the Ethiopian Biodiversity Institute (EBI), Addis Ababa, Ethiopia, which has been extensively collecting and documenting several landraces from the widely producing corridors of the country. So far, the institute has documented more than 17,000 collections from the major barleygrowing regions and zones of the country, including the ArsiBale highlands [14]. Most of those collections were kept at the Ethiopian Biodiversity Institute (EBI) or former Gene Bank for several decades as an exsitu conservation effort and are completely absent from the farmers’ field. In addition, the landraces have not yet been taken back to the farmers’ fields to be evaluated for their important agronomic traits, an important step for further largescale utilization and conservation efforts of the available genetic resources [15]. As a consequence, there is a paucity of valuable uptodate information on the performance of the Ethiopian barley landraces under the currently changing climatic and soil conditions. However, such information is essential for parental selection in order to develop bestperforming, highly productive and good quality varieties, as well as for production [16,17,18,19], and for planning efficient germplasm conservation and utilization strategies [20, 21].
Therefore, the present study was initiated to assess the overall performance of ex situ conserved barley landraces deposited at the Ethiopian Biodiversity Institute (EBI) for a long period so as to generate comprehensive and wellorganized information on the extent of their germination rate and genetic diversity using agromorphological traits and to study the relationship between yieldrelated traits. The information generated would be used as an important input for further improvement and conservation programs. In addition, it could be used as baseline information to develop national and international plans regarding the characterization and utilization of ex situ conserved landraces.
Results
Germination tests
Germination performance, including germination percentage, mean germination time (MGT), and germination rate per day for the evaluated exsitu conserved barley accessions is presented under Table 1. Most of the accessions showed promising germination percentage, where the highest score of 100% was recorded in three accessions, namely EBL025 (3246), EBL052 (3273), and EBL133 (212844), all of which were collected 39–44 years ago (the first two in 1979 and the third in 1984). The smallest germination percentage (46.2%) was recorded in EBL089 (4462) that was of recent collection (1980) as compared to several other accessions. The range and mean germination time across the accessions showed very minimum variation, the smallest being 9.31 in EBL023 (1723) and the largest being 13.99 in EBL136 (212847). These two extreme accessions had been exsitu conserved for nearly 40 years at the EBI. Likewise, the germination rate showed a minor variation across the accessions where the highest (0.15) was exhibited by EBL124 or 64245 and the smallest (0.07) in EBL083 or 3833 and EBL146 or 215373.
Qualitative morphological traits
Performance of the accessions assessed in terms of the eleven qualitative traits is presented under Table 2. In this regard, most of the accessions (51.34%) had six kernel row numbers (KRN). All the remaining accessions had two rows (28.7%) and irregular KRN with variable lateral florets (20%). Long awns and awn less accessions accounted for nearly equal proportions of the sixKRN accessions at 26.7% and 24.7%, respectively. Of the tworowed barley accessions, those with lateral florets accounting for 19% and were the most common type as compared to the tworowed deficient type, which accounted for 9%. In terms of kernel covering (KC), the majority (65%) had covered grains, while 20% and 15% of the accessions studied, respectively represented naked and semicovered grains. Among all lema color (LC) types, yellow (47%) was found to be the most common, and this was followed by black/gray (25%) lema color, while tan/red (13%) and purple (15%) were fewer. The distribution of lemma awn barb (LAB) was nearly proportionate across the accessions, with specific distributions containing smooth (35%), intermediate (32%), and rough (33%). However, the distribution of lemma type varied across the accessions, where most had lemma teeth (38%) and lemma hair (37%), while the remaining (25%) had no lemma teeth.
With regards to growth habit (GH), most of the accessions (65%) were erect, followed by intermediate (21%) and prostate (14%). Accessions with white (35%) awn color were more frequent as compared to those with yellow (21%), brown (16%), black (15%), and reddish colors (13%). Regarding spike density (SD), a large number of accessions had an intermediate density (45%), which was followed by dense (37.3%), and lax (18.0%). Most of the accessions (61%) had long rachilar hair (LH), while some 39% had short lachilar hair. A large number of the accessions had green stem pigmentation (58%), which was followed by purple (basal only) (22%), and purple half or many (20%). The frequency distribution of the accessions with regard to glumes color (GLC) indicated that white (47%) was the most frequent, and was followed by yellow (24%), brown (16%), and black (13%), which was the least frequent.
Quantitative morphological traits
Descriptive statistics
Descriptive statistics of the quantitative traits measured in the barley accessions are presented in Table 3. In general, the exsitu conserved barley landraces showed a wide range of variability, as evidenced by the wide range for most of the quantitative traits. Accordingly, seed yield per hectare (SYPH) showed the widest range (3,210.8 kg/ha) with an average mean performance value of 2,249.9 ± 9.81 kg/ha. This was followed by plant height (PH), which had a combined average performance of 112.0 ± 0.34 cm with a range of 112.0. The accessions performed well for disease traits such as recovery rate per stand (RPS), with a combined mean performance of 82.0 ± 0.23 and range of 50.00 and net blotch (NB) with mean performance of 75.9 ± 0.14 and range unit of 44.
Analysis of variance (ANOVA)
The analysis of variance for the quantitative traits computed using the data combined over the three locations and two seasons (years) is presented in Table 4. Most of the quantitative traits showed highly significant (P < 0.01) variation over years, locations, and yearbylocation interactions. Likewise, the mean square values in most traits (eleven out of the total eighteen) showed a highly significant (P < 0.01) or significant (P < 0.05) variation among the accessions, accessionbyyear and accessionbylocation interactions. On the other hand, mean square values of only some traits (six out of the total eighteen) showed a significant (P < 0.05) variation for yearbylocationbyaccession interactions.
All of the quantitative traits showed a coefficient of genetic determination (R^{2}) nearly greater than or equal to 0.5 with greater scores in days to maturity (DTM) and days to emergence (DTE) (each 0.8 and 0.7, respectively). Similarly, the coefficients of variation (CV) were moderate for most of the traits and within the acceptable range except in effective fertile tiller (EFT) and rate of infestation (Inf) that had 47.6 and 55.6, respectively (Table 4).
Analysis of components of variance
Estimates of variance components of the quantitative traits computed using combined data is presented under Table 5. Estimate of both phenotypic (δ^{2}_{p}) and genotypic (δ^{2}_{g}) variances showed a wide range of variation (0.0 and 0.1 in leaf width or LW to 1709.9 and 174529.2 in seed yield per hectare or SYPH). Similarly, estimate of variance due to genotypeyear (σ^{2}gy), genotypelocation (σ^{2}gl), and genotypeyearlocation (σ^{2}gyl) interactions showed a wide range of variation when the traits are considered all together. However, the variations in all the traits showed narrow range when disregarding seed yield per hectare (SYPH). Estimate of error (environment) variance (δ^{2}_{e}) also revealed a wide range (0.1 in leaf width or LW to 166971.3 in seed yield per hectare or SYPH).
Likewise, estimates of both phenotypic (PCV) and genotypic coefficients of variations (GCV) showed a wide range of variations (PCV = 6.0 in days to maturity or DTM to 58.82 in rate of infestation or Inf; GCV = 0.0 in leaf width or LW to 5.9 in infestation or Inf). All the traits considered showed a lower estimate of GCV, the highest score being 5.9 in rate of infestation (Inf). On the other hand, four traits such as rate of infestation (Inf), leaf rust (LR), effective fertile tiller (EFT), and single leaf area (SLA) had higher (> 20%) PCV estimates. In addition, PCV estimate is by far greater than the corresponding GCV values in all the traits considered. Similarly, estimate of genotype environment coefficients of variation (GECV) showed a slightly wide range of variations among the traits considered (Table 5).
Estimates of heritability in broad sense and genetic advance
Estimate of heritability in broad sense (H^{2}%) in the quantitative traits considered revealed a medium (44.82 in thousand seed weight or TSW) to lowe (38.51 in DTH and less in others) and a wide range of variation (0.00% in leaf width (LW) to 44.82% in thousand seed weight (TSW). Similarly, estimates of genetic advance (GA) revealed a wide range of variation (0.00 in leaf width LW to 26.62 in seed yield per hectare of land or SYPH). A similar trend of wide range has been shown in genetic advance as a per cent of traits mean (GA as % mean) (0.00 in leaf width or LW to 11.01 in thousand seed weight or TSW) (Table 5).
Analysis of correlation coefficients
Result of the pairwise correlation coefficients between the quantitative traits studied is presented in Table 6. In this regard, considerable number of the traits showed significant (P < 0.05) phenotypic (below diagonal) pairwise correlations. However, only few traits showed significant genotypic pairwise correlations. Seed yield per hectare of land (SYPH), one of the important quantitative traits, is among the traits that showed a significant (P < 0.05) and positive phenotypic correlation with traits such as single leaf area (SLA) (0.06), thousand seed weight (TSW) (0.84), recovery rate per stand (RPS) (0.12), number of grains per plant (NGPP) (0.08) and a significant to highly significant negative correlation with days to emergence (DTE) (0.05) (P < 0.05), and rate of infestation (Inf) (0.07) (P < 0.001). Similarly, it showed a highly significant (P < 0.001) and positive genotypic correlation with only two traits such as single leaf area (SLA) (0.21), and thousand seed weight (TSW) (0.84). Despite the significance level, the extents or magnitude of correlation sounds smaller (r < 0.5) in all those traits and most of the remaining except SYPH vs TSW (r = 0.84), and RPS vs Inf (r = 0.81 for phenotypic and r = 0.80 for genotypic).
Principal Components Analysis (PCA)
PC analysis, conducted using the 18 standardized quantitative traits revealed that the first eight principal axes (eigenvalue ≥ 1.06) accounted for 69% of the total variation (Table 7). The first principal component (PC1) accounted for 14.00% of the total variation and had high contributing factor loadings from thousand seed weight (TSW) (0.50), single leaf area (SLA) (0.48), and seed yield per hectare (SYPH) (0.45). The second PC axis accounted for 12.00% of the total variation and differentiated the accessions largely on the bases of recovery rate per stand (RPS) (0.52), rate of infestation (Inf) (0.49), and days to maturity (DTM) (0.42). The third PC axis contributed 9.00% of the total variation and had greater contributing factor loadings from seed yield per hectare (SYPH) (0.34), and flag leaf length (FLL) (0.31). The fourth and fifth PC axes each accounted for 8.00% and 7.00% (in that order) of the total variation and differentiated the accessions largely on the bases of leaf rust (LR) (0.41), number of grains per plant (NGPP) (0.35), net blotch (NB) (0.33), spike length (SL) (0.44). The sixth, seventh and eighth pcs axes accounted for 7.00%, 6.00% and 6.00%, respectively, of the total variation and had greater contributing factors from flag leaf length (FLL) (0.42) (sixth), leaf rust (LR) (0.36) (seventh), and leaf number (LN) (0.54), days to heading (DTH) (0.44) (eighth) (Table 7).
PCA loading plot showed a loose positive and negative correlation among the traits considered (Fig. 1). For example, single leaf area (SLA), leaf width (LW), flag leaf length (FLL) and seed yield per hectare (SYPH) had a weak positive association with each other. Similarly, rate of infestation (Inf) and recovery per stand (RPS) showed a strong negative association as did effective fertile tiller (EFT) and days to maturity (DTM). On the other hand, three traits such as seed yield per hectare (SYPH), thousand seed weight (TSW), and single leaf area (SLA) showed a strong and positive association.
PCA score plot grouped the accessions nearly into five clusters (Fig. 2). The grouping pattern was not parallel to the specific locations of collections. PCA biplot also revealed a weak association or contribution of most of the traits to the grouping patterns of the accessions (Fig. 3).
Cluster analysis
Cluster analysis was conducted using standardized data in order to have a good picture of the genetic association between the accessions studied. Accordingly, the patterns of grouping of individual accessions revealed eight clusters in which larger number of accessions were grouped under clusters 1 and 3 (C1 and C3; each contained 24 accessions or 16% of the total). Clusters 2 and 7 (C2 and C7) were the second largest (each contained 20 accessions or 13.33%), followed by clusters 8 (C8) (19), 5 (C5) (18) and 4 (C4) (17). The pattern seems poor in revealing specific geographic regions of origin (collection) of the accessions and hence, accessions from different specific localities appeared on the same cluster and vice versa (Fig. 4; Table 8).
The estimate of pairwise generalized square distance between the clusters revealed moderate value and the range is nearly the same in several of the clusters. For example, clusters 1 and 3 (24.01, the largest distance), 6 and 7 (23.21), 5 and 8 (22.67), 3 and 6 (22.62), and 2 and 5 (21.83) revealed nearly closer pairwise distances. Clusters 1 and 3 showed the smallest pairwise distance (8.61). Similarly, intra (within) cluster distance among the accessions in each cluster showed a similar pattern (3.58 in cluster 3 to 4.24 in clusters 4 and 5) except for those accessions on cluster 6 that showed a relatively greater within cluster distance (6.13). Regarding estimate of mean distance, all the clusters are nearly equidistant from each other (14.17 in cluster 3 to 16.87 in cluster 1) except cluster 6 which is a bit distant from the others (mean cluster distance of 20.30) (Table 9).
Discussions
Germination efficiency of the exsitu conserved landraces
While there is no clear distinction between specific storage conditions and the termination of life processes in different species, the storage conditions of different seeds have a significant impact on the termination of their life processes. In this regard [22], reported that there is no clearly specified duration for seed death because it is a gradual and cumulative process in which more and more cells die until certain critical parts of the seed become unable to perform their essential function. However, there is a general understanding regarding significant disparities in storage duration between orthodox and recalcitrant seeds where the first could maintain its moisture content for an extended period and remain viable as compared to the second type,which can only persist for a relatively shorter spanunless special storage practices are considered. Moore [23] reported that seed becomes less vigorous as the neverceasing aging process moves onward toward death; but even long before death, the seed becomes questionable or worthless for planting purpose, especially under field conditions that are not highly favorable for germination and seedling development. In this regard, all the exsitu conserved barley landraces considered exhibited a promising performance and revealed a good germination rate regardless of their long storage duration (36 to 57 years) at the EBI coldroom or longer years of collections (1964—1985) from the field. The result suggests good handling from the EBI which could be further enhanced as it is the only source for some germplasms which are totally missing from the farmer’s field. On the other hand, the present result is partly contradicting the general report by [24] suggesting the negative effects of long storage duration on germination rate and related issues and this could be partly attributed to the different environmental and genetic factors including storage temperature, seed moisture content, and genetic variability. The slight variation among the tested landraces seems normal since seed longevity vary among different genotypes, cultivars, and accessions because of genetic variations [25].
Patterns of genetic variability in the landraces
The landraces studied showed a different pattern of genetic variability with respect to most of the qualitative morphological traits assessed, indicating the significantly stronger effects of selection pressure for various enduse qualities. For example, dominance of sixrowed types in the evaluated landraces partly suggests artificial selection pressure by farmers with the intention of obtaining more yield per plot of land. Similarly, larger number of the accessions with being long kernels and yellow or white color indicate preference of the farmers for different enduse qualities. Dominance of the accessions with erect growth habit indicates preference by the farmers because of their suitability for traditional and/or mechanized harvesting. Similarly, dominance of the accessions with white seed and glumes colors, as compared to the reddish and black types, once again qualify preference by the local community for different end use qualities. There are similar reports regarding the existence of wide variability in terms of frequency of those qualitative traits in Ethiopian barley accessions [13, 26, 27].
In general, the patterns of genetic variability in qualitative morphological traits suggest that the Ethiopian barley landraces, although conserved ex situ over a longer period of time, have higher morphological diversity, which is due to the country’s being the center of origin and having high ecological heterogeneity. This higher genetic variation and the good performance of exsitu conserved landraces are key for selection breeding.
Together with qualitative traits, standard quantitative traits are among the important phenotypic markers that are widely used in breeding and conservation of plants, animals and other organisms. In this regard, most of the quantitative traits used to evaluate the exsitu conserved Ethiopian barley landraces revealed a wide range of variability and wide ranges between the maximum and minimum mean values suggesting the ample variability in the landraces for further breeding work. Moreover, quantitative traits such as seed yield per hectare of land (SYPH), thousand seed weight (TSW), number of grains per plant (NGPP), number of effective fertile tillers (EFT) and disease related traits such as rate of infestation (Inf) and recovery per stand (RPS) are important targets of selection. In addition, analysis of variance (ANOVA) supported importance of those traits and others in targeted selective breeding as they showed a significant variation among the tested accessions. Similar results have been reported on the significant variations of several of the quantitative traits in Ethiopian barley landraces [28,29,30,31].
However, the existing variation could not be fully exploited as intended because of several genetic and environment related factors which are largely explained in terms of interactions between or among genes and interaction between genes and environmental factors. Consequently, the highly significant variation in mean square values for most of the traits considered over test years, locations, and yearlocation interactions observed along with interactions between accessions and test years, and locations signals the necessity of exercising the utmost caution and emphasizes the significance of testing the landraces at multiple locations over a couple of years to clearly indicate the amount of genetic based available variations for further use. There had been similar reports by [7, 13, 21, 32, 33] on different barley landraces from Ethiopia.
Similarly, the estimation of variance components is used to disclose the extent of genetic variation and the effects of interaction, especially interaction with the environment, for further uses. In this regard, the tested landraces showed wide range of variation in both phenotypic (δ^{2}p) and genotypic (δ^{2}g) variance estimates which is in agreement with results reported by [13, 34] on Ethiopian barley landraces collected from farmer’s field. Such significant and wide range of variations show the existence of large variability among the tested exsitu conserved barley landraces in particular and barley landraces in Bale and Arsi zones in general. The significant and wide range of variations due to genotypeyear (σ^{2}gy), genotypelocation (σ^{2}gl), genotypeyearlocation (σ^{2}gyl) interactions and error variance (σ^{2}e) in some of the traits especially in seed yield per hectare of land (SYPH) revealed detectable impact of environmental factors to the variations. Thus, care should be taken while selecting the landraces for yield improvement though it has polygenic inheritance pattern. The high coefficient of genetic determination (R^{2}) detected, particularly, in most of the yield contributing traits, also suggest possibility of identifying superior accessions with respect to the traits.
The insightful effect of environmental factors on the detected variation was further explicated by estimates of genotypeenvironment interaction (δ^{2}gl) which was greater than zero and PCV values which were by far greater than the corresponding GCV values in all the traits. Such huge environment effect could be attributed to the current un predictable climate change which is quite different from the years back conditions when the landraces were collected. As a result, the exsitu conserved landraces require some sort of multiple tastings to either develop adaptive potential though they are expected to have excelled adaptive and unique genes as compared to germplasms under cultivation.
Four traits such as rate of infestation (Inf), leaf rust (LR), effective fertile tiller (EFT), and single leaf area (SLA) had scored higher (> 20%) PCV and lower GCV estimate following [35]. This observation is in agreement with the result of [13, 36]. Similarly, the higher difference between the GCV and PCV estimates in these traits suggest the highly pronounced environmental influence. However, this finding is contradicting the reports of [13, 34] on barley landraces from farmers’ fields and thus, there might be differing environmental conditions from the situation when the landraces had been collected.
Traits heritability and seed yield improvement
The concept of heritability pertains to the proportion of phenotypic variability that can be ascribed to genetic variability.. Its estimation is helpful in predicting the expected progress to be achieved through selection process since it indicates the heritable portion of the total variations which is the point of interest in morphological traitsbased genetic performance analysis.
Its value could be very high (≥ 80%), moderately high (60–79%), medium (40–59%), or low (< 40%) following [37] benchmark. In this sense, all the traits considered except thousand seed weight (TSW) with a medium heritability, had low broad sense heritability (H^{2}) estimates and eventually low GCV. The result implies that selection for the characters require special attention due to high environmental effects. Moreover, most of the traits had a direct link with seed yield per hectare of land (SYPH), an important trait for breeding, though it follows polygenic inheritance and thus, environmental factors and other quantitative traits should be seriously managed under the current unpredictable environmental conditions. Similar result has been reported by [13] on barley landraces from collected from farmers field in Ethiopia stressing the profound effect of environmental condition on traits’ heritability.
Genetic advance (GA) can be delineated as the enhancement of traits genotypic values for the new population that result from selection relative to the base population, under one cycle of selection at a given selection intensity [37]. To this end, estimates of GA for seed yield per hectare of land (SYPH) was 8.43 with the associated expected genetic advance values expressed as a percentage of the genotypes mean (GAM) of 0.37 which is low. GAM signposts the gain that could be expected from selection of the top 5% of the populations. Johnson et al. [38] categorized GAM as low (< 10%), moderate (10–20%), and high (> 20%). In this regard, all of the traits considered including seed yield per hectare of land (SYPH) an important target trait for improvement had low estimate.
Performance of the landraces in terms of pairwise correlation coefficients
The pairwise correlation coefficient analysis determines the extent and degree of the relationship between two characteristics. The association could be attributed to genotypic component, linkage between genes or gene effect [39], or to environmental effects (phenotypic), or both [40]. When determining how strongly traits are correlated, both correlation coefficients are crucial in determining whether selection for a given trait leads to either progress or retrogression, especially when it comes to quantitative traits like yield [41]. In this case, substantial number of the traits considered showed significant phenotypic and genotypic correlations. For example, seed yield per hectare of land (SYPH) is among those traits that showed significant positive and negative phenotypic and genotypic correlations with several traits. Those correlated traits, depending up on the magnitude of correlation, are important in improving the yield of barley landraces.
Patterns of genetic relationship in the landraces
The extent and trends of genetic relatedness in any population or target sample can be revealed using cluster analysis methods, including PCA and cluster analysis. The practical application of PCA lies in its ability to identify the traits that have most contributed to the observed variation within a group of samples or populations. This makes it useful for selecting parental lines during breeding. With regards to its magnitude, traits with coefficients of the eigenvector close to one show a strong influence on a given trait and vice versa [42]. Therefore, traits with higher coefficients, typically 0.6 and above, on the PC axes should be considered more important [43]. Similarly, characters with higher factor loadings contribute more to the lumping together or scattering apart of accessions and thus, are given much attention on choosing the clusters for any desired purpose of breeding or conservation [44]. With this fact, the first eight principal axes (eigen value ≥ 1.06) accounted for 69.00% of the total variation. However, most of the traits in these PC axes had lower factor loadings (< 0.50) except four traits such as leaf number (LN), recovery per stand (RPN), thousand seed weight (TSW) and number of grains per plant (NGPP). These three traits are supposed to play a great role for the divergence and exhibited great influence on the phenotype of the accessions and could be targets of selection breeding. The result is concordant with the previous reports by [9, 25] using barley landrace collections from Ethiopia.
Likewise, cluster analysis is also used to show genetic relatedness between or among the subjects of study. In this regard, the exsitu conserved barley landraces considered formed eight major genetic clusters with a weak trend of association between or among accessions from the same geographic location of collection and vice versa. Falconer [45] reported that variation in origin (geographical separation), ancestral relationship, gene frequency and morphology are the probable sources of genetic diversity. However, it is evidenced that, though genetic diversity is associated with geographical diversity, they are not necessarily directly related. To this end, the weak tendency of association between geographical proximity and genetic diversity of the accessions revealed a moderate divergence among the tested accessions.
Conclusions
In recent days, Ethiopian barley landraces are diminishing at faster rate owing to environmental constraints and research focus of the country. Larger number of the landraces is found exsitu conserved (one of the effective ways of preserving germplasms for longer duration) at the Ethiopian Biodiversity Institute (EBI), Ethiopia. In the present study, the landraces showed a promising germination rate and percentage regardless of their very long storage duration at the Institute. Moreover, the study generated basic information on the extents of their genetic variability that promotes the potentiality and high economic values of ex situ conserved barley landraces and promoted detailed studies at more locations over a couple of years to clearly exploit the actual genetic based variability in the current scenario of climate change. In addition, more robust molecular markers are mandatory to clearly reveal the genetic based variability for further utilizations.
Materials and methods
Experimental materials
The study involved a total of 150 ex situ conserved cultivated barley landrace accessions. The samples were obtained from the Ethiopian Biodiversity Institute (EBI), Addis Ababa, Ethiopia. The accessions had been collected from the ArsiBale highlands before 1986 and were assigned full passport data (Additional file 1).
Seed viability test
Seed germination test was conducted at the Biology Department Laboratory of Madda Walabu University (BaleRobe, Ethiopia). Ten clean barley seeds of each accession were soaked in 75 mL distilled water in a separate sterile petri dish for 24 h. The seeds were then allowed to germinate on a Whatman filter paper at 20 °C for 7 days, following the procedures specified by the International Seed Testing Agency (ISTA) [46]. The experiment was conducted using a complete randomized design (CRD) in two replications.
Field experiment
The field experiment was conducted under rainfed conditions over two years (2021 and 2022) at three locations, namely Madda Walabu University (MWU) Integrated Research Field (Bale Robe, Ethiopia), Sinana Agricultural Research Institute (SARI) (Bale, Ethiopia), and Agarfa Agricultural Technical, Vocational and Education Training (TVET) College (Agarfa, Ethiopia). Description of those areas is presented below (Table 10).
First season seed planting were done from 5—15 June 2021. Similarly, second round seed planting were done from 8—16 June 2022. Harvesting was done from November 10 – 17, 2021 for first year planting and from November 14 – 21 for second year planting.
Experimental design and procedures
The experiment was set up in an incomplete block design called alpha lattice, with two replicates per site. Each accession was grown in four rows forming a plot 2 m long and 1.2 m wide. The distances between the blocks, plots, rows and plants were 1.5 m, 0.6 m, 0.4 m and 10 cm respectively. Sowing was done by hand at the correct depth (3–6 cm) in the moist soil to cover the seeds evenly and thus maintain moist conditions for vigorous and healthy germination and growth.
Data collection and statistical analysis
Physiological seed qualities such as standard germination and seed vigor tests were conducted following the International Seed Testing Agency (ISTA) [46] and Maguire [47]. Accordingly, the germination percentage was given as:
Seed vigor was evaluated by using seed germination rate, which relies on the quantity of normal seedlings and average germination time. Hence, seed germination rate was calculated following [47] which is given as Seed germination rate = ∑n/∑D, where n is the number of seeds germinated on day D out of 100 seeds sown, and D is the number of days counted from the beginning of the test.
Mean germination time (MGT) was calculated following [48]. Accordingly, \(\mathrm{MGT}= \sum_{\mathrm{i}=1}^{\mathrm{k}}{\mathrm{n}}_{\mathrm{i}}{\mathrm{t}}_{\mathrm{i}}/ \sum_{\mathrm{i}=1}^{\mathrm{k}}{\mathrm{n}}_{\mathrm{i}}\), where \({\mathrm{n}}_{\mathrm{i}}\) is the number of seeds germinated at the time i; \({\mathrm{t}}_{\mathrm{i}}\) is the time from the start of the experiment to the i^{th} observation, and k is the time of last germination.
Field performance data was collected for a total of 29 traits, including 14 quantitative and 11 qualitative as well as four diseaserelated traits adopted from the standard barley descriptors [49].
The frequency analysis for the qualitative traits was carried out using MINITAB® Release 19 [50] statistical software. Following Hartley’s Fmax based error variance homogeneity test [51], analyses of variance (ANOVA) were calculated for each site and combined across sites for quantitative and disease traits using the GLM procedure of the SAS software based on the following statistical model:
Where, b = effect of block i, g = effect of genotype (accession) j, l = effect of location k, y = effect of year l, gl = effect of interaction of genotype j by location k, gy = the effect of interaction of genotype j by year l, gyl = effect of interaction of genotype j, by year l, and location k, and e = effect of interaction of genotype by block i, genotype j, location k, and year l.
Locations, years, accessions and their interactions were considered as random variables in the analysis according to [52]. The variance components assigned to the accessions and their interactions were calculated using the VARCOMP procedure of SAS.
The estimation of environmental, genotypic and phenotypic variance components and their coefficients of variation per site and combined across sites was computed based on the methods of [53, 54]. Accordingly,
Phenotypic variance (\({\delta }_{p}^{2}\)) per location \(={\delta }_{g}^{2}+{\delta }_{e}^{2}\) where, \({\delta }_{p}^{2}\) = phenotypic variance; \({\delta }_{g}^{2}\) = genotypic variance and \({\delta }_{e}^{2}\) = environmental variance = error variance
Genotypic variance (\({\updelta }_{\mathrm{g}}^{2}\)) per location \(=({\mathrm{MS}}_{\mathrm{g}}{\mathrm{MS}}_{\mathrm{e}})/\mathrm{r}\) where, \({\mathrm{MS}}_{\mathrm{g}}\) = mean square of genotype; \({\mathrm{MS}}_{\mathrm{e}}\) is mean square of error and \(\mathrm{r}\) is the number of replications
Phenotypic Coefficient of Variation (\(\mathrm{PCV }(\mathrm{\%})\)) per location \(=(\sqrt{{\updelta }_{\mathrm{p}}^{2}}/\mathrm{m})\mathrm\times100\), where, PCV = phenotypic coefficient of variation; m = population mean for the trait considered
Genotypic Coefficient of Variation (\(\mathrm{GCV }(\mathrm{\%})\)) per location = \((\sqrt{{\updelta }_{\mathrm{g}}^{2}}/\mathrm{m})\mathrm{\times}100\) where, \(\mathrm{GCV}\) = genotypic coefficient of variation
Genotypic variance (\({\updelta }_{\mathrm{g}}^{2}\)) combined over location = \({\updelta }_{\mathrm{g}}^{2}=({\mathrm{MS}}_{\mathrm{g}}{\mathrm{MS}}_{\mathrm{gl}})/\mathrm{rl}\), where, \({\mathrm{MS}}_{\mathrm{g}}\) = mean square of genotype; \({\mathrm{MS}}_{\mathrm{gl}}\) is mean square due to genotype by environment interaction; \(\mathrm{l}\) = number of locations; \(\mathrm{r}\) = number of replications
G x E interaction variance (\({\updelta }_{\mathrm{gl}}^{2}\)) combined over location = \(({\mathrm{MS}}_{\mathrm{gl}}{\mathrm{MS}}_{\mathrm{e}})/\mathrm{r}\), where, \({\mathrm{MS}}_{\mathrm{gl}}\) = mean square due to genotype by environment interaction; \({\mathrm{MS}}_{\mathrm{e}}\) = combined error means square (\({\updelta }_{\mathrm{e}}^{2}\))
Phenotypic variance (\({\delta }_{p}^{2}\)) combined over location = \({\updelta }_{\mathrm{g}}^{2}+({\updelta }_{\mathrm{gl}}^{2}/\mathrm{l})+({\updelta }_{\mathrm{e}}^{2}/\mathrm{rl})\),
Phenotypic coefficient of variance (\(\mathrm{PCV }(\mathrm{\%})\)) combined over location = (√δ^{2}p/m) × 100, \(\mathrm{PCV }(\mathrm{\%})=(\sqrt{{\updelta }_{\mathrm{p}}^{2}}/\mathrm{m})\mathrm\times100\) where, PCV = phenotypic coefficient of variation; \({\updelta }_{\mathrm{p}}^{2}\) = phenotypic variance and \(\mathrm{m}\) = population mean for the trait considered
Genotypic coefficient of variation (\(\mathrm{GCV }(\mathrm{\%})\)) combined over locations = (√δ^{2} g/m) × 100, \(\mathrm{GCV}=(\sqrt{{\updelta }_{\mathrm{g}}^{2}}/\mathrm{m})\mathrm\times100\) where, GCV = genotypic coefficient of variation; \({\updelta }_{\mathrm{g}}^{2}\) = genotypic variance; \(m\) = population mean for the trait considered
G x E interaction coefficient of variation (\(\mathrm{GECV}\)) = \((\sqrt{{\updelta }_{\mathrm{gl}}^{2}}/\mathrm{m})\mathrm\times100\) where \({\updelta }_{\mathrm{gl}}^{2}\), = genotypic x environment variance; \(\mathrm{m}\) = population mean for the trait considered
Heritability in broad sense (H^{2} or h^{2}) was estimated according to [54] as: \({\mathrm{H}}^{2}=({\updelta }_{\mathrm{g}}^{2}/{\updelta }_{\mathrm{p}}^{2})\mathrm\times100\) where, \({\updelta }_{\mathrm{p}}^{2}={\updelta }_{\mathrm{g}}^{2}+({\updelta }_{\mathrm{gl}}^{2}/\mathrm{l})+({\updelta }_{\mathrm{e}}^{2}/\mathrm{rl})\)
Expected genetic advance under selection assuming the selection intensity at 5% was also computed following [54] as: \(\mathrm{GA}=(\mathrm{K})({\updelta }_{\mathrm{p}})({\mathrm{H}}^{2})\) where \(\mathrm{GA}\) = expected genetic advance; \(\mathrm{K}\) = selection differential that varies depending up on the selection intensity and stands at 2.056 for selecting 5% of the genotypes. \({\updelta }_{\mathrm{p}}\) = phenotypic standard deviation and, \({\mathrm{H}}^{2}\) = heritability in broad sense
Genetic advance as percent of mean was obtained by the formula of [55] as; GA (% of mean) = (GA/m) × 100, where, GA = genetic advance; \(m\) = population mean for the trait considered
Phenotypic and genotypic correlation coefficients between two traits were determined by using PROC CANDISC procedure of SAS software following the variance and covariance components [53, 55].
\({\mathrm{r}}_{\mathrm{p}(\mathrm{xy})}={\mathrm{COV}}_{\mathrm{P}(\mathrm{x},\mathrm{y})}/\sqrt{({\updelta }_{\mathrm{px}}^{2})({\updelta }_{\mathrm{py}}^{2})}\), where \({\mathrm{COV}}_{\mathrm{P}(\mathrm{x},\mathrm{y})}\) = phenotypic covariance between traits X and Y, \({\mathrm{r}}_{\mathrm{p}(\mathrm{xy})}\) = phenotypic correlation coefficient between traits X and Y, \({\updelta }_{\mathrm{px}}^{2}\)= phenotypic variance of trait X; \({\updelta }_{\mathrm{py}}^{2}\)= phenotypic variance of trait Y.
\({\mathrm{r}}_{\mathrm{g}(\mathrm{xy})}={\mathrm{COV}}_{\mathrm{g}(\mathrm{x},\mathrm{y})}/\sqrt{({\updelta }_{\mathrm{gx}}^{2})({\updelta }_{\mathrm{gy}}^{2})}\), where \({\mathrm{COV}}_{\mathrm{g}(\mathrm{x},\mathrm{y})}\) = genotypic covariance between traits X and Y, \({\mathrm{r}}_{\mathrm{g}(\mathrm{xy})}\) = genotypic correlation coefficient between traits X and Y, \({\updelta }_{\mathrm{gx}}^{2}\)= genotypic variance of trait X; \({\delta }_{gy}^{2}\)= genotypic variance of trait Y.
Phenotypic and genotypic correlation coefficients were tested for significance using the formula proposed by [53, 56], using the ttable with (g2) degrees of freedom at 5% and 1% significance levels; g is the number of genotypes (treatments) used in the study.
\({\mathrm{t}}_{\mathrm{p}}={\mathrm{r}}_{\mathrm{p}(\mathrm{xy})}/{\mathrm{SE}}_{\mathrm{p}(\mathrm{xy})}\) and \({\mathrm{t}}_{\mathrm{g}}={\mathrm{r}}_{\mathrm{g}(\mathrm{xy})}/{\mathrm{SE}}_{\mathrm{g}(\mathrm{xy})}\), respectively, where \({\mathrm{SE}}_{\mathrm{p}(\mathrm{xy})}\) and \({\mathrm{SE}}_{\mathrm{g}(\mathrm{xy})}\), represent standard error for phenotypic and genotypic correlation, and were computed as:
\({\mathrm{SE}}_{\mathrm{p}(\mathrm{xy})}=\sqrt{{(1{\mathrm{r}}_{\mathrm{p}\left(\mathrm{xy}\right)})}^{2}/(2{\mathrm{H}}_{\mathrm{x}}{\mathrm{H}}_{\mathrm{y}})}\), and \({\mathrm{SE}}_{\mathrm{g}(\mathrm{xy})}=\sqrt{{(1{\mathrm{r}}_{\mathrm{g}\left(\mathrm{xy}\right)})}^{2}/(2{\mathrm{H}}_{\mathrm{x}}{\mathrm{H}}_{\mathrm{y}})}\), where, \({\mathrm{H}}_{\mathrm{x}}\) and \({\mathrm{H}}_{\mathrm{y}}\) are heritability estimate for traits \(\mathrm{x}\) and \(\mathrm{y}\).
Multivariate analysis that includes clustering and PCA was also performed using SAS and Minitab software. Using SAS software version 9.0, the pseudoF and pseudot^{2} statistics were used to calculate the number of clusters. The generalized Mahalanobis D^{2} statistic was used to determine the genetic distance between clusters:
\({\mathrm{D}}_{\mathrm{p}}^{2}={({\overline{\mathrm{X}} }_{\mathrm{i}}{\overline{\mathrm{X}} }_{\mathrm{j}})}^{\mathrm{T}}{\mathrm{S}}^{1}({\overline{\mathrm{X}} }_{\mathrm{i}}{\overline{\mathrm{X}} }_{\mathrm{j}})\) where, \({\mathrm{D}}_{\mathrm{p}}^{2}\) = total generalized distance based on p characters, \({\overline{\mathrm{X}} }_{\mathrm{i}}\) and \({\overline{\mathrm{X}} }_{\mathrm{j}}\) are the p (sample) mean vectors of accessions i and j, respectively, the superscript T denotes matrix transpose, and S denotes the (biascorrected) sample covariance matrix of the n observations in the observed sample which is given as;
\(\mathrm{S}=(\left({\mathrm{n}}_{1}1\right){\mathrm{S}}_{1}+\left({\mathrm{n}}_{2}1\right){\mathrm{S}}_{2})/\mathrm{N}\), [57] where, \({\mathrm{n}}_{1}\) and \({\mathrm{n}}_{2}\) represent random samples of sizes drawn from groups G1 and G2, \({\mathrm{S}}_{1}\) and \({\mathrm{S}}_{2}\) designate (biascorrected) sample covariance matrices, and \(\mathrm{N}={\mathrm{n}}_{1}+{\mathrm{n}}_{2}2\).
The D^{2} value for the cluster pairs was considered as a calculated chisquared value (χ^{2}) and tested for significance at the required probability level against the tabulated values of χ^{2} for p degrees of freedom (d.f = n1), where p is the number of characteristics considered [53].
Principal component analysis (PCA) was performed for the combined and standardized accession mean using MINITAB® Release 19 statistical software [50].
Availability of data and materials
The list of accessions used during the current study is included as “Additional file 1”. Further requests for the raw data and other essential materials can be accessed from the corresponding author; Fekadu Gadissa, email fikega2000@gmail.com, Phone + 251 911 909582.
Abbreviations
 ANOVA:

Analysis of variance
 EBI:

Ethiopian Biodiversity Institute
 GCV:

Genotypic Coefficient of variance
 NMSA:

National Meteorology Service Authority
 PCA:

Principal Components Analysis
 PCV:

Phenotypic Coefficient of variance
References
Vavilov NI. The origin, variation, immunity and breeding of cultivated plants. Chro Bot. 1951;13:1–366.
Harlan JR, Zohary D. Distribution of wild wheats and barley. Science. 1966;153:1074–80.
Bunder J, Loeber A, Broers JEW, Havertkort B. An integrated approach to biotechnology development. In: Bunders J, Haverkort B, Hiemstra W, editors. Biotechnology; building on farmers’ knowledge. London and Basingstoke: Macmillan; 1996. p. 201–27.
Hawkes JG. The diversity of crop plants. London, England: Harvard University Press; 1983.
Chahal GS, Gosal SS. Principles and procedures of plant breeding: biotechnological and conventional approaches. New Delhi: Narosa Publishing House; 2002.
Bekele E. A differential rate of regional distribution of barley flavonoid patterns in Ethiopia, and a view on the centre of origin of barley. Hereditas. 1983;98:269–80.
Negassa M. Patterns of phenotypic diversity in an Ethiopian barley collection, and the ArsiBale Highland as a centre of origin of barley. Hereditas. 1985;102:139–50.
Assamere A, Mamuye A, Amare K, Fiseha M. Diversity, distribution, agronomic and post harvest management of local barley (Hordeum vulgare L.) variety in South Wollo, Ethiopia. PLoS One. 2021;16:e0250899.
Asfaw Z. The barleys of Ethiopia. In: Brush SB, editor. Genes in the field: on farm conservation of crop diversity. Boca Raton: Lewis Publisher; 2000. p. 77–108.
Piffanelli P, Ramsay L, Waugh R, Benabdelmouna A, D’Hont A, Hollricher K, … Panstruga R. A barley cultivationassociated polymorphism conveys resistance to powdery mildew. Nature. 2004;430(7002):887–891.
Mulatu B, Grando S (eds). Barley Research and Development in Ethiopia. Proceedings of the 2^{nd} National Barley Research and Development Review Workshop. 28–30 November 2006, HARC, Holetta, Ethiopia. ICARDA, PO Box 5466, Aleppo, Syria. 2011. pp xiv + 391.
Asfaw Z. Home gardens in Ethiopia: some observations and generalizations. Home gardens and in situ conservation of plant genetic resources in farming systems. 2022. p. 118125.
Gadissa F, Abebe M, Bekele T. Agromorphological traitsbased genetic diversity assessment in Ethiopian barley (Hordeum vulgare L.) landrace collections from Bale highlands, Southeast Ethiopia. Agric Food Secur. 2021;67(10):1–51.
Angassa D, Mohammed J. Agromorphological variability study of Ethiopian barley (Hordeum vulgare L.) accessions for their important agronomical traits at Hadiya zone, Southern Ethiopia. J Plant Sci. 2022;10(1):19–25.
Abtew WG. Characterization of genetic variation among Ethiopian barley (Hoerdeum vulgare L.) genotypes. Stuttgart: PhD dissertation. Institute of Plant Breeding, Seed Science and Population Genetics, Hohenheim University; 2020.
Singh RP, HuertaEspino J, William HM. Genetics and breeding for durable resistance to leaf and stripe rusts in wheat. Turk J Agric For. 2005;29:121–7.
Bartual R, Carbonell EA, Green DE. Multivariate analysis of a collection of soybean cultivars for southeastern Spain. Euphytica. 1985;34:113–23.
Dale MFB, FordLloyd BV, Arnold MH. Variation in some agronomically important characters in a germplasm collection of beet (Beta vulgaris L.). Euphytica. 1985;34:449–55.
Jaradat AA. Phenotypic divergence for morphological and yield related traits among landrace genotypes of durum wheat from Jordan. Euphytica. 1991;52:155–64.
Rezai A, Frey KJ. Multivariate analysis of variation among wild oat accessionsseed traits. Euphytica. 1990;49:111–9.
Demissie A, Bjørnstad A. Geographical, altitude and agroecological differentiation of isozyme and hordein genotypes of landrace barleys from Ethiopia: implications to germplasm conservation. Genet Resour Crop Evol. 1997;44:43–55.
Woodstock LW. Physiological and biochemical tests for seed vigor. Seed Sei and Technol. 1973;1:127–215.
Moore RP. Life alone is not enough how alive are seeds? Seedmen’s Digest. 1955;6(9):12–21.
Walter M, Mertens C, Rhein M. Mixing estimates from a largescale hydrographic survey in the North Atlantic. Geophys Res Lett L13605. 2005;32:13.
James E, Bass LN, Clark DC. Varietal differences in longevity of vegetable seeds and their response to various storage conditions. Am Soc Hort Sci. 1967;92:521–8.
Tenaw S, Tanto T. Genetic diversity of qualitative traits of barley (Hordeum Vulgare L.) landrace populations collected from Gamo Highlands of Ethiopia. Int J Biodiver Conserv. 2014;6(9):663–73.
Tsehaye T, Kebebew F. Morphological diversity and geographic distribution of adaptive traits in finger Millet (EleusineCoracana (L.) Gaertn. [Poaceae]) populations from Ethiopia. Ethiop J Biol Sci. 2002;1:37–62.
Alemayehu F, Parlevliet JE. Variation between and within Ethiopian barley landraces. Euphytica. 1997;94:183–9.
Assefa A, Labuschagne MT. Phenotypic variation in barley (Hordeum vulgare L.) landraces from north Shewa in Ethiopia. Biodiv Conserv. 2004;13:1441–51.
Lakew B, Semeane Y, Alemayehu F, et al. Exploiting the diversity of barley landraces in Ethiopia. Genet Resour Crop Evol. 1997;44:109–16.
Semeane Y, Lakew B, Alemayehu F, van Leur JAG, Grando S, Ceccarelli S. Variation in Ethiopian barley landrace populations for resistance to barley leaf scald and net blotch. Plant Breed. 1998;117:419–23.
Degu HD, Tehelku TF, Kalousova M, Sato K. Genetic diversity and population structure of barley landraces from Southern Ethiopia’s Gumer district: Utilization for breeding and conservation. PLoS One. 2023;18(1):e0279737.
Kassahun D, Asefa D, G/Kidan D. Synergistic effect of a combined application of organic and inorganic fertilizers on yield and yield components of Tef (Eragrostis tef (Zucc) Trotter) under terminal drought at Adiha, Northern Ethiopia. J Drylands. 2010;3(1):158–64.
Jalata Z, Ayana A, Zeleke H. Variability, heritability and genetic advance for some yield and yield related traits in Ethiopian barley (Hordeum vulgare L.) landraces and crosses. Int J Plant Breed Genet. 2011;5:44–52.
Deshmukh SN, Basu MS, Reddy PS. Genetic variability, character association and path coefficient analysis of quantitative traits in Virginia bunch varieties of groundnut. Indian J Agric Sci. 1986;56:816–21.
Andonov KL, Sariev BS, Zhundibaev LP. Structure of phenotypic variability in traits of spring barley. Acta Agril Shanghai. 1979;22:187–8.
Singh BD. Plant breeding: principles and methods. New Delhi: Kalyani Publishers; 2001. p. 896.
Johnson HW, Robinson HW, Comstock R. Estimates of genetic and environmental variability in soybeans. Agron J. 1955;47:314–8.
Shafique MS, Ahsan M, Mehmood Z, Abdullah M, Shakoor A, Ahmad MI. Genetic variability and interrelationship of various agronomic traits using correlation and path analysis in Chickpea (Cicer arietinum L.). Acad J Agric Res. 2016;4(2):082–5.
Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. Longman Group Limited: Harlow; 1996. p. 438.
Akinwale MG, Gregorio G, Nwilene F, Akinyele BO, Ogunbayo SA, Odiyi AC. Heritability and correlation coefficient analysis for yield and its components in rice (Oryza sativa L.). Afr J Plant Sci. 2011;5:207–12.
De Lacy IH, Cooper M. Pattern analysis for the analysis of regional variety trials. In: Kang MS, editor. Genotypebyenvironment interaction and plant breeding. Baton Rouge: Louisiana State University; 1990. p. 301–34.
Balkaya A, Ozbakir M, Kurtar ES. The phenotypic diversity and fruit characterization of winter squash (Cucurbita maxima) populations from the Black Sea region of Turkey. Afr J Biotech. 2010;9:152–62.
Jagadev PN, Samal KM, Lenka L. Genetic divergence in rape mustard. Indian J Genet. 1991;51:465–6.
Falconer IS. Introduction to quantitative genetics. 2nd ed. New York: Longman Inc.; 1981.
ISTA [International Seed Testing Association]. International rules for seed testing. Seed Sci Technol. 1996;21(Suppl):1–288.
Maguire ID. Speed of germinationaid in selection and evaluation for seedling emergence and vigour. Crop Sci. 1962;2:176–7.
Ranal MA, Santana DG. How and why to measure the germination process? Rev Brasil Botân. 2006;29:1–11.
IPGRI. Description of Barley (Hordeum vulgare L). Rome: International Plant Genetic Resource Institute; 1994.
Minitab LLC. Minitab statistical package for windows. State college: Minitab Inc.; 2019.
Gomez KA, Gomez A. Statistical procedures for agricultural research. 2nd ed. New York: Wiley; 1984. p. 680.
Peterson DM, Wesenberg DM, Burrup DE, Erickson CA. Relationships among agronomic traits and grain composition in oat genotypes grown in different environments. Crop Sci. 2005;45(4):1249–55.
Singh RK, Chaudhary BD. Biometrical methods in quantitative genetic analysis. New Delhi: Kalyani Publishers; 1985. p. 13–55.
Allard RW. Principle of plant breeding. New York: Wiley; 1960.
Robertson GR. The sampling variances of the genetic correlation coefficients. Biometrics. 1959;15:469–85.
Sharma JR. Statistical and biometrical techniques in plant breeding. Calcutta: New Age International Limited Publishers; 1998.
Mahalanobis PC. On the generalized distance in statistics. In: Proceedings of National Institute of Science of the India, New Delhi. 1936;2(1):4955.
Acknowledgements
The authors would like to thank Mada Walabu University for material and technical supports for this research. We would also like to thank the Ethiopian Biodiversity Institute (EBI) for offering us the germplasm.
Funding
The study was partly supported by Madda Walabu University research directorate through annual research budget which has no budget or grant number and code.
Author information
Authors and Affiliations
Contributions
FG and TBG did all the field experiment, data collection, data clearing; FG data analysis, data interpretation and write up; TBG edited the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All the procedures in accordance with the relevant national, international and institutional guidelines have been followed. The barley accessions used were obtained from the Ethiopian Biodiversity Institute (EBI).
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
List of the exsitu conserved barley landraces, along with their passport data, considered in the present study.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gadissa, F., Gudeta, T.B. Phenotypic characterization and seed viability test in exsitu conserved Ethiopian cultivated barley (Hordeum vulgare L.) landraces. BMC Plant Biol 23, 613 (2023). https://0doiorg.brum.beds.ac.uk/10.1186/s12870023046287
Received:
Accepted:
Published:
DOI: https://0doiorg.brum.beds.ac.uk/10.1186/s12870023046287