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

Description of a New Species of the Pareas hamptoni Complex from Yunnan, China, with Confirmation of P. hamptoni Sensu Stricto in China (Squamata, Pareidae)

We describe a new species of the genus Pareas, based on three specimens collected from Guanyinshan Provincial Nature Reserve in Yuanyang County, Honghe Prefecture, Yunnan Province, China. The new species is distinguished from its congeners by one preocular, one postocular or postocular fused with subocular, loreal not bordering the orbit, one row enlarged vertebral scales, five rows keeled mid-dorsal scales at the middle of the body, 189-192 ventral scales and 72-89 subcaudal scales. The dorsal surfaces of the head and body are yellowish red or yellowish brown, and the belly and ventral surfaces of the head and tail are pinkish yellow or yellow with more or less small black spots. Phylogenetic analyses of mitochondrial DNA recovered the new species being the sister taxon to P. hamptoni sensu stricto. The genetic divergences between the new species and P. hamptoni sensu stricto were 4.2% in the Cyt b sequences and 5.0% in the ND4 sequences. In addition, based on specimens collected from Honghe and Wenshan prefectures, we confirmed that P. hamptoni sensu stricto is distributed in China.

A multidisciplinary approach to tackling invasive species: barcoding, morphology, and metataxonomy of the leafhopper Arboridia adanae

The leafhopper genus Arboridia includes several species that feed on Vitis vinifera and cause leaf chlorosis. We report the first alien Arboridia infestation in Italy in 2021 in an Apulian vineyard. To confirm the taxonomic status of the species responsible for crop damage, and reconstruct its demographic history, we barcoded individuals from Apulia together with Arboridia spp. from Crete (Greece), A. adanae from Central Turkey and other specimens of the presumed sister species, A. dalmatina from Dalmatia (Croatia). Molecular phylogenies and barcoding gap analysis identified clades not associated with sampling locations. This result is incongruent with classical specimen assignment and is further supported by morphological analyses, which did not reveal significant differences among the populations. Therefore, we propose A. dalmatina as a junior synonym of A. adanae, which would become the only grapevine-related Arboridia species in the eastern Mediterranean. To further characterise A. adanae evolution, we performed a molecular clock analysis that suggested a radiation during the Pleistocene glaciations. Finally, to assess whether the Apulian individuals carried microorganisms of agricultural relevance, we sequenced their bacterial microbiota using 16S rRNA amplicon sequencing identifying three phytopathogens not generally associated with Arboridia activities as well as Wolbachia in one Apulian haplogroup. We discuss the agricultural implications of this infestation.

Climbing on the La Canna Volcanic Sea Stack to Obtain First-Hand Data on the Tiniest Population of the Critically Endangered Aeolian Wall Lizard Podarcis raffonei

Among the extant populations of the critically endangered Aeolian wall lizard, the most vulnerable is the one surviving on La Canna, a columnar volcanic stack off the Filicudi Island. Here, I report the results of the first climbing expedition by a biologist on La Canna, that contributed direct observations and updated information on the size, morphology, and genetic variability of this population. Lizard density at the sampling site (a small terrace at 50 m of elevation) was 1.7 m^-2, twice of a previous estimate. Standard methods for estimating population size are unsuitable for La Canna. An educated guess of about a hundred individuals can be drawn, considering the extent of habitat available on the stack and the number of observed lizards. Lizards on La Canna were not fearless, despite what was reported by alpinists, possibly because of aggressive intraspecific interactions or high environmental temperatures during sampling. Biometric data significantly extend the body size of La Canna's lizards and indicate that it is not smaller than other P. raffonei populations. A complete lack of genetic diversity was found at the mitochondrial nd4 gene, in line with previous allozyme data and with estimates on other microinsular Podarcis populations. The small size of the La Canna population implies severe genetic drift and an extremely high level of inbreeding, as supported by low heterozygosity found across the genome. Detrimental effects of inbreeding depression are evident as cephalic malformations observed in all captured lizards of La Canna and might represent the more immediate threat to the persistence of this population.

