Multilocus phylogeography of the common midwife toad, Alytes obstetricans (Anura, Alytidae): Contrasting patterns of lineage diversification and genetic structure in the Iberian refugium.

Recent investigations on the evolutionary history of the common midwife toad (Alytes obstetricans) revealed high levels of geographically structured genetic diversity but also a situation where delineation of major historical lineages and resolution of their relationships are much more complex than previously thought. We studied sequence variation in one mitochondrial and four nuclear genes throughout the entire distribution range of all recognized A. obstetricans subspecies to infer the evolutionary processes that shaped current patterns of genetic diversity and population subdivision. We found six divergent, geographically structured mtDNA haplogroups diagnosing population lineages, and varying levels of admixture in nuclear markers. Given the timeframe inferred for the splits between major lineages, the climatic and environmental changes that occurred during the Pleistocene seem to have shaped the diversification history of A. obstetricans. Survival of populations in allopatric refugia through the Ice Ages supports the generality of the "refugia-within-refugia" scenario for the Iberian Peninsula. However, lineages corresponding to subspecies A. o. almogavarii, A. o. pertinax, A. o. obstetricans, and A. o. boscai responded differently to Pleistocene climatic oscillations after diverging from a common ancestor. Alytes o. obstetricans expanded northward from a northern Iberian refugium through the western Pyrenees, leaving a signal of contrasting patterns of genetic diversity, with a single mtDNA haplotype north of the Pyrenees from SW France to Germany. Both A. o. pertinax and A. o. boscai are widespread and genetically diverse in Iberia, the latter comprising two divergent lineages with a long independent history. Finally, A. o. almogavarii is mostly restricted to the north-eastern corner of Iberia north of the Ebro river, with additional populations in a small region in south-eastern France. This taxon exhibits unparalleled levels of genetic diversity and little haplotype sharing with other lineages, suggesting a process of incipient speciation.

Copy number polymorphism in the α-globin gene cluster of European rabbit (Oryctolagus cuniculus).

Comparative genomic studies have revealed that mammals typically possess two or more tandemly duplicated copies of the α-globin (HBA) gene. The domestic rabbit represents an exception to this general rule, as this species was found to possess a single HBA gene. Previous electrophoretic surveys of HBA polymorphism in natural populations of the European rabbit (Oryctolagus cuniculus) revealed extensive geographic variation in the frequencies of three main electromorphs. The variation in frequency of two electromorphs is mainly partitioned between two distinct subspecies of European rabbit, and a third is restricted to the hybrid zone between the two rabbit subspecies in Iberia. Here we report the results of a survey of nucleotide polymorphism, which revealed HBA copy number polymorphism in Iberian populations of the European rabbit. By characterizing patterns of HBA polymorphism in populations from the native range of the European rabbit, we were able to identify the specific amino-acid substitutions that distinguish the previously characterized electromorphs. Within the hybrid zone, we observed the existence of a second HBA gene duplicate, named HBA2, that mostly represents a novel sequence haplotype, which occurs in higher frequency within the hybrid zone, and thus appears to have arisen in hybrids of the two distinct subspecies. Although this novel gene is also present in other wild Iberian populations, it is almost absent from French populations, which suggest a recent ancestry, associated with the establishment of the post-Pleistocene contact zone between the two European rabbit subspecies.

Permanent genetic resources added to Molecular Ecology Resources Database 1 December 2010-31 January 2011.

This article documents the addition of 238 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Alytes dickhilleni, Arapaima gigas, Austropotamobius italicus, Blumeria graminis f. sp. tritici, Cobitis lutheri, Dendroctonus ponderosae, Glossina morsitans morsitans, Haplophilus subterraneus, Kirengeshoma palmata, Lysimachia japonica, Macrolophus pygmaeus, Microtus cabrerae, Mytilus galloprovincialis, Pallisentis (Neosentis) celatus, Pulmonaria officinalis, Salminus franciscanus, Thais chocolata and Zootoca vivipara. These loci were cross-tested on the following species: Acanthina monodon, Alytes cisternasii, Alytes maurus, Alytes muletensis, Alytes obstetricans almogavarii, Alytes obstetricans boscai, Alytes obstetricans obstetricans, Alytes obstetricans pertinax, Cambarellus montezumae, Cambarellus zempoalensis, Chorus giganteus, Cobitis tetralineata, Glossina fuscipes fuscipes, Glossina pallidipes, Lysimachia japonica var. japonica, Lysimachia japonica var. minutissima, Orconectes virilis, Pacifastacus leniusculus, Procambarus clarkii, Salminus brasiliensis and Salminus hilarii.

