Rapid genetic assimilation of native wall lizard populations (Podarcis muralis) through extensive hybridization with introduced lineages

Unravelling the origin of the common wall lizards (Podarcismuralis) in south-eastern Europe using mitochondrial evidence

The origin of the common wall lizards (Podarcismuralis) populations in south-eastern Europe (namely in Bulgaria and Romania), representing the north-eastern range border of this species, was addressed using mitochondrial DNA. We compared cytochrome b sequences from Bulgaria and Romania with those from the contiguous range in Central Europe that are available from previous studies. We recorded five main haplogroups in Bulgaria and Romania, belonging to the Central Balkan clade. However, haplogroup III was recorded in more localities than previously found. Additionally, signs of haplotype admixture were identified in several populations along the Danube River. The presence of the Southern Alps haplotype in one population from Otopeni, Bucharest (Romania) and its close phylogenetic relationships to north Italy populations suggests human-mediated introductions of this wall lizard clade in Romania. Our results confirm that P.muralis can have non-native lineages and admixture through active human-mediated transport.

Population genomics of wall lizards reflects the dynamic history of the Mediterranean Basin

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an “Out-of-Italy” expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

The genic view of hybridization in the Anthropocene

Human impact is noticeable around the globe, indicating that a new era might have begun: the Anthropocene. Continuing human activities, including land-use changes, introduction of non-native species and rapid climate change, are altering the distributions of countless species, often giving rise to human-mediated hybridization events. While the interbreeding of different populations or species can have detrimental effects, such as genetic extinction, it can be beneficial in terms of adaptive introgression or an increase in genetic diversity. In this paper, I first review the different mechanisms and outcomes of anthropogenic hybridization based on literature from the last five years (2016-2020). The most common mechanisms leading to the interbreeding of previously isolated taxa include habitat change (51% of the studies) and introduction of non-native species (34% intentional and 19% unintentional). These human-induced hybridization events most often result in introgression (80%). The high incidence of genetic exchange between the hybridizing taxa indicates that the application of a genic view of speciation (and introgression) can provide crucial insights on how to address hybridization events in the Anthropocene. This perspective considers the genome as a dynamic collection of genetic loci with distinct evolutionary histories, giving rise to a heterogenous genomic landscape in terms of genetic differentiation and introgression. First, understanding this genomic landscape can lead to a better selection of diagnostic genetic markers to characterize hybrid populations. Second, describing how introgression patterns vary across the genome can help to predict the likelihood of negative processes, such as demographic and genetic swamping, as well as positive outcomes, such as adaptive introgression. It is especially important to not only quantify how much genetic material introgressed, but also what has been exchanged. Third, comparing introgression patterns in pre-Anthropocene hybridization events with current human-induced cases might provide novel insights into the likelihood of genetic swamping or species collapse during an anthropogenic hybridization event. However, this comparative approach remains to be tested before it can be applied in practice. Finally, the genic view of introgression can be combined with conservation genomic studies to determine the legal status of hybrids and take appropriate measures to manage anthropogenic hybridization events. The interplay between evolutionary and conservation genomics will result in the constant exchange of ideas between these fields which will not only improve our knowledge on the origin of species, but also how to conserve and protect them.

Genome-Wide SNPs Detect Hybridisation of Marsupial Gliders (Petaurus breviceps breviceps × Petaurus norfolcensis) in the Wild

Petaurus breviceps and Petaurus norfolcensis have produced hybrids in captivity, however there are no reported cases of Petaurus hybridisation in the wild. This study uses morphological data, mitochondrial DNA, and nuclear genome-wide SNP markers to confirm P. breviceps breviceps × P. norfolcensis hybridisation within their natural range on the central coast of New South Wales, Australia. Morphological data identified a potential hybrid that was confirmed with next-generation sequencing technology and 10,111 genome-wide SNPs. Both STRUCTURE and NewHybrid analyses identified the hybrid as a P. norfolcensis backcross, which suggests an initial F1 hybrid was fertile. The mitochondrial DNA matched that of a P. b. breviceps, indicating that a P. b. breviceps female initially mated with a P. norfolcensis male to produce a fertile female offspring. Our study is an important example of how genome-wide SNPs can be used to identify hybrids where the distribution of congeners overlaps. Hybridisation between congeners is likely to become more frequent as climate changes and habitats fragment, resulting in increased interactions and competition for resources and mates.