A New Species of Pareas (Squamata, Pareidae) from Guangxi Province, China

We described a new species of genus Pareas from Baise City, Guangxi Zhuang Autonomous Region, China, based on morphological and molecular evidence. Pareas baiseensis sp. nov. is distinguished from its congeners by the combination of (1) Yellowish-brown body colouration; (2) Frontal subhexagonal to diamond-shaped with its lateral sides converging posteriorly; (3) The anterior pair of chin shields is longer than it is broad; (4) Loreal not in contact with the eye, prefrontal in contact with the eye, two or three suboculars; (5) Rows of 15-15-15 dorsal scales, five rows of mid-dorsal scales keeled at the middle of the body, one vertebral scale row enlarged; (6) 187-191 ventrals, 89-97 subcaudals, all divided, cloacal plate single; (7) Two postocular stripes, the nuchal area forming a dark black four-pointed fork collar with the middle tines shorter than the outside tines. The genetic divergence (uncorrected p-distance) between the new species and other representatives of Pareas ranged from 13.9% to 24.4% for Cytochrome b (Cyt b) and 12.1% to 25.5% for NADH dehydrogenase subunit 4 (ND4). Phylogenetic analyses of mitochondrial DNA gene data recovered the new species from being the sister taxon to (P. boulengeri + P. chinensis) from China.

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.

Shared haemogregarine infections in competing lacertids

In parasite–host interactions host species may differ in their ability to fight parasitic infections, while other ecological interactions, including competition, may differentially alter their physiological state, making them even more susceptible to parasites. In this study, we analyse the haemogregarine blood parasites infecting two competing lizard species, Iberolacerta horvathi and Podarcis muralis, and explore host–parasite relationships under different host competition scenarios. Both species were infected with haemogregarine parasites belonging to the genus Karyolysus. Using the 18S rRNA gene, six new Karyolysus haplotypes were identified clustering with other Central and Eastern European samples, and widely shared between both lizard hosts. Haemogregarine infections were detected at all sampled sites with over 50% of individuals parasitized. Overall, I. horvathi was more frequently and also more intensely parasitized than P. muralis, with higher infection rates observed in syntopy. Males of both species tended to be more frequently infected and showed a higher infection intensity than conspecific females. The results suggest that parasitisation by haemogregarines may be relevant in the dynamics of the competitive relationship between these lizard species. More studies, including immunological response analysis, and the identification of the vectors are needed to better understand host–parasite relationships and competition.

An integrative taxonomic revision of slug-eating snakes (Squamata: Pareidae: Pareineae) reveals unprecedented diversity in Indochina

Slug-eating snakes of the subfamily Pareinae are an insufficiently studied group of snakes specialized in feeding on terrestrial mollusks. Currently Pareinae encompass three genera with 34 species distributed across the Oriental biogeographic region. Despite the recent significant progress in understanding of Pareinae diversity, the subfamily remains taxonomically challenging. Here we present an updated phylogeny of the subfamily with a comprehensive taxon sampling including 30 currently recognized Pareinae species and several previously unknown candidate species and lineages. Phylogenetic analyses of mtDNA and nuDNA data supported the monophyly of the three genera Asthenodipsas, Aplopeltura, and Pareas. Within both Asthenodipsas and Pareas our analyses recovered deep differentiation with each genus being represented by two morphologically diagnosable clades, which we treat as subgenera. We further apply an integrative taxonomic approach, including analyses of molecular and morphological data, along with examination of available type materials, to address the longstanding taxonomic questions of the subgenus Pareas, and reveal the high level of hidden diversity of these snakes in Indochina. We restrict the distribution of P. carinatus to southern Southeast Asia, and recognize two subspecies within it, including one new subspecies proposed for the populations from Thailand and Myanmar. We further revalidate P. berdmorei, synonymize P. menglaensis with P. berdmorei, and recognize three subspecies within this taxon, including the new subspecies erected for the populations from Laos and Vietnam. Furthermore, we describe two new species of Pareas from Vietnam: one belonging to the P. carinatus group from southern Vietnam, and a new member of the P. nuchalis group from the central Vietnam. We provide new data on P. temporalis, and report on a significant range extension for P. nuchalis. Our phylogeny, along with molecular clock and ancestral area analyses, reveal a complex diversification pattern of Pareinae involving a high degree of sympatry of widespread and endemic species. Our analyses support the "upstream" colonization hypothesis and, thus, the Pareinae appears to have originated in Sundaland during the middle Eocene and then colonized mainland Asia in early Oligocene. Sundaland and Eastern Indochina appear to have played the key roles as the centers of Pareinae diversification. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns of Pareinae, with tectonic movements, orogeny and paleoclimatic shifts being the probable drivers of diversification. Our study brings the total number of Pareidae species to 41 and further highlights the importance of comprehensive taxonomic revisions not only for the better understanding of biodiversity and its evolution, but also for the elaboration of adequate conservation actions.

