Skip to main content

Adaptation and diversity along an altitudinal gradient in Ethiopian barley (Hordeum vulgare L.) landraces revealed by molecular analysis



Among the cereal crops, barley is the species with the greatest adaptability to a wide range of environments. To determine the level and structure of genetic diversity in barley (Hordeum vulgare L.) landraces from the central highlands of Ethiopia, we have examined the molecular variation at seven nuclear microsatellite loci.


A total of 106 landrace populations were sampled in the two growing seasons (Meher and Belg; the long and short rainy seasons, respectively), across three districts (Ankober, Mojanawadera and Tarmaber), and within each district along an altitudinal gradient (from 1,798 to 3,324 m a.s.l). Overall, although significant, the divergence (e.g. FST) is very low between seasons and geographical districts, while it is high between different classes of altitude. Selection for adaptation to different altitudes appears to be the main factor that has determined the observed clinal variation, along with population-size effects.


Our data show that barley landraces from Ethiopia are constituted by highly variable local populations (farmer's fields) that have large within-population diversity. These landraces are also shown to be locally adapted, with the major driving force that has shaped their population structure being consistent with selection for adaptation along an altitudinal gradient. Overall, our study highlights the potential of such landraces as a source of useful alleles. Furthermore, these landraces also represent an ideal system to study the processes of adaptation and for the identification of genes and genomic regions that have adaptive roles in crop species.


Among the cereal crops, barley is the species that presents the highest adaptability to a wide range of environments. It is cultivated from arctic latitudes to tropical areas, and it is grown at the highest altitudes. In Tibet, Nepal, Ethiopia and the Andes, farmers cultivate barley on mountain slopes at altitudes higher than any other cereals [1, 2].

Ethiopia is probably the region of barley cultivation that presents the highest variability for climatic and edaphic conditions. It is cultivated from 1,400 m above sea level (a.s.l.) to over 4,000 m a.s.l., and it has adapted to specific sets of agro-ecological and microclimatic regimes throughout the country [3]. Landraces represent over 90% of the barley cultivated in Ethiopia. In contrast to the genetic uniformity of modern cultivars, landraces show variation both between and within populations. This within-population diversity of these barley landraces might allow them to cope with environmental stresses, which is very important for achieving yield stability [4]. For this reason, landraces represent a very interesting model to study the processes of adaptation and for identification of genes and genomic regions that have adaptive roles in a crop species, through association mapping [58] and through scanning for signatures of selection at the molecular level [913].

Knowledge of the population structure of Ethiopian barley landraces, together with a deeper understanding of the nature and extent of their variation, is an important prerequisite for efficient conservation and use of the existing plant materials. Several studies have reported a high level of genetic diversity in barley populations from Ethiopia, such as those based on morphological traits [1416], and on biochemical data [17, 18]. Moreover, Ethiopian barley is a precious source of genes that control important agronomic traits, such as resistance to disease (e.g. powdery mildew, barley yellow dwarf virus, net blotch, scald and loose smut) and to insect attack [1925], high lysine, and protein quality and content [26], and malting and brewing quality [27]. Furthermore, the great variability of environmental conditions in Ethiopia that promote adaptive divergence, and the cultivation of barley in two growing seasons per year [16], have probably driven the structure of variation of these landraces.

Numerous studies on Ethiopian barley have detected correlations between altitude and the frequency of morphological [1416, 28], agronomic [24, 29] and biochemical [30] traits, including pathogen resistance [20, 24, 25, 29]. Altitudinal gradients provide substantial changes in numerous environmental variables, which also include atmospheric pressure, temperature, clear-sky turbidity, UV-B radiation, and humidity. Furthermore, the combined environmental changes across altitudes can influence various biological processes in plants, and can thus result in adaptive changes and constraints on the genetic diversity of plant populations. Association between environmental variables and allele frequencies can be maintained by a balance between selection and gene flow [3133]. However, in contrast to latitudinal clines, altitudinal gradients involve dramatic ecological transitions over relatively short linear distances, and so they require especially strong selection to counterbalance the homogenising effects of gene flow [34]. Nevertheless, the interactions between drift and spatially restricted gene flow or historical patterns of colonisation (admixture between previously isolated populations) might also explain clinal variation.

The present study is based on a large collection of Ethiopian barley landraces (3,170 individual plants, for a total of 106 landrace populations) that was established in 2005 in the North Shewa zone (Amhara region, Ethiopia), from farmers who have used their own seed for generations [16]. The collection was obtained through visiting the same three districts (Ankober, Mojanawadera and Tarmaber) in both the Meher and Belg growing seasons [16] (see [16], for further details). Here, we have analysed the molecular diversity and genetic structure of this collection, using simple sequence repeats (SSRs), while taking into account not only the environmental adaptation (geographical and altitude factors), but also, for the first time, the two growing seasons per year.


Level of polymorphism

Seven mapped nuclear SSRs were used to examine the levels and patterns of genetic variation of the barley landraces collected in North Shewa, in the central highlands of Ethiopia [16]. Table 1 provides a summary of the statistics computed considering each of the two seasons (Belg and Meher), the three districts (Ankober, Mojanawadera, and Tarmaber) and the three altitude classes (<2,300, 2,300-2,800, >2,800 m a.s.l.). The same statistics were also computed for each locus, and are given in Additional file 1. A total of 66 alleles were detected, with the number of alleles per locus (n a ) ranging from four (HVM20 and HVM67) to 23 (Bmac0156), with an average of 9.4 alleles per locus.

Table 1 Summary statistics for the seasons, districts and altitude classes, and for the whole sample.

The variation between seasons was not significant across any of the statistics considered (Wilcoxon signed-paired rank test), while the only significant difference among the districts (Wilcoxon signed-paired rank test, after Bonferroni correction, P = 0.03) was for the allelic richness between the Ankober (R S = 8.22) and Tarmaber districts (R S = 7.11; Table 1). For the altitude classes, more marked differences were seen (Table 1). The number of alleles (n a ) and the allelic richness (R S ) of the low altitude class (<2,300 m a.s.l.: n a = 39 and R S = 5.57) were lower than those of the intermediate altitude class (2,300-2,800 m a.s.l.: n a = 50 and R S = 7.10; difference marginally non-significant, P = 0.06, Wilcoxon signed-paired rank test, after Bonferroni correction) and of the high altitude class (>2800 m a.s.l.: n a = 57, P = 0.06, and R S = 7.40, P = 0.09). The genetic diversity (He) of the intermediate altitude class (He = 0.65) was significantly higher (Wilcoxon signed-paired rank test, after Bonferroni correction, P = 0.03) than that of the high altitude class (He = 0.57).

Considering the pattern of private alleles, the two seasons did not share 14 (21%) out of the 66 SSR alleles, with each of these showing a frequency lower than 5%. The three districts showed a total of 15 private alleles (23%), three of which had a frequency higher than 5%. Finally, compared to seasons and districts, the altitude classes had a higher number of private alleles (17; 26%), four of which (one for <2,300 m a.s.l., and three for >2,800 m a.s.l.) had a frequency higher than 5%. Moreover, the average frequency of private alleles that showed a frequency higher than 5% was 0.00 for seasons, 0.063 for districts, and 0.085 for altitude classes.

The overall level of observed heterozygosity was very low (0.003). Only three individuals, all located at high altitude, showed one or two heterozygous loci. Overall, the average genetic diversity (He, [35]) in the Ethiopian barley landraces was 83.5% of the diversity found in the Syrian and Jordanian barley landraces (Table 2). This difference was significant (Wilcoxon signed-paired rank test, P = 0.05). A significant gap (-31.2%) was also seen for the number of alleles, n a (Wilcoxon signed-paired rank test, P = 0.01).

Table 2 Diversity of landraces from Ethiopia compared with those from Syria, Jordan and Nepal, and with modern varieties.

Although it is important to consider the differences in the sample sizes, the Ethiopian barley landraces showed a higher number of alleles (n a = 62, not considering the HVM20 locus) than seen for both the Nepalese barley landraces (n a = 44) [36] and the modern varieties (n a = 41) [37] (Table 2).

The divergence (FST and RST) estimates between the seasons, districts, and altitude classes are given in Table 3. The genetic differentiation between the two seasons was low for both FST (0.02) and RST (0.01), even if in contrast to RST, the FST was significantly different from zero (P < 0.001). A slightly higher level of differentiation was seen for districts (FST = 0.02, P < 0.01; RST = 0.03, P < 0.05). Among the altitude classes, the levels of differentiation were higher compared to seasons and districts, both for FST (0.10, P < 0.001) and RST (0.09, P < 0.001). The divergence (FST and RST) estimates between seasons, districts and altitude classes for each of the SSR loci analysed are given in Additional file 2. The individual SSR loci divergence (FST,) between seasons was significant for four loci (Bmag0013, HVM67, Bmac0040, and Bmac0134), even if the values were low, and varied from 0.02 to 0.05 (Additional file 2). The Bmag0013 locus was the only one for which the divergence was higher between seasons (FST = 0.05) than among districts or altitude classes (FST = 0.01 and 0.03 respectively; Additional file 2). A similar trend was seen for districts, with three low significance values ranking from 0.02 (Bmac0156) to 0.04 (Bmac0134). For the altitude classes, the situation was different; indeed, all of the SSR loci were significantly differentiated among the altitude classes, with values varying from 0.03 (Bmag0013) to 0.28 (HVM67).

Table 3 Divergence (FST and RST) estimates between seasons, districts and altitude classes.

Table 4 reports the divergence estimates (FST) among the three districts within the same altitude classes (<2,300 m a.s.l.; 2,300-2,800 m a.s.l., and >2,800 m a.s.l.). The differentiation among the low altitude classes of the three districts was higher (FST = 0.10, P < 0.01) than that among the three districts within the intermediate (FST = 0.04, P > 0.05) and high (FST = 0.01, P > 0.05) altitude classes.

Table 4 Divergence (FST) estimates between the three different districts within the same altitude classes.

Population structure and altitude cline of genetic variation

Recently, various studies (mainly in the field of human genetics [38, 39]) have questioned the reliability of the STRUCTURE results in cases of complex genetic structures. In particular, this debate has considered whether Bayesian clustering algorithms are appropriate tools for studying genetic structures in populations with continuous variation of allele frequencies. To address this problem, the population structure analysis of our collection was carried out using two different methods: the first was the Pritchard et al. [40] non-spatial Bayesian clustering method implemented in the STRUCTURE software, version 2.1 [41], and the second was based on a Bayesian clustering algorithm that incorporates spatial information when identifying clusters of individuals, and was implemented in the TESS software, version 2.0 [42].

The plot of the average ln likelihood values over 100 runs for K values ranging from 1 to 10 shows that the ln likelihood estimates increase progressively as K increases. Thus, we used the Evanno et al. [43] method to provide a better estimate of the 'true' number of clusters (K) that characterised our sample. The highest ΔK value was found at K = 2, although high values of ΔK were found also for K = 3, 4 and 10. To minimise the effects of outlier runs, we also computed the ΔK value considering five random re-samplings of 80 of the 100 runs for each K value. In all cases, the highest ΔK corresponded to K = 2. Thus, hereafter, we will refer to cluster S1 and cluster S2 for the K = 2 populations identified by STRUCTURE.

The percentages of membership (q), or inferred group ancestries, of the individuals in each of the two clusters were computed. One hundred and ninety-seven individuals (93%) had a q higher than 0.70, in 192 individuals (91%), q was higher than 0.80, and in 183 (86%), it was higher than 0.90. Considering the threshold for membership of q ≥ 0.70, cluster S1 included 110 individuals, and cluster S2, 87 individuals.

Figure 1 shows the collection site coordinates and the average q values for clusters S1 and S2 of each barley landrace.

Figure 1

Map of the collection sites. Collection site coordinates and average q values for clusters S1 (green) and S2 (red) of each of the barley landraces. Landraces that showed an average q > 0.70 were considered completely assigned to one of the two clusters identified by STRUCTURE.

When the individuals were sorted according to altitude, the STRUCTURE analysis revealed a well-defined cline of variation; indeed, at low altitudes, cluster S1 was mainly seen, and this tended to be substituted by cluster S2 at high altitudes. The relationships between altitude and q for cluster S2 were assessed using a simple linear regression model. This regression was significant (R2= 0.42, F = 152.5 and P < 0.0001; Figure 2), and it was consistent within the Meher (R2= 0.45, F = 85.0 and P < 0.0001) and the Belg (R2= 0.36, F = 59.4 and P < 0.0001) seasons (Figure 3), and within the three districts (P < 0.0001, R2 = 0.54, 0.21, and 0.40, F = 83.6, 17.8, and 50.0; for Ankober, Mojanawadera and Tarmaber, respectively) (Figure 4).

Figure 2

Linear regression analysis for altitude and membership coefficient ( q ) for cluster S2 identified by STRUCTURE.

Figure 3

Population structure across seasons and altitudes using STRUCTURE. STRUCTURE assignment at K = 2 for the 108 individuals collected during the Belg growing season and the 104 individuals collected during the Meher season, ordered according to crescent altitude levels. Green, cluster S1; red, cluster S2.

Figure 4

Population structure across districts and altitudes using STRUCTURE. STRUCTURE assignment at K = 2 for the individuals collected in the three districts (Ankober, Mojanawadera and Tarmaber), ordered according to crescent altitude levels. Green, cluster S1; red, cluster S2.

The second approach used for the population structure analysis was the algorithm implemented in the TESS software that allows for prior spatial information when identifying clusters of individuals. This TESS software allows the testing of different values of the spatial dependent parameter (Ψ); this parameter weights the relative importance given to spatial connectivities. We tested three different Ψ: 0.0 (no spatial correlation between individuals), and 0.6 and 1.0 (moderate and strong spatial correlation between individuals, respectively). For each Ψ value considered, we performed 100 runs for K from 2 to 10. Then the DIC values of the 10 runs with the smallest DIC were averaged. In plotting the average DIC values versus K (from 2 to 10), we obtained an inflection point at K = 6. This was the same for all of the different interaction parameter (Ψ) values considered. Thus the best model was obtained for K = 6 and the lowest DIC value corresponded to Ψ = 0.6, indicating a moderate spatial correlation between individuals. The average over 10 runs with the smallest DIC value q was computed. We defined the six TESS clusters as T1, T2, T3, T4, T5 and T6. One hundred and twenty-five individuals (59%) were assigned to one of the six clusters with a q higher than 0.60, with 104 (49%) higher than 0.70, and 70 (33%) higher than 0.80. Figure 5 shows the assignment results for all of the individuals when ordered according to altitude. Averaging the q for the six TESS clusters for the three altitude classes, we observed that the previous cline of variation along altitude identified by the STRUCTURE assignment was also maintained for TESS; indeed, at low altitudes, the T1 and T2 clusters were mainly present, and these tended to be substituted by the T5 and T6 clusters at high altitudes (Figure 6). The T3 cluster was mostly present at intermediate and high altitudes, while T4 was present in all of the altitude classes (Figure 6). However, the genotype assignment for T4 is more scattered at intermediate and high altitudes than at low altitudes (Figure 5). Indeed, at low altitudes, this genotype assignment for T4 is characteristic of two populations collected during the Meher season. It is important to note that the percentage of individuals that had a q to one of the six TESS clusters higher than 0.70 is highest (76%) at low altitudes (<2,300 m a.s.l.), and decreases at intermediate (35%) and high (45%) altitudes.

Figure 5

Population structure along the altitude gradient using TESS. TESS assignment at K = 6 for the 212 individuals, ordered according to crescent altitude levels. Blue-purple, T1; red, T2; pink, T3; yellow, T4; grey, T5; green, T6.

Figure 6

Average membership coefficient ( q ) for the six TESS clusters for the three altitude classes. See legend to Figure 5 for colour key.

The correlation analysis using Spearman's rho (ρ) for non-parametric tests showed that cluster S1 identified by STRUCTURE was significantly correlated with TESS clusters T1 and T2, and cluster S2 with T5 and T6 (Additional file 3). The step-wise multiple regression analysis showed that mainly TESS clusters T1, T2, T3 and T4 provide a better prediction for cluster S1 (R2 = 0.84), and T4, T5 and T6 are the major contributors for cluster S2 (R2 = 0.84) (Additional file 4); TESS cluster T4 contributes to both of the S1 and S2 clusters identified by STRUCTURE. These results are consistent with the relationships between the altitude and the q for the TESS groups, which was tested by performing a simple regression analysis; indeed this relationship was highly significant (P < 0.0001) for the TESS clusters T1, T2, T5 and T6 (R2 = 0.25, 0.14, 0.10, 0.22 respectively), significant (P = 0.05) for T3, but not significant for T4 (P = 0.08).

The ANOVA analysis revealed significant correlations between all of the loci considered and the altitude (Additional file 5).

Considering the morphological data, 82% of the character states considered and studied by Tanto Hadado et al. [16] were significantly correlated with altitude (P < 0.05), 79% of which had a probability lower than 0.01; moreover, at least one character state for each of the morphological traits considered was significantly correlated with altitude at P < 0.01 (Table 5).

Table 5 Relationships between frequencies of the character states for each population and altitude (Spearman's rho).

Association between morphological and molecular characterisations

To investigate the correspondence between the morphological and molecular characterisations, a contingency analysis was performed. The genetic groups defined by TESS analysis were considered, and only the genotypes that were assigned to one of the groups with a q higher than 0.70 were included in the analysis (104 genotypes in all). No inference was carried out for TESS cluster T3 because no genotypes were assigned to this group according to this threshold. The analysis was conducted considering the seven morphological traits: kernel-row number, spike density, lemma awn barbs, glume colour, lemma type, length of rachilla hair, and lemma colour [16].

The results show that all of the traits except lemma awn barbs (P = 0.066), were significantly associated with the TESS genetic clusters. The association was higher for length of rachilla hairs, glume colour and spike density (P < 0.0001), followed by lemma colour (P < 0.001), row number (P = 0.003) and lemma type (P = 0.011). The proportion of variance of each morphological trait explained by genetic groups (R2) varied from 36.6% (length of rachilla hairs) to 13.0% (lemma type), with an average of 22.5%. The types and frequencies of each of these seven qualitative traits in each of the TESS clusters (with the exclusion of T3, as indicated above) are given in Additional file 6.

The relationships between the genetic groups according to morphological characterisation were assessed using a non-parametric correlation analysis (Spearman's rho) (Additional file 7). The clusters T1 and T6 (the former mainly present at low altitudes, and the latter at high altitudes) showed the highest morphological differentiation. T1 included mainly six-rowed (65.2%) and irregular (21.7%) types, and was the only group that included two-rowed types (4.4%), and it showed the highest frequency among all of the groups for the two-rowed deficient types (8.7%). Both the regular and deficient two-rowed types were absent in T6, which included mainly irregular types (74.2%). T1 was characterised by a high frequency of dense spike types (47.8%), while T6 was characterised by lax types (41.9%), and both of these groups showed high frequencies of intermediate spike density types (43.5% and 54.9%, respectively). T1 was monomorphic for the lemma awn barbs (rough) and glume colour (white) traits; in contrast, T6 also showed intermediate lemma awn barb types (12.9%), and brown (35.5%) and black (6.4%) glume colour types. All of the individuals of T6 had short rachilla hairs, while those of T1 had both short (52.2%) and long (47.8%) rachilla hair types. Finally, differences were seen also for the lemma colour trait: T1 was mostly represented by white types (82.6%), while T6 was mostly represented by black types (74.2%).

Spatial structure and landscape analysis

Mantel tests between the genetic and geographical distance matrices were significant (r = 0.12, P < 0.001). However, the correlation between molecular and geographical distance was about 2-fold lower than that between the molecular distances and differences in altitude (r = 0.28, P < 0.001), and also lower than that between geography and altitude (r = 0.14, P < 0.001).

The same pattern was seen for the spatial autocorrelation analysis, which showed that within a distance class of 5 km, positive and significant Moran's I values were detected. This indicates that the proximal individuals (included those of the same population/field) were more genetically similar (i.e. related) than expected from a random distribution (Figure 7a). In this class, the individuals at similar altitudes within each district were compared. The degrees of correlation, however, rapidly decreased, and there were already negative Moran's I values for the 10-km distance class. In this class, individuals from relatively different altitudes were mostly compared, mainly within each district. A similar situation was seen in the 20-km and 30-km distance classes. The degrees of correlation, however, tended to increase, and interestingly, positive and significant Moran's I values were found for the 40-km and 50-km distance classes: when individuals of different districts but at similar altitudes (either low or high) were compared, they tended to be genetically more similar than for a random distribution (Figure 7a). This is particularly evident considering the relationships between the genetic and altitude distances; indeed, the individuals that were collected at similar altitudes showed higher similarities than those collected at different altitudes (Figure 7b).

Figure 7

Autocorrelation analysis. Results of the autocorrelation analysis using all 212 genotypes and all loci. The 99% probability envelopes are indicated: light blue, upper limit; green, lower limit. Dark blue, observed data: for genetic distance vs geographical distance). Inset: the first 20 km from panel a) (40 classes, 500 m each).

Finally, it is important to note that based on the distribution of the collection sites, when individuals from similar altitudes were compared, they also tended to be geographically more distant. In contrast, the individuals from different altitudes were geographically closer (Figure 7c).

Spatial autocorrelation between geographical and genetic distances was also performed separately for the three altitude classes (Additional file 8), which confirmed the previous trends seen, showing at low altitudes (<2,300 m a.s.l.) a clear geographical effect (isolation by distance relationship among plants), while at high altitudes (>2,800 m a.s.l.), this effect disappeared, as all of the similarity values were not significantly different from random values, and thus no geographical structure was evident. A non-significant, intermediate (but as for the >2,800 m a.s.l. class) trend was seen at moderate altitudes (2,300-2,800 m a.s.l).

Sliding-window analysis

Figure 8 shows the results of the 'sliding-windows' analysis. In Figure 8, each point represents: for the y axis, the He, R S and LD estimates of a group of 40 individuals from 20 landrace populations (window); and for the x axis, the mean altitude over all of the 20 populations. Thus, the diversity and LD estimates were computed along the altitudinal cline. This analysis was carried out separately for each of the two seasons. The Belg season was not represented at low altitudes, while the Meher season was characterised by the whole altitude range considered.

Figure 8

'Sliding-window' analysis for the two seasons: Belg (blue) and Meher (orange). a) Genetic diversity (He); b) allelic richness (RS); c) multilocus linkage disequilibrium (LD) along the altitudinal cline, with P values also shown. The altitude of each window (20 populations, 40 individuals) was determined by averaging the altitude of the populations (farmers' fields); hence, the differences in the average altitudes between the windows are not constant.

To test the significance of the trends seen, we grouped the individuals into different and non-overlapping altitude ranges: three for Belg: B1, B2 and B3 (average altitudes of 2,397, 2,859 and 3,022 m a.s.l., respectively), and four for Meher: M1, M2, M3, M4 (average altitudes of 1,962, 2,416, 2,909 and 3,138 m a.s.l., respectively). The significance was tested between adjacent altitude ranges, using the Wilcoxon signed-paired rank test and the Bonferroni correction. The choice of the altitude ranges was designed to have the average altitudes more uniform between the Belg (B1, B2, and B3) and the corresponding Meher (M2, M3, M4) classes (the M1 class did not have an equivalent in Belg, which was not represented at low altitudes in this analysis). The number of individuals per altitude range varied from 34 to 40 for Belg, and from 16 to 36 for Meher.

Considering the Meher season, the diversity (He) showed a significant increase according to altitude, to about 2,500 m a.s.l.. Indeed, M2 showed a higher level of diversity than M1 (Wilcoxon signed-paired rank test, after Bonferroni correction: P = 0.018); although no other comparisons were significant in the Meher season, a reduction at high altitudes was seen. A significant reduction of He at high altitudes was seen for the Belg season (B2-B3: P = 0.046) (Figure 8a).

The allelic richness, R S , showed a marked increase from low to intermediate altitudes in the Meher season (Wilcoxon signed-paired rank test, after Bonferroni correction: M1 < M2: P = 0.017), and reached a plateau at intermediate and high altitudes (M2 vs M3; M3 vs M4) (Figure 8b). For the Belg season, the trend of allelic richness resulted in a slight increase at intermediate altitudes (B1 < B2: P = 0.046) and a decrease at high altitudes (B2 > B3, P = 0.018).

The multilocus LD among our unlinked loci (Figure 8c) decreased with altitude for the season, and interestingly, only that for the Belg season reached non-significant levels in the interval between 2,685 and 3,010 m a.s.l. (Figure 8c).


Several studies have analysed the diversity of Ethiopian barley landraces. However, along with our previous study on the same barley landrace collection that was based on morphological data [16], for the first time the present study has considered an ad-hoc collection of barley landraces from Ethiopia that was obtained on a single-plant basis, and that also included the two main growing seasons.

Level of polymorphism in Ethiopian barley landraces from North Shewa

In this study, we analysed a collection of barley landraces from an area of about 2,000 km2 in the central highlands of Ethiopia, in North Shewa. All of the seven microsatellites used for the analyses are in common with the study of Russel et al. [44] which included barley landraces from several different agro-ecological areas of Syria and Jordan, and six are also shared with the studies of both Pandey et al. [36], which analysed barley landraces from the highlands of Nepal, and Macaulay et al. [37], which analysed a sample of modern barley varieties. Interestingly, as the different barley samples were analysed with the same microsatellites, this has allowed us to make comparisons between the results of the present study and those obtained in these other studies. We have shown here a relevant level of genetic diversity in the Ethiopian barley landrace collection, which was higher than both the Nepalese barley landraces and the sample of modern varieties. The substantial amount of genetic variation detected in the present study was also in agreement with the high level of morphological diversity seen in the same sample of landrace populations by Tanto Hadado et al. [16].

The Ethiopian landraces from North Shewa showed a lower genetic variation than those of the Syrian and Jordanian landraces. However, it is important to consider that Russel et al. [44] worked at a country level (i.e. on a larger geographical scale), with their study of barley landraces collected from several different agro-ecological areas of Syria and Jordan. For this reason, the level of diversity found in the much smaller area in our survey in the present study is undoubtedly remarkable.

Population structure

Overall, the population structure appears to be very low between seasons and geographical districts, while it is highly associated with altitude differences. Indeed, we show significant, but very low, divergence (FST = 0.02) with no private alleles (p < 0.05) between the two different planting seasons. Moreover, the two seasons show very similar levels of diversity. These data are very consistent to those for the morphological traits on the same populations [16], and they confirm that seed flow occurs between seasons, as was also suggested by 5% of the farmers' fields during the collection, where a sequential cultivation of barley was seen in the two seasons. In other words, barley is often planted in one season using the seeds obtained from the harvest in the previous season. A slightly higher divergence was seen between districts, with Ankober having a higher diversity compared to the other two districts.

In contrast, the divergence estimators show strong effects due to altitude. In the present study, we used two different Bayesian approaches (the STRUCTURE and TESS programmes) to make inferences on the population structure of our sample. STRUCTURE identified two main clusters: rather than a subdivision among seasons or districts, these two clusters highlighted a clear altitudinal cline, and this tendency was maintained both within seasons and districts. Even though the analysis conducted with TESS identified more genetic groups (six), these were in line with the STRUCTURE results. Indeed, five TESS clusters were highly associated to the two clusters identified by STRUCTURE, with only one exception: TESS cluster T4. This indicates that the data obtained with the two different approaches are not completely redundant, with a more detailed and refined picture being provided by TESS. However, independent of a clear or a less clear definition of the "real" numbers of the populations that characterised our sample of genotypes, it is important to stress that both of these tools were very useful and that they were in agreement in the defining of the major aspects of the genetic structure of these barley landrace populations. For instance, as for the STRUCTURE clusters, the TESS groups reflect the altitudinal cline of variation, with five of them significantly correlated with altitude. The same effect of altitude was also seen for all of the SSR loci used in the present study. Moreover, most (84%) of the character-state frequencies obtained from Tanto Hadado et al. [16] (at least one for each trait) were found to be significantly associated with altitude. Thus, in line with other published works cited in the Introduction, the present study has identified a strong association between altitude and genetic composition of barley landraces from Ethiopia, with a clear altitudinal cline shown.

Two possible evolutionary scenarios can be considered to explain the altitudinal cline of variation that we have detected here: a clinal variation might have originated by various patterns of migration, or by different selections along the geographical or environmental gradients [45]. The results in the present study for the spatial autocorrelation between altitude and genetic distances reflects a continuous and constant reduction in the Moran I index, from positive to negative values; such a trend is characteristic of a typical strong selection gradient [45]. However, when there is a large difference between the immigrant and the resident gene frequencies, this trend can also be generated by a specific case of unidirectional migration [45]. The hypothesis of unidirectional migration in our case would be in a vertical direction, from low to high altitudes, although even if this is theoretically possible, under our conditions it is not realistic. This appears evident from the geographical structure of the collection area, with the three districts. Moreover, when we analysed the correlogram related to the comparison between geographical and genetic distances, for the first distance class (between genotypes at a distance lower than 5 km), there were positive Moran I values that became significantly negative just between the genotypes with a distance of 5 km to 10 km, where the differences in altitude are already substantial (on average, 125-250 m). Furthermore, the Moran I values became positive only at high geographical distances, when comparisons were mainly made between genotypes at similar altitudes, but in different districts. These data confirm the hypothesis that selection for adaptation at different altitudes is the main factor that determines the clinal variation seen. This divergence between different districts but within the same altitude class is consistent with the selection hypothesis, and it is suggestive of different selection intensities at low and high altitudes. Indeed, when we compared individuals at low altitude across the three districts, there was a significant and relatively high FST (0.10, P < 0.01), while there were very low and non-significant FST values between individuals at intermediate (0.04) and high (0.01) altitudes from different districts. These data are confirmed by a separate spatial autocorrelation at different altitude classes, which shows a correlation between geographical distances and genetic distances only at low altitudes. This suggests that the homogeneous selection for adaptation to high altitudes is the major factor that explains these data, along with isolation by distance (IBD) at low altitudes. Alternatively, we should consider the hypothesis of different migration patterns at low vs high altitudes: a higher level of seed flow between districts at high altitudes than at low altitudes.

In support of the role of adaptive selection in shaping the allele frequencies of our collection, we found that the two morphological traits analysed that relate to colour (kernel and glume colour) are significantly associated with altitude. With their colour linked to the presence of anthocyanins, the coloured types in particular were positively related to altitude, while the white types showed a negative relationship with altitude. These anthocyanins are a subgroup of the flavonoids, for which an important role is well documented in responses to both biotic and abiotic stresses (see [46], for review). Thus, the higher presence of coloured types at high as opposed to low altitudes might be an adaptive response to the increased UV-B radiation characteristics at high altitudes.

However, even though there is a large amount of coherent evidence that suggests that selection is a crucial factor in shaping the diversity of Ethiopian barley landrace populations, particularly at high elevations, selection alone is not enough to explain the diversity patterns seen. Indeed, the diversity significantly increases from low to moderate altitudes (2,500 m a.s.l.), both in its richness and its expected heterozygosity. This occurred only for the Meher season, because barley is grown in the Belg season mainly at higher altitudes. This pattern can be explained by the increasing number of populations (fields) and by the higher population sizes (larger plots), as shown by Tanto Hadado et al. [16]. At low altitudes, the environmental conditions are more favourable on average, and many crops can be grown easily; thus the farmers tend to grow many different crops (e.g. tef, maize and sorghum), and in general, barley tends to be cultivated in smaller plots, as compared to at higher elevations. Thus, the increased diversity might be related to the greater population size of the barley landraces, which is not only related to the barley plot size increase, but also to the two seasons of barley growth per year at intermediate and high altitudes; thus, the effective number of generations per year is higher than one [16]. The maximum diversity is seen at moderate altitude (2,460 m a.s.l.), where due to the spring rains, the barley is grown in both seasons. When the altitude further increases, for both seasons, the diversity measured as expected heterozygosity decreases significantly with altitude, while the richness remains almost constant, even if at very high altitudes a small but significant reduction is seen. This trend is in agreement with the hypothesis of an increasing selection intensity for adaptation at high altitudes. Selection is expected to increase the frequency of favourable alleles, which would reduce the expected heterozygosity, while the richness is less affected by selection, compared to the expected heterozygosity; this would be favoured by the greater population size. However, at very high altitudes, a reduction in the number of alleles is seen. Thus, we explain the diversity pattern seen by the combined effects of selection and population size. While the plot (field) sizes constantly increase with altitude in a linear fashion, the number of fields increases up to 2,800-3,000 m. a.s.l., followed by a reduction at higher altitudes [16]. Thus, the overall population size increases exponentially according to altitudes up to 2,800-2,900 m a.s.l., with a reduction at higher altitudes.

Our interpretation is further supported by the reduction in the multilocus LD that was seen from low to high altitudes, which was even non-significant in the Belg season, while at the extreme altitudes, the multilocus LD increased slightly. The reduction in the LD might also be explained by both the greater population size, which would increase the effective recombinations, and the selection, which will favour recombination that will produce novel multilocus genotypes.

The level of observed heterozygosity was very low (0.003), in agreement with the strict selfing nature of barley where the level of outcrossing is lower than 1-2% [4750]. Only three individuals, all located at high altitude, showed one or two heterozygous loci. Clearly, our study does not have the power to discriminate between different altitudes for the level of heterozygosity. However, a higher level of heterozygosity at high altitude might partially explain the pattern of LD in our study, because of a higher effective recombination rate at high altitude. To test this hypothesis a further study should be conducted using a larger sampling.


Landraces are a key component of agro-biodiversity, and they represent a crucial reservoir of genetic diversity for plant breeding. Moreover, in-situ conservation of landraces can provide a number of advantages, including the potential for adaptation to environmental changes because of an ongoing evolutionary process. Nevertheless, few studies have described the population structures of landraces or have analysed the roles of different evolutionary forces in the shaping of their genetic diversity. Here, we show that barley landraces from Ethiopia have highly variable local populations (farmer's fields) with large within-population diversity. These landraces also appear to be locally adapted, and the major driving force that has shaped their population structure is selection for adaptation along an altitudinal gradient. Moreover, our data suggest that the two-season system (which characterises barley cultivation in Ethiopia) and its effects on landrace population size have important roles in counterbalancing the homogenising effects of selection. Our study highlights the potential of barley landraces from Ethiopia as a source of useful alleles. They are also an ideal system to study the processes of adaptation and for identification of genes and genomic regions that have adaptive roles in a crop species, which can be achieved through association mapping and scanning for signatures of selection for molecular-diversity structure.


Plant materials

The plant materials were derived from a barley landrace collection that was carried out in 2005, in the North Shewa zone in the central highlands of Ethiopia [16]. This collection [16] was obtained by visiting the same three districts (Ankober, Mojanawadera and Tarmaber) in both the long (Meher) and the short (Belg) growing seasons. Within each farm visited, and according to the information from the farmers, all of the different landraces of barley that were grown and kept separated by the farmers during the seed selection process were collected separately from different fields and considered as different landrace populations (including those sampled in different seasons from the same farmer). In most cases (80 farmers), only one landrace was grown, while 10 of the farmers grew two landraces, and only two grew three landraces. For this reason, per farmer, more than one landrace population was sampled (on average, 1.15). Overall, the collection includes 106 barley landrace populations (fields), and Additional file 9 shows the collection site coordinates of these barley landrace populations. Within each field, 100 spikes (one spike per plant) were randomly sampled all along a diagonal of the field, with the plants sampled from 5-10 m apart. The geographical position of each field (as latitude, longitude and altitude) was determined using the Global Positioning System (GPS). Before threshing, 30 spikes per population were randomly sampled from the 100 spikes collected from each field. The morphological evaluation of these materials was based on eight morphological traits of the mature spikes (kernel-row number, spike density, lemma awn barbs, glume colour, lemma type, length of rachilla hair, kernel cover, and lemma colour), as reported by the Bioversity International Barley descriptors [51]; the details of the sampling strategy and morphological evaluation are given in [16].

From this collection, we randomly sampled two individual spikes for each population (farmer's field) collected from different individual plants. From each spike, a single seed was grown to the three-leaf stage and used for DNA extraction. Overall we analysed 212 genotypes from 106 barley landrace populations.

Genotypic data

Seven SSR markers (HVM20, Bmac0134, Bmag0013, HVM67, Bmac0113, Bmac0040, Bmac0156), as one marker per chromosome (see Additional file 10), were selected from Ramsay et al. [52] and used for the genetic characterisation of the 212 genotypes considered. The DNA was obtained from young leaves of single plants, using the miniprep extraction method of Doyle and Doyle [53]. The amplification conditions are reported in Additional file 10. The genotyping of the seven SSR markers was carried out with the ABI Prism 3100-Avant Genetic Analyser automatic sequencer, with GENESCAN 7.0 analysis software (PE Applied Biosystems, Foster City, CA, USA).

Statistical analysis

Level of polymorphism

The number of alleles (n a ), the average number of observed alleles per locus (n o ), the effective number of alleles per locus (n e , [54]), the Levene [55] observed heterozygosity (H o ) and the Nei's [35] unbiased genetic diversity (He) estimates based on allele frequencies, were calculated for each SSR locus. These were analysed as averages for each of the seasons (Belg and Meher), the districts (Ankober, Mojanawadera, Tarmaber), and the three altitude classes (<2,300, 2,300-2,800, >2,800 m a.s.l.), and for the whole sample, using the POPGENE software, version 1.31 [56]. As the number of alleles observed was highly dependent on the sample size, we also computed the allelic richness (R S , [57]), using the FSTAT software [58], a methodology that estimates the number of alleles independent of the sample size. The number of private alleles and their average frequencies for each of the above-mentioned groups were determined by inspection of the allele distributions. This computation was carried out also considering a minimum threshold frequency of 5%, to reduce the effects of sampling error [59]. The differences between seasons, districts, and altitude classes for the genetic diversity estimates (n a , n e , He and R S ) were tested using the Wilcoxon signed-ranks non-parametric test for two groups, arranged for paired observations (i.e. one pair of estimates for each locus) [60].

The partitioning of the genetic diversity was obtained using an analysis of molecular variance framework (AMOVA, [61]). This AMOVA analysis was performed using the Arlequin software, version 3.1 [62]. The divergence between seasons, districts, and altitude classes for each of the SSR loci was analysed, and the averages were quantified with FST [63] and RST (R-statistics, [64]) estimators. FST and RST differ in sensitivity when estimated on SSRs [65]; indeed FST can underestimate the magnitude of differentiation when populations are highly structured or are in a situation where the SSRs show high mutation rates; in contrast, RST is independent of the mutation rate, even if it has a high associated variance. We also considered the three altitude classes separately within each district, and we performed AMOVA analysis to compute the average unweighted FST estimates between different districts within the same altitude classes (<2,300, 2,300-2,800, >2,800 m a.s.l.). Similarly, the pairwise FST estimates between districts within the same altitude classes were computed.

All of the seven SSRs used were the same as those for the study of Russel et al. [44], who worked on barley landraces from several different agro-ecological areas of Syria and Jordan. Similarly, six (excluding HVM20) were common to the studies of Pandey et al. [36], who worked on Nepalese barley landraces, and Macaulay et al. [37], who worked on a set of modern barley varieties. This has allowed us to make comparisons among and between the diversity levels detected previously in Ethiopian, Syrian, Jordanian and Nepalese landraces, and with modern varieties of Hordeum vulgare.

Population structure

To further investigate the population structure of our sample, a Bayesian-model-based approach was used, as proposed by Pritchard et al. [40] and implemented in the software STRUCTURE, version 2.1 [41], to assign the genotypes into genetically structured groups. The software was run for presumed populations (K) from 1 to 10, following the admixture ancestry model. The run length was 100,000 MCMC repetitions and 100,000 burn-in periods, with 100 independent replicates for each K, to achieve consistent results. An ad-hoc statistic (ΔK, [43]) was used to estimate the "true" K number. The percentages of membership (q) of the individuals in each of the inferred K clusters were computed by an additional run using 1,000,000 MCMC repetitions and 1,000,000 burn-in periods.

An additional cluster analysis was performed with TESS, version 2.0 [42], a programme that introduces spatial correlation between individuals, in contrast to STRUCTURE, which assumes that all of the individuals are equally unrelated. The incorporation of a spatial component into the clustering model has the potential to determine if the clines provide a sensible description of the underlying pattern of variation [6668]. TESS implements a Bayesian clustering algorithm that uses a hidden Markov random field (HMRF) model to compute the proportion of individual genomes originating in K populations. The HMRF represents spatial connectivities as links in a network of individuals. Furthermore, it incorporates decay of the membership coefficient (q) correlation with distance, a property similar to isolation-by-distance. The network was automatically generated by the TESS programme using a Dirichlet tessellation obtained from the spatial coordinates of the samples. Considering the spatial distribution of our samples, we modified the network by removing several links, to account for potential geographical barriers. Runs were based on a burn-in period of 20,000 cycles followed by 30,000 iterations. One hundred replicates were performed for K values from 2 to 10, and for all of the runs the admixture model was used. In the analysis, we considered three values of the spatial dependence parameter (Ψ), 0.0, 0.6 and 1.0. This parameter weights the relative importance given to the spatial connectivities (Ψ = 0 recovers the model underlying STRUCTURE, while Ψ = 0.6 and 1.0 indicate moderate and strong values, respectively). For each run, the programme computed the Deviance Information Criterion (DIC), which is a model-complexity-penalised measure of how well the model fits the data. The lower DIC values represent the better fits for the data. Thus, we selected the 10 runs for each K that corresponded to the 10 lowest values of the DIC, and we averaged these values to determine the 'true' K number.

The TESS algorithm incorporates an additional regularisation feature that generally leads to a less ambiguous determination of K [42, 68]. Moreover, it can achieve an accuracy similar to that obtained with non-spatial methods, while using a smaller number of genetic markers [68]. Thus, we averaged the estimated q over the 10 runs with the smallest values of the DIC for the K value identified. We used the software CLUMPP, version 1.1 [69], which implements the Greedy algorithm, to allow for label switching and to decide which of the clusters of each run corresponded to a specific label.

The results of both of the STRUCTURE and TESS Bayesian clustering programmes were visualised using DISTRUCT, version 1.1 [70]. The output represents each individual as a vertical line, partitioned into K coloured segments, which represent the individually estimated membership fractions in the K clusters.

A simple linear regression model was used to test the relationships between the altitude and the q for the STRUCTURE and TESS clusters. The q values for STRUCTURE and TESS were used to investigate the relationships between the clusters identified by these two different programmes. The analysis was performed using Spearman's rho (ρ) for non-parametric correlation. Moreover, a step-wise multiple regression model was also performed, with the STRUCTURE clusters as dependent variables and the TESS clusters as independent variables. These analyses were carried out using the JMP 7 software (SAS Institute, Cary, USA).

Association between the genetic and morphological characterisations

To investigate the associations between the genetic and morphological characterisations, a contingency analysis was performed with the likelihood ratio chi-squared test, using the JMP 7 software. The proportion of the total uncertainty attributed to the model fit (R2) was also calculated. For this analysis, we considered the TESS genetic groups and only the genotypes assigned to one of the groups with a q higher than 0.70. No naked (hulless) barley was present among the 212 individuals analysed, thus the analysis was conducted considering seven morphological traits: kernel-row number, spike density, lemma awn barbs, glume colour, lemma type, length of rachilla hair, and lemma colour [16].

Spatial structure, clinal variation and landscape analysis

Spatial autocorrelation between spatial (geographical and altitude) and genetic distances were computed separately using Spatial Genetic Software (SGS), version 1.0 d [71]. These calculations were carried out using Moran's I [72, 73] for spatial distance classes (in metres), the dimension of which was 5,000 m for geographical distances (9 classes), and 102 m for altitude differences (9 classes). The sizes and numbers of the classes were fixed, to retain biological meaning and to guarantee at least 1,000 pairs of data points in each class. The significances of the observed average Moran's I values were assessed by comparing them with the corresponding values derived by randomly permuting the multilocus genotypes over the spatial coordinates of the samplings (500 times). The 99% confidence envelopes were estimated.

We separately obtain the Moran's I correlograms over all of the seven SSRs and for each locus. Moreover, we also designed the distogram between geographical distance and altitude difference by considering the former as a quantitative trait and the latter as a linear geographical distance, and using the city-block distance.

The AIS software [74] was used to calculate the associations between genetic distances among individuals and geographical and altitude distances, using a Mantel test [75]. The genetic distances between individuals were calculated following Equation of Miller [74], which is similar to Nei's [76] measure of genetic distance, except that allelic similarities are measured between individuals rather than populations.

To test for correlation between the SSR alleles and altitude, an analysis of variance (ANOVA) was performed with the Wilcoxon non-parametric test. Moreover, using the morphological data of the Tanto Hadado et al. [16] study (thus 106 populations, 30 individuals per population), we computed the frequencies of the characters for each population (character states with a frequency lower than 0.05 and higher than 0.95 were not included in the analysis), and we tested their relationships with altitude using the Spearman's rho (ρ) for non-parametric correlations. These analyses were carried out using the JMP 7 software.

Sliding-window analysis

To obtain a representation of the total genetic structure of the Ethiopian barley landraces from North Shewa that is as accurate as possible, we performed a sliding windows analysis. This approach is commonly used in ecology [77, 78], but, to the best of our knowledge, has not used in molecular analyses. This allowed us to represent the variations of the genetic diversity (He), the allelic richness (RS), and the linkage disequilibrium (LD) along the altitudinal cline.

We reconstructed populations of 20 landraces (fields), which corresponded to 40 individuals (windows); each window differed by only one landrace population (two individuals), that represented the step size, moving forward from the lowest to the highest 20 landrace populations along the altitudinal gradient. Thus, considering the step size of one population, the windows partially overlapped, sharing 38 individuals (19 landrace populations). This analysis was carried out separately for the two seasons. The altitudes for each window were obtained by averaging the altitudes of the populations (farmers' fields) included in the window. Thus the differences in altitudes between the windows varied along the altitudinal cline.

To examine the LD (non-random association between alleles at different loci) between our unlinked loci, a multilocus summary statistic of association between loci, rd, was calculated [79]. The significance of its deviation from zero was obtained using the MultiLocus software, version 1.2 [80], which performs shuffling of genotypes with each locus and the overall loci. The expectation of this statistic is that it is independent of the number of polymorphic loci within the population examined, allowing comparisons between different samples.

To test the significance of the trends observed for He and R S , we grouped the individuals following three altitudinal ranges for Belg, and four for Meher. The differences were tested between adjacent altitudinal ranges within each season using the Wilcoxon signed-ranks non-parametric test [60].

Authors' Information

This study is a component of the PhD thesis of Tesema Tanto Hadado, which focused on an analysis of phenotypic and molecular diversity in barley landraces from the central highlands of Ethiopia. All of the contributing authors are broadly interested in molecular evolution and ecological genetics, with particular interests in the conservation and use of biodiversity and crop genetic resources.


  1. 1.

    Matz SA: The Chemistry and Technology of Cereals as Food and Feed. Van Nostrand Reinhold: New York; 1991.

    Google Scholar 

  2. 2.

    von Bothmer R, Sato K, Komatsudam T, Yasuda S, Fischbeck G: The domestication of cultivated barley. Diversity in barley (Hordeum vulgare). Edited by: von Bothmer R, van Hintum T, Knüpffer H, Sato K. Amsterdam: Elsevier Science BV; 2003. 9-27. full_text.

    Google Scholar 

  3. 3.

    Asfaw Z: The barleys of Ethiopia. Genes in the Field: On-farm Conservation of Crop Diversity. Edited by: Brush SB. Boca Raton, Florida, USA: Lewis Publisher; 2000.77-108.

    Google Scholar 

  4. 4.

    Zhu Y, Chen H, Fan J, Wang Y, Li Y, Chen J, Fan J, Yang S, Hu L, Leung H, Mew TW, Teng PS, Wang Z, Mundt CC: Genetic diversity and disease control in rice. Nature. 2000, 406: 718-722. 10.1038/35021046.

    PubMed  CAS  Article  Google Scholar 

  5. 5.

    Lynch M, Walsh B: Genetics and Analysis of Quantitative Traits. Sunderland, MA, USA: Sinauer Associates; 1997.

    Google Scholar 

  6. 6.

    Thornsberry JM, Goodman MM, Doebley J, Kresovich S, Nielsen D, Buckler ES: Dwarf8 polymorphisms associate with variation in flowering time. Nature Genet. 2001, 28: 286-289. 10.1038/90135.

    PubMed  CAS  Article  Google Scholar 

  7. 7.

    Flint-Garcia SA, Thuillet AC, Yu J, Pressoir G, Romero SM, Mitchell SE, Doebley J, Kresovich S, Goodman MM, Buckler ES: Maize association population: a high-resolution platform for quantitative trait locus dissection. Plant J. 2005, 44: 1054-1064. 10.1111/j.1365-313X.2005.02591.x.

    PubMed  CAS  Article  Google Scholar 

  8. 8.

    Mazzucato A, Papa R, Bitocchi E, Mosconi P, Nanni L, Negri V, Picarella ME, Siligato F, Soressi G, Tiranti B, Veronesi F: Genetic diversity, structure and marker-trait associations in a collection of Italian tomato (Solanum lycopersicum L.) landraces. Theor Appl Genet. 2008, 116: 657-669. 10.1007/s00122-007-0699-6.

    PubMed  Article  Google Scholar 

  9. 9.

    Kohn MH, Pelz HJ, Wayne RK: Natural selection mapping of the warfarin-resistance gene. Proc Natl Acad Sci USA. 2000, 97: 7911-7915. 10.1073/pnas.97.14.7911.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  10. 10.

    Vigouroux Y, McMullen M, Hittinger CT, Houchins K, Schulz L, Kresovich S, Matsuoka Y, Doebley J: Identifying genes of agronomic importance in maize by screening microsatellites for evidence of selection during domestication. Proc Natl Acad Sci USA. 2002, 99: 9650-9655. 10.1073/pnas.112324299.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  11. 11.

    Luikart G, England PR, Tallmon D, Jordon S, Taberlet P: The power and promise of population genomics: from genotyping to genome typing. Nature Rev Genet. 2003, 4: 981-994. 10.1038/nrg1226.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Papa R, Bellucci E, Rossi M, Leonardi S, Rau D, Gepts P, Nanni L, Attene G: Tagging the signatures of domestication in common bean (Phaseolus vulgaris) by means of pooled DNA samples. Ann Botany. 2007, 100: 1039-1051. 10.1093/aob/mcm151.

    CAS  Article  Google Scholar 

  13. 13.

    Bitocchi E, Nanni L, Rossi M, Rau D, Bellucci E, Giardini A, Buonamici A, Vendramin GG, Papa R: Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy. Mol Ecol. 2009, 18: 603-621. 10.1111/j.1365-294X.2008.04064.x.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    Engels JMM: Genetic diversity in Ethiopia in relation to altitude. Genet Resour Crop Evol. 1994, 41: 67-73. 10.1007/BF00053050.

    Article  Google Scholar 

  15. 15.

    Demissie A, Bjørnstad A: Phenotypic diversity of Ethiopian barley in relation to geographical regions, altitudinal range, and agro-ecological zones as an aid to germplasm collection and conservation strategy. Hereditas. 1996, 124: 17-29. 10.1111/j.1601-5223.1996.00017.x.

    Article  Google Scholar 

  16. 16.

    Tanto Hadado T, Rau D, Bitocchi E, Papa R: Genetic diversity of barley (Hordeum vulgare L.) landraces from the central highlands of Ethiopia: comparison between the 'Belg' and 'Meher' growing seasons using morphological traits. Genet Resour Crop Evol. 2009, 56: 1131-1148. 10.1007/s10722-009-9437-z.

    Article  Google Scholar 

  17. 17.

    Bekele E: A differential rate of regional distribution of barley flavonoid patterns in Ethiopia, and a view on the center of origin of barley. Hereditas. 1983, 98: 269-280. 10.1111/j.1601-5223.1983.tb00605.x.

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    Demissie A, Bjørnstad A: Geographical, altitude and agro-ecological differentiation of isozyme and hordein genotypes of landrace barleys from Ethiopia: implications to germplasm conservation. Genet Resour Crop Evol. 1997, 44: 43-55. 10.1023/A:1008686009866.

    Article  Google Scholar 

  19. 19.

    Wiberg A: Sources of resistance to powdery mildew in barley. Hereditas. 1974, 78: 1-40. 10.1111/j.1601-5223.1974.tb01426.x.

    PubMed  CAS  Article  Google Scholar 

  20. 20.

    Qualset CO: Sampling germplasm in a centre of diversity: an example of disease resistance in Ethiopian barley. Crop Genetic Resources for Today and Tomorrow. Edited by: Frankel OH, Hawkes JG. Cambridge University Press; 1975. 81-96.

    Google Scholar 

  21. 21.

    Zhang Q, Webster RK, Allard RW: Geographical distribution and association between resistance to four races of Rhynchosporium secalis. Phytopathology. 1987, 77: 352-357. 10.1094/Phyto-77-352.

    Article  Google Scholar 

  22. 22.

    Fukuyama T, Takeda H: Survey of resistance to Scald in world collection of barley. Japan J Breed. 1992, 42: 761-768.

    Article  Google Scholar 

  23. 23.

    Jørgensen JH: Discovery, characterization and exploitation of Mlo powdery mildew resistance in barley. Euphytica. 1992, 63: 141-152. 10.1007/BF00023919.

    Article  Google Scholar 

  24. 24.

    Alemayehu F: Genetic variation between and within Ethiopian barley landraces with emphasis on durable resistance. PhD Thesis. Landbouw Universiteit Wageningen, Holland; 1995.

    Google Scholar 

  25. 25.

    Yitbarek S, Berhane L, Fikadu A, 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-423. 10.1111/j.1439-0523.1998.tb01966.x.

    Article  Google Scholar 

  26. 26.

    Munck L, Karisson KE, Hagberg A, Eggum BO: Gene for improved nutritional value in barley seed protein. Science. 1970, 168: 985-987. 10.1126/science.168.3934.985.

    PubMed  CAS  Article  Google Scholar 

  27. 27.

    Lance RCM, Nilan RA: Screening for low acid soluble β-glucan barleys. Barley Genet Newsl. 1980, 10: 41.

    Google Scholar 

  28. 28.

    Asfaw Z: Relationships between spike morphology, hordeins and altitude within Ethiopian barley. Hereditas. 1989, 110: 203-209. 10.1111/j.1601-5223.1989.tb00782.x.

    Article  Google Scholar 

  29. 29.

    Alemayehu F, Parlevliet JE: Variation between and within Ethiopian barley landraces. Euphytica. 1997, 94: 183-189. 10.1023/A:1002946004268.

    Article  Google Scholar 

  30. 30.

    Asfaw Z: Variation in hordein polypeptide pattern within Ethiopian barley, Hordeum vulgare L (Poaceae). Hereditas. 1989, 110: 185-191. 10.1111/j.1601-5223.1989.tb00779.x.

    Article  Google Scholar 

  31. 31.

    Haldane JBS: The theory of a cline. J Genet. 1948, 48: 277-284. 10.1007/BF02986626.

    PubMed  CAS  Article  Google Scholar 

  32. 32.

    Slatkin M: Gene flow and selection in a cline. Genetics. 1973, 75: 733-756.

    PubMed  CAS  PubMed Central  Google Scholar 

  33. 33.

    Endler JA: Geographic Variation, Speciation, and Clines. Princeton, New Jork: Princeton University Press; 1977.

    Google Scholar 

  34. 34.

    Stortz FJ, Dubach JM: Natural selection drives altitudinal divergence at the albumin locus in deer mice (Peromyscus maniculatus). Evolution. 2004, 58: 1342-1352.

    Article  Google Scholar 

  35. 35.

    Nei M: Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978, 89: 583-590.

    PubMed  CAS  PubMed Central  Google Scholar 

  36. 36.

    Pandey M, Wagner C, Friedt W, Ordon F: Genetic relatedness and population differentiation of Himalayan hulless barley (Hordeum vulgare L) landraces inferred with SSRs. Theor Appl Genet. 2006, 113: 715-729. 10.1007/s00122-006-0340-0.

    PubMed  CAS  Article  Google Scholar 

  37. 37.

    Macaulay M, Ramsay L, Powel W, Waugh R: A representative highly informative genotyping set of barley SSRs. Theor and Appl Genet. 2001, 102: 801-809. 10.1007/s001220000487.

    CAS  Article  Google Scholar 

  38. 38.

    Serre D, Pääbo S: Evidence for gradients of human genetic diversity within and among continents. Genome Res. 2004, 14: 1679-1685. 10.1101/gr.2529604.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  39. 39.

    Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, Feldman MW: Clines, clusters, and the effect of study design on the inference of human population structure. PLoS Genet. 2005, 1: e70-10.1371/journal.pgen.0010070.

    PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Pritchard JK, Stephens M, Donnelly P: Inference of population structure from multilocus genotype data. Genetics. 2000, 155: 945-959.

    PubMed  CAS  PubMed Central  Google Scholar 

  41. 41.

    Falush D, Stephens M, and Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.

    PubMed  CAS  PubMed Central  Google Scholar 

  42. 42.

    Chen C, Durand E, Forbes F, François O: Bayesian clustering algorithms ascertaining spatial population structure: a new computer program and a comparison study. Mol Ecol Notes. 2007, 7: 747-756. 10.1111/j.1471-8286.2007.01769.x.

    Article  Google Scholar 

  43. 43.

    Evanno G, Reganut E, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.

    PubMed  CAS  Article  Google Scholar 

  44. 44.

    Russell JR, Booth A, Fuller JD, Baum M, Ceccarelli S, Grando S, Powell W: Patterns of polymorphism detected in the chloroplast and nuclear genomes of barley landraces sampled from Syria and Jordan. Theor Appl Genet. 2003, 107: 413-421. 10.1007/s00122-003-1261-9.

    PubMed  CAS  Article  Google Scholar 

  45. 45.

    Sokal RR, Jacquez GM, Wooten MC: Spatial autocorrelation analysis of migration and selection. Genetics. 1989, 121: 845-855.

    PubMed  CAS  PubMed Central  Google Scholar 

  46. 46.

    Pourcel L, Routaboul JM, Cheynier V, Lepiniec L, Debeaujon I: Flavonoid oxidation in plants: from biochemical properties to physiological functions. Trends Plant Sci. 2007, 12: 29-36. 10.1016/j.tplants.2006.11.006.

    PubMed  CAS  Article  Google Scholar 

  47. 47.

    Brown AHD, Zohary D, Nevo E: Outcrossing rates and heterozygosity in natural populations of Hordeum spontaneum Koch in Israel. Heredity. 1978, 41: 49-62. 10.1038/hdy.1978.63.

    Article  Google Scholar 

  48. 48.

    Chaudhary HR, Jana S, Acharya SN: Outcrossing rates in barley populations in the Canadian prairies. Can J Genet Cytol. 1980, 22: 353-360.

    Article  Google Scholar 

  49. 49.

    Parzies HK, Spoor W, Ennos RA: Outcrossing rates of barley landraces from Syria. Plant Breed. 2000, 119: 520-522. 10.1046/j.1439-0523.2000.00532.x.

    Article  Google Scholar 

  50. 50.

    Abdel-Ghani AH, Parzies HK, Omary A, Geiger HH: Estimating the outcrossing rate of barley landraces and wild barley populations collected from ecologically different regions of Jordan. Theor Appl Genet. 2004, 109: 588-595. 10.1007/s00122-004-1657-1.

    PubMed  Article  Google Scholar 

  51. 51.

    International Plant Genetic Resources Institute: Bioversity International Descriptor for Barley (Hordeum vulgare L.). Rome, Italy. 1994

    Google Scholar 

  52. 52.

    Ramsay L, Macaulay M, Ivanissevich S, MacLean K, Cardle L, Fuller J, Edwards KJ, Tuvesson S, Morgante M, Massari A, Maestri E, Marmiroli N, Sjakste T, Ganal M, Powell W, Waugh R: A simple sequence repeat-based linkage map of barley. Genetics. 2000, 156: 1997-2005.

    PubMed  CAS  PubMed Central  Google Scholar 

  53. 53.

    Doyle JJ, Doyle JL: A rapid DNA isolation procedure from small quantities of fresh leaf tissues. Phytochem Bull. 1987, 19: 11-15.

    Google Scholar 

  54. 54.

    Kimura M, Crow JF: The number of alleles that can be maintained in a finite population. Genetics. 1964, 49: 725-738.

    PubMed  CAS  PubMed Central  Google Scholar 

  55. 55.

    Levene H: On a matching problem arising in genetics. Ann Math Stat. 1949, 20: 91-94. 10.1214/aoms/1177730093.

    Article  Google Scholar 

  56. 56.

    Yeh FC, Boyle TJB: POPGENE, Version 1.1. Department of Renewable Resources, University of Alberta, Edmonton, Alberta; 1997.

    Google Scholar 

  57. 57.

    El Mousadik A, Petit RJ: High level of genetic differentiation for allelic richness among populations of the argan tree (Argania spinosa L. Skeels) endemic to Morocco. Theor Appl Genet. 1996, 92: 832-839. 10.1007/BF00221895.

    PubMed  CAS  Article  Google Scholar 

  58. 58.

    Goudet J: FSTAT a program to estimate and test gene diversities and fixation indices (Version 2.9.3). 2001, []

    Google Scholar 

  59. 59.

    Mohammad A, Al-Rabab' AH, Williams CG: An ancient bottleneck in the Lost Pines of central Texas. Mol Ecol. 2004, 13: 1075-1084. 10.1111/j.1365-294X.2004.02142.x.

    Article  Google Scholar 

  60. 60.

    Sokal RR, Rohlf FJ: Biometry: the Principles and Practice of Statistics in Biological Research 3rd edition. Freeman WH and Co. New York; 1995.

    Google Scholar 

  61. 61.

    Excoffier L, Smouse P, Quattro J: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mithocondrial DNA restriction data. Genetics. 1992, 131: 479-491.

    PubMed  CAS  PubMed Central  Google Scholar 

  62. 62.

    Laval LG, Schneider S: Arlequin ver 30: An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.

    PubMed Central  Google Scholar 

  63. 63.

    Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.

    Article  Google Scholar 

  64. 64.

    Slatkin M: A measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995, 139: 457-462.

    PubMed  CAS  PubMed Central  Google Scholar 

  65. 65.

    Balloux F, Moulin NL: The estimation of population differentiation with microsatellite markers. Mol Ecol. 2002, 11: 155-165. 10.1046/j.0962-1083.2001.01436.x.

    PubMed  Article  Google Scholar 

  66. 66.

    Wasser S, Shedlock A, Comstock K, Ostrander E, Mutayoba B, Stephens M: Assigning African elephants DNA to geographic region of origin: applications to the ivory trade. Proc Natl Acad Sci USA. 2004, 101: 14847-14852. 10.1073/pnas.0403170101.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  67. 67.

    Guillot G, Estoup A, Mortier F, Cosson J: A spatial statistical model for landscape genetics. Genetics. 2005, 170: 1261-1280. 10.1534/genetics.104.033803.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  68. 68.

    François O, Ancelet S, Guillot G: Bayesian clustering using hidden Markov random fields in spatial population genetics. Genetics. 2006, 174: 805-816. 10.1534/genetics.106.059923.

    PubMed  PubMed Central  Article  Google Scholar 

  69. 69.

    Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806. 10.1093/bioinformatics/btm233.

    PubMed  CAS  Article  Google Scholar 

  70. 70.

    Rosenberg NA: DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004, 4: 137-138. 10.1046/j.1471-8286.2003.00566.x.

    Article  Google Scholar 

  71. 71.

    Degen B, Petit R, Kremer A: SGS -- Spatial Genetic Software: a computer program for analysis of spatial genetic and phenotypic structures of individuals and populations. J Hered. 2001, 92: 447-448. 10.1093/jhered/92.5.447.

    PubMed  CAS  Article  Google Scholar 

  72. 72.

    Moran PAP: Notes on continuous stochastic phenomena. Biometrika. 1950, 37: 17-23.

    PubMed  CAS  Article  Google Scholar 

  73. 73.

    Sokal RR, Oden NL: Spatial autocorrelation in biology: I Methodology. Biol J Linnean Soc. 1978, 10: 199-228. 10.1111/j.1095-8312.1978.tb00013.x.

    Article  Google Scholar 

  74. 74.

    Miller MP: Alleles in space (AIS): computer software for the joint analysis of interindividual spatial and genetic information. J Hered. 2005, 96: 722-724. 10.1093/jhered/esi119.

    PubMed  CAS  Article  Google Scholar 

  75. 75.

    Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Research. 1967, 27: 209-220.

    PubMed  CAS  Google Scholar 

  76. 76.

    Nei M: Genetic diversity and the neutral mutation theory. Heredity. 1983, 51: 531.

    Google Scholar 

  77. 77.

    Craig P, Haslett J, Unwin A, Wills G: Moving statistics - an extension of brushing for spatial data. Computing Science and Statistics. Proceedings of the 21st Symposium on the Interface. Berlin: Springer Verlag, 1989:170-174.

    Google Scholar 

  78. 78.

    Haining R: Spatial Data Analysis Theory and Practice. Cambridge University Press: Cambridge, UK 2003.

    Google Scholar 

  79. 79.

    Burt A, Koufopanou V, Taylor JW: Population genetics of human-pathogenic fungi. The molecular epidemiology of infectious diseases. Edited by: Thompson RCA. London: Chapman and Hall, 2000:229-244.

    Google Scholar 

  80. 80.

    Agapow PM, Burt A: Indices of multilocus linkage disequilibrium. Mol Ecol Notes. 2001, 1: 101-102. 10.1046/j.1471-8278.2000.00014.x.

    CAS  Article  Google Scholar 

Download references


We are grateful to the Ethiopian farmers for their kind help and to whom we dedicate this study. The authors would like to particularly thank Stefano Leonardi, A. H. D. Brown, T. Hodgkin and D. I. Jarvis, for their valuable advice. We thank the guides, Asefa Mekonnen, Mengesha Ergeta, Getachew Adere, and Legesse Bejiga, for their support during the survey and collection work. The authors also thank Bioversity International and the Institute of Biodiversity Conservation of Ethiopia for their financial and material support that were provided during the field work in Ethiopia. T. Tanto Hadado thanks the Bioversity International for the support for his stay in Ancona, Italy.

Author information



Corresponding author

Correspondence to Roberto Papa.

Additional information

Authors' contributions

RP conceived, designed and coordinated the study; TTH and DR carried out the barley landrace collection in Ethiopia; EB performed the DNA extraction and molecular analysis; TTH DR EB and RP analysed the data and contributed to the drafting and writing of the manuscript. All of the authors have read and approved this version of the manuscript.

Electronic supplementary material

Additional file 1: Summary statistics computed for each locus considering the two seasons, three districts, and three altitude classes, and for the whole sample. (DOC 136 KB)

Additional file 2: Divergence (FST and RST) estimates for each of the SSR loci analysed, computed considering the seasons, districts, and altitude classes. (DOC 47 KB)

Additional file 3: Non-parametric correlation (Spearman's rho) between TESS and STRUCTURE clusters. (DOC 25 KB)

Additional file 4: Step-wise multiple regression analysis for STRUCTURE and TESS. The model was performed considering the STRUCTURE clusters as dependent variables and the TESS clusters as independent variables, from the data illustrated in Figures 4 and 5. (DOC 31 KB)

Additional file 5: Relationships between the molecular data and altitude. (DOC 43 KB)

Additional file 6: Types and frequencies of each of the seven qualitative traits, computed considering the genotypes that were assigned to one of the TESS clusters with a coefficient of membership ( q ) higher than 0.70 (from the data illustrated in Figure 5). (DOC 88 KB)

Additional file 7: Non-parametric correlation (Spearman's rho) between the TESS clusters for their morphological traits (from the data illustrated in Figure 5). (DOC 46 KB)

Additional file 8: Spatial autocorrelation analysis between the geographical and genetic distances, performed separately for the three altitude classes. See legend to Figure 7 for colour key. (PDF 24 KB)

Additional file 9: Collection site coordinates of the barley landraces analysed. (DOC 3 MB)

Additional file 10: List of SSRs used in the present study. (DOC 50 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Tanto Hadado, T., Rau, D., Bitocchi, E. et al. Adaptation and diversity along an altitudinal gradient in Ethiopian barley (Hordeum vulgare L.) landraces revealed by molecular analysis. BMC Plant Biol 10, 121 (2010).

Download citation


  • Altitudinal Gradient
  • Deviance Information Criterion
  • Spike Density
  • Barley Landrace
  • Landrace Population