Allozyme variation in cork oak (Quercus suber L.): the role of phylogeography and genetic introgression by other Mediterranean oak species and human activities

 Genetic variation in the cork oak (Quercus suber L.) was investigated using 11 loci from seven enzyme systems in 40 populations sampled over the entire distribution of this species in the western Mediterranean Basin. Mean heterozygosity values over the polymorphic loci (Ho=0.283), the percentage of polymorphic populations (M=0.76), and the total genetic diversity (Ht=0.31) from which 11% was accounted for among-population variation, are among the highest recorded in oak species. In contrast to previous results in holm oak (Q. ilex L.), another evergreen species in the same area, cork oak possessed a smaller allele pool and a lower average number of alleles per locus and per population (A=2.0). More particularly, very few low-frequency alleles were observed in cork oak except for eight populations in which allozyme polymorphism at locus Pgi 1, diagnostic between both species, indicates that these low-frequency alleles are introgressed from holm oak. On the basis of the genetic distance estimated from allozyme frequencies, 32 of the 40 cork oak populations studied were classified into two very distinct sets which also corresponded to distinct geographic areas. One set gathered together the 18 populations from the Iberian peninsula and two adjacent areas in France, i.e. the centre of origin of cork oak, according to paleobotanical data. This set was characterized by a larger allele pool, a higher within-population genetic diversity and a lower differentiation between populations than was observed in the other set, which comprised the populations from North Africa, Sicily, Sardinia, Corsica, continental Italy and the region of Provence (southeastern France). In these more southern and eastern disjunct areas, cork oak migration from Iberia may have occurred at different periods since the end of the Tertiary. The possible effect of human activity on cork oak genetic structure, i.e. the selection of good-quality cork, acorn over-use for animal food, and even human nutrition, is discussed. Received: 3 March 1998 / Accepted: 19 March 1998     

Phylogeographical lineages in brown trout (Salmo trutta): investigating microgeographical differentiation between native populations from Northern Spain

SUMMARY 1. The large microgeographical differentiation revealed by allozyme studies in brown trout ( Salmo trutta) populations is one of the most striking features of this species. Additionally, allozymes showed great genetic differences between Atlantic and Mediterranean populations on a macrogeographical scale.
2. This study was carried out in order to assess whether the great differences observed between Atlantic and Mediterranean populations persisted where the two are geographically close (the 'microgeographical scale'). Sixteen populations of brown trout, S. trutta , were screened for genetic variation at 25 allozyme loci. The sampling sites, which occupied a relatively small geographical area, were distributed across Cantabrian (Atlantic) and Mediterranean drainages in Northern Spain.
3. The neighbour-joining tree, inferred from Nei's genetic distance, showed that brown trout populations clustered into two different groups. These groups corresponded to the Cantabrian and the Mediterranean groups of populations, although no clear geographical pattern emerged within each of the groups. This geographical pattern is basically caused by significant differences in the frequency distribution of the CK-A1 * locus, with a higher frequency of * 115 in Cantabrian samples (0.586 ± 0.091) while allele * 100 was more frequent in Mediterranean samples (0.931 ± 0.038). In addition, this study revealed alleles exclusive to the Mediterranean and Cantabrian populations, agreeing with previous findings.
4. Genetic differentiation between Cantabrian and Mediterranean regions (14.19%) was similar to that estimated in Spain at a larger scale (13%), showing that most of the differences between the regions can be observed even in a small geographical area.     

Genetic variation at allozyme and RAPD loci in sessile oak Quercus petraea (Matt.) Liebl.: the role of history and geography

The nuclear genetic variation within and among 21 populations of sessile oak was estimated at 31 RAPD loci in conjunction with previous estimates of variation at eight allozyme loci. The aim of the study was to assess the relative role of isolation-by-distance and postglacial history on patterns of nuclear variation. Because of its small effective population size and maternal transmission, the chloroplast genome is a good marker of population history. Both kinds of nuclear variation (RAPD and allozyme) were therefore compared, first, to the geographical distances among populations and, secondly, to chloroplast DNA restriction polymorphism in the same populations. Multiple Mantel tests were used for this purpose. Although RAPDs revealed less genetic diversity than allozymes, levels of genetic differentiation ( G _ST) were identical. The standard genetic distance calculated at all RAPD loci was correlated with geographical distances but not with the genetic distance calculated from chloroplast DNA data. Conversely, allozyme variation was correlated with chloroplast DNA variation, but not with geography. Possibly, divergent selection at two allozyme loci during the glacial period could explain this pattern. Because of its greater number of loci assayed, RAPDs probably provided a less biased picture of the relative role of geography and history.     

