Spatial patterns of human gene frequencies in Europe

The aims of this study of spatial patterns of human gene frequencies in Europe are twofold. One is to present new methodology developed for the analysis of such data. The other is to report on the diversity of spatial patterns observed in Europe and their interpretation as evidence of population processes. Spatial variation in 59 allele and haplotype frequencies (26 genetic systems) for polymorphisms in blood antigens, enzymes, and proteins is analyzed for an aggregate of 3,384 localities, using homogeneity tests, one-dimensional and directional spatial correlograms, and SYMAP interpolated surfaces. The data matrices are reduced to reveal the principal patterns by clustering techniques. The findings of this study can be summarized as follows: 1) There is significant heterogeneity in allele frequencies among the localities for all but one genetic system. 2) There are significant spatial patterns for most allele frequencies. 3) There is a substantial minority of clinal patterns in these populations. Clinal trends are found more frequently in HLA alleles than for other variables. North-south and northwest-southwest gradients predominate. 4) There is a strong decline in overall genetic similarity with geographic distance for most variables. 5) There are few, if any, appreciable correlations in pairs of allele frequencies over the continent, and there is little interesting correlation structure in the resulting correlation matrix. 6) Few spatial correlograms are markedly similar to each other, yet they form well-defined clusters. Spatial variation patterns, therefore, differ among allele frequencies. Patterns of human gene frequencies in modern Europe are diverse and complex. No single model suffices for interpretation of the observed genetic structure. Some clinal patterns reported here support the Neolithic demic-expansion hypothesis, others suggest latitudinal selection. Most of the clinal patterns are in HLA alleles, but there is also evidence from ABO for east-west migration diffusion. The majority of patterns are patchy, consistent with hypotheses of isolation by distance or of settlement of genetically differing, subsequently expanding ethnic groups. While undoubtedly there has been an ongoing stochastic process of differentiation consistent with the isolation-by-distance model, this has not obscured the directional patterns caused by migration (demic diffusion), and has perhaps only reinforced the contribution from settlement of ethnic units to patterns of genetic variation. However, the impact of the latter is most difficult to discern and requires further methodological developments.     

Diversity of some gene frequencies in European and Asian populations. VI. Geographic patterns of PGM and ACP

Genetic variation at the phosphoglucomutase-1 (PGM) and acid phosphatase (ACP) loci, in 160 and 120 European and Asian populations, respectively, is described by spatial autocorrelation statistics and directional variograms. Short-distance patterns of gene frequencies correspond to those predicted by models of isolation by distance, but long-range differentiation of populations is observed as well. A possible role of climatic selection in maintaining the PGM polymorphism is supported by the north-south orientation of the gradient for that locus, but not by biochemical evidence. By and large, the observed patterns of gene frequencies seem to reflect a combination of demographic processes, subdivision and isolation of local populations among them.     

HOW SEDENTARY ARE GREYWING FRANCOLINS (FRANCOLINUS AFRICANUS)?

Phasianids are considered to be sedentary birds with limited dispersal so that populations may be expected to show genetic isolation by distance. To test this, we examined genetic variability in 618 greywing francolins (Francolinus africanus) at 24 localities over a 1,500 km² area. We subdivided the samples to measure genetic population structure among localities separated by 6–60 km, and among coveys separated by 0.1–6 km. Thirteen of 30 (43%) allozyme loci were polymorphic, and heterozygosity ranged from 5.3 to 8.5% over 24 localities and averaged 7.0%, a value much larger than that found for other phasianids. Significant allele-frequency heterogeneity was detected among localities and among coveys at several localities for several loci. Mantel's test, however, showed that there was no correlation between geographical distance and the allele-frequency difference between localities for all but one allele. Although spatial autocorrelation was detected with Moran's I and Geary's c for two alleles, the geographical patterns of I in correlograms of 18 independent alleles showed a “crazy-quilt” pattern of allele-frequency patches. This shows that the isolation-by-distance model of subpopulation structure is inappropriate for these birds. Individuals, therefore, appear to disperse far beyond neighboring populations. “Private-allele” and F_ST estimates of migration under the island model were 8–9 individuals between localities of each generation. Allele-frequency heterogeneity, large amounts of gene flow, and the general lack of spatial autocorrelation imply that the small, socially-structured populations of greywing are subject to high rates of turnover, founder effects, and random drift.     