The possible origin of the common wall lizard, Podarcis muralis (Laurenti, 1768) in Ukraine

The phylogenetic relationships and possible origin of a putative non-native population of Podarcis muralis in Ukraine were assessed based on sequences of the mitochondrial gene cytochrome b. Phylogenetic analysis showed that the Ukrainian lizards belong to two distinct mitochondrial lineages (haplogroups), both occurring within the Central Balkan clade, which includes most of central and south-eastern European populations. From overall three detected Ukrainian haplotypes, one haplotype share same genetic signal with the hyplotype from the locality Bjala (Bulgaria), the other two are unique for Ukrainian population. Two of haplotypes correspond with haplogroup covering large geographic region of Bulgaria, Serbia, and Romania. These results reinforce previous findings that the species has the ability to establish new populations out of its native range. While most introductions to Germany and Britain have been deliberate, it appears likely that human transport of goods via the Danube river of goods is responsible for the range expansion into Ukraine.

A reference DNA barcode library for Austrian amphibians and reptiles

In the last few years, DNA barcoding became an established method for species identification in biodiversity inventories and monitoring studies. Such studies depend on the access to a comprehensive reference data base, covering all relevant taxa. Here we present a comprehensive DNA barcode inventory of all amphibian and reptile species native to Austria, except for the putatively extinct Vipera ursinii rakosiensis and Lissotriton helveticus, which has been only recently reported for the very western edge of Austria. A total of 194 DNA barcodes were generated in the framework of the Austrian Barcode of Life (ABOL) initiative. Species identification via DNA barcodes was successful for most species, except for the hybridogenetic species complex of water frogs (Pelophylax spp.) and the crested newts (Triturus spp.), in areas of sympatry. However, DNA barcoding also proved powerful in detecting deep conspecific lineages, e.g. within Natrix natrix or the wall lizard (Podarcis muralis), resulting in more than one Barcode Index Number (BIN) per species. Moreover, DNA barcodes revealed the presence of Natrix helvetica, which has been elevated to species level only recently, and genetic signatures of the Italian water frog Pelophylax bergeri in Western Austria for the first time. Comparison to previously published DNA barcoding data of European amphibians and reptiles corroborated the results of the Austrian data but also revealed certain peculiarities, underlining the particular strengths and in the case of the genus Pelophylax also the limitations of DNA barcoding. Consequently, DNA barcoding is not only powerful for species identification of all life stages of most Austrian amphibian and reptile species, but also for the detection of new species, the monitoring of gene flow or the presence of alien populations and/or species. Thus, DNA barcoding and the data generated in this study may serve both scientific and national or even transnational conservation purposes.

Hybridization and polyploidy in the weeping lizard Liolaemus chiliensis (Squamata: Liolaemidae)

In reptiles, polyploidy is an unusual phenomenon that can originate from interspecific hybridization. In Chile, the lizard Liolaemus chiliensis provides a unique model with which to study the origin of polyploidy because it occurs in populations of diploid individuals and a few populations that also contain triploid and diploid–triploid mosaic lizards. To test whether L. chiliensis polyploids were hybrids between lineages within the species, we analysed the karyotype, the mitochondrial gene Cytb, seven microsatellite loci, and the linear and geometric morphometry of individuals throughout the species distribution. We found polyploidy at different localities throughout the distribution. Hybrids were detected with the microsatellite loci and morphological data. A clear relationship between hybridization and polyploidy could not be established because triploid and mosaic lizards were hybrids and purebreds. However, our results support the occurrence of both phenomena, and further research is needed to clarify how hybridization could be involved in the origin of polyploidy.

Validating buccal swabbing as a minimal-invasive method to detect pesticide exposure in squamate reptiles

