Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations

Intraspecific phylogeographical patterns largely depend on the life history traits of a species. Especially species with a high degree of cold tolerance, limited requirements towards habitat preferences, and relatively low active dispersal capacities may have responded in a different way to the Pleistocene climatological fluctuations than the majority of taxa studied so far. To evaluate this possibility, we studied Arion fuscus (Muller, 1774), a common and widespread European terrestrial slug, from 88 locations (N = 964). Sequence variation was assessed for fragments of the mitochondrial 16S rDNA and COI genes by means of single-strand conformation polymorphisms (SSCP) and subsequent DNA sequencing. Additionally, eight allozyme loci were scored in 843 individuals. Phylogenetic analysis revealed the presence of two major evolutionary lineages, one in the Balkan region and another in the Alps and the rest of Europe. The sequence divergence between the two lineages was limited (3.3%), but gene flow between the regions was absent, suggesting that the two regions have been isolated since the late Pliocene or early Pleistocene. Allozyme differentiation among geographical regions and mitochondrial DNA (mtDNA) lineages was low. The geographical patterns observed in our data showed that (i) haplotype and nucleotide diversities are very low in northern Europe, suggesting that single haplotypes rapidly colonized large areas; (ii) recently expanded haplotype clades have restricted distribution ranges, suggesting that current gene flow is low; and (iii) genetic diversity in the Alps is much higher than in other regions and estimated past gene flow from the Eastern Alps to other regions was high, suggesting that this was a refugial zone during the Pleistocene. This full-range phylogeography suggests the existence of an alternative refugial zone, situated north of the refugial areas currently recognized in most other taxa.     

Integrating molecular ecology and predictive modelling: implications for the conservation of the barbastelle bat (Barbastella barbastellus) in Portugal

Here we employ multiple approaches (population genetics, past and present predictive modelling) to determine the population structure and the potential distribution in Portugal of the barbastelle, Barbastella barbastellus. Population structure was analysed by sequencing two mitochondrial DNA fragments: cytochrome-b and D-loop. The latter fragment showed almost no variation with only three haplotypes found. In contrast, 13 haplotypes were identified in cytochrome-b sequences although without a significant population structure despite the existence of unique haplotypes in the northern and southern regions. Predictive modelling showed that, under current conditions, the distribution of B. barbastellus is fragmented and highly dependent on the existence of native woodland, with the majority of the suitable areas located in the centre and north of the country, while the predicted occurrence in the south was very restricted. Our bioclimatic models for the Last Glacial Maximum also exhibited a similar pattern which, in agreement with the location of the unique haplotypes, suggests that some degree of isolation has existed between northern and southern populations. Nevertheless, the separation of these populations is not significant in the genetic analysis; hence, Portuguese populations should be managed as a single unit, where the maintenance of the existing mature native woodland is paramount for conservation. In conclusion, this study illustrates how combining genetic analysis with predictive modelling can help in the conservation management of a poorly documented species in a rapid and low-cost manner.     

Postglacial colonization of Northern Europe by distinct phylogeographic lineages of the bullhead, Cottus gobio

Three major phylogeographic lineages of the cottid fish Cottus gobio (bullhead) were identified in northern Europe from mitochondrial DNA sequences and allozyme data. The largely separate freshwater distributions of the lineages demonstrate distinct postglacial colonization histories. West of the Baltic Sea, Swedish lakes were invaded from the southwest (Germany). Another, eastern lineage has colonized the inland waters northeast and east of the Baltic, from refugia in northwest Russia; this lineage comprises a distinct subgroup found only from Estonia. The third lineage, found south and southeast of the Baltic, probably descended from rivers draining to the Black Sea from the north (e.g. Dnepr). In coastal waters of the Baltic Sea, and in near-coast inland waters, the lineages are now found intermixed in various combinations. The alternating fresh- and saltwater phases of the Baltic basin have variously enabled and disabled the use of coastal waters as colonization routes. Hypotheses on the chronology of dispersal and lineage mixing can be based on the distribution of the marker genes and the paleohydrographical record. The diversity of the Fennoscandian bullhead thus comprises anciently diverged (probably mid-Pleistocene) refugial lineages that in their freshwater range constitute distinct evolutionarily significant units. The thorough mixing of the various genomic origins in and around the Baltic, however, refutes the controversial view of distinct species status for the western and eastern ('Cottus koshewnikowi') bullheads. The postglacial contact of the lineages has created new diversity that cannot be interpreted in a conventional hierarchical framework of taxonomic or conservation units.     