High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent, multiple refugia: evidence from mtDNA, microsatellites and nuclear genealogies.

Pleistocene glaciations often resulted in differentiation of taxa in southern European peninsulas, producing the high levels of endemism characteristic of these regions (e.g. the Iberian Peninsula). Despite their small ranges, endemic species often exhibit high levels of intraspecific differentiation as a result of a complex evolutionary history dominated by successive cycles of fragmentation, expansion and subsequent admixture of populations. Most evidence so far has come from the study of species with an Atlantic distribution in northwestern Iberia, and taxa restricted to Mediterranean-type habitats remain poorly studied. The Iberian Midwife toad (Alytes cisternasii) is a morphologically conserved species endemic to southwestern and central Iberia and a typical inhabitant of Mediterranean habitats. Applying highly variable genetic markers from both mitochondrial and nuclear genomes to samples collected across the species' range, we found evidence of high population subdivision within A. cisternasii. Mitochondrial haplotypes and microsatellites show geographically concordant patterns of genetic diversity, suggesting population fragmentation into several refugia during Pleistocene glaciations followed by subsequent events of geographical and demographic expansions with secondary contact. In addition, the absence of variation at the nuclear beta-fibint7 and Ppp3caint4 gene fragments suggests that populations of A. cisternasii have been recurrently affected by episodes of extinction and recolonization, and that documented patterns of population subdivision are the outcome of recent and multiple refugia. We discuss the evolutionary history of the species with particular interest in the increasing relevance of Mediterranean refugia for the survival of genetically differentiated populations during the Pleistocene glaciations as revealed by studies in co-distributed taxa.

The genomic legacy from the extinct Lepus timidus to the three hare species of Iberia: contrast between mtDNA, sex chromosomes and autosomes.

Extensive interspecific genetic introgression is often reported, and appraising its genomic impact can serve to determine whether it results from selection on specific loci or from demographic processes affecting the whole genome. The three species of hares present in the Iberian Peninsula harbour high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species now extinct in the region. This could result from the invasive replacement of L. timidus by the temperate species during deglaciation but should then have left traces in the nuclear genome. We typed single nucleotide polymorphisms (SNPs) discovered by sequencing 10 autosomal loci, two X-linked and one Y-linked in species-wide samples of the four taxa. Based on lineage-diagnostic SNPs, we detected no trace of L. timidus sex chromosomes in Iberia. From the frequencies of inferred haplotypes, autosomal introgression into L. granatensis appeared mostly sporadic but always widespread instead of restricted to the north as mtDNA. Autosomal introgression into Iberian L. europaeus, inhabiting the Pyrenean foothills, was hardly detectable, despite quasi-fixation of L. timidus mtDNA. L. castroviejoi, endemic to the Cantabrian Mountains and fixed for L. timidus mtDNA, showed little traces of autosomal introgression. The absence of sex-chromosome introgression presumably resulted from X-linked hybrid male unfitness. The contrasting patterns between the autosomes and mtDNA could reflect general gender asymmetric processes such as frequency-dependent female assortative mating, lower mtDNA migration and higher male dispersal, but adaptive mtDNA introgression cannot be dismissed. Additionally, we document reciprocal introgression between L. europaeus and both L. granatensis in Iberia and L. timidus outside Iberia.

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula.

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.

The limits of mtDNA phylogeography: complex patterns of population history in a highly structured Iberian lizard are only revealed by the use of nuclear markers.