The use of enzymatic biomarkers constitutes a widely used approach in ecotoxicology. However, standard sampling procedures are invasive, requiring tissue, organ or blood extraction. This leads to concerns regarding conservation practice, animal welfare and restrictions in study design. New techniques are needed to avoid these problems, but still generate reliable data. Here, we aimed at validating the use of buccal swabs as a reliable method to detect pesticide exposure in reptiles. Common wall lizards (Podarcis muralis) were divided into control, dermal and oral treatment groups and exposed to different pesticide formulations. Subsequently, buccal swabs were taken and enzymatic activity was analyzed. We were able to confirm the suitability of the method to detect effects of pesticide exposure on the enzymatic level. While exposure to the formulation Roundup Ultramax® didn't match when compared to effects previously observed in situ when compared to other glyphosate based formulations, effects could still be detected. This can be seen as a strong indicator that the active ingredient of a formulation may not always be the mian driver for ecotoxicological effectsat the enzymatic level. At the same time, exposure towards the single formulation Vivando® didn't result in any effects. However, individuals residing in agricultural landscapes will mostly be exposed to pesticide mixes containing different formulations. Our results strongly advocate that buccal swabbing is a reliable minimal invasive method to generate samples for detecting effects of pesticide exposure in reptiles. Due to its easy handling, we believe it will provide new opportunities concerning study designs.

Phylogeographic evidence for multiple long-distance introductions of the common wall lizard associated with human trade and transport

The common wall lizard has been widely introduced across Europe and overseas. We investigated the origin of putatively introduced Podarcis muralis populations from two southern Europe localities: (i) Ljubljana (Slovenia), where uncommon phenotypes were observed near the railway tracks and (ii) the port of Vigo (Spain), where the species was recently found 150 km far from its previously known range. We compared cytochrome-b mtDNA sequences of lizards from these populations with published sequences across the native range. Our results support the allochthonous status and multiple, long-distance origins in both populations. In Ljubljana, results support two different origins, Serbia and Italy. In Vigo, at least two separate origins are inferred, from western and eastern France. Such results confirm that human-mediated transport is promoting biological invasion and lineage admixture in this species. Solid knowledge of the origin and invasion pathways, as well as population monitoring, is crucial for management strategies to be successful.

Herps without borders: a new newt case and a review of transalpine alien introductions in western Europe

Biogeographic processes have led to different evolutionary taxa occurring in the northern and southern edges of the Alpine Mountains in Western Europe. The integrity of this diversity is being challenged by frequent human-mediated trans-alpine translocations, sometimes leading to biological invasions. Several alien terrestrial vertebrates of south Alpine origins (Italy, Swiss Ticino) are causing damages to native north Alpine fauna. In this paper, we used molecular tools to characterize the understudied case of the Mediterranean smooth newt (Lissotriton vulgaris meridionalis) expanding in the outskirts of Geneva since its introduction before 1975. Mitochondrial and nuclear DNA sequencing suggest that these exotic populations are a mixture between two diverged L. v. meridionalis lineages from central Italy, and traces of potential hybridization with the native L. v. vulgaris was detected. This situation echoes many other trans-alpine alien introductions. We review all comparable cases of southern to northern Alps introductions in vertebrates, including seven reptiles and four amphibians. The majority of south alpine alien lineages were presumably imported voluntarily by enthusiasts and appear to perform better in the disturbed habitats found in the anthropogenic landscapes of Western Europe compared to their native north Alpine counterparts. Most pose serious threats to related species of similar ecology, through direct competition, predation and introgressive hybridization. Difficulties to detect alien species on time lead to significant conservation costs. Better education together with more appropriate and reactive management plans will be necessary to limit the impact of future alien introductions.

Admixture of hybrid swarms of native and introduced lizards in cities is determined by the cityscape structure and invasion history

Introductions of non-native lineages increase opportunities for hybridization. Non-native lineages of the common wall lizard, Podarcis muralis, are frequently introduced in cities where they hybridize with native populations. We aimed at unravelling the invasion history and admixture of native and non-native wall lizards in four German cities using citywide, comprehensive sampling. We barcoded and genotyped 826 lizards and tested if gene flow in populations composed of admixed native and introduced lineages is facilitated by similar environmental factors to those in native populations by comparing fine-scale landscape genetic patterns. In cities with non-native lineages, lizards commonly occurred in numerous clusters of hybrid swarms, which showed variable lineage composition, consisting of up to four distinct evolutionary lineages. Hybrid swarms held vast genetic diversity and showed recent admixture with other hybrid swarms. Landscape genetic analyses showed differential effects of cityscape structures across cities, but identified water bodies as strong barriers to gene flow in both native and admixed populations. By contrast, railway tracks facilitated gene flow of admixed populations only. Our study shows that cities represent unique settings for hybridization, caused by multiple introductions of non-native taxa. Cityscape structure and invasion histories of cities will determine future evolutionary pathways at these novel hybrid zones.

