BACKGROUND
cultivated barley , the domesticated form of hordeum spontaneum c. koch, was one of the founder crops of agriculture in western asia, first appearing in the archaeological record in the 8th and 7th millennia bc  <cit> .
barley was also one of the principal crops that accompanied the spread of agriculture into europe during the 6th and 5th millennia bc.
today it is grown throughout the continent, mainly for animal feed and malt for brewing  <cit> .
until the introduction of modern cultivars in the mid-20th century, european barley comprised a large number of landraces, each of these a 'dynamic population or populations of a cultivated plant that has historical origin, distinct identity and lacks formal crop improvement, as well as often being genetically diverse, locally adapted and associated with traditional farming systems'  <cit> .
many landraces died out during the last century, but seeds representing several thousand types from all parts of europe are available from germplasm collections  <cit> .
increasingly, these landrace collections are being looked on as important sources of germplasm with which to enrich the genepool of modern barley cultivars  <cit> .
exploitation of these landraces in modern crop breeding requires understanding not only of their phenotypic attributes and environmental adaptations, but also their evolutionary relationships.
this information would enable landraces with similar valuable traits but different genetic backgrounds to be identified, providing alternative strategies for the introduction of the traits into modern germplasm.
an understanding of the relationships between different landraces might also provide information relating to the spread and subsequent development of barley cultivation in europe.
this possibility has been recognized for some time  <cit> , but has not been extensively explored.
recent studies have suggested, however, that european barley landraces are differentiated into genetically defined populations.
for example, dna sequencing of the photoperiod response gene ppd-h <dig> in european barley landraces has revealed three distinct groups of haplotypes, two  with the daylight responsive phenotype found in most wild barleys, and one  with the derived nonresponsive phenotype that confers adaptation to the long annual growth season of northern europe  <cit> .
multilocus studies have also been carried out, avoiding the risks of inferring population history from single gene data.
population structure, linked at least in part with ecogeography and/or agronomic factors, has been revealed by analysis of microsatellites in barley landraces from iberia  <cit> , the levant  <cit>  and the himalayas  <cit> .
a relationship between microsatellite genotypes and ecogeography has also been demonstrated at the microscale for wild barley  <cit> .
in this paper we show that european barley landraces can be divided into populations based on their microsatellite genotypes.
we assess the extent to which the population structure can be explained by human selection, environmental adaptation, geographical partitioning occurring during the initial spread of barley cultivation into europe, and/or the later introduction of daylength-nonresponsive landraces.
RESULTS
microsatellite genotypes
we studied  <dig> accessions of cultivated barley  and typed  <dig> microsatellites .
each of these microsatellites displayed variability among the landraces that were tested, the number of alleles observed per locus ranging between 2- <dig> with a mean of  <dig>  .
pic values varied between  <dig> - <dig>  , in broad agreement with values previously obtained for spanish barley landraces  <cit> .
there were no significant differences between the overall diversities of two- and six-rowed barley landraces, though several individual microsatellites did show significant differences when these two sets of landraces were compared.
missing data  varied from  <dig> - <dig> % .
missing data, which usually arise when a landrace has a sequence polymorphism within the annealing site for one of the two primers used to amplify the microsatellite, will lead to underestimates of overall diversity, but are not an issue for subsequent data analysis because structure is able to compensate for gaps in the overall dataset.
identification of populations
to assess whether the barley landraces could be divided into populations, structure was used to calculate the probability distribution pr for values of k from 2- <dig>  where x is the genotypes of the sampled landraces and k is the number of populations.
multiple runs were carried out and population assignments tested for reproducibility by mutual correlation of the q-matrix outputs  <cit> .
q-matrices were accepted as reproducible if the allocations to each population could be identified as similar in a pair of q-matrices and were highly correlated .
for k =  <dig>   <dig>   <dig>   <dig>   <dig>   <dig> and  <dig> reproducibility was achieved after duplicate runs, while for k =  <dig>   <dig>   <dig>   <dig>   <dig>   <dig> and  <dig> two reproducible results were obtained after a third run.
the population assignments at k =  <dig> were symmetrical, each accession being given an equal allocation  to every population, and reproducible results could not be obtained for values of k from 17- <dig>  these results indicate that the dataset displays population structure and that the most likely value of k is ≤ <dig>
three methods were used to identify the most likely value of k more accurately.
the first method was based on the prediction that the upper limit of k is indicated by the smallest value that captures the major structure in the data  <cit> .
a plot of lnpr against k was therefore made for k =  <dig> to  <dig> .
this plot shows lnpr increasing until k = 12- <dig>  after which the line reaches a plateau.
this observation suggests that the most likely value of k is ≤ <dig>  in the second method the rate of change of lnpr for successive values of k was plotted against k .
the position of a plateau in a plot of this type indicates the most likely value of k  <cit> .
examination of figure 1b therefore suggests that k = 9- <dig>  the third method used to identify the most likely number of populations assumes that genetic structure should relate to phenotype  <cit> .
data on spring or winter growth habit, two-row or six-row ear morphology and caryopsis structure were converted into binary characters and logistic regressions of these phenotypes against the q-matrices for k = 2- <dig> carried out with the r statistics package.
in this type of analysis, the lowest values for akaike's information criterion  identify the range of k where population structure is best correlated with phenotype.
for growth habit, the lowest aic values were found for k = 9- <dig>  for row number at k =  <dig>  and for caryopsis structure at k = 8- <dig> .
taken together, the three analyses summarized in figures  <dig>   <dig> suggest that the most likely value of k for the barley landraces is between  <dig> and  <dig>  while we do not attempt to assign a 'correct' value for k, we chose a value of  <dig> as a starting point for examination of population structure.
a graphical representation of the population structure at k =  <dig> revealed that each population was partially admixed with other populations, the overall degree of admixture being similar for each one .
half the landraces  had a proportional membership of ≥ <dig>  for their primary population, and only  <dig> had a primary proportional membership of <  <dig> .
the relationships between these populations were studied in two ways.
first, the microsatellite data were used to construct a neighbour-joining tree and those accessions with a proportional membership of ≥ <dig>  in their primary population marked on the tree .
accessions belonging to a single population clustered together, except for population  <dig>  and to a lesser extent  <dig>  which appeared in two distinct regions of the tree.
populations 1- <dig> and some members of  <dig> and  <dig> grouped close to one another, as did population  <dig> alongside the remainder of  <dig> and  <dig>  populations  <dig>   <dig> and  <dig> each formed a separate group, although population  <dig> was split into two parts each with its own deep root towards the base of the tree.
to explore further the population relationships, the groups resulting from structure analyses with preset k values of 4- <dig> were investigated.
again, only those landraces that displayed a proportional membership of ≥ <dig>  in their primary populations were considered.
comparing the allocations for each of these landraces as k increased revealed a hierarchical structure to the landrace populations .
landraces that were grouped together at k =  <dig> were still grouped together at k = 7- <dig>  a subset of landraces from two of the populations at k =  <dig>  a and b, were placed in a shared population at k = 5- <dig>  then at higher k values re-sorted into distinct lineages that corresponded to their original allocation between populations a and b. these landraces made up populations 1- <dig> at k =  <dig>  the remaining two populations at k =  <dig>  c and d, both gave rise to lineages that remained distinct as k increased to  <dig>
phenotypic features of the populations
for those accessions for which data were available, ear row number, caryopsis structure, growth habit and daylength responsiveness were compared with the population structure .
two-rowed barley predominated in populations 1- <dig>  and six-row barley predominated in populations 5- <dig> , these differences being significant .
naked barleys were predominant in population  <dig> , and formed a sizeable minority of population  <dig> , but were virtually absent from the other seven populations .
winter growth habit was predominant only in population  <dig> .
spring growth habit accounted for > 90% of the accessions in populations 1- <dig> and  <dig> according to passport data, and in population  <dig> also according to vrn genotype.
the remaining two populations  comprised a mixture of spring and winter barleys.
the discrepancies within populations  <dig> and  <dig> between the passport data for growth habit and the habit predicted from the vrn genotype are due to many of these accessions coming from the mediterranean region, where barley with a facultative growth habit is favoured.
these barleys are autumn sown, able to survive mild winter conditions, and flower in the late spring.
they are therefore recorded as having a winter growth habit although they display the spring genotype when the vrn genes are typed.
of
a top row, number of landraces with each phenotype, with the total that could be phenotyped in brackets; bottom row, percentage of landraces with each phenotype
for  <dig> accessions, daylength responsiveness was deduced from the ppd-hi genotype  <cit>  .
all of the accessions of populations  <dig>   <dig> and  <dig> that were typed, as well as 95% of those in population  <dig>  were daylength nonresponsive, and all those typed from population  <dig> and 91% of those from population  <dig> were responsive.
the other four populations were mixed, ranging from 20-91% responsive.
where possible, those accessions with the daylength responsive genotype were placed in group b or c  <cit> .
group b predominated  only in population  <dig> and ranged from 0-13% in the other populations.
geographical distributions of the populations
all  <dig> accessions were included in the geographical analyses.
the mean centre and standard distance for each population  were calculated, along with the mean pairwise geographical distance within each population to assess the degree of clustering .
all of the landraces were included.
landraces of each population were placed on a map of europe and standard deviation ellipses plotted .
visual inspection of the maps indicated that the populations had non-identical distribution patterns.
three populations  were distributed over a relatively broad area of central and western europe.
population  <dig> had a similar but more western distribution, largely due to a preponderance of british landraces.
populations  <dig> and  <dig> centred on switzerland and the carpathian mountains, with population  <dig> tightly clustered in a small area of west-central switzerland.
population  <dig> was mainly located in the balkan region of southeast europe, population  <dig> predominated in north central europe including scandinavia and the baltic states, and population  <dig> clustered in the mediterranean region.
to assess the apparent differences between these distribution patterns, the great circle distances between pairs of populations were tested for significance using student's t-test  against a null hypothesis that the mean centres were  <dig> km apart.
this analysis suggested that  <dig> of the  <dig> possible population pairs had a highly significant  difference in their mean centres, and a further four  had a significant  difference.
the only non-significant differences were between paired populations  <dig> and  <dig>   <dig> and  <dig>   <dig> and  <dig>  and  <dig> and  <dig>
the relationship between the population structure and environment was investigated by comparing the between-population variance and within-population variance for a range of climatic factors.
this analysis was carried out separately for the spring and winter accessions  as these different growth habits might be expected to result in different patterns of environmental adaptation.
for randomized populations, the f ratio for between-population variance compared with within-population variance was close to one for each of the climate variables at all months of the year.
when all the climate variables were combined, the between-population variance for spring barleys was significantly higher than within-population variance for each month, with the highest values occurring during the growing season for these accessions .
between-population variance was also significantly greater than within-population variance for winter barleys when all climate variables were combined  but with no obvious seasonal trend.
discussion
identification of populations
structure was used to determine whether european barley landraces can be divided into populations based on their microsatellite genotypes.
in this context, a population is defined as a group of individuals that share a characteristic set of allele frequencies at the loci that are studied  <cit> .
structure places individuals in populations in such a way as to minimize within-population deviations from hardy-weinberg equilibrium and linkage equilibrium.
it therefore assumes that the individuals are fully outcrossing, and modelling studies have suggested that with partially inbreeding taxa the results of structure analysis can be spurious  <cit> .
although cultivated barley has an outcrossing rate of less than 2%  <cit> , the outcomes of previous structure analyses of genetic data from barley have been supported by other analyses of the same data and are in agreement with the conclusions of earlier and later work .
in these projects it therefore appears that structure identified authentic populations despite barley's low outcrossing rate.
an explanation might be provided by other simulations which have shown that, over hundreds of generations, the pattern of multilocus marker inheritance in a population of plants displaying 2% outbreeding is indistinguishable from that displayed by a panmictic population  <cit> .
the results of our structure analysis were reproducible for values of k up to  <dig> and, using standard methods, we concluded that the most likely number of populations was between 8- <dig>  the population structure was hierarchical, and from k = 7- <dig> the only effect of each incremental increase in k was to subdivide an existing population with the memberships of the other populations remaining unchanged .
this consistent pattern of population assignment indicates that the results of the structure analysis were not spurious.
to further assess the validity of the analysis, we constructed a neighbour-joining tree for all  <dig> accessions and marked the positions of those accessions with a proportional membership of ≥ <dig>  in their primary population at k =  <dig> .
accessions belonging to a single population clustered together, with the exception of populations  <dig> and  <dig>  whose members are distributed in two parts of the tree.
the tree topology therefore provides independent support for the structure results, and also suggests that it is reasonable to use k =  <dig> as the basis for interpretation of the population structure.
association between phenotype and population structure
for a domesticated plant such as cultivated barley, one possible way in which population structure could arise is as a result of selection for particular phenotypic traits.
the phenotypes of greatest agronomic importance in modern farming are ear row number, caryopsis structure , growth habit and flowering time.
wild barley has a two-rowed ear, each spikelet having a fertile central floret flanked by two infertile laterals which, when combined with a long awn, takes on an arrowlike form that is an effective aid to seed dispersal and burial  <cit> .
many cultivated barleys retain this ancestral head structure but in the derived six-rowed form the two lateral florets are fertile.
six-rowed barley is more often used as an animal or human feed, whereas two-rowed barley is favoured for malting and brewing.
wild and most cultivated barleys have hulled grains where the outer lemma and inner palea adhere to the pericarp epidermis at maturity.
this form is favoured by brewers because the hull debris aids wort filtration, whereas the free-threshing 'naked' varieties are preferred when barley is grown for direct human consumption  <cit> .
wild barley has a winter growth habit, meaning that it requires vernalization - exposure to a prolonged period of cold - in order to promote subsequent flowering.
the majority of european landraces lack this requirement and have a spring growth habit, where plants avoid periods of cold weather by completing their growth cycle during a single season, rather than over wintering as plants in a vegetative state  <cit> .
finally, most wild barleys display a photoperiod response that triggers flowering early in the season, before the conditions become too dry for further vegetative growth.
many cultivated barleys, especially landraces from northern europe, are daylength nonresponsive, and so continue vegetative growth until flowering later in the summer  <cit> , allowing them to take advantage of the longer growing season in northern europe.
we compared phenotypic data for ear row number, caryopsis structure, growth habit and flowering time with the population memberships .
populations 1- <dig> display a similar set of phenotypic features, most of these accessions being two-rowed , hulled , spring habit , and daylength nonresponsive .
populations  <dig> and  <dig>  which contained a high proportion of accessions with the facultative growth habit, also show some similarities when other phenotypes are considered, being largely six-row , hulled  and daylength responsive .
in population  <dig>  however, the majority of the responsive accessions were members of group c , whereas in population  <dig> the majority were group b , these two groups having distinct evolutionary histories  <cit> .
no other similarities between the range of phenotypes displayed by different populations were apparent.
these comparisons indicate that there is a degree of association between phenotype and population structure, suggesting that selection may have played a role in the origin and/or maintenance of one or more populations.
association between climatic factors and population structure
the spread of agriculture involved the dispersal of barley well beyond the native range of the wild species into the variety of environments found in europe.
adaptation to these new conditions is reflected in a north-south clinal distribution of landraces with the daylength responsive and nonresponsive genotypes of the photoperiod gene ppd-h <dig>  nonresponsive forms more common in the cooler northern latitudes  <cit> .
with wild barley, there is a strong correlation between population structure and temperature and precipitation  <cit> .
it might therefore be anticipated that similar climatic correlations may be discernable in the population structure of cultivated barley.
analysis of a series of climate variables supported these expectations .
between-population variance was significantly higher than within-population variance for both spring and winter barleys.
this trend was apparent at all months of the year, but for spring barleys was strongest during the growing season.
for winter barleys the seasonal trend was less clear.
the results indicate that the accessions in each population are adapted, at least to some extent, to their environment, but do not reveal whether this adaptation was a factor in the origin of individual populations, or merely reflects the more recent evolution of landraces to the environments in which they are being grown.
origins of the populations
the relationships inferred from the groupings revealed by neighbour-joining , along with the phenotypic and geographic data, enable possible origins for the populations to be deduced.
populations 1- <dig> are closely related, forming a distinct group in the neighbour-joining tree, and have identical phenotypes, virtually all of their members being two-rowed and hulled with spring growth habit and daylength nonresponsiveness .
we have previously shown that the nonresponsive phenotype of european barleys originated in the eastern fertile crescent and that the first nonresponsive plants probably entered europe after the initial spread of agriculture  <cit> .
this population of nonresponsive plants would almost certainly have had a distinct genetic makeup compared with the barley already present in europe, which originated in the western fertile crescent.
populations 1- <dig> are almost exclusively nonresponsive  and could be the descendents of this original population of nonresponsive plants.
these three populations possess the wild phenotypes for ear row number and caryopsis structure, but have acquired a spring growth habit, whereas their wild progenitors would have been winter types.
the presence of some members of population  <dig> in the same region of the neighbour-joining tree as populations 1- <dig> is indicative of past cross-hybridisation between these populations, which we discuss below.
population  <dig> is also made up entirely of daylength nonresponsive accessions.
this population is located some distance from populations 1- <dig> in the tree topology.
population  <dig> has a narrow geographical distribution in switzerland and the carpathian mountains  and is the only population in which the majority of accessions have naked rather than hulled grains.
the apparent lack of a close relationship between population  <dig> and populations 1- <dig> might indicate that the former is not directly descended from the latter.
instead, population  <dig> could have become homogeneous for daylength nonresponsiveness via a founder effect operating on a population that contained a mixture of responsive and nonresponsive types.
the tree topology suggests that this progenitor of population  <dig> might have been related to the modern populations  <dig> and/or  <dig>
population  <dig> forms a separate cluster in the neighbour-joining tree, but has a mixture of phenotypes, including two- and six-row barleys, hulled and naked forms, spring and winter habits and both daylight responsive and nonresponsive.
there is little uniformity to the combination of phenotypes possessed by individual accessions, and the two deeply rooted groups within the population  <dig> cluster are equally mixed.
these features, along with the broad geographical distribution, suggests that this population has not been subject to selection.
with a crop such as barley, one way in which a distinct genetic population might arise is by geographical partitioning during or soon after the initial spread of agriculture.
populations might be expected to arise in this way if the process of spread involves two or more trajectories that isolate different parts of the crop so that cross-hybridization between the nascent populations is restricted.
the original spread of agriculture into europe is thought to have followed at least two trajectories, one along a northern route through the balkans, hungary and danube and rhine valleys, and the other through the mediterranean basin to italy and iberia  <cit> .
the lack of evidence for human or environmental selection might therefore indicate that population  <dig> is a relict of a population that originated from the geographical partitioning that occurred during this initial period of spread along the northern trajectory.
another candidate as a relict is population  <dig>  as the core area of distribution of this population lies within those regions of mediterranean europe where crops are thought to have spread via the southern trajectory.
if the spread of cultivation along this trajectory resulted in evolution of a distinct population of barley then that population, at least initially, would have had a geographical distribution very similar to that displayed today by population  <dig>
population  <dig> is predominantly six-rowed, hulled and daylight responsive, with a mixture of winter and spring types.
population  <dig> has similar phenotypic features to population  <dig> but contains a greater proportion of landraces with the winter growth habit and is exclusively daylight responsive, whereas population  <dig> includes some nonresponsive types.
the possibility that the two populations might have an evolutionary relationship is supported by the structure analysis, the two populations being grouped as one at k =  <dig>  not splitting into separate populations until k =  <dig> , but the topology of the neighbour-joining tree gives less evidence for a close relationship.
the final two populations,  <dig> and  <dig>  are grouped as one by structure at k ≤  <dig>  and their accessions are located together in the neighbour-joining tree, albeit in three separate parts of the topology.
their geographical distributions are largely non-overlapping, with population  <dig> centering on the northern balkans, hungary and romania, and population  <dig> in northern europe, scandinavia and the baltic states.
this suggests that originally they formed a single population spanning most of the eastern half of europe, subsequently splitting into two, possibly by geographical partitioning.
they are largely six-row, entirely hulled and predominantly spring growth habit, but they contain a mixture of daylength responsive and nonresponsive forms.
the latter are located almost exclusively within the lower part of the tree shown in figure  <dig>  alongside populations 1- <dig>  the implication is that cross-hybridization resulted in transfer of the daylength nonresponsive phenotype from populations 1- <dig> to some members of populations  <dig> and  <dig>  daylight nonresponsiveness and spring growth habit can be advantageous for the successful growth of barley in the more northerly regions of europe.
acquisition of daylength nonresponsiveness by a group of early barley landraces that had already evolved a spring growth habit might therefore have been one of the evolutionary adaptations that enabled cultivation of those plants to be extended further north into the regions now occupied by populations  <dig> and  <dig>  it might therefore be hypothesized that these populations represent a derived form of barley that evolved during the spread of agriculture into central and northern europe.
we explore these and other archaeological interpretations of the population structure in more detail elsewhere .
CONCLUSIONS
we have shown that barley landraces can be divided into populations based on their microsatellite genotypes, and that these populations have different core distributions in europe.
the population structure is partly associated with phenotype, suggesting that human selection and/or environmental adaptation may have played a role in the origin and/or maintenance of one or more populations, but there is also evidence that at least some of the population structure originated during the initial spread of barley cultivation into europe, or reflects the later introduction of daylength-nonresponsive varieties.
the dissection of population structure, combined with examination of their phenotypic attributes and environmental adaptations, enables a rational approach to the identification of landraces that might be used as sources of valuable germplasm for modern breeding programmes.

 cultivated barley, the domesticated form of hordeum spontaneum c. koch, was one of the founder crops of agriculture in western asia, first appearing in the archaeological record in the 8th and 7th millennia bc. barley was also one of the principal crops that accompanied the spread of agriculture into europe during the 6th and 5th millennia bc. today it is grown throughout the continent, mainly for animal feed and malt for brewing. until the introduction of modern cultivars in the mid-20th century, european barley comprised a large number of landraces, each of these a'dynamic population or populations of a cultivated plant that has historical origin, distinct identity and lacks formal crop improvement, as well as often being genetically diverse, locally adapted and associated with traditional farming systems'. many landraces died out during the last century, but seeds representing several thousand types from all parts of europe are available from germplasm collections. increasingly, these landrace collections are being looked on as important sources of germplasm with which to enrich the genepool of modern barley cultivars. exploitation of these landraces in modern crop breeding requires understanding not only of their phenotypic attributes and environmental adaptations, but also their evolutionary relationships. this information would enable landraces with similar valuable traits but different genetic backgrounds to be identified, providing alternative strategies for the introduction of the traits into modern germplasm. an understanding of the relationships between different landraces might also provide information relating to the spread and subsequent development of barley cultivation in europe. this possibility has been recognized for some time, but has not been extensively explored. recent studies have suggested, however, that european barley landraces are differentiated into genetically defined populations. for example, dna sequencing of the photoperiod response gene ppd-h in european barley landraces has revealed three distinct groups of haplotypes, two with the daylight responsive phenotype found in most wild barleys, and one with the derived nonresponsive phenotype that confers adaptation to the long annual growth season of northern europe. multilocus studies have also been carried out, avoiding the risks of inferring population history from single gene data. population structure, linked at least in part with ecogeography and/or agronomic factors, has been revealed by analysis of microsatellites in barley landraces from iberia, the levant and the himalayas. a relationship between microsatellite genotypes and ecogeography has also been demonstrated at the microscale for wild barley. in this paper we show that european barley landraces can be divided into populations based on their microsatellite genotypes. we assess the extent to which the population structure can be explained by human selection, environmental adaptation, geographical partitioning occurring during the initial spread of barley cultivation into europe, and/or the later introduction of daylength-nonresponsive landraces.  microsatellite genotypes we studied accessions of cultivated barley and typed microsatellites. each of these microsatellites displayed variability among the landraces that were tested, the number of alleles observed per locus ranging between 2- with a mean of. pic values varied between -, in broad agreement with values previously obtained for spanish barley landraces. there were no significant differences between the overall diversities of two- and six-rowed barley landraces, though several individual microsatellites did show significant differences when these two sets of landraces were compared. missing data varied from - %. missing data, which usually arise when a landrace has a sequence polymorphism within the annealing site for one of the two primers used to amplify the microsatellite, will lead to underestimates of overall diversity, but are not an issue for subsequent data analysis because structure is able to compensate for gaps in the overall dataset. identification of populations to assess whether the barley landraces could be divided into populations, structure was used to calculate the probability distribution pr for values of k from 2- where x is the genotypes of the sampled landraces and k is the number of populations. multiple runs were carried out and population assignments tested for reproducibility by mutual correlation of the q-matrix outputs. q-matrices were accepted as reproducible if the allocations to each population could be identified as similar in a pair of q-matrices and were highly correlated. for k = and reproducibility was achieved after duplicate runs, while for k = and two reproducible results were obtained after a third run. the population assignments at k = were symmetrical, each accession being given an equal allocation to every population, and reproducible results could not be obtained for values of k from 17- these results indicate that the dataset displays population structure and that the most likely value of k is ≤ three methods were used to identify the most likely value of k more accurately. the first method was based on the prediction that the upper limit of k is indicated by the smallest value that captures the major structure in the data. a plot of lnpr against k was therefore made for k = to. this plot shows lnpr increasing until k = 12- after which the line reaches a plateau. this observation suggests that the most likely value of k is ≤ in the second method the rate of change of lnpr for successive values of k was plotted against k. the position of a plateau in a plot of this type indicates the most likely value of k. examination of figure 1b therefore suggests that k = 9- the third method used to identify the most likely number of populations assumes that genetic structure should relate to phenotype. data on spring or winter growth habit, two-row or six-row ear morphology and caryopsis structure were converted into binary characters and logistic regressions of these phenotypes against the q-matrices for k = 2- carried out with the r statistics package. in this type of analysis, the lowest values for akaike's information criterion identify the range of k where population structure is best correlated with phenotype. for growth habit, the lowest aic values were found for k = 9- for row number at k = and for caryopsis structure at k = 8-. taken together, the three analyses summarized in figures suggest that the most likely value of k for the barley landraces is between and while we do not attempt to assign a'correct' value for k, we chose a value of as a starting point for examination of population structure. a graphical representation of the population structure at k = revealed that each population was partially admixed with other populations, the overall degree of admixture being similar for each one. half the landraces had a proportional membership of ≥ for their primary population, and only had a primary proportional membership of <. the relationships between these populations were studied in two ways. first, the microsatellite data were used to construct a neighbour-joining tree and those accessions with a proportional membership of ≥ in their primary population marked on the tree. accessions belonging to a single population clustered together, except for population and to a lesser extent which appeared in two distinct regions of the tree. populations 1- and some members of and grouped close to one another, as did population alongside the remainder of and populations and each formed a separate group, although population was split into two parts each with its own deep root towards the base of the tree. to explore further the population relationships, the groups resulting from structure analyses with preset k values of 4- were investigated. again, only those landraces that displayed a proportional membership of ≥ in their primary populations were considered. comparing the allocations for each of these landraces as k increased revealed a hierarchical structure to the landrace populations. landraces that were grouped together at k = were still grouped together at k = 7- a subset of landraces from two of the populations at k = a and b, were placed in a shared population at k = 5- then at higher k values re-sorted into distinct lineages that corresponded to their original allocation between populations a and b. these landraces made up populations 1- at k = the remaining two populations at k = c and d, both gave rise to lineages that remained distinct as k increased to phenotypic features of the populations for those accessions for which data were available, ear row number, caryopsis structure, growth habit and daylength responsiveness were compared with the population structure. two-rowed barley predominated in populations 1- and six-row barley predominated in populations 5-, these differences being significant. naked barleys were predominant in population, and formed a sizeable minority of population, but were virtually absent from the other seven populations. winter growth habit was predominant only in population. spring growth habit accounted for > 90% of the accessions in populations 1- and according to passport data, and in population also according to vrn genotype. the remaining two populations comprised a mixture of spring and winter barleys. the discrepancies within populations and between the passport data for growth habit and the habit predicted from the vrn genotype are due to many of these accessions coming from the mediterranean region, where barley with a facultative growth habit is favoured. these barleys are autumn sown, able to survive mild winter conditions, and flower in the late spring. they are therefore recorded as having a winter growth habit although they display the spring genotype when the vrn genes are typed. of a top row, number of landraces with each phenotype, with the total that could be phenotyped in brackets; bottom row, percentage of landraces with each phenotype for accessions, daylength responsiveness was deduced from the ppd-hi genotype. all of the accessions of populations and that were typed, as well as 95% of those in population were daylength nonresponsive, and all those typed from population and 91% of those from population were responsive. the other four populations were mixed, ranging from 20-91% responsive. where possible, those accessions with the daylength responsive genotype were placed in group b or c. group b predominated only in population and ranged from 0-13% in the other populations. geographical distributions of the populations all accessions were included in the geographical analyses. the mean centre and standard distance for each population were calculated, along with the mean pairwise geographical distance within each population to assess the degree of clustering. all of the landraces were included. landraces of each population were placed on a map of europe and standard deviation ellipses plotted. visual inspection of the maps indicated that the populations had non-identical distribution patterns. three populations were distributed over a relatively broad area of central and western europe. population had a similar but more western distribution, largely due to a preponderance of british landraces. populations and centred on switzerland and the carpathian mountains, with population tightly clustered in a small area of west-central switzerland. population was mainly located in the balkan region of southeast europe, population predominated in north central europe including scandinavia and the baltic states, and population clustered in the mediterranean region. to assess the apparent differences between these distribution patterns, the great circle distances between pairs of populations were tested for significance using student's t-test against a null hypothesis that the mean centres were km apart. this analysis suggested that of the possible population pairs had a highly significant difference in their mean centres, and a further four had a significant difference. the only non-significant differences were between paired populations and and and and and the relationship between the population structure and environment was investigated by comparing the between-population variance and within-population variance for a range of climatic factors. this analysis was carried out separately for the spring and winter accessions as these different growth habits might be expected to result in different patterns of environmental adaptation. for randomized populations, the f ratio for between-population variance compared with within-population variance was close to one for each of the climate variables at all months of the year. when all the climate variables were combined, the between-population variance for spring barleys was significantly higher than within-population variance for each month, with the highest values occurring during the growing season for these accessions. between-population variance was also significantly greater than within-population variance for winter barleys when all climate variables were combined but with no obvious seasonal trend. discussion identification of populations structure was used to determine whether european barley landraces can be divided into populations based on their microsatellite genotypes. in this context, a population is defined as a group of individuals that share a characteristic set of allele frequencies at the loci that are studied. structure places individuals in populations in such a way as to minimize within-population deviations from hardy-weinberg equilibrium and linkage equilibrium. it therefore assumes that the individuals are fully outcrossing, and modelling studies have suggested that with partially inbreeding taxa the results of structure analysis can be spurious. although cultivated barley has an outcrossing rate of less than 2%, the outcomes of previous structure analyses of genetic data from barley have been supported by other analyses of the same data and are in agreement with the conclusions of earlier and later work. in these projects it therefore appears that structure identified authentic populations despite barley's low outcrossing rate. an explanation might be provided by other simulations which have shown that, over hundreds of generations, the pattern of multilocus marker inheritance in a population of plants displaying 2% outbreeding is indistinguishable from that displayed by a panmictic population. the results of our structure analysis were reproducible for values of k up to and, using standard methods, we concluded that the most likely number of populations was between 8- the population structure was hierarchical, and from k = 7- the only effect of each incremental increase in k was to subdivide an existing population with the memberships of the other populations remaining unchanged. this consistent pattern of population assignment indicates that the results of the structure analysis were not spurious. to further assess the validity of the analysis, we constructed a neighbour-joining tree for all accessions and marked the positions of those accessions with a proportional membership of ≥ in their primary population at k =. accessions belonging to a single population clustered together, with the exception of populations and whose members are distributed in two parts of the tree. the tree topology therefore provides independent support for the structure results, and also suggests that it is reasonable to use k = as the basis for interpretation of the population structure. association between phenotype and population structure for a domesticated plant such as cultivated barley, one possible way in which population structure could arise is as a result of selection for particular phenotypic traits. the phenotypes of greatest agronomic importance in modern farming are ear row number, caryopsis structure, growth habit and flowering time. wild barley has a two-rowed ear, each spikelet having a fertile central floret flanked by two infertile laterals which, when combined with a long awn, takes on an arrowlike form that is an effective aid to seed dispersal and burial. many cultivated barleys retain this ancestral head structure but in the derived six-rowed form the two lateral florets are fertile. six-rowed barley is more often used as an animal or human feed, whereas two-rowed barley is favoured for malting and brewing. wild and most cultivated barleys have hulled grains where the outer lemma and inner palea adhere to the pericarp epidermis at maturity. this form is favoured by brewers because the hull debris aids wort filtration, whereas the free-threshing'naked' varieties are preferred when barley is grown for direct human consumption. wild barley has a winter growth habit, meaning that it requires vernalization - exposure to a prolonged period of cold - in order to promote subsequent flowering. the majority of european landraces lack this requirement and have a spring growth habit, where plants avoid periods of cold weather by completing their growth cycle during a single season, rather than over wintering as plants in a vegetative state. finally, most wild barleys display a photoperiod response that triggers flowering early in the season, before the conditions become too dry for further vegetative growth. many cultivated barleys, especially landraces from northern europe, are daylength nonresponsive, and so continue vegetative growth until flowering later in the summer, allowing them to take advantage of the longer growing season in northern europe. we compared phenotypic data for ear row number, caryopsis structure, growth habit and flowering time with the population memberships. populations 1- display a similar set of phenotypic features, most of these accessions being two-rowed, hulled, spring habit, and daylength nonresponsive. populations and which contained a high proportion of accessions with the facultative growth habit, also show some similarities when other phenotypes are considered, being largely six-row, hulled and daylength responsive. in population however, the majority of the responsive accessions were members of group c, whereas in population the majority were group b, these two groups having distinct evolutionary histories. no other similarities between the range of phenotypes displayed by different populations were apparent. these comparisons indicate that there is a degree of association between phenotype and population structure, suggesting that selection may have played a role in the origin and/or maintenance of one or more populations. association between climatic factors and population structure the spread of agriculture involved the dispersal of barley well beyond the native range of the wild species into the variety of environments found in europe. adaptation to these new conditions is reflected in a north-south clinal distribution of landraces with the daylength responsive and nonresponsive genotypes of the photoperiod gene ppd-h nonresponsive forms more common in the cooler northern latitudes. with wild barley, there is a strong correlation between population structure and temperature and precipitation. it might therefore be anticipated that similar climatic correlations may be discernable in the population structure of cultivated barley. analysis of a series of climate variables supported these expectations. between-population variance was significantly higher than within-population variance for both spring and winter barleys. this trend was apparent at all months of the year, but for spring barleys was strongest during the growing season. for winter barleys the seasonal trend was less clear. the results indicate that the accessions in each population are adapted, at least to some extent, to their environment, but do not reveal whether this adaptation was a factor in the origin of individual populations, or merely reflects the more recent evolution of landraces to the environments in which they are being grown. origins of the populations the relationships inferred from the groupings revealed by neighbour-joining, along with the phenotypic and geographic data, enable possible origins for the populations to be deduced. populations 1- are closely related, forming a distinct group in the neighbour-joining tree, and have identical phenotypes, virtually all of their members being two-rowed and hulled with spring growth habit and daylength nonresponsiveness. we have previously shown that the nonresponsive phenotype of european barleys originated in the eastern fertile crescent and that the first nonresponsive plants probably entered europe after the initial spread of agriculture. this population of nonresponsive plants would almost certainly have had a distinct genetic makeup compared with the barley already present in europe, which originated in the western fertile crescent. populations 1- are almost exclusively nonresponsive and could be the descendents of this original population of nonresponsive plants. these three populations possess the wild phenotypes for ear row number and caryopsis structure, but have acquired a spring growth habit, whereas their wild progenitors would have been winter types. the presence of some members of population in the same region of the neighbour-joining tree as populations 1- is indicative of past cross-hybridisation between these populations, which we discuss below. population is also made up entirely of daylength nonresponsive accessions. this population is located some distance from populations 1- in the tree topology. population has a narrow geographical distribution in switzerland and the carpathian mountains and is the only population in which the majority of accessions have naked rather than hulled grains. the apparent lack of a close relationship between population and populations 1- might indicate that the former is not directly descended from the latter. instead, population could have become homogeneous for daylength nonresponsiveness via a founder effect operating on a population that contained a mixture of responsive and nonresponsive types. the tree topology suggests that this progenitor of population might have been related to the modern populations and/or population forms a separate cluster in the neighbour-joining tree, but has a mixture of phenotypes, including two- and six-row barleys, hulled and naked forms, spring and winter habits and both daylight responsive and nonresponsive. there is little uniformity to the combination of phenotypes possessed by individual accessions, and the two deeply rooted groups within the population cluster are equally mixed. these features, along with the broad geographical distribution, suggests that this population has not been subject to selection. with a crop such as barley, one way in which a distinct genetic population might arise is by geographical partitioning during or soon after the initial spread of agriculture. populations might be expected to arise in this way if the process of spread involves two or more trajectories that isolate different parts of the crop so that cross-hybridization between the nascent populations is restricted. the original spread of agriculture into europe is thought to have followed at least two trajectories, one along a northern route through the balkans, hungary and danube and rhine valleys, and the other through the mediterranean basin to italy and iberia. the lack of evidence for human or environmental selection might therefore indicate that population is a relict of a population that originated from the geographical partitioning that occurred during this initial period of spread along the northern trajectory. another candidate as a relict is population as the core area of distribution of this population lies within those regions of mediterranean europe where crops are thought to have spread via the southern trajectory. if the spread of cultivation along this trajectory resulted in evolution of a distinct population of barley then that population, at least initially, would have had a geographical distribution very similar to that displayed today by population population is predominantly six-rowed, hulled and daylight responsive, with a mixture of winter and spring types. population has similar phenotypic features to population but contains a greater proportion of landraces with the winter growth habit and is exclusively daylight responsive, whereas population includes some nonresponsive types. the possibility that the two populations might have an evolutionary relationship is supported by the structure analysis, the two populations being grouped as one at k = not splitting into separate populations until k =, but the topology of the neighbour-joining tree gives less evidence for a close relationship. the final two populations, and are grouped as one by structure at k ≤ and their accessions are located together in the neighbour-joining tree, albeit in three separate parts of the topology. their geographical distributions are largely non-overlapping, with population centering on the northern balkans, hungary and romania, and population in northern europe, scandinavia and the baltic states. this suggests that originally they formed a single population spanning most of the eastern half of europe, subsequently splitting into two, possibly by geographical partitioning. they are largely six-row, entirely hulled and predominantly spring growth habit, but they contain a mixture of daylength responsive and nonresponsive forms. the latter are located almost exclusively within the lower part of the tree shown in figure alongside populations 1- the implication is that cross-hybridization resulted in transfer of the daylength nonresponsive phenotype from populations 1- to some members of populations and daylight nonresponsiveness and spring growth habit can be advantageous for the successful growth of barley in the more northerly regions of europe. acquisition of daylength nonresponsiveness by a group of early barley landraces that had already evolved a spring growth habit might therefore have been one of the evolutionary adaptations that enabled cultivation of those plants to be extended further north into the regions now occupied by populations and it might therefore be hypothesized that these populations represent a derived form of barley that evolved during the spread of agriculture into central and northern europe. we explore these and other archaeological interpretations of the population structure in more detail elsewhere.  we have shown that barley landraces can be divided into populations based on their microsatellite genotypes, and that these populations have different core distributions in europe. the population structure is partly associated with phenotype, suggesting that human selection and/or environmental adaptation may have played a role in the origin and/or maintenance of one or more populations, but there is also evidence that at least some of the population structure originated during the initial spread of barley cultivation into europe, or reflects the later introduction of daylength-nonresponsive varieties. the dissection of population structure, combined with examination of their phenotypic attributes and environmental adaptations, enables a rational approach to the identification of landraces that might be used as sources of valuable germplasm for modern breeding programmes. 