Biogeographic and demographic history of the Mediterranean snakes Malpolon monspessulanus and Hemorrhois hippocrepis across the Strait of Gibraltar

BACKGROUND

The contribution of North Africa to the assembly of biodiversity within the Western Palaearctic is still poorly documented. Since the Miocene, multiple biotic exchanges occurred across the Strait of Gibraltar, underlying the high biogeographic affinity between the western European and African sides of the Mediterranean basin. We investigated the biogeographic and demographic dynamics of two large Mediterranean-adapted snakes across the Strait and assess their relevance to the origin and diversity patterns of current European and North African populations.

RESULTS

We inferred phylogeographic patterns and demographic history of M. monspessulanus and H. hippocrepis, based on range-wide multilocus data, combined with fossil data and species distribution modelling, under present and past bioclimatic envelopes. For both species we identified endemic lineages in the High Atlas Mountains (Morocco) and in eastern Iberia, suggesting their persistence in Europe during the Pleistocene. One lineage is shared between North Africa and southern Iberia and likely spread from the former to the latter during the sea-level low stand of the last glacial stage. During this period M. monspessulanus shows a sudden demographic expansion, associated with increased habitat suitability in North Africa. Lower habitat suitability is predicted for both species during interglacial stages, with suitable areas restricted to coastal and mountain ranges of Iberia and Morocco. Compiled fossil data for M. monspessulanus show a continuous fossil record in Iberia at least since the Pliocene and throughout the Pleistocene.

CONCLUSIONS

The previously proposed hypothesis of Pleistocene glacial extinction of both species in Europe is not supported based on genetic data, bioclimatic envelopes models, and the available fossil record. A model of range retraction to mountain refugia during arid periods and of glacial expansion (demographic and spatial) associated to an increase of Mediterranean habitats during glacial epochs emerges as a general pattern for mesic vertebrates in North Africa. Moreover, the phylogeographic pattern of H. hippocrepis conforms to a well-established biogeographic partition between western and eastern Maghreb.

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.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae)

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Viability, behavior, and color expression in the offspring of matings between common wall lizard Podarcis muralis color morphs