Tracing the maternal origin of the common wall lizard (Podarcis muralis) on the northern range margin in Central Europe

The maternal origin of isolated populations of the common wall lizard (Podracis muralis) in the Czech Republic, representing the north-eastern range border of the species, was addressed. We compared mitochondrial DNA sequences of the cytochrome b gene of samples from these populations with those from within the continuous range in Slovakia, the northern Balkan region, and those available from previous studies. We recorded five main haplogroups in the studied region, with all available Central European samples belonging to the same haplogroup. The star-like structure of this haplogroup suggests a scenario of relatively recent, post-glacial population expansion, which is further supported by a coalescent-based demographic analysis. The presence of unique haplotypes in two of the three isolated Czech populations together with close phylogenetic relationships to adjacent Slovak populations suggests either autochthonous origin or human-mediated introductions from geographically and genetically closest populations. We therefore support conservation programs for all three isolated Czech populations.

A simple SNP genotyping method reveals extreme invasions of non-native haplotypes in pale chub Opsariichthys platypus, a common cyprinid fish in Japan

Biological invasion by non-native subspecies or populations is one of the most serious threats to ecosystems, because these species might be easily established in the introduced area and can negatively affect native populations through competition and hybridization. Pale chub Opsariichthys platypus, one of the most common fish in East Asia, exhibits clear genetic differentiation among regional populations; however, introgression and subsequent loss of genetic integrity have been occurring throughout Japan due to the artificial introduction of non-native conspecifics. In this study, we developed a simple SNP genotyping method to discriminate between native and non-native mitochondrial DNA (mtDNA) haplotypes in pale chub using real-time PCR assay. We then investigated the distribution patterns of non-native pale chub in Tokai region, located in the center of Honshu Island, Japan and developed a predictive model of the occurrence of non-natives to reveal the factors influencing their invasion. The specificity and accuracy of the genotyping method were confirmed by using samples whose haplotypes were determined previously. Extensive occurrence of non-native haplotypes in Tokai region was detected by this method. In addition, our models suggested that the presence of non-natives varied greatly depending on the river system, and was positively influenced by the impounded water areas. Our method could accurately distinguish between native and non-native haplotypes of pale chub in Japan and suggested key environmental factors associated with the presence of non-natives. This approach can greatly reduce experimental costs be a great contribution for quantitative investigation.

Male behaviour drives assortative reproduction during the initial stage of secondary contact

Phenotypic divergence in allopatry can facilitate speciation by reducing the likelihood that individuals of different lineages hybridize during secondary contact. However, few studies have established the causes of reproductive isolation in the crucial early stages of secondary contact. Here, we establish behavioural causes of assortative reproduction between two phenotypically divergent lineages of the European wall lizard (Podarcis muralis), which have recently come into secondary contact. Parentage was highly assortative in experimental contact zones. However, despite pronounced divergence in male phenotypes, including chemical and visual sexual signals, there was no evidence that females discriminated between males of the two lineages in staged interactions or under naturalistic free-ranging conditions. Instead, assortative reproduction was driven by male mate preferences and, to a lesser extent, male-male competition. The effects were more pronounced when the habitat structure promoted high lizard densities. These results emphasize that assortative reproduction can occur in the absence of female choice and that male behaviour may play an important role in limiting hybridization during the initial stages of secondary contact.

Wide prevalence of hybridization in two sympatric grasshopper species may be shaped by their relative abundances

BACKGROUND

Hybridization between species is of conservation concern as it might threaten the genetic integrity of species. Anthropogenic factors can alter hybridization dynamics by introducing new potentially hybridizing species or by diminishing barriers to hybridization. This may even affect sympatric species pairs through environmental change, which so far has received little attention. We studied hybridization prevalence and the underlying behavioral mechanisms in two sympatric grasshopper species, a rare specialist (Chorthippus montanus) and a common generalist (Chorthippus parallelus). We conducted a mate choice experiment with constant intraspecific density and varying heterospecific density, i.e. varying relative frequency of both species.