Phylogeographic analyses based on the sole use of the mitochondrial DNA (mtDNA) molecule reveal only a small part of the evolutionary history of a species or a set of related species. In this study, we have combined the application of slow- and fast-evolving nuclear markers (proteins and microsatellites, respectively) together with the analysis of two-gene genealogies to further understand the history of the Iberian endemic Schreiber's green lizard, Lacerta schreiberi, a species for which a well established phylogeographical scenario is available. In sharp contrast with the observation of four divergent and almost allopatric mtDNA clades, our nuclear data revealed how two groups of populations diverged, persisted and began to admix along the mountains of the Iberian Central System. In addition, the combination of mtDNA and nuclear data showed how the core area of the species distribution responded to ice ages, first by relatively old processes of population expansion to the south followed by episodes of contraction that are at the origin of present-day isolates, and more recently by a postglacial expansion to the Iberian Northwest where new habitats were made available after climatic amelioration. Taken together with recently published results for a variety of other organisms, our results suggest that complex processes of fragmentation, expansion and admixture can only be properly addressed through the use of several and complementary types of molecular markers. Finally, we also suggest that southern European refugia are both hotspots and melting pots of genetic diversity.

Evidence for contrasting modes of selection at interacting globin genes in the European rabbit (Oryctolagus cuniculus).

In hybrid zones between genetically differentiated populations, variation in locus-specific rates of introgression may reflect adaptation to different environments or adaptation to different genetic backgrounds. The European rabbit, Oryctolagus cuniculus, is well-suited to studies of such hybrid zone dynamics because it is composed of two genetically divergent subspecies that hybridize in a zone of secondary contact in central Iberia. A species-wide survey of allozyme variation revealed a broad range of locus-specific divergence levels (F(ST) ranged from 0 to 0.54, mean F(ST)=0.16). Interestingly, the two loci that fell at opposite ends of the distribution of F(ST) values, haemoglobin alpha-chain (HBA) and haemoglobin beta-chain (HBB), encode interacting subunits of the haemoglobin protein. The contrasting patterns of spatial variation at these two loci could not be reconciled under a neutral model of population structure. The HBA gene exhibited higher-than-expected levels of population differentiation, consistent with a history of spatially varying selection. The HBB gene exhibited lower-than-expected levels of population differentiation, consistent with some form of spatially uniform selection. Patterns of linkage disequilibrium and allele frequency variation do not appear to fit any simple model of two-locus epistatic selection.

Participatory dam systems modelling: a case study of the transboundary Guadiana River in the Iberian Peninsula.

Modelling tools have been widely used to investigate best management practices. But in contrast to the plethora of modelling studies, the practical implementation of outcomes is comparatively small. There is an urgent need to implement results and to show the practical validation of the concept developed, especially against the context of water stress mitigation. The participative development of modelling studies as a joint effort of stakeholders and modellers is seen as a key to achieve a wider identification, acceptance, trust and applicability of results. Participatory planning in the water sector is also increasingly requested in water management, where tasks have been for clarified decades through different institutional arrangements and national laws. Stakeholder involvement in water resources management have been limited to what was long time seen as participation, merely information on action to be taken. In the last decade the need for participation has been reflected in different ways. In Europe, the implementation of the Water Framework Directive (WFD) under Art. 14, requires all the European countries to involve stakeholders in decision making processes on water resource management. The aim of this paper is to present and discuss the research framework and possible results of investigating dam modelling through participatory systems modelling. We introduce a structured approach to use participatory modelling (PM) for stimulating the integration of modelling and decision making, also as a way of implementation of some articles of the WFD. The results and the framing of this paper are part of the AQUASTRESS Project. The conceptual modelling has been developed by a multidisciplinary research team, local stakeholders and local experts. Some results are discussed and recommendations made.

Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation.