Color polymorphisms are widely studied to identify the mechanisms responsible for the origin and maintenance of phenotypic variability in nature. Two of the mechanisms of balancing selection currently thought to explain the long-term persistence of polymorphisms are the evolution of alternative phenotypic optima through correlational selection on suites of traits including color and heterosis. Both of these mechanisms can generate differences in offspring viability and fitness arising from different morph combinations. Here, we examined the effect of parental morph combination on fertilization success, embryonic viability, newborn quality, antipredator, and foraging behavior, as well as inter-annual survival by conducting controlled matings in a polymorphic lacertid Podarcis muralis, where color morphs are frequently assumed to reflect alternative phenotypic optima (e.g., alternative reproductive strategies). Juveniles were kept in outdoor tubs for a year in order to study inter-annual growth, survival, and morph inheritance. In agreement with a previous genome-wide association analysis, morph frequencies in the year-old juveniles matched the frequencies expected if orange and yellow expressions depended on recessive homozygosity at 2 separate loci. Our findings also agree with previous literature reporting higher reproductive output of heavy females and the higher overall viability of heavy newborn lizards, but we found no evidence for the existence of alternative breeding investment strategies in female morphs, or morph-combination effects on offspring viability and behavior. We conclude that inter-morph breeding remains entirely viable and genetic incompatibilities are of little significance for the maintenance of discrete color morphs in P. muralis from the Pyrenees.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

No evidence for differential sociosexual behavior and space use in the color morphs of the European common wall lizard (Podarcis muralis)

Explaining the evolutionary origin and maintenance of color polymorphisms is a major challenge in evolutionary biology. Such polymorphisms are commonly thought to reflect the existence of alternative behavioral or life-history strategies under negative frequency-dependent selection. The European common wall lizard Podarcis muralis exhibits a striking ventral color polymorphism that has been intensely studied and is often assumed to reflect alternative reproductive strategies, similar to the iconic "rock-paper-scissors" system described in the North American lizard Uta stansburiana. However, available studies so far have ignored central aspects in the behavioral ecology of this species that are crucial to assess the existence of alternative reproductive strategies. Here, we try to fill this gap by studying the social behavior, space use, and reproductive performance of lizards showing different color morphs, both in a free-ranging population from the eastern Pyrenees and in ten experimental mesocosm enclosures. In the natural population, we found no differences between morphs in site fidelity, space use, or male-female spatial overlap. Likewise, color morph was irrelevant to sociosexual behavior, space use, and reproductive success within experimental enclosures. Our results contradict the commonly held hypothesis that P. muralis morphs reflect alternative behavioral strategies, and suggest that we should instead turn our attention to alternative functional explanations.

Spatial variation in gene flow across a hybrid zone reveals causes of reproductive isolation and asymmetric introgression in wall lizards

Hybrid zones provide insights into the evolution of reproductive isolation. Sexual selection can contribute to the evolution of reproductive barriers, but it remains poorly understood how sexual traits impact gene flow in secondary contact. Here, we show that a recently evolved suite of sexual traits that function in male-male competition mediates gene flow between two lineages of wall lizards (Podarcis muralis). Gene flow was relatively low and asymmetric in the presence of exaggerated male morphology and coloration compared to when the lineages share the ancestral phenotype. Putative barrier loci were enriched in genomic regions that were highly differentiated between the two lineages and showed low concordance between the transects. The exception was a consistently low genetic exchange around ATXN1, a gene that modulates social behavior. We suggest that this gene may contribute to the male mate preferences that are known to cause lineage-assortative mating in this species. Although female choice modulates the degree of reproductive isolation in a variety of taxa, wall lizards demonstrate that both male-male competition and male mate choice can contribute to the extent of gene flow between lineages.

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.

Genetic diversity and inferred dispersal history of the Schlegel's Japanese Gecko (Gekko japonicus) in Northeast Asia based on population genetic analyses and paleo-species distribution modelling

To understand the genetic diversity and dispersal history of Schlegel's Japanese gecko (Gekko japonicus), we performed genetic analyses and paleo-species distributional modelling. For the genetic analysis, we analyzed mitochondrial DNA (mtDNA) (cytochrome b [Cytb] and NADH dehydrogenase 2 [ND2]) and seven microsatellite loci of 353 individuals from 11 populations (2 east coast China, 4 west and central coast Japan and 5 Korea). For the paleo-species distribution modelling, we used 432 occurrence data points (125 China, 291 Japan and 16 Korea) over the Pleistocene and Holocene. China is inferred to be the source population, which had higher genetic diversity (mtDNA) and more private alleles (mtDNA) compared to Japanese and Korean populations. Differences between the three counties were very small in the mtDNA haplotype network despite some genetic structure among the three countries. Microsatellite analysis inferred that genetic exchange has actively occurred among the Chinese, Japanese and Korean populations. Suitable habitats in Japan should have been plentiful by the mid-Holocene, but have only recently become available in Korea. These results suggest that dispersal of G. japonicus occurred after the Holocene warming from the east coast of China to the west and central coasts of Japan and Korea, and gene flow is actively occurring among the three countries.