RESULTS

Mate choice was frequency-dependent in both species with a higher risk of cross-mating with increasing heterospecific frequency, while conspecific mating increased linearly with increasing conspecific density. This illustrates that reproductive barriers could be altered by environmental change, if the relative frequency of species pairs is affected. Moreover, we performed a microsatellite analysis to detect hybridization in twelve syntopic populations (and four allotopic populations). Hybrids were detected in nearly all syntopic populations with hybridization rates reaching up to 8.9 %. Genetic diversity increased for both species when hybrids were included in the data set, but only in the common species a positive correlation between hybridization rate and genetic diversity was detected.

CONCLUSION

Our study illustrates that the relative frequency of the two species strongly determines the effectiveness of reproductive barriers and that even the more choosy species (Ch. montanus) may face a higher risk of hybridization if population size decreases and its relative frequency becomes low compared to its sister species. The asymmetric mate preferences of both species may lead to quasi-unidirectional gene flow caused by unidirectional backcrossing. This might explain why genetic diversity increased only in the common species, but not in the rare one. Altogether, the hybridization rate was much higher than expected for a widely sympatric species pair.

Widespread primary, but geographically restricted secondary, human introductions of wall lizards, Podarcis muralis

Establishing the introduction pathways of alien species is a fundamental task in invasion biology. The common wall lizard, Podarcis muralis, has been widely introduced outside of its native range in both Europe and North America, primarily through escaped pets or deliberate release of animals from captive or wild populations. Here, we use Bayesian clustering, approximate Bayesian computation (ABC) methods and network analyses to reconstruct the origin and colonization history of 23 non-native populations of wall lizards in England. Our analyses show that established populations in southern England originate from at least nine separate sources of animals from native populations in France and Italy. Secondary introductions from previously established non-native populations were supported for eleven (47%) populations. In contrast to the primary introductions, secondary introductions were highly restricted geographically and appear to have occurred within a limited time frame rather than being increasingly common. Together, these data suggest that extant wall lizard populations in England are the result of isolated accidental and deliberate releases of imported animals since the 1970s, with only local translocation of animals from established non-native populations. Given that populations introduced as recently as 25 years ago show evidence of having adapted to cool climate, discouraging further translocations may be important to prevent more extensive establishment on the south coast of England.

Historical Invasion Records Can Be Misleading: Genetic Evidence for Multiple Introductions of Invasive Raccoons (Procyon lotor) in Germany

Biological invasions provide excellent study systems to understand evolutionary, genetic and ecological processes during range expansions. There is strong evidence for positive effects of high propagule pressure and the associated higher genetic diversity on invasion success, but some species have become invasive despite small founder numbers. The raccoon (Procyon lotor) is often considered as a typical example for such a successful invasion resulting from a small number of founders. The species' largest non-native population in Germany is commonly assumed to stem from a small number of founders and two separate founding events in the 1930s and 1940s. In the present study we analyzed 407 raccoons at 20 microsatellite loci sampled from the invasive range in Western Europe to test if these assumptions are correct. Contrary to the expectations, different genetic clustering methods detected evidence for at least four independent introduction events that gave rise to genetically differentiated subpopulations. Further smaller clusters were either artifacts or resulted from founder events at the range margin and recent release of captive individuals. We also found genetic evidence for on-going introductions of individuals. Furthermore a novel randomization process was used to determine the potential range of founder population size that would suffice to capture all the alleles present in a cluster. Our results falsify the assumption that this species has become widespread and abundant despite being genetically depauperate and show that historical records of species introductions may be misleading.

The genetic integrity of the ex situ population of the European wildcat (Felis silvestris silvestris) is seriously threatened by introgression from domestic cats (Felis silvestris catus)