Alleles of sAAT-1,2* isoloci in brown trout: potential diagnostic marker for tracking routes of post-glacial colonization in northern Europe

Distribution alleles of sAAT-1,2 * isoloci in brown trout Salmo trutta populations of northern Europe revealed two glacial refugia, 'Baltic' and 'Atlantic'.     

Reconstructing the demise of Tethyan plants: climate‐driven range dynamics of Laurus since the Pliocene

Aim Climate changes are thought to be responsible for the retreat and eventual extinction of subtropical lauroid species that covered much of Europe and North Africa during the Palaeogene and early Neogene; little is known, however, of the spatial and temporal patterns of this demise. Herein we calibrate ecological niche models to assess the climatic requirements of Laurus L. (Lauraceae), an emblematic relic from the Tethyan subtropical flora, subsequently using these models to infer how the range dynamics of Laurus were affected by Plio‐Pleistocene climate changes. We also provide predictions of likely range changes resulting from future climatic scenarios. Location The Mediterranean Basin and Macaronesian islands (Canaries, Madeira, Azores). Methods We used a maximum‐entropy algorithm (Maxent) to model the relationship between climate and Laurus distribution over time. The models were fitted both to the present and to the middle Pliocene, based on fossil records. We employed climatic reconstructions for the mid‐Pliocene (3 Ma), the Last Glacial Maximum (21 ka) and a CO₂‐doubling future scenario to project putative species distribution in each period. We validated the model projections with Laurus fossil and present occurrences. Results Laurus preferentially occupied warm and moist areas with low seasonality, showing a marked stasis of its climatic niche. Models fitted to Pliocene conditions successfully predicted the current species distribution. Large suitable areas existed during the Pliocene, which were strongly reduced during the Pleistocene, but humid refugia within the Mediterranean Basin and Macaronesian islands enabled long‐term persistence. Future climate conditions are likely to re‐open areas suitable for colonization north of the current range. Main conclusions The climatic requirements of Laurus remained virtually unchanged over the last 3 Myr. This marked niche conservatism imposed largely deterministic range dynamics driven by climate conditions. This species's relatively high drought tolerance might account for the survival of Laurus in continental Europe throughout the Quaternary whilst other Lauraceae became extinct. Climatic scenarios for the end of this century would favour an expansion of the species's range towards northern latitudes, while severely limiting southern populations due to increased water stress.     

Species delimitation and phylogeography of the Abies chensiensis complex inferred from morphological and molecular data

Species delimitation in the Abies chensiensis complex, consisting of A. chensiensis, A. recurvata, A. ernestii and A. ernestii var. salouenensis, has been the subject of a long‐term dispute, as various taxonomic combinations have been proposed. We combined different lines of evidence, including morphological characters and mitochondrial and plastid DNA sequence data, to assess species delimitation. Plastid DNA sequence variation was highly consistent with the morphological result and separated A. recurvata from the other three taxa, but the mitochondrial DNA genealogy was more strongly related to geographical distribution rather than species delimitation, as proposed in previous studies. On the basis of morphological characters and plastid DNA sequence variation, the current species status of A. recurvata was well supported, and we also accepted A. ernestii and A. ernestii var. salouenensis as varieties of A. chensiensis. Furthermore, the phylogeographical history of A. chensiensis was surveyed by molecular data and ecological niche modelling. On the basis of molecular data, the current phylogeographical pattern of the A. chensiensis complex was probably shaped by recent rapid expansions from multiple refugia in the Hengduan Mountains and habitat contraction after the last glaciation. This hypothesis was supported by the results of ecological niche modelling and previous fossil records. Our findings suggest that these fir taxa might have continued to expand their range after the largest glaciation until the end of the last glaciation, but have contracted since then with the increasing temperature and habitat changes. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 175–188.     