Allozyme, mitochondrial-DNA, and morphometrie variability indicate cryptic species of anchovy (Engraulis encrasicolus)

Previous surveys of population structure in the Atlantic-Mediterranean anchovy Engraulis encrasicolus L. have reported heterogeneity in morphology, allozyme frequencies, and mitochondrial DNA haplotype frequencies at a regional scale. In particular, two stocks of anchovy have been detected in the Adriatic Sea. In this paper, the available data is reviewed with the aim to relate genetic variation to geography at the widest possible geographical scale, for investigating the evolutionary mechanisms underlying stock structure in anchovy. Correspondence analysis of allozyme frequencies (24 samples, three polymorphic loci) compiled from the literature indicates three distinct entities in the Mediterranean Sea. Open-sea or oceanic anchovy populations are genetically different from inshore-water populations within a region (Nei's ^ G _ST = 0.035–0.067), while broadscale geographical variation is weak for each of these two habitat-specific forms (^ G _ST = 0.005–0.006). Mitochondrial-DNA haplotype frequencies support the distinction between an inshore form and an oceanic form (^ G _ST = 0.067–0.107), with virtually no genetic differences among oceanic populations across the Gulf of Biscay, the western Mediterranean and the Ionian Sea (^ G _ST = −0.001). If natural selection on marker loci is unimportant, these results indicate the occurrence of two parapatric, genetically distinct, habitat-specific forms that are widely distributed throughout the Mediterranean Sea. Persistent allele and haplotype-frequency differences between these forms indicate reproductive isolation and the presence of an E. encrasicolus species complex in the Mediterranean. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society 75 : 261–269.     

High variability of chloroplast DNA in three Mediterranean evergreen oaks indicates complex evolutionary history

Chloroplast DNA variation was studied in three evergreen Quercus species (Q. suber L., Q. ilex L. and Q. coccifera L.) from the Western Mediterranean Basin using PCR-RFLP. We studied five primer pair/enzyme combinations, four of them previously used in other European Quercus, obtaining a large number of haplotypes (81) grouped in three main types (suber type, ilex-coccifera I type and ilex-coccifera II type). Such level of haplotype diversity is higher than previously reported for the genus. Remarkable differences in haplotype richness between species have been found. Q. ilex and Q. coccifera usually share the same haplotypes, while a number of Q. suber populations possesses variants of the ilex-coccifera I type. This fact is interpreted as a result of genetic introgression between Q. suber and Q. ilex. Reproductive factors that could determine this exchange are discussed, as well as the influence of different species histories on the present structure of evergreen Quercus in the Western Mediterranean Basin.     

Phenotypic and allozyme variation in Mediterranean and desert populations of wild barley,Hordeum spontaneum Koch

Populations of wild barley, Hordeum spontaneum Koch, were collected in two distinct climatic regions, desert and Mediterranean. Plants from five desert and five Mediterranean populations were compared and contrasted for extent and structure of phenotypic variation. These same 10 and one other population from each region were analyzed for allozyme variation. In a field trial of phenotypic diversity, two phenological and 14 morphological traits were examined. Study of allozyme variation was performed using eight enzyme systems encoding for 13 loci. Plants from the desert and Mediterranean regions were significantly different in seven of 16 phenotypic traits, exhibited a high (30%) interregional component of phenotypic variation, and showed a high degree of segregation on a principal component scattergram indicating ecotypic differentiation. Mediterranean populations were twice as variable as desert populations in reproductive growth parameters (stem and spike length) and grain filling (spikelet weight), but half as variable for onset of reproduction. The extent and structure of phenotypic and allozyme variation did not match. The Mediterranean and desert populations did not differ in amount of allozyme variation as estimated by mean number of alleles per locus, effective number of alleles, polymorphism, and gene diversity (n(a), n(e), P, and H(e)), did not segregate on the basis of population genetic distances, and exhibited a low proportion of interregion allozyme diversity (2%). No effect of selection on allozyme distribution was detected. Our results suggest that the adaptation of plants originating from desert and Mediterranean environments is reflected in phenotypic but not in allozyme variation.     

The distribution of Quercus suber chloroplast haplotypes matches the palaeogeographical history of the western Mediterranean

Combining molecular analyses with geological and palaeontological data may reveal timing and modes for the divergence of lineages within species. The Mediterranean Basin is particularly appropriate for this kind of multidisciplinary studies, because of its complex geological history and biological diversity. Here, we investigated chloroplast DNA of Quercus suber populations in order to detect possible relationships between their geographical distribution and the palaeogeographical history of the western Mediterranean domain. We analysed 110 cork oak populations, covering the whole distribution range of the species, by 14 chloroplast microsatellite markers, among which eight displayed variation among populations. We identified five haplotypes whose distribution is clearly geographically structured. Results demonstrated that cork oak populations have undergone a genetic drift geographically consistent with the Oligocene and Miocene break-up events of the European-Iberian continental margin and suggested that they have persisted in a number of separate microplates, currently found in Tunisia, Sardinia, Corsica, and Provence, without detectable chloroplast DNA modifications for a time span of over 15 million years. A similar distribution pattern of mitochondrial DNA of Pinus pinaster supports the hypothesis of such long-term persistence, in spite of Quaternary climate oscillations and of isolation due to insularity, and suggests that part of the modern geographical structure of Mediterranean populations may be traced back to the Tertiary history of taxa.     

Allozyme diversity in Japanese populations of Hemerocallis thunbergii, H. middendorffii, and H. exaltata (Liliaceae)

Eight populations of Hemerocallis thunbergii, H. middendorffii , and H. exaltata from Japan were analyzed by starch gel electrophoresis to measure genetic variation at ten allozyme loci. The former two species are widely distributed in eastern Asia, whereas H. exaltata is restricted to the Tobishima and Sado Islands in western Japan. Populations of the three species maintain high levels of allozyme variation: mean expected heterozygosities were 0.288, 0.209, and 0.169, respectively. The UPGMA phenogram obtained from Nei's genetic distance values indicated that H. exaltata is allozymically the most distinct species, whereas populations of H. thunbergii and H. middendorffii are closely related. However, based on gross morphology H. exaltata and H. middendorffii are closely related. This discrepancy and the observed high levels of allozyme variation within populations of the three Hemerocallis species in Japan may in part be due to genetic change among populations of the two Japanese archipelago species.     

Genetic structure of the widespread and common Mediterranean bryophyte Pleurochaete squarrosa (Brid.) Lindb. (Pottiaceae) - evidence from nuclear and plastidic DNA sequence variation and allozymes

The Mediterranean Basin as one the world's most biologically diverse regions provides an interesting area for the study of plant evolution and spatial structure in plant populations. The dioecious moss Pleurochaete squarrosa is a widespread and common bryophyte in the Mediterranean Basin. Thirty populations were sampled for a study on molecular diversity and genetic structure, covering most major islands and mainland populations from Europe and Africa. A significant decline in nuclear and chloroplast sequence and allozyme variation within populations from west to east was observed. While DNA sequence data showed patterns of isolation by distance, allozyme markers did not. Instead, their considerable interpopulation genetic differentiation appeared to be unrelated to geographic distance. Similar high values for coefficients of gene diversity (G(ST)) in all data sets provided evidence of geographic isolation and limited gene flow among populations (i) within islands, (ii) within mainland areas, and (iii) between islands and mainland. Notably, populations in continental Spain are strongly genetically isolated from all other investigated areas. Surprisingly, there was no difference in gene diversity and G(ST) between islands and mainland areas. Thus, we conclude that large Mediterranean islands may function as 'mainland' for bryophytes. This hypothesis and its implication for conservation biology of cryptogamic plants warrant further investigation. While sexually reproducing populations were found all over the Mediterranean Basin, high levels of multilocus linkage disequilibrium provide evidence of mainly vegetative propagation even in populations where sexual reproduction was observed.     

Genetic structure of Araucaria angustifolia (Araucariaceae) populations in Brazil: implications for the in situ conservation of genetic resources

The distribution of the genetic variation within and among natural populations of A. ANGUSTIFOLIA growing in different regions in Brazil was assessed at microsatellite and AFLP markers. Both markers revealed high gene diversity ( H = 0.65; AR = 9.1 for microsatellites and H = 0.27; P = 77.8 % for AFLPs), moderate overall differentiation ( RST = 0.13 for microsatellites and FST = 0.10 for AFLPs), but high divergence of the northernmost, geographically isolated population. In a Bayesian analysis, microsatellite data suggested population structure at two levels: at K = 2 and at K = 3 in agreement to the geographical distribution of populations. This result was confirmed by the UPGMA dendrogram based on microsatellite data (bootstrap support > 95 %). Non-hierarchical AMOVA revealed high variation among populations from different A POSTERIORI defined geographical groups. The genetic distance between sample locations increased with geographical distance for microsatellites ( R = 0.62; P = 0.003) and AFLPs ( R = 0.32; P = 0.09). This pattern of population differentiation may be correlated with population history such as geographical isolation and postglacial colonization of highlands. Implications of the population genetic structure for the conservation of genetic resources are discussed.     

Natural hybridisation between kermes (Quercus coccifera L.) and holm oaks (Q. ilex L.) revealed by microsatellite markers

Hybridisation between species of the genus Quercus is a common phenomenon as a result of weak reproductive isolation mechanisms between phylogenetically close species that frequently co-occur in mixed stands. In this study, we use microsatellite markers to analyse introgression between kermes ( Quercus coccifera L.) and holm ( Q. ilex L.) oak, two closely related taxa that frequently dominate the landscape in extensive areas in the Mediterranean region. All tested microsatellites amplified and were polymorphic in both kermes and holm oaks. Bayesian admixture analyses showed a good correspondence between each species and one of the two inferred genetic clusters. Five sampled individuals were a priori tentatively identified as hybrids on the basis of intermediate morphological characteristics, and it was confirmed that they also presented mixed genotypes. However, we also detected different levels of genetic introgression among morphologically pure individuals, suggesting that successful backcrossing and/or reduced phenotypic expression of genetic variance in certain individuals may have resulted in strong convergence towards a single species phenotype.     

Geographical variability of allozyme variation and genetic divergence between Tunisian lizards of the sand-dwelling lacertid genus Acanthodactylus

Acanthodactylus is a widespread lacertid genus including about thirthy species whose systematics remains much debated.The present study is aimed at improving the taxonomic status of the Tunisian species whose geographical distribution and ecological requirements have been studied elsewhere. Allozymic variation at 22 loci was scored in 14 populations for a total of 249 individuals. High genetic variability is confirmed in all samples; it is higher in species displaying a wide bioclimatic habitat range. No reduction of genetic variability was observed in two insular populations.Any interspecific hybrid could be detected from our genetic data, even in sympatric areas arguing in favor of specific status of four taxons in Northern, Central and Southern Tunisia.Phylogenetic relationships are discussed on the basis of UPGMA and Fitch-Margoliash procedures applied to the genetic distances matrix and centered data analysis of allozyme frequencies.     

Population genetics, phylogeography and speciation of Cystodytes (Ascidiacea) in the western Mediterranean Sea

Several morphotypes that so far have been attributed to the allegedly cosmopolitan ascidian Cystodytes dellechiajei occur in the Mediterranean Sea. Colour variation is the difference most frequently reported. In this study, we addressed the genetic structure of this ascidian in relation to geographical location and colour morph. Partial sequences of the gene cytochrome  c oxidase subunit 1 (COI) were obtained from seven populations of the western Mediterranean, encompassing eight colour varieties. All population genetic analyses (exact test, pairwise F _ST, hierarchical analysis of molecular variance, multidimensional scaling, nested clade analysis) indicated clearly that differences between colour morphs are large enough to obscure any geographical differentiation when colours are combined within localities. When variance due to colour divergence was removed, however, a significant geographical variability between localities remained. The genetic divergence between the colour morphs analysed was significant in comparisons of the brown and purple forms with the others, but not among the green, blue, and white morphs. Phylogeographic analyses suggest that population fragmentation and range expansions have shaped the present-day distribution of the haplotypes. Taken together with existing chemotype information, our results indicate that several species are present in the area, and that a thorough revision of the genus is necessary.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 203–214.     

Genetic structure of Atriplex halimus populations in the Mediterranean Basin

BACKGROUND AND AIMS: The saltbush Atriplex halimus is a chenopodiaceous plant well adapted to dry saline habitats and widely distributed in the Mediterranean Basin. A study was carried out to analyse the genetic diversity of A. halimus at the level of the Mediterranean Basin. METHODS: To assess the intra- and interpopulational variation of A. halimus a total of 51 populations and six plants per populations was analysed with the RAPD-PCR technique. For the study of the phylogeny of the populations, 21 samples of A. halimus and seven samples of other species of Atriplex were analysed by the sequencing of the ITS (internal transcribed spacer) region of the ribosomal DNA. KEY RESULTS: The AMOVA analysis of the RAPD results showed that populations were divided into two discrete genetic groups, as the variation among groups accounted for 54.36 % of the total variance of the collection. At the same time, the intrapopulational diversity was high, as 301 out of 306 plants analysed constituted an individual RAPD haplotype. The sequencing of the ITS region also showed a significant separation of the two genetic groups, with a genetic distance of 0.023 nucleotide substitutions per site. Using A. breweri, A. canescens, A. glauca and A. prostrata as outgroups in the phylogenetic analysis, A. breweri and A. canescens are the species closest to A. halimus from this group, while A. prostrata is the most distant. CONCLUSIONS: The present work indicates that two genetic groups of A. halimus can be distinguished after analysing the genetic diversity of 51 populations from ten countries in the Mediterranean Basin.     

Comparative anatomical analysis of the cotyledonary region in three Mediterranean Basin Quercus (Fagaceae)

Anatomical changes at the cotyledonary node from the embryo to the seedling stage in Quercus coccifera, Q. ilex, and Q. humilis were investigated by light and scanning electron microscopy techniques. Mature embryos of Q. humilis possess 2-3 pairs of leaf primordia and a pair of cotyledonary buds, whereas in Q. coccifera and Q. ilex there are two incipient primordia, and cotyledonary buds are not observed until 1 wk after germination. In all three species the cotyledonary buds multiply, forming bud clusters, and a vascular connection is well established within 5-6 wk after germination. As development proceeds, the cotyledonary region becomes woody, but buds, which are exogenous in origin, never become embedded in the periderm. In comparison with Q. suber, another native Mediterranean Basin oak, the cotyledonary node is short and axillary buds are not present below the insertion of cotyledons. In addition, starch accumulation in the cotyledonary region is not observed from histological analysis in the three oaks. Therefore, in Q. coccifera, Q. ilex, and Q. humilis seedlings the cotyledonary node can be considered to be an important regenerative structure enabling them to resprout after the elimination of the shoot above the cotyledons, despite the absence of a lignotuberous structure.     

Dactylorhiza majalis s.l. (Orchidaceae) in acid habitats: variation patterns, taxonomy, and evolution

Morphometric and allozyme data were collected from 11 populations of Dactylorhiza majalis s.l. in the Ardennes (Belgium, France) and North Jutland (Denmark). All study populations were growing in obviously acid habitats. Genetic distances were calculated and illustrated by UPGMA dendrograms, while principal components analyses (PCAs) were used to summarize morphological variation patterns. Based on the PCAs, the probability of each character to distinguish correctly between the main groups of populations was estimated. The study populations from the Ardennes and from the Danish island of Læsø seem to agree with populations from Sweden previously studied by Hedrén—and to belong to D. majalis subsp. sphagnicola (comb. nov.). It should be noted that allozyme data suggest a polytopic origin of this taxon. The study populations from the Danish region of Thy differ morphologically from the others and exhibit much less genetic variation. They are described as a new subspecies, D. majalis subsp. calciñigiens. The evolution and biogeography of the two taxa are discussed, and brief taxonomic accounts are provided.     

Micro-evolution in the weevil genus Larinus: the formation of host biotypes and speciation

Data are presented on allozyme variation between 15 populations of the stenophagous capitulum weevil, Larinus cynarae , and three populations of its congener, L. latus , that had been collected throughout the northern mediterranean range of these species. A phenetic analysis of these data revealed no direct relationship between genetic variation and host-plant association within L. cynarae , but there was a strong geographical structuring of allozyme patterns. Most of the genetic variation was due to differences between geographical regions and variation within these was small. Wright's F _ST values showed that Italian and Greek populations of L. cynarae were most distinct from L. latus , with southern Iberian, northern Spanish and French populations increasingly less so. This pattern was associated with a cline in the frequencies of certain alleles along this geographical arc from France to Greece. A phenogram of Nei's genetic distances indicated the close genetic relationship between the two species of Larinus and separated the populations of L. cynarae into three allopatric groups. These groups have different host-plant spectra - dominated by Cynara cardunculus in Italy and Greece, Cynara humilis/Onopordum in southern Iberia and Onopordum spp. in France/Northern Spain - and can be considered to be host biotypes of L. cynarae. L. latus , which occurs in Greece and further east is also an Onopordum specialist. An analysis of the phylogeny of this group of Larinus indicates a primary separation into eastern ( L. latus ) and western ( L. cynarae ) taxa, with further branching of the L. cynarae lineage into the putative host-biotypes. An hypothesis for the evolution of these taxa is given, based on the evolutionary history of host-plant taxa and geographical constraints.     

Population genetic variation and structure of the invasive weed Mikania micrantha in southern China: consequences of rapid range expansion

Invasive plants such as Mikania micrantha provide valuable opportunities for studying population genetic consequences of rapid range expansion. Twenty-eight populations of M. micrantha throughout its introduced range in southern China were examined by using intersimple sequence repeat markers. Population genetic parameters were estimated by Bayesian approaches as well as conventional methods. Bottleneck signature, multilocus linkage disequilibrium, character compatibility, and cluster analyses were conducted to assay the factors that may act to shape population variability. High levels of genetic variation and differentiation were detected in the introduced populations of M. micrantha. All populations experienced severe bottlenecks. Most of them demonstrated significant linkage disequilibrium and matrix compatibility. Populations were mainly clustered into 2 groups, and those from different regions intermingled in the unweighted pair group method with arithmetic mean (UPGMA) dendrogram. No geographical signature was found in the pattern of population genetic variation. This research indicates that during M. micrantha invasion, multiple introductions mitigated the loss of genetic variation associated with bottlenecks. Nonetheless, bottlenecks enhanced the population differentiation. Human-mediated long-distance dispersal events of seeds or propagules explain the lack of geographic structure in genetic variation. Although asexual reproduction is the predominant mating mode in M. micrantha, it has little effect on the population genetic composition.     

Discordant patterns of nuclear and mitochondrial introgression in Iberian populations of the European common frog (Rana temporaria)

Amphibians often show complex histories of intraspecific and interspecific genetic introgression, which might differ in mitochondrial and nuclear genes. In our study of the genetic differentiation of the European common frog, Rana temporaria (159 specimens from 23 populations were analyzed for 24 presumptive allozyme loci; 82 specimens were sequenced for a 540-bp fragment of the mitochondrial 16S rRNA gene), multilocus correspondence analysis (CA) and Bayesian assignment tests of the nuclear data were concordant in identifying 2 population groups corresponding to 1) the Pyrenees in the east and 2) the Galicia and Asturias regions in the west, the latter corresponding to the subspecies R. temporaria parvipalmata. Geographically intermediate populations were genetically intermediate in the allozyme CA and, less clearly in the Bayesian assignment, with mitochondrial haplotypes exclusively belonging to the parvipalmata group. This indicates different degrees of introgression in the mitochondrial and nuclear genomes. Although Pyrenean high-altitude populations are morphologically distinct from low-altitude populations, these 2 groups were not separate clusters in any analysis. This suggests that the morphological differences may be due to fast adaptation to elevational gradients, likely under maintenance of gene flow, and that the underlying genetic changes are not detectable by the analyzed markers. We argue that a parsimonious explanation for the observed pattern along the east-west axis in northern Spain may be competition between invading and resident populations, with no need to invoke selection. However, in order to conclusively rule out selective processes, additional and finer scale data are required to test for asymmetric mating preference/behaviour, sex-biased gene flow, or sex-biased survival of potential hybrids.     