Genetic differentiation predicts body size divergence between island and mainland populations of common wall lizards (Podarcis muralis)

Small-bodied vertebrates sometimes evolve gigantism on islands, but there is a lack of consistent association with ecological factors or island characteristics. One possible reason is that, even if the ecological conditions are right, body size might fail to diverge on islands that were isolated recently or if there is gene flow between islands and the mainland. We studied body size, ventral colour polymorphism and genetic structure across nine islands and adjacent mainland populations of common wall lizards (Podarcis muralis) off the western coast of France. Population genetic data suggested that island populations might have maintained gene flow after their geographical isolation from the mainland. Island lizards were larger and heavier than mainland lizards on average, but the extent of gigantism varied substantially between islands. Island size and distance from the mainland were poor predictors of body size, but lizards from populations that were highly genetically differentiated from the mainland were larger than lizards from less differentiated populations. Colour morphs that were rare on the mainland tended to be more common on islands. We propose that genetic isolation or bottlenecks promote body size evolution in island lizards, which makes it challenging to identify ecological causes of island gigantism without complementary genetic information.

Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard (Podarcis muralis), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [sepiapterin reductase (SPR)] and carotenoid [beta-carotene oxygenase 2 (BCO2)] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard

Reptiles show an amazing color diversity based on variation in melanins, carotenoids, and pterins. This study reveals genes controlling differences between three color morphs (white, orange, and yellow) in the common wall lizard. Orange pigmentation, due to high levels of orange/red pterins in skin, is caused by genetic changes in the sepiapterin reductase gene. Yellow skin, showing high levels of yellow carotenoids, is controlled by the beta-carotene oxygenase 2 locus. Thus, the color polymorphism in the common wall lizard is associated with changes in two small regions of the genome containing genes with crucial roles in pterin and carotenoid metabolism. These genes are likely to have pleiotropic effects on behavior and other traits associated with the different color morphs.

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard (Podarcis muralis), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [sepiapterin reductase (SPR)] and carotenoid [beta-carotene oxygenase 2 (BCO2)] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

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.

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.

Extra-Mediterranean glacial refuges in barred and common grass snakes (Natrix helvetica, N. natrix)

Extra-Mediterranean glacial refugia of thermophilic biota, in particular in northern latitudes, are controversial. In the present study we provide genetic evidence for extra-Mediterranean refugia in two species of grass snake. The refuge of a widely distributed western European lineage of the barred grass snake (Natrix helvetica) was most likely located in southern France, outside the classical refuges in the southern European peninsulas. One genetic lineage of the common grass snake (N. natrix), distributed in Scandinavia, Central Europe and the Balkan Peninsula, had two distinct glacial refuges. We show that one was located in the southern Balkan Peninsula. However, Central Europe and Scandinavia were not colonized from there, but from a second refuge in Central Europe. This refuge was located in between the northern ice sheet and the Alpine glaciers of the last glaciation and most likely in a permafrost region. Another co-distributed genetic lineage of N. natrix, now massively hybridizing with the aforementioned lineage, survived the last glaciation in a structured refuge in the southern Balkan Peninsula, according to the idea of 'refugia-within-refugia'. It reached Central Europe only very recently. This study reports for the first time the glacial survival of a thermophilic egg-laying reptile species in Central Europe.

Unraveling climate influences on the distribution of the parapatric newts Lissotriton vulgaris meridionalis and L. italicus