Patterns of gene flow inferred from genetic distances in the Mediterranean region.

The analysis of population structure may lead to inferences about demographic phenomena. In particular, regions of sharp genetic differentiation suggest the existence of factors that impaired gene flow and increased the evolutionary role of genetic drift. Here, we present an analysis of a data set of 10 allele frequencies in 39 populations of the Mediterranean region. As a preliminary step, we describe spatial patterns of allele frequencies using spatial autocorrelation analysis. We then construct a network connecting localities and estimate genetic distances along the edges of the network. By applying specific algorithms, we locate on the map the areas of sharpest genetic differentiation, or genetic boundaries. The main boundaries separate the northern and the southern coasts, especially in their western portions; in addition, several localities appear genetically isolated. The comparatively high genetic differentiation across the western Mediterranean, where the sea distances between localities are shorter, strongly suggests that the sea distance by itself can hardly be regarded as a major isolating factor among these populations. On the contrary, the decrease in genetic resemblance between populations of the 2 coasts as one proceeds westward may reflect an increased genetic exchange in the eastern Mediterranean basin or independent human dispersal along the 2 coasts or both.     

Evaluating the utility of microsatellites for investigations of autopolyploid taxa

Autopolyploid taxa present numerous challenges for population genetic analyses due to difficulties determining allele dosage. Dosage ambiguity hinders accurate assessment of allele frequencies, multilocus genotypes (MLGTs), as well as levels and patterns of clonality. The pervasiveness of polyploidy in the evolutionary history of plant taxa makes this a recurring problem. Whereas diploidization of loci may occur over time, duplication of at least some loci is still frequently evident. Fortunately, with high-quality allozyme gels, it is possible to accurately infer allele dosage and, thus, determine exact MLGTs. However, accurately assessing dosage of microsatellite peaks is nearly impossible when studying wild populations with a large number of alleles per locus. Even if precise knowledge of genotypes is not required, for comparable numbers of alleles per locus and loci, the number of "phenotypes" is always lower with microsatellites than allozymes due to the inability to assess allele dosage. Microsatellite loci typically have more alleles per locus relative to allozymes although fewer loci are generally employed. Here, we present a mathematical model for comparing the relative utility of simple sequence repeat (SSR) versus allozyme markers to discriminate MLGTs. For example, the average plant allozyme study (2.6 alleles per locus, 10 polymorphic loci) has better discriminating power than SSR markers with 10 alleles at each of 3 loci, 9 alleles at 4 loci, 6 alleles at 5 loci, 5 alleles at 6 loci, and 4 alleles at 8 loci, demonstrating the value of assessing the relative discriminating power of these markers.     

Identification, Inheritance, and Variation of Microsatellite Markers in the Black Scallop Mimachlamys varia

Five polymorphic microsatellite loci were identified in the black scallop Mimachlamys varia after construction of a genomic library enriched for (GT)n. To examine the transmission pattern of microsatellite alleles, several families were created and genotypes scored for three loci. The expected Mendelian ratios were found in 12 of 14 segregations examined. Unexpected segregations may be explained by a genotyping error (allelic dropout), given that when a specific allele was treated as dominant, the phenotypic ratios conformed to Mendelian expectations. The five loci were also examined in two samples from the Spanish coast. The two localities displayed similar mean values for the number of alleles per locus (7.2-8.4), allelic richness (7.2-7.9), and observed (0.389-0.484) and expected heterozygosity (0.545-0.618). Significant Hardy-Weinberg deviations were observed at three loci, with heterozygote deficiency occurring in all cases. Global multilocus θ value and allele frequencies at one locus revealed significant differentiation between the two localities.     

Impact of negative frequency-dependent selection on mating pattern and genetic structure: a comparative analysis of the S-locus and nuclear SSR loci in Prunus lannesiana var. speciosa