Phylogeography and ecological niche modelling of the New Zealand stick insect Clitarchus hookeri (White) support survival in multiple coastal refugia

Aim Increasing our understanding of the effects of the Last Glacial Maximum (LGM) and determining the location of refugia requires studies on widely distributed species with dense sampling of populations. We have reconstructed the biogeographic history of Clitarchus hookeri (White), a widespread species of New Zealand stick insect that exhibits geographic parthenogenesis, using phylogeographic analysis and ecological niche modelling. Location New Zealand. Methods We used DNA sequence data from the mitochondrial cytochrome c oxidase subunit I gene to reconstruct phylogenetic relationships among haplotypes from C. hookeri and two undescribed Clitarchus species. We also used distribution data from our own field surveys and museum records to reconstruct the geographic distribution of C. hookeri during the present and the LGM, using ecological niche modelling. Results The ecological niche models showed that the geographic distribution of C. hookeri has expanded dramatically since the LGM. Our model predicted large areas of suitable LGM habitat in upper North Island, and small patches along the east coast of South Island. The phylogeographic analysis shows that populations in the northern half of North Island contain much higher levels of genetic variation than those from southern North Island and South Island, and is congruent with the ecological niche model. The distribution of bisexual populations is also non-random, with males completely absent from South Island and very rare in southern North Island. Main conclusions During the LGM C. hookeri was most likely restricted to several refugia in upper North Island and one or more smaller refugia along the east coast of South Island. The unisexual populations predominate in post-glacial landscapes and are clearly favoured in the recolonization of such areas. Our study exemplifies the utility of integrating ecological niche modelling and phylogeographic analysis.     

Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen

Norway spruce (Picea abies [L.] Karst.) is a broadly distributed European conifer tree whose history has been intensively studied by means of fossil records to infer the location of full‐glacial refugia and the main routes of postglacial colonization. Here we use recently compiled fossil pollen data as a template to examine how past demographic events have influenced the species’ modern genetic diversity. Variation was assessed in the mitochondrial nad1 gene containing two minisatellite regions. Among the 369 populations (4876 trees) assayed, 28 mitochondrial variants were identified. The patterns of population subdivision superimposed on interpolated fossil pollen distributions indicate that survival in separate refugia and postglacial colonization has led to significant structuring of genetic variation in the southern range of the species. The populations in the northern range, on the other hand, showed a shallow genetic structure consistent with the fossil pollen data, suggesting that the vast northern range was colonized from a single refugium. Although the genetic diversity decreased away from the putative refugia, there were large differences between different colonization routes. In the Alps, the diversity decreased over short distances, probably as a result of population bottlenecks caused by the presence of competing tree species. In northern Europe, the diversity was maintained across large areas, corroborating fossil pollen data in suggesting that colonization took place at high population densities. The genetic diversity increased north of the Carpathians, probably as a result of admixture of expanding populations from two separate refugia.     

Species' diversity in the New Caledonian endemic genera Cephalidiosus and Nobarnus (Insecta: Heteroptera: Tingidae), an approach using phylogeny and species' distribution modelling

The patterns of local endemism in New Caledonia were analysed in two endemic genera of Tingidae (Insecta, Heteroptera), Cephalidiosus and Nobarnus , through a phylogenetic analysis and species' distribution modelling. The aim was to determine the possible causes of diversification and endemism in New Caledonia. Our results show that environmental conditions are probably important for the distribution of the genus Cephalidiosus , in conjunction with other factors such as resource (host plant) distribution, but suggest that the same environmental conditions have not influenced the speciation processes and diversification in the genus.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 177–184.     

Mitochondrial haplotype diversity among Portuguese brown trout Salmo trutta L. populations: relevance to the post-pleistocene recolonization of northern Europe