Background

Climate is often considered as a key ecological factor limiting the capability of expansion of most species and the extent of suitable habitats. In this contribution, we implement Species Distribution Models (SDMs) to study two parapatric amphibians, Lissotriton vulgaris meridionalis and L. italicus, investigating if and how climate has influenced their present and past (Last Glacial Maximum and Holocene) distributions. A database of 901 GPS presence records was generated for the two newts. SDMs were built through Boosted Regression Trees and Maxent, using the Worldclim bioclimatic variables as predictors.

Results

Precipitation-linked variables and the temperature annual range strongly influence the current occurrence patterns of the two Lissotriton species analyzed. The two newts show opposite responses to the most contributing variables, such as BIO7 (temperature annual range), BIO12 (annual precipitation), BIO17 (precipitation of the driest quarter) and BIO19 (precipitation of the coldest quarter). The hypothesis of climate influencing the distributions of these species is also supported by the fact that the co-occurrences within the sympatric area fall in localities characterized by intermediate values of these predictors. Projections to the Last Glacial Maximum and Holocene scenarios provided a coherent representation of climate influences on the past distributions of the target species. Computation of pairwise variables interactions and the discriminant analysis allowed a deeper interpretation of SDMs’ outputs. Further, we propose a multivariate environmental dissimilarity index (MEDI), derived through a transformation of the multivariate environmental similarity surface (MESS), to deal with extrapolation-linked uncertainties in model projections to past climate. Finally, the niche equivalency and niche similarity tests confirmed the link between SDMs outputs and actual differences in the ecological niches of the two species.

Conclusions

The different responses of the two species to climatic factors have significantly contributed to shape their current distribution, through contractions, expansions and shifts over time, allowing to maintain two wide allopatric areas with an area of sympatry in Central Italy. Moreover, our SDMs hindcasting shows many concordances with previous phylogeographic studies carried out on the same species, thus corroborating the scenarios of potential distribution during the Last Glacial Maximum and the Holocene emerging from the models obtained.

Electronic supplementary material

The online version of this article (10.1186/s12983-017-0239-4) contains supplementary material, which is available to authorized users.

Inferring the demographic history of an oligophagous grasshopper: Effects of climatic niche stability and host-plant distribution

Understanding the consequences of past environmental changes on the abiotic and biotic components of the landscape and deciphering their impacts on the demographic trajectories of species is a major issue in evolutionary biogeography. In this study, we combine nuclear and mitochondrial genetic data to study the phylogeographical structure and lineage-specific demographic histories of the scrub-legume grasshopper (Chorthippus binotatus binotatus), a montane taxon distributed in the Iberian Peninsula and France that exclusively feeds on certain scrub-legume species. Genetic data and paleo-distribution modelling indicate the presence of four main lineages that seem to have diverged in allopatry and long-term persisted in Iberian and French refugia since the Mid Pleistocene. Comparisons of different demographic hypotheses in an Approximate Bayesian Computation (ABC) framework supported a population bottleneck in the northwestern French clade and paleo-distribution modelling indicate that the populations of this lineage have experienced more severe environmental fluctuations during the last 21 000 years than those from the Iberian Peninsula. Accordingly, we found that nuclear genetic diversity of the populations of scrub-legume grasshopper is positively associated with local stability of suitable habitats defined by both Pleistocene climate changes and historical distributional shifts of host-plant species. Overall, our study highlights the importance of integrating the potential effects of abiotic (i.e. climate and geography) and biotic components (i.e. inter-specific interactions) into the study of the evolutionary and demographic history of specialist taxa with narrow ecological requirements.

Phylogenetic relationships of Podarcis siculus (Rafinesque-Schmaltz, 1810) and Podarcis tauricus (Pallas, 1814) in Turkey, based on mitochondrial DNA