Mating processes of local demes and spatial genetic structure of island populations at the self-incompatibility (S-) locus under negative frequency-dependent selection (NFDS) were evaluated in Prunus lannesiana var. speciosa in comparison with nuclear simple sequence repeat (SSR) loci that seemed to be evolutionarily neutral. Our observations of local mating patterns indicated that male-female pair fecundity was influenced by not only self-incompatibility, but also various factors, such as kinship, pollen production and flowering synchrony. In spite of the mating bias caused by these factors, the NFDS effect on changes in allele frequencies from potential mates to mating pollen was detected at the S-locus but not at the SSR loci, although the changes from adult to juvenile cohorts were not apparent at any loci. Genetic differentiation and isolation-by-distance over various spatial scales were smaller at the S-locus than at the SSR loci, as expected under the NFDS. Allele-sharing distributions among the populations also had a unimodal pattern at the S-locus, indicating the NFDS effect except for alleles unique to individual populations probably due to isolation among islands, although this pattern was not exhibited by the SSR loci. Our results suggest that the NFDS at the S-locus has an impact on both the mating patterns and the genetic structure in the P. lannesiana populations studied.     

BIOCHEMICAL GENETICS OF MOSQUITOFISH. IV. CHANGES OF ALLELE FREQUENCIES THROUGH TIME AND SPACE

Gene frequency data from samples of Gambusia affinis populations at 76 localities across the Savannah River drainage were used to investigate temporal and spatial patterns in population genetic structure. Localities in the Par Pond system on the Savannah River Plant were sampled in 1971, 1977, and 1979. Allelic frequencies in these populations were generally stable through time, although significant temporal changes were observed among samples from Pond C, an impoundment receiving thermal effluent. Significant spatial heterogeneity in allele frequencies was observed on both microgeographic and regional scales. Populations within the Par Pond system were spatially subdivided at four of the five loci surveyed (mean F_ST = 0.051). Subdivision was even more pronounced when samples from across the Savannah River drainage were compared (mean F_ST = 0.196). A hierarchial analysis of gene diversity (G_ST) demonstrated that most of the genic diversity across the drainage exists as within-subdivision diversity. Even when populations from such contrasting habitats as rivers, creeks, ponds, and reservoirs are compared, an average of only 13% of the total gene diversity was attributed to between-group diversity. Greatest between-group gene diversity was observed when reservoirs were compared with one another. This general pattern of low between-habitat diversity suggests that differential selection pressures are not playing a major role in producing the observed levels of subdivision. In the Par Pond system, neither single locus nor multilocus genetic distances were significantly associated with geographic distance or with its reciprocal. For samples from over the Savannah River drainage, significant correlations between genetic and geographic distance were observed only for the Gpi-2 and Pgm-2 loci. Thus, there was a general lack of concordance between genetic and geographic distances. Spatial autocorrelation demonstrated patterns consistent with Wright's isolation by distance model. Significant positive correlations in allelic frequencies among neighboring populations were observed for five of six alleles; allelic frequencies in more distantly separated populations were typically not correlated.     

意大利蜜蜂产浆性状10个微卫星位点的遗传分析

我国培育成功的世界上蜂王浆产量最高的蜂种（Ea）是从20世纪30年代引进我国的意大利蜜蜂（Eb）中选育的。采用10个微卫星位点对蜂王浆高产蜜蜂、原种意大利蜜蜂（Ee）和本地意大利蜜蜂进行研究,以探明由于人工选择和地理隔离造成的其分子进化上的一些特征。结果,10个微卫星位点在3个蜂种中共扩增到96个等位基因,其中有48个等位基因是不同的,表明10个微卫星位点在3个蜂种中的高度多态性,而且由于人工选择和地理隔离已造成3蜂种在遗传基础上的一些分化。3蜂种的多态信息含量（PIC）分别为0．57、0．50、0．57,杂合度分别为0．60、0．57、0．61,二者均无显著差异。遗传距离的分析结果显示,Ee和Eb（0．14）、Eb和Ea（0．16）之间的遗传距离较近,而Ea和Ee（0.25）之间的遗传距离较远。等位基因频率的分析结果表明,6个位点的7个等位基因的频率（A29的159bp、A24的100bp和104bp、A7的110bp、A43的126bp、A14的221bp和A113的221bp）以Ee、Eb、Ea的顺序递增,Ea的这7个等位基因频率分别显著高于Ee和Eb。同时4个位点的4个等位基因的频率（A24的106bp、A43的140bp;A113的215bp和A14的219bp）以Ea、Eb、Ee的Ⅲ页序下降,Ea在这4个位点的频率分别显著低于Eb和Ee。这些位点的等位基因可能与蜂王浆产量有关。     

VNTR loci reveal differentiation between and structure within populations of the eastern barred bandicoot Perameles gunnii