Studies on the genetic diversity and relatedness of zoo populations are crucial for implementing successful breeding programmes. The European wildcat, Felis s. silvestris, is subject to intensive conservation measures, including captive breeding and reintroduction. We here present the first systematic genetic analysis of the captive population of Felis s. silvestris in comparison with a natural wild population. We used microsatellites and mtDNA sequencing to assess genetic diversity, structure and integrity of the ex situ population. Our results show that the ex situ population of the European wildcat is highly structured and that it has a higher genetic diversity than the studied wild population. Some genetic clusters matched the breeding lines of certain zoos or groups of zoos that often exchanged individuals. Two mitochondrial haplotype groups were detected in the in situ populations, one of which was closely related to the most common haplotype found in domestic cats, suggesting past introgression in the wild. Although native haplotypes were also found in the captive population, the majority (68%) of captive individuals shared a common mtDNA haplotype with the domestic cat (Felis s. catus). Only six captive individuals (7.7%) were assigned as wildcats in the STRUCTURE analysis (at K = 2), two of which had domestic cat mtDNA haplotypes and only two captive individuals were assigned as purebred wildcats by NewHybrids. These results suggest that the high genetic diversity of the captive population has been caused by admixture with domestic cats. Therefore, the captive population cannot be recommended for further breeding and reintroduction.

Sexual selection against natural hybrids may contribute to reinforcement in a house mouse hybrid zone

Sexual selection may hinder gene flow across contact zones when hybrid recognition signals are discriminated against. We tested this hypothesis in a unimodal hybrid zone between Mus musculus musculus and Mus musculus domesticus where a pattern of reinforcement was described and lower hybrid fitness documented. We presented mice from the border of the hybrid zone with a choice between opposite sex urine from the same subspecies versus hybrids sampled in different locations across the zone. While no preference was evidenced in domesticus mice, musculus males discriminated in favour of musculus signals and against hybrid signals. Remarkably, the pattern of hybrid unattractiveness did not vary across the hybrid zone. Moreover, allopatric populations tested in the same conditions did not discriminate against hybrid signals, indicating character displacement for signal perception or preference. Finally, habituation-discrimination tests assessing similarities between signals pointed out that hybrid signals differed from the parental ones. Overall, our results suggest that perception of hybrids as unattractive has evolved in border populations of musculus after the secondary contact with domesticus. We discuss the mechanisms involved in hybrid unattractiveness, and the potential impact of asymmetric sexual selection on the hybrid zone dynamics and gene flow between the two subspecies.

Persistence across Pleistocene ice ages in Mediterranean and extra-Mediterranean refugia: phylogeographic insights from the common wall lizard

BACKGROUND

Pleistocene climatic oscillations have played a major role in structuring present-day biodiversity. The southern Mediterranean peninsulas have long been recognized as major glacial refugia, from where Northern Europe was post-glacially colonized. However, recent studies have unravelled numerous additional refugia also in northern regions. We investigated the phylogeographic pattern of the widespread Western Palaearctic lizard Podarcis muralis, using a range-wide multilocus approach, to evaluate whether it is concordant with a recent expansion from southern glacial refugia or alternatively from a combination of Mediterranean and northern refugia.

RESULTS

We analyzed DNA sequences of two mitochondrial (cytb and nd4) and three nuclear (acm4, mc1r, and pdc) gene fragments in individuals from 52 localities across the species range, using phylogenetic and phylogeographic methods. The complex phylogeographic pattern observed, with 23 reciprocally monophyletic allo- parapatric lineages having a Pleistocene divergence, suggests a scenario of long-term isolation in multiple ice-age refugia across the species distribution range. Multiple lineages were identified within the three Mediterranean peninsulas - Iberia, Italy and the Balkans - where the highest genetic diversity was observed. Such an unprecedented phylogeographic pattern - here called "refugia within all refugia" - compasses the classical scenario of multiple southern refugia. However, unlike the southern refugia model, various distinct lineages were also found in northern regions, suggesting that additional refugia in France, Northern Italy, Eastern Alps and Central Balkans allowed the long-term persistence of this species throughout Pleistocene glaciations.

CONCLUSIONS

The phylogeography of Podarcis muralis provides a paradigm of temperate species survival in Mediterranean and extra-Mediterranean glacial refugia. Such refugia acted as independent biogeographic compartments for the long-term persistence of this species, for the differentiation of its genetic lineages, and for the short-distance post-glacial re-colonization of neighbouring areas. This finding echoes previous findings from recent phylogeographic studies on species from temperate ecoregions, thus suggesting the need for a reappraisal of the role of northern refugia for glacial persistence and post-glacial assembly of Holarctic biota.