Pleistocene climatic oscillations were a major force shaping genetic variability in many taxa. We analyse the relative effects of the ice ages across a latitudinal gradient in the Western Mediterranean region, testing two main predictions: (i) species with historical distributions in northern latitudes should have experienced greater loss of suitable habitat, resulting in higher extinction of historical lineages than species distributed in southern latitudes, where the effects of the ice ages were not as drastic. This would be reflected in the observation of lower diversity and number of differentiated lineages in northern areas. (ii) a signature of demographic expansion following the climate amelioration should be obvious in northern species, whereas in the south evidence of long-term effective population size stability should be observed. We used as models three species of wall lizards (Podarcis bocagei, Podarcis carbonelli and Podarcis vaucheri) that replace each other along the study area. We investigated the patterns of mitochondrial DNA diversity and subdivision and obtained demographic parameter estimates for each species. Our results suggest that P. bocagei, the northernmost species, bears low genetic diversity, a shallow coalescent history and marks of a demographic expansion. In contrast, P. vaucheri, the species with a southernmost distribution, shows deeper coalescence events, complex geographical substructure and no evidence for population growth. The species with an intermediate distribution, P. carbonelli, shows average levels of diversity, substructure and population growth. Taken together, these results conform to our main predictions and are explained by a differential influence of the ice ages on distinct latitudes.

The rise and fall of the mountain hare (Lepus timidus) during Pleistocene glaciations: expansion and retreat with hybridization in the Iberian Peninsula.

The climatic fluctuations during glaciations have affected differently arctic and temperate species. In the northern hemisphere, cooling periods induced the expansion of many arctic species to the south, while temperate species were forced to retract in southern refugia. Consequently, in some areas the alternation of these species set the conditions for competition and eventually hybridization. Hares in the Iberian Peninsula appear to illustrate this phenomenon. Populations of Iberian hare (Lepus granatensis), brown hare (Lepus europaeus) and broom hare (Lepus castroviejoi) in Northern Iberia harbour mitochondrial haplotypes from the mountain hare (Lepus timidus), a mainly boreal and arctic species presently absent from the peninsula. To understand the history of this past introgression we analysed sequence variation and geographical distribution of mitochondrial control region and cytochrome b haplotypes of L. timidus origin found in 378 specimens of these four species. Among 124 L. timidus from the Northern Palaearctic and the Alps we found substantial nucleotide diversity (2.3%) but little differentiation between populations. Based on the mismatch distribution of the L. timidus sequences, this could result from an expansion at a time of temperature decrease favourable to this arctic species. The nucleotide diversity of L. timidus mtDNA found in Iberian L. granatensis, L. europaeus and L. castroviejoi (183, 70 and 1 specimens, respectively) was of the same order as that in L. timidus over its range (1.9%), suggesting repeated introgression of multiple lineages. The structure of the coalescent of L. granatensis sequences indicates that hybridization with L. timidus was followed by expansion of the introgressed haplotypes, as expected during a replacement with competition, and occurred when temperatures started to rise, favouring the temperate species. Whether a similar scenario explains the introgression into Iberian L. europaeus remains unclear but it is possible that it hybridized with already introgressed L. granatensis.

Genetic variation at chemokine receptor CCR5 in leporids: alteration at the 2nd extracellular domain by gene conversion with CCR2 in Oryctolagus, but not in Sylvilagus and Lepus species.

Whereas in its natural host (Sylvilagus sps.) the effects of myxoma virus infections are benign, in European rabbit (Oryctolagus cuniculus), it causes a highly infectious disease with very high mortality rate, known as myxomatosis. There is evidence that, as with HIV-1 virus in human, myxoma virus may use chemokine receptors such as CCR5 of the host target cell for entry and activation of pathways of immune avoidance. We have characterized and compared CCR5 genes of leporid species with different susceptibility levels to myxomatosis. The CCR5 protein of O. cuniculus differs markedly from all those known from other species. The most striking was the replacement of a specific peptide motif of the second extracellular loop (ECL2) by a motif, which in other species characterizes the CCR2 molecules. While absent in Sylvilagus and Lepus species, this CCR2 imposed CCR5-ECL2 alteration was observed in all genomes of 25 European rabbits, representing the subspecies O. cuniculus algirus and O. cuniculus cuniculus. Allelic variation at the rabbit CCR5 locus confirmed that the gene conversion predates the subspecies split (1-2 Ma).