The Eastern Barred Bandicoot Perameles gunnii has declined in abundance within mainland south-eastern Australia, to a relict wild population of less than 100 individuals in Hamilton, Victoria. It is more common, but is also declining in Tasmania. Genomic DN A variability was compared within and between surviving populations of P. gunnii using variable number of tandem repeat (VNTR) markers in one of two ways. First, average percentage differences (APDs) were determined between profiles for two VNTR probe—endonuclease combinations. Secondly, because one of these combinations revealed two multiallelic VNTR loci, genotypes were assigned and analysed for homogeneity of allele frequencies among subpopulations, for deviation of heterozygosity from Hardy-Weinberg equilibrium within populations and for genetic structuring among individuals from different subpopulations. The results of both the APD and defined locus approaches showed consistent trends within and between populations. Genetic variability was higher among mainland P. gunnii than in Tasmanian populations (higher APDs, number of alleles, and heterozygosity at one locus), despite the known decline and subdivision of the Hamilton population. Eleven per cent of the variability detected in Hamilton was attributed to genetic differentiation between east and west subdivisions of the population. Departure from random mating indicating local inbreeding within collecting localities was evident for one locus in both north and south Tasmania, particularly at one locality. AH alleles at both loci were unique to either Hamilton or Tasmanian P. gunnii. The initial captive colony contains high heterozygosity for these loci. It is concluded that VNTR markers can be of benefit for use in studies of population differentiation and for conservation management.     

Temporal and Microgeographic Variation in Allozyme Frequencies in a Natural Population of DROSOPHILA BUZZATII

Temporal variation in allozyme frequencies at six loci was studied by making monthly collections over 4 yr in one population of the cactophilic species Drosophila buzzatii. Ten sites were defined within the study locality, and for all temporal samples, separate collections were made at each of these sites. Population structure over microgeographic space and changes in population structure over time were analyzed using F-statistic estimators, and multivariate analyses of allele and genotype frequencies with environmental variables were carried out. Allele frequencies showed significant variation over time, although there were no clear cyclical or seasonal patterns. A biplot analysis of allele frequencies over seasons within years and over years showed clear discrimination among years by alleles at four loci. During the 4 yr, three alleles showed directional changes which were associated with directional changes in environmental variables. Significant associations with one or more environmental variables were found for allele frequencies at every locus and for both expected and observed heterozygosities (except those for Est-1 and Est-2). Thus, variation in allele frequencies over time cannot be attributed solely to drift. Significant linkage disequilibria were detected among three loci (Est-2, Hex and Aldox), but there was no evidence for spatial or temporal patterns. The F-statistic analyses showed significant differentiation among months within years for all loci, but the statistic used (coancestry) was heterogeneous among loci. Estimates of F (inbreeding) for all loci were significantly different from zero, with the loci in four groups, Adh-1 (negative), Pgm(small positive), Est-2 and Hex (intermediate) and Est-1 and Aldox (high positive). The correlation of genes within individuals within populations (f) for each locus in each month by site sample differed among loci, as did the (f) for each locus in each month by site sample differed among loci, as did the patterns of change in f over time (seasons). Heterogeneity in the F-statistic estimates indicates that natural selection is directly or indirectly affecting allele and genotype frequencies at some loci. However, the F-statistic analyses showed essentially no microgeographic structure (i.e., among sites), although there was significant heterogeneity in allele frequencies among flies emerging from individual rots. Thus, microspatial heterogeneity probably is most important at the level of individual rots, and coupled with habitat selection, it could be a major factor promoting diversifying selection and the maintenance of polymorphism.(ABSTRACT TRUNCATED AT 400 WORDS)     

Concordance of allozyme and microsatellite differentiation in a marine fish, but evidence of selection at a microsatellite locus

Previous studies have reported higher levels of divergence for microsatellites than for allozymes in several species, suggested to reflect stabilizing selection on the allozymes. We compared the differentiation patterns of 11 allozyme and nine microsatellite loci using 679 spawning Atlantic herring (Clupea harengus) collected in the Baltic and North Seas to test for differential natural selection on these markers. Observed distributions of F statistics for the two types of markers are conspicuously dissimilar, but we show that these differences can largely be explained by sampling phenomena caused by different allele frequency distributions and degrees of variability. The results show consistently low levels of differentiation for both marker types, with the exception of one outlier microsatellite locus with a notably high F(ST). The aberrant pattern at this locus is primarily due to two alleles occurring at markedly high frequencies in the Baltic, suggesting selection at this locus, or a closely linked one. When excluding this locus, the two marker types show similar, weak differentiation patterns with F(ST) values between the Baltic and the North Seas of 0.001 and 0.002 for allozymes and microsatellites, respectively. This small heterogeneity, and weak isolation by distance, is easier to distinguish statistically with microsatellites than with allozymes that have fewer alleles and skewed frequency distributions. The allozymes, however, also detect surprisingly low levels of divergence. Our results support suggestions that previously described differences between marker types are primarily caused by a small number of outlier loci.     