Mitochondrial haplotype diversity in seven Portuguese populations of brown trout, Salmo trutta L., was investigated by sequencing the 5' end of the mitochondrial DNA (mtDNA) control region. Five new haplotypes were described for this species, each two to three mutational steps distant from the common north Atlantic haplotype. Significant population subdivision of mtDNA haplotypes was also apparent. Based on these results, as well as on published data describing the distribution of both mtDNA haplotypes and allozyme alleles throughout Europe, the postglacial recolonization of northern Europe was re-evaluated. It is argued that the available data do not support the contribution of two major glacial refugia (southwest Atlantic and Ponto-Caspian Basin) to this postglacial recolonization, as proposed in a recently published model. The unique genetic architecture of Portuguese brown trout within the Atlantic-basin clade of this species represents a highly valuable genetic resource that should be protected from introgression with nonendemic strains of hatchery fish.     

Complex patterns of colonization and refugia revealed for European grayling Thymallus thymallus,based on complete sequencing of the mitochondrial DNA control region

The complete mitochondrial DNA (mtDNA) control region (1043 base pairs) and 162-bp of flanking transfer RNA genes were sequenced in 316 European grayling, Thymallus thymallus, from 44 populations throughout the Western European range of the species. A total of 58 haplotypes were revealed with pairwise divergence ranging from 0.001 to 0.038. An inferred intraspecific phylogenetic tree revealed two well-supported clades within the Danube basin, one highly divergent clade in the Adriatic basin, and one large, diverse group representing most other populations. A deeply divergent haplotype fixed in the Loire basin in central France, more groups of haplotypes from distinct Danubian tributaries, and a relatively ancestral haplotype fixed in former tributaries of the Elbe in Denmark all suggest a complex pattern of interglacial and postglacial expansions originating from disjunct refugia throughout central Europe. Despite some evidence of human-mediated stock transfers, parsimony-network-based nested-clade analysis (NCA) supported specific inferences relating to corridors of postglacial expansion such as the lower Rhine (Moselle) and Elbe systems (Danish populations) serving as sources for expansion into the Baltic to the north as well as the upper Rhine and Danube to the south; and specific Rhine populations (Doller, Orbe and Reuss) serving as sources for colonization of the Rhone. The multiple divergent clades representing populations in the upper Danube, as well as the deeply divergent haplotypes found in the Adriatic and Loire basins (> 5% divergence from Asian outgroups) support the theory that European grayling have had a long history in Western Europe, pre-dating Pleistocene glacial cycles. The patterns of mtDNA divergence shown here support a perspective of rich inter- and intrabasin genetic diversity that should be protected from current trends to translocate brood stocks for rearing and release in response to declining populations, especially in southern European basins.     

Phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: contrasting results from mitochondrial and microsatellite data

Phylogeographical studies are typically based on haplotype data, occasionally on nuclear markers such as microsatellites, but rarely combine both. This is unfortunate because the use of markers with contrasting modes of inheritance and rates of evolution might provide a more accurate and comprehensive understanding of a species' history. Here we present a detailed study of the phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum , using 1098 bp of the mitochondrial ND2 gene from 45 localities from across its Palaearctic range to infer population history. In addition, we re-analysed a large microsatellite data set available for this species and compared the results of both markers to infer population relationships and the historical processes influencing them. We show that mtDNA, the most popular marker in phylogeography studies, yielded a misleading result, and would have led us to conclude erroneously that a single expansion had taken place in Europe. Only by combining the mitochondrial and microsatellite data sets are we able to reconstruct the species' history and show two colonization events in Europe, one before the Last Glacial Maximum (LGM) and one after it. Combining markers also revealed the importance of Asia Minor as an ancient refugium for this species and a source population for the expansion of the greater horseshoe bat into Europe before the LGM.     

Past tree range dynamics in the Iberian Peninsula inferred through phylogeography and palaeodistribution modelling: A review