The Italian wall lizard and the Balkan wall lizard have a series of taxonomic revisions. However, their phylogenetic relationships still remain uncertain in Turkey. In the present study, we have assessed taxonomic relationships, both of Podarcis siculus and Podarcis tauricus through estimation of phylogenetic relationships among 43 and 42 specimens, respectively, using mtDNA (16 S rRNA and cytb) from great main populations in Turkey. The genetic distances among the populations of P. siculus in Turkey were very low and they were ranged from 0.2 to 1.6% in 16 S rRNA while they were ranged from 0.0% to 3.3% in cytb. On the other hand, the p-distances among the populations of P. tauricus were ranged from 0.0 to 0.6% in 16 S rRNA while they were 0.2% cytb in Turkey. Finally, most of the topologically identical trees of phylogenetic analyses and p-distances showed that monophyly was found in extant populations of P. siculus and P. tauricus. The nominate subspecies, P. s. siculus and P. t. tauricus are representatives of these lizards in Turkey.

Microevolution of the noble crayfish (Astacus astacus) in the Southern Balkan Peninsula

Background

The noble crayfish (Astacus astacus) displays a complex historical and contemporary genetic status in Europe. The species divergence has been shaped by geological events (i.e. Pleistocene glaciations) and humanly induced impacts (i.e. translocations, pollution, etc.) on its populations due to species commercial value and its niche degradation. Until now, limited genetic information has been procured for the Balkan area and especially for the southernmost distribution of this species (i.e. Greece). It is well known that the rich habitat diversity of the Balkan Peninsula offers suitable conditions for genetically diversified populations. Thus, the present manuscript revisits the phylogenetic relationships of the noble crayfish in Europe and identifies the genetic make-up and the biogeographical patterns of the species in its southern range limit.

Results

Mitochondrial markers (i.e. COI and 16S) were used in order to elucidate the genetic structure and diversity of the noble crayfish in Europe. Two of the six European haplotypic lineages, were found exclusively in Greece. These two lineages exhibited greater haplotypic richness when compared with the rest four (of “Central European” origin) while they showed high genetic diversity. Divergence time analysis identified that the majority of this divergence was captured through Pleistocene, suggesting a southern glacial refugium (Greece, southern Balkans). Furthermore, six microsatellite markers were used in order to define the factors affecting the genetic structure and demographic history of the species in Greece. The population structure analysis revealed six to nine genetic clusters and eight putative genetic barriers. Evidence of bottleneck effects in the last ~5000 years (due to climatic and geological events and human activities) is also afforded. Findings from several other research fields (e.g. life sciences, geology or even archaeology) have been utilized to perceive the genetic make-up of the noble crayfish.

Conclusions

The southernmost part of Balkans has played a major role as a glacial refugium for A. astacus. Such refugia have served as centres of expansion to northern regions. Recent history of the noble crayfish in southern Balkans reveals the influence of environmental (climate, geology and/or topology) and anthropogenic factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0971-6) contains supplementary material, which is available to authorized users.

Digging up the roots of an insular hotspot of genetic diversity: decoupled mito-nuclear histories in the evolution of the Corsican-Sardinian endemic lizard Podarcis tiliguerta

BACKGROUND

Mediterranean islands host a disproportionately high level of biodiversity and endemisms. Growing phylogeographic evidence on island endemics has unveiled unexpectedly complex patterns of intra-island diversification, which originated at diverse spatial and temporal scales. We investigated multilocus genetic variation of the Corsican-Sardinian endemic lizard Podarcis tiliguerta with the aim of shedding more light on the evolutionary processes underlying the origin of Mediterranean island biodiversity.

RESULTS