second, evolutionary information would enable landraces with similar traits but different genetic backgrounds to be identified, providing alternative strategies for the introduction of these traits into modern germplasm. first, the evolutionary relationships between different landraces might provide information on the spread and subsequent development of barley cultivation, including the adaptation of the crop to new environments and its response to human selection. there was also evidence that at least some of the population structure derived from geographical partitioning set up during the initial spread of barley cultivation into europe, or reflected the later introduction of novel varieties. these three populations probably originated in the eastern fertile crescent and entered europe after the initial spread of agriculture. the evolutionary relationships between barley landraces were inferred from the genotypes for microsatellites. human selection and/or environmental adaptation may therefore have played a role in the origin and/or maintenance of one or more of the barley landrace populations. the discovery of population structure, combined with knowledge of associated phenotypes and environmental adaptations, enables a rational approach to identification of landraces that might be used as sources of germplasm for breeding programs. comparisons with ear row number, caryopsis structure, seasonal growth habit and flowering time revealed a degree of association between population structure and phenotype, and analysis of climate variables indicated that the landraces are adapted, at least to some extent, to their environment. the population structure also enables hypotheses concerning the prehistoric spread and development of agriculture to be addressed. in particular, three closely-related populations were made up almost entirely of plants with the daylength nonresponsive version of the photoperiod response gene ppd-h conferring adaptation to the long annual growth season of northern europe.
 understanding the evolution of cultivated barley is important for two reasons. first, the evolutionary relationships between different landraces might provide information on the spread and subsequent development of barley cultivation, including the adaptation of the crop to new environments and its response to human selection. second, evolutionary information would enable landraces with similar traits but different genetic backgrounds to be identified, providing alternative strategies for the introduction of these traits into modern germplasm.  the evolutionary relationships between barley landraces were inferred from the genotypes for microsatellites. the landraces could be divided into nine populations, each with a different geographical distribution. comparisons with ear row number, caryopsis structure, seasonal growth habit and flowering time revealed a degree of association between population structure and phenotype, and analysis of climate variables indicated that the landraces are adapted, at least to some extent, to their environment. human selection and/or environmental adaptation may therefore have played a role in the origin and/or maintenance of one or more of the barley landrace populations. there was also evidence that at least some of the population structure derived from geographical partitioning set up during the initial spread of barley cultivation into europe, or reflected the later introduction of novel varieties. in particular, three closely-related populations were made up almost entirely of plants with the daylength nonresponsive version of the photoperiod response gene ppd-h conferring adaptation to the long annual growth season of northern europe. these three populations probably originated in the eastern fertile crescent and entered europe after the initial spread of agriculture.  the discovery of population structure, combined with knowledge of associated phenotypes and environmental adaptations, enables a rational approach to identification of landraces that might be used as sources of germplasm for breeding programs. the population structure also enables hypotheses concerning the prehistoric spread and development of agriculture to be addressed. 