Cultural barriers associated with large gene frequency differences among Italian populations.

Analysis of geographic variation for eight red cell markers in Italy shows significant spatial structure for most alleles. Effective population sizes estimated from FST values at these loci are much smaller than those predicted from data on consanguineous marriage, suggesting the presence of factors (presumably barriers) that have reduced gene flow and enhanced the evolutionary weight of genetic drift. Most regions of sharp gene frequency change correspond to geographic and linguistic barriers. Two allele frequencies are significantly correlated with measures of linguistic differentiation but not with indexes describing broad religious and social attitudes. The similarity between patterns of genetic and linguistic variation in Italy, also observed in a previous study, suggests that in specific areas linguistic diversity has acted as a biological barrier constraining mating, dispersal, or both. There is no evidence for a similar role of other extent cultural barriers.     

Microsatellite Variation in Two Populations of Free-Ranging Yellow Baboons (Papio hamadryas cynocephalus)

We investigated genetic variation at six microsatellite (simple sequence repeat) loci in yellow baboons (Papio hamadryas cynocephalus) at two localities: the Tana River Primate Reserve in eastern Kenya and Mikumi National Park, central Tanzania. The six loci (D1S158, D2S144, D4S243, D5S1466, D16S508, and D17S804) were all originally cloned from and characterized in the human genome. These microsatellites are polymorphic in both baboon populations, with the average heterozygosity across loci equal to 0.731 in the Tana River sample and 0.787 in the Mikumi sample. The genetic differentiation between the two populations is substantial. Kolmogornov–Smirnov tests indicate that five of the six loci are significantly different in allele frequencies in the two populations. The mean F _ST across loci is 0.069, and Shriver's measure of genetic distance, which was developed for microsatellite loci (Shriver et al., 1995), is 0.255. This genetic distance is larger than corresponding distances among human populations residing in different continents. We conclude that (a) the arrays of alleles present at these six microsatellite loci in two geographically separated populations of yellow baboons are quite similar, but (b) the two populations exhibit significant differences in allele frequencies. This study illustrates the potential value of human microsatellite loci for analyses of population genetic structure in baboons and suggests that this approach will be useful in studies of other Old World monkeys.     

Measures of spatial structure in samples of genotypes for multiallelic loci.

Various spatial autocorrelation statistics have been widely used both in theoretical population genetics and to study the spatial distribution of diploid genotypes in many plant and animal populations. However, previous simulation studies have considered only diallelic loci. In this paper, we use a large number of space-time simulations to characterize for the first time the parametric and statistical values of Moran's I-statistics for converted individual genotypes as well as for join-count statistics. A wide range of levels of dispersal and numbers of alleles and allele frequencies are modelled and the results reveal the different general effects of each of these factors on these statistics. We also examine the range of appropriate sampling designs and sizes for which predicted values can be interpolated for specific sampling schemes for any given population genetic field survey. Numbers of alleles and allele frequencies each affect some statistics but not others. The results indicate generally low standard deviations. The results also develop precise and efficient methods of estimating gene dispersal, based on the various autocorrelation measures of standing spatial patterns of genetic variation within populations. The results also extend these methods to loci with multiple alleles, typical of those studied through modern molecular methods.     

The Genetic Structure of a Tribal Population, the Yanomama Indians. Xv. Patterns Inferred by Autocorrelation Analysis