The forests in the Iberian Peninsula have been strongly influenced by past climatic changes, but reconstructing their historical distributions and dynamics is very difficult due to the complex climatic characteristics and relief of the region. Research disciplines such as phylogeography and species distribution modelling can describe the past range dynamics of individual tree species in relatively great detail and help elucidate how these species have reacted to climatic changes. Here we review phylogeographical and modelling studies from species representative of the major Iberian forest types and attempt to extract general trends from the diversity of individual species histories in the Peninsula. To date most studies focus on geographical range dynamics during the Pleistocene, but an increasing body of evidence shows that some species have also retained the genetic imprints of much more ancient processes. Many widespread Iberian species show a deep and often remarkably clear-cut divide between populations from the Mediterranean and from the Atlantic regions of the Iberian Peninsula, suggesting that both areas have independently sustained viable populations over extended periods. In fact, phylogeographical studies commonly find that species had several glacial refugia across the Iberian Peninsula. On the other hand, distribution models help identifying further suitable areas that could have sustained so far undetected refugia. Such studies are of interest for species conservation, because refugium populations are high-priority targets due to their long-term persistence and unique evolutionary trajectory. Overall, we conclude that palaeoecology, phylogeography and species distribution modelling have a great potential to inform each other because of their complementary perspectives and results. A true integration of these approaches is therefore fundamental for further progress in our understanding of past Iberian environments and the organisms they harboured.     

The shaping of genetic variation in edge‐of‐range populations under past and future climate change

With rates of climate change exceeding the rate at which many species are able to shift their range or adapt, it is important to understand how future changes are likely to affect biodiversity at all levels of organisation. Understanding past responses and extent of niche conservatism in climatic tolerance can help predict future consequences. We use an integrated approach to determine the genetic consequences of past and future climate changes on a bat species, Plecotus austriacus. Glacial refugia predicted by palaeo‐modelling match those identified from analyses of extant genetic diversity and model‐based inference of demographic history. Former refugial populations currently contain disproportionately high genetic diversity, but niche conservatism, shifts in suitable areas and barriers to migration mean that these hotspots of genetic diversity are under threat from future climate change. Evidence of population decline despite recent northward migration highlights the need to conserve leading‐edge populations for spearheading future range shifts.     

An assessment of the effect of data partitioning on the performance of modelling algorithms for habitat suitability for ticks