The Balearic Islands: a reservoir of cpDNA genetic variation for evergreen oaks

Aim To analyse the role of the Balearic Islands as a refuge area for evergreen Quercus (cork oak: Quercus suber L., holm oak: Q. ilex L., kermes oak: Q. coccifera L.), by using molecular, historical and palaeobotanical data. Location The Western Mediterranean Basin (Balearic Islands, eastern Iberia, Provence, Sardinia, Corsica, Sicily, Malta, Italy, Northern Africa). Methods We sampled 108 populations and used the PCR‐RFLP technique with five universal cpDNA primers to define haplotypes in the sampled populations. Diversity, differentiation parameters and spatial analysis of the populations, using a spatial version of amova , were linked to the geological history of the Western Mediterranean Basin in order to explain the present spatial pattern of the evergreen Quercus populations in the Balearics. Results Evergreen Quercus cpDNA shows a complex structure, with remnants of ancient diversity in the Balearics. Balearic populations of holm oak are related to Iberian populations, while for cork and kermes oaks, we found both Tyrrhenian and Iberian haplotypes. Main conclusions The complex spatial patterns of cpDNA in Balearic evergreen Quercus appears explicable in terms of a combination of physical (vicariance and long distance dispersal) and biological (introgressive hybridization) factors. The Balearics constitute a glacial refuge area and a reservoir of genetic variation with traces of ancient diversity from Messinian–Pliocene stages.