We analysed DNA sequences of mitochondrial (12S and nd4) and nuclear (acm4 and mc1r) gene fragments in 174 individuals of P. tiliguerta from 81 localities across the full range of the species in a geographic and genealogical framework. We found surprisingly high genetic diversity both at mitochondrial and nuclear loci. Seventeen reciprocally monophyletic allopatric mitochondrial haplogroups were sharply divided into four main mitochondrial lineages (two in Corsica and two in Sardinia) of Miocene origin. In contrast, shallow divergence and shared diversity within and between islands was observed at the nuclear loci. We evaluated alternative biogeographic and evolutionary scenarios to explain such profound discordance in spatial and phylogenetic patterning between mitochondrial and nuclear genomes. While neutral models provided unparsimonious explanations for the observed pattern, the hypothesis of environmental selection driving mitochondrial divergence in the presence of nuclear gene flow is favoured.

CONCLUSIONS

Our study on the genetic variation of P. tiliguerta reveals surprising levels of diversity underlining a complex phylogeographic pattern with a striking example of mito-nuclear discordance. These findings have profound implications, not only for the taxonomy and conservation of P. tiliguerta. Growing evidence on deep mitochondrial breaks in absence of geographic barriers and of climatic factors associated to genetic variation of Corsican-Sardinian endemics warrants additional investigation on the potential role of environmental selection driving the evolution of diversity hotspots within Mediterranean islands.

A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus

BACKGROUND

The current distribution of genetic diversity is the result of a vast array of microevolutionary processes, including short-term demographic and ecological mechanisms and long-term allopatric isolation in response to Quaternary climatic fluctuations. We investigated past processes that drove the population differentiation and spatial genetic distribution of the Italian wall lizard Podarcis siculus by means of sequences of mitochondrial cytb (n = 277 from 115 localities) and nuclear mc1r and β-fibint7genes (n = 262 and n = 91, respectively) from all its distribution range. The pattern emerging from the genetic data was compared with current and past (last glacial maximum) species distribution modeling (SDM).

RESULTS

We identified seven deeply divergent parapatric clades which presumably remained isolated in different refugia scattered mainly throughout the Tyrrhenian coast. Conversely, the Adriatic coast showed only two haplogroups with low genetic variability. These results appear to agree with the SDM prediction at the last glacial maximum (LGM) indicating a narrow area of habitat suitability along the Tyrrhenian coast and much lower suitability along the Adriatic one. However, the considerable land exposure of the Adriatic coastline favored a glacial colonization of the Balkan Peninsula.

CONCLUSIONS

Our population-level historical demography showed a common trend consistent with glacial expansions and regional persistence during the last glacial maximum. This complex genetic signature appears to be inconsistent with the expectation of the expansion-contraction model and post-LGM (re)colonizations from southern refugia. Hence it is one of an increasing number of cases in which these assumptions are not met, indicating that long-term fragmentation and pre-LGM events such as glacial persistence were more prominent in shaping genetic variation in this temperate species.

Inter- and intra-population variability of the protein content of femoral gland secretions from a lacertid lizard

Femoral glands of male lizards produce waxy secretions that are involved in inter- and intraspecific chemical communication. The main components of these secretions are proteins and lipids, the latter having been extensively studied and already associated to male quality. On the opposite, the composition and role of proteins are nearly unknown, the only available information coming from few studies on iguanids. These studies got the conclusion that proteins might have a communicative function, notably they could signal individual identity. A generalization of these findings requires the extension of protein analysis to other lizard families, and the primary detection of some patterns of individual variability. Using the common wall lizard Podarcis muralis as a model species, the protein fraction of the femoral pore secretions was investigated to provide the first characterization of this component in a lacertid lizard and to explore its source of variability, as a first step to support the hypothesized communicative role. Samples of proteins from femoral secretions were collected from 6 Italian populations and subjected to 1-dimensional electrophoresis. The binary vector of the band presence/absence was used to define the individual profiles. Protein fraction is found to have a structured pattern, with both an individual and a population component. Although the former supports the potential communicative role of proteins, the latter offers a double interpretation, phylogenetic or environmental, even though the phylogenetic effect seems more likely given the climatic resemblance of the considered sites. Further studies are necessary to shed light on both these issues.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0669-1) contains supplementary material, which is available to authorized users.