A comparison of the performance of five modelling methods using presence/absence (generalized additive models, discriminant analysis) or presence-only (genetic algorithm for rule-set prediction, ecological niche factor analysis, Gower distance) data for modelling the distribution of the tick species Boophilus decoloratus (Koch, 1844) (Acarina: Ixodidae) at a continental scale (Africa) using climate data was conducted. This work explicitly addressed the usefulness of clustering using the normalized difference vegetation index (NDVI) to split original records and build partial models for each region (cluster) as a method of improving model performance. Models without clustering have a consistently lower performance (as measured by sensitivity and area under the curve [AUC]), although presence/absence models perform better than presence-only models. Two cluster-related variables, namely, prevalence (commonness of tick records in the cluster) and marginality (the relative position of the climate niche occupied by the tick in relation to that available in the cluster) greatly affect the performance of each model (P < 0.05). Both sensitivity and AUC are better for NDVI-derived clusters where the tick is more prevalent or its marginality is low. However, the total size of the cluster or its fragmentation (measured by Shannon's evenness index) did not affect the performance of models. Models derived separately for each cluster produced the best output but resulted in a patchy distribution of predicted occurrence. The use of such a method together with weighting procedures based on prevalence and marginality as derived from populations at each cluster produced a slightly lower predictive performance but a better estimation of the continental distribution of the tick. Therefore, cluster-derived models are able to effectively capture restricting conditions for different tick populations at a regional level. It is concluded that data partitioning is a powerful method with which to describe the climate niche of populations of a tick species, as adapted to local conditions. The use of this methodology greatly improves the performance of climate suitability models.     

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central–marginal hypothesis

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central–marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold‐adapted species likely used two isolated glacial refugia in southern France, in permafrost‐free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.     

Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia

Amborella trichopoda Baill. (Amborellaceae, Amborellales), the sole living member of the sister group to all other extant angiosperms, is endemic to New Caledonia. We addressed the intraspecific phylogeography of Amborella by investigating whether its present population genetic structure could be related to its current and past habitats. We found moderate range‐wide genetic diversity based on nuclear microsatellite data and detected four well‐differentiated, geographically distinct genetic groups using Bayesian clustering analyses. We modelled the ecological niche of Amborella based on the current climatic and environmental conditions. The predictive ability of the model was very good throughout the Central East mainland zone, but Amborella was predicted in the northern part of the island where this plant has not been reported. Furthermore, no significant barrier was detected based on habitat suitability that could explain the genetic differentiation across the area. Conversely, we found that the main genetic clusters could be related to the distribution of the suitable habitat at the last glacial maximum (LGM, c. 21 000 years BP), when Amborella experienced a dramatic 96.5% reduction in suitable area. At least two lineages survived in distinct putative refugia located in the Massif des Lèvres and in the vicinity of Mount Aoupinié. Our findings finally confirmed the importance of LGM rainforest refugia in shaping the current intra‐ and interspecific diversity in New Caledonian plants and revealed the possibility of an as yet unreported refugium. The combination of niche modelling and population genetics thereby offered novel insight into the biogeographical history of an emblematic taxon.     

Refugial races and postglacial colonization history of the freshwater amphipod Gammarus lacustris in Northern Europe

The systematic structure and postglacial population history of the freshwater amphipod Gammarus lacustris were explored in an allozyme survey of 65 populations across Northern Europe. A strong multilocus pattern of differentiation discriminated populations of the north‐east (north‐eastern Norway, northern Finland) from those in the west and the south (southern and central Scandinavia, Denmark, Poland). This principal division is attributed to postglacial colonization of the area by two main refugial races or lineages, one from the east (Russia), the other from the south (north‐western European continent). The strongly diverged Eastern and Western races (Nei's D= 0.3, from 22 loci) now meet in a secondary contact zone across a narrow sector of northernmost Norway. Genetic population compositions in this zone vary in a mosaic pattern, and show no evidence of reproductive incompatibility. Similar contacts of eastern and western lineages, far older than the latest glaciation, are now known from a number of taxa and they constitute a general pattern in Fennoscandian phylogeography. Within the Western Gammarus race, the populations through coastal north‐western Norway are further distinguished from those in southern Scandinavia and Denmark by a set of unique alleles at high frequencies (D = 0.12). This suggests an independent early colonization of the coastal region by another distinct stock, either along an early deglaciated coastal corridor from the south‐west, or directly from the ice‐free continental shelf off the Norwegian coast – a hypothesis that has also previously been presented for G. lacustris, and parallels controversial suggestions of local refugia for other taxa in Scandinavia. The coastal population type only later could come into contact with Gammarus invading over the mountains from the south; these two population types now smoothly intergrade. © 2003 The Linnean Society of London, Biological Journal of the Linnean Society, 2003, 79, 523–542.     

Analysis of molecular differentiation in a hybrid zone between chromosomally distinct races of the common shrew Sorex araneus (Insectivora: Soricidae) suggests their common ancestry

During postglacial colonization, populations that diverged in different refugia produced a patchwork of genomes, often delimited with sharp hybrid zones. The outcome of hybridization following the secondary contact of two genetically distinct populations is hard to predict. In this context, the present study investigated the genetic structure of the hybrid zone between the Drnholec and Białowieża chromosome races of the common shrew ( Sorex araneus ) in Poland using biparentally inherited (seven autosomal microsatellites) and uniparentally inherited (Y-linked microsatellite and mtDNA) molecular markers. On the basis of diagnostic chromosomes, the Drnholec and Białowieża races were classified to different karyotypic groups, which were believed to have independent glacial histories. It was found that genetic differentiation between the Drnholec and Białowieża races was weak and nonsignificant with respect to all molecular markers. However, these results are in contrast with the chromosomal structure of this hybrid zone. The very sharp frequency clines of the diagnostic chromosomes strongly suggest that gene flow between the Drnholec and Białowieża races was reduced. Nonsignificant correlations between genetic differentiation and both the presence of an environmental barrier and geographical distance reveal that only differences in karyotypes might be a reason for limited gene exchange between the races. It is assumed that a lack of molecular differences between the Drnholec and Białowieża races results from a shared ancestral variation.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 79–90.