Genealogy of the nuclear beta-fibrinogen locus in a highly structured lizard species: comparison with mtDNA and evidence for intragenic recombination in the hybrid zone.

The study of nuclear genealogies in natural populations of nonmodel organisms is expected to provide novel insights into the evolutionary history of populations, especially when developed in the framework of well-established mtDNA phylogeographical scenarios. In the Iberian Peninsula, the endemic Schreiber's green lizard Lacerta schreiberi exhibits two highly divergent and allopatric mtDNA lineages that started to split during the late Pliocene. In this work, we performed a fine-scale analysis of the putative mtDNA contact zone together with a global analysis of the patterns of variation observed at the nuclear beta-fibrinogen intron 7 (beta-fibint7). Using a combination of DNA sequencing with single-strand conformational polymorphism (SSCP) analysis, we show that the observed genealogy at the beta-fibint7 locus reveals extensive admixture between two formerly isolated lizard populations while the two mtDNA lineages remain essentially allopatric. In addition, a private beta-fibint7 haplotype detected in the single population where both mtDNA lineages were found in sympatry is probably the result of intragenic recombination between the two more common and divergent beta-fibint7 haplotypes. Our results suggest that the progressive incorporation of nuclear genealogies in investigating the ancient demography and admixture dynamics of divergent genomes will be necessary to obtain a more comprehensive picture of the evolutionary history of organisms.

The evolution of the immunoglobulin heavy chain variable region (IgVH) in Leporids: an unusual case of transspecies polymorphism.

In domestic rabbit (Oryctolagus cuniculus), three serological types have been distinguished at the variable domain of the antibody H chain, the so-called V(H) a allotypes a1, a2, and a3. They correspond to highly divergent allelic lineages of the V(H) 1 gene, which is the gene rabbit utilizes in more than 80% of VDJ rearrangements. The sharing of serological V(H) a markers between rabbit and snowshoe hare (Lepus americanus) has suggested that the large genetic distances between rabbit V(H) 1 alleles (9-14% nucleotide differences) can be explained by unusually long lineage persistence times (transspecies polymorphism). Because this interpretation of the serological data is uncertain, we have determined the nucleotide sequences of V(H) genes expressed in specimens of Lepus species. Two sequence groups were distinguished, one of which occurred only in hare specimen displaying serological motifs of the rabbit V(H) a-a2 allotype. Sequences of this group are part of a monophyletic cluster containing the V(H) 1 sequences of the rabbit a2 allotype. The fact that this "transspecies a2 cluster" did not include genes of other rabbit V(H) a allotypes (a1, a3, and a4) is incompatible with the existence of a common V(H) a ancestor gene within the species, and suggests that the divergence of the V(H) a lineages preceded the Lepus vs Oryctolagus split. The sequence data are furthermore compatible with the hypothesis that the V(H)a polymorphism can be two times older than the divergence time between the Lepus and Oryctolagus lineages, which was estimated at 16-24 million years.

Epileptiform activity triggers long-term plasticity of GABA(B) receptor signalling in the developing rat hippocampus.