Fifteen allele frequencies have previously been determined for 50 villages of the Yanomama, an Amerindian tribe from southern Venezuela and northern Brazil. These frequencies were subjected to spatial autocorrelation analysis to investigate their population structure. There are significant spatial patterns for most allele frequencies. Clinical patterns, investigated by one-dimensional and directional spatial correlograms, were relatively few in number and were moderate in strength. Overall, however, there is a marked decline in genetic similarity with geographic distance. The results are compatible with a hierarchic population structure superimposed on the geography, and generated by a stochastic fission-fusion model of village propagation, followed by localized gene flow. Strong temporal autocorrelations of allele frequencies based on linguistic-historical distances representing time since divergence were also found. There appears to be a stronger relation between geography and linguistic-historical hierarchic subdivisions than between either feature and genetic distances. These findings confirm by different approaches the results of earlier analyses concerning the important roles of both stochastic and social factors in determining village allele frequencies and the occurrence within this tribe of some allele frequency clines most likely due to the operation of chance historical processes.     

Nucleotide variation at Msn and Alas2, two genes flanking the centromere of the X chromosome in humans

The centromeric region of the X chromosome in humans experiences low rates of recombination over a considerable physical distance. In such a region, the effects of selection may extend to linked sites that are far away. To investigate the effects of this recombinational environment on patterns of nucleotide variability, we sequenced 4581 bp at Msn and 4697 bp at Alas2, two genes situated on either side of the X chromosome centromere, in a worldwide sample of 41 men, as well as in one common chimpanzee and one orangutan. To investigate patterns of linkage disequilibrium (LD) across the centromere, we also genotyped several informative sites from each gene in 120 men from sub-Saharan Africa. By studying X-linked loci in males, we were able to recover haplotypes and study long-range patterns of LD directly. Overall patterns of variability were remarkably similar at these two loci. Both loci exhibited (i) very low levels of nucleotide diversity (among the lowest seen in the human genome); (ii) a strong skew in the distribution of allele frequencies, with an excess of both very-low and very-high-frequency derived alleles in non-African populations; (iii) much less variation in the non-African than in the African samples; (iv) very high levels of population differentiation; and (v) complete LD among all sites within loci. We also observed significant LD between Msn and Alas2 in Africa, despite the fact that they are separated by approximately 10 Mb. These observations are difficult to reconcile with a simple demographic model but may be consistent with positive and/or purifying selection acting on loci within this large region of low recombination.     

广西石灰岩地区蜈蚣蕨居群的遗传多样性研究

采用等位酶分析方法 ,研究了广西石灰岩地区蜈蚣蕨居群的遗传多样性 ,分析了其空间变化趋势。检测了 8个酶系统 ,1 5个酶位点。分析结果表明 :广西石灰岩地区蜈蚣蕨居群遗传多样性程度较高 ,每个位点的等位基因平均数为 1 .6 7,多态位点为 5 7.78% ,平均期望杂合度为 0 .2 4 9。蜈蚣蕨居群的遗传组成在居群间有一定的差异 ,但差异的程度并不与空间距离成正比     

THE RELATIVE IMPORTANCE OF HISTORICAL EVENTS AND GENE FLOW ON THE POPULATION STRUCTURE OF A MEDITERRANEAN RAGWORT,SENECIO GALLICUS (ASTERACEAE)

Comparisons of cytoplasmic and nuclear diversity within and among natural plant populations have the potential to distinguish the relative influences of seed and pollen dispersal on contemporary gene flow, or alternatively, may permit inferences of the colonization history of a species via seed. We examined patterns of cpDNA and allozyme variation in Senecio gallicus, a diploid, annual plant that occurs in both coastal and ruderal inland areas of the Iberian Peninsula and southern France. The species appears to have a strong propensity for long-distance seed dispersal. Five cpDNA haplotypes were found by RFLP analysis among a sample of 111 individuals derived from 11 populations. Differences in haplotype frequencies across populations were most evident with respect to a dramatic increase in the frequency of a derived haplotype from coastal to inland localities. The level of cpDNA differentiation among populations within the inland group (θ₀ = 0.07) was significantly less than that seen within the coastal group (θ₀ = 0.41). In contrast, for allozymes, no significant difference in population structure was evident between collections from coastal and inland habitats. At the rangewide geographic scale, there was only a very weak association between inferred levels of gene flow and geographic distance for cpDNA, and no such association was found for allozymes. It appears that while seed movement in the species might be sufficiently great to disturb the pattern of isolation by distance for cpDNA, it cannot fully account for the nearly randomized spatial structure at polymorphic allozyme loci. It is suggested that isolation of populations in Atlantic-Mediterranean coastal refugia during previous glacial maxima, and the effects of subsequent colonization events in inland areas, have had an important effect on molding the present genetic structure of the species.