GABA(B) receptor (GABA(B)R)-mediated presynaptic inhibition regulates neurotransmitter release from synaptic terminals. In the neonatal hippocampus, GABA(B)R activation reduces GABA release and terminates spontaneous network discharges called giant depolarizing potentials (GDPs). Blocking GABA(B)Rs transforms GDPs into longer epileptiform discharges. Thus, GABA(B)R-mediated presynaptic inhibition of GABA release (GABA auto-inhibition) controls both spontaneous network activity and excitability in the developing hippocampus. Here we show that extensive release of endogenous GABA during epileptiform activity impairs GABA auto-inhibition, but not GABA(B)R-mediated inhibition of glutamate release, leading to hyperexcitability of the neonatal hippocampal network. Paired-pulse depression of GABA release (PPD) and heterosynaptic depression of glutamate release were used to monitor the efficacy of presynaptic GABA(B)R-mediated inhibition in slices. PPD, but not heterosynaptic depression, was dramatically reduced after potassium (K+)-induced ictal-like discharges (ILDs), suggesting a selective impairment of GABA(B)R-dependent presynaptic inhibition of GABAergic terminals. Impairing GABA auto-inhibition induced a 44% increase in GDP width and the appearance of pathological network discharges. Preventing GABA-induced activation of GABA(B)Rs during ILDs avoided PPD loss and most modifications of the network activity. In contrast, a partial block of GABA(B)Rs induced network discharges strikingly similar to those observed after K+-driven ILDs. Finally, neither loss of GABA auto-inhibition nor network hyperexcitability could be observed following synchronous release of endogenous GABA in physiological conditions (during GDPs at 1 Hz). Thus, epileptiform activity was instrumental to impair GABA(B)R-dependent presynaptic inhibition of GABAergic terminals. In conclusion, our results indicate that endogenous GABA released during epileptiform activity can reduce GABA auto-inhibition and trigger pathological network discharges in the newborn rat hippocampus. Such functional impairment may play a role in acute post-seizure plasticity.

High levels of nucleotide diversity in the European rabbit (Oryctolagus cuniculus) SRY gene.

We have sequenced 2,388 bp of the European rabbit sex determining region Y (SRY) gene. These data provide a 10-fold increase in the coverage of the Y chromosome in this species, including the entire open reading frame of the SRY, the polyadenylation signal, and two repetitive sequences in the 5' -region. A survey of 2021 bp of this gene in eight domestic breeds and four wild individuals revealed a total of nine single nucleotide polymorphisms and one indel, defining two deeply divergent lineages. The resulting estimation of nucleotide diversity (pi=1.34 x10(-3)) is very high when compared with other species, but no variability was detected among the domestic breeds. This study represents a first step in the characterization of the European rabbit Y chromosome and its variability. These sequences can be used in additional phylogeographical analyses of the European rabbit and other Leporid species, as well as in evolutionary studies of sex determination and the Y chromosome in wild species.

Invasion from the cold past: extensive introgression of mountain hare (Lepus timidus) mitochondrial DNA into three other hare species in northern Iberia.

Mitochondrial DNA introgression from Lepus timidus into Lepus granatensis and Lepus europaeus was recently reported in Iberia, although L. timidus presumably retreated from this region at the end of the last ice age. Here we assess the extent of this ancient mtDNA introgression by RFLP analysis of 695 specimens representing the three hare species present in Iberia. The introgressed L. timidus lineage was found in 23 of the 37 populations sampled. It is almost fixed in L. europaeus across its Iberian range in the Pyrenean foothills, and in L. granatensis, which occupies the rest of the peninsula, it is predominant in the north and gradually disappears further south. We also found it in Lepus castroviejoi, a species endemic to Cantabria. Multiple hybridizations and, potentially, a selective advantage for the L. timidus lineage can explain the remarkable taxonomic and geographical range of this mitochondrial introgression.

Genetic exchange across a hybrid zone within the Iberian endemic golden-striped salamander, Chioglossa lusitanica.

The study of hybrid zones resulting from Pleistocene vicariance is central in examining the potential of genetically diverged evolutionary units either to introgress and merge or to proceed with further isolation. The hybrid zone between two mitochondrial lineages of Chioglossa lusitanica is located near the Mondego River in Central Portugal. We used mitochondrial and nuclear diagnostic markers to conduct a formal statistical analysis of the Chioglossa hybrid zone in the context of tension zone theory. Key results are: (i) cline centres are not coincident for all markers, with average widths of ca. 2-15 km; (ii) heterozygote deficit was not observed across loci near the transect centre; (iii) associations of parental allele combinations ('linkage disequilibrium'R) were not detected either across loci or across the transect. These observations suggest that the Chioglossa hybrid zone is not a tension zone with strong selection against hybrids but instead one shaped mostly by neutral mixing. The patterns uncovered suggest a complex history of populations over a small scale that may be common in southern Pleistocene refugia.