Evolutionary history of Ichthyosaura alpestris (Caudata, Salamandridae) inferred from the combined analysis of nuclear and mitochondrial markers

Incorporating microbiome analyses can enhance conservation of threatened species and ecosystem functions

Conservation genomics is a rapidly growing subdiscipline of conservation biology that uses genome-wide information to inform management of biodiversity at all levels. Such efforts typically focus on species or systems of conservation interest, but rarely consider associated microbes. At least three major approaches have been used to study how microorganisms broadly contribute to conservation areas: (1) diversity surveys map out microbial species distribution patterns in a variety of hosts, natural environments or regions; (2) functional surveys associate microbial communities with factors of interest, such as host health, symbiotic interactions, environmental characteristics, ecosystem processes, and biological invasions; and (3) manipulative experiments examine the response of changes to microbial communities or determine the functional roles of specific microbes within hosts or communities by adding, removing, or genetically modifying microbes. In practice, multiple approaches are often applied simultaneously. The results from all three conservation genomics approaches can be used to help design practical interventions and improve management actions, some of which we highlight below. However, experimental manipulations allow for more robust causal inferences and should be the ultimate goal of future work. Here we discuss how further integration of microbial research of a host's microbiome and of free living microbes into conservation biology will be an essential advancement for conservation of charismatic organisms and ecosystem functions in light of ongoing global environmental change.

Predicting the invasiveness of alpine newts in the UK

Predicting invasion risk to novel environments is essential for risk management and conservation decision making but the evolutionary lineage at which to make these predictions is often unclear. Here we predict the current suitability across the United Kingdom (UK) for the alpine newt Ichthyosaura alpestris, a species with a complex evolutionary history, a broad native range, a growing number of introduced populations and anecdotal reports of ecological consequences to native amphibian communities. We use species distribution and ecological niche modelling to predict environmental suitability of the alpine newt in the UK at both the species-level and lineage-level and to quantify evolutionary lineage niche overlap. We show good model transferability at the species-level and parts of the UK-especially central and eastern England and parts of central and northern Scotland-to be highly environmentally suitable for the alpine newt. Our findings provide evidence of environmental niche differences at the lineage-level, with the Greek lineage being distinct from most other lineages, but with low confidence in maxent predictions for the Greek, Balkan and Italian lineages due to high levels of extrapolation. In contrast, the niche of the UK records appear to share the same niche as the Central lineage. We find 66% of currently known alpine newt records to fall within areas predicted to be environmentally suitable at the species-level, providing a series of testable hypotheses to better understand the invasion ecology of this species in the UK.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10530-025-03543-2.

Phylogenomics yields new systematic and taxonomical insights for Lissotriton newts, a genus with a strong legacy of introgressive hybridization

The ease with which genome-wide data can nowadays be collected allows complicated phylogenetic questions to be re-evaluated. Phylogenetic relationships among newts have often proven difficult to resolve due to the prevalence of incomplete lineage sorting and introgressive hybridization. For the newt genus Lissotriton, phylogenetic relationships are not settled and there is controversy surrounding the species status of several taxa. We obtain c. 7 k nuclear DNA markers with target enrichment by sequence capture and conduct a concatenated analysis with RAxML, gene-tree summarization with ASTRAL, and species tree estimation with SNAPPER. We explore introgression between evolutionary lineages with TreeMix and Dsuite and compare how introgression events influence the different phylogenetic tools employed. We retrieve tree topologies that are discordant with previous mtDNA-based attempts, particularly concerning the phylogenetic placement of L. italicus and the L. vulgaris species complex. Yet, we also observe deviations between the phylogenetic hypotheses resulting from the different analyses. We interpret the placement of L. montandoni deep within the L. vulgaris species complex by SNAPPER, rather than as the sister taxon to the remainder of the L. vulgaris species complex according to RAxML and ASTRAL, as an artifact of introgression - well-documented in previous work and backed up by TreeMix and Dsuite analyses. Our analyses allow us to make some taxonomical recommendations: we confirm the recently proposed species status of L. kosswigi and L. graecus and propose that L. v. lantzi and L. v. schmidtleri had better be treated as subspecies. Our work also highlights areas for further taxonomic research: range-wide phylogenomic data are required to disentangle the L. boscai - L. maltzani species complex and the northern and southern lineages of L. v. vulgaris. Our study illustrates the power of target enrichment by sequence capture in tackling longstanding questions in taxa with an extensive history of hybridization and introgression.

Combining Historical and Molecular Data to Study Nearly Extinct Native Italian Grey Partridge (Perdix perdix) at the Turn of the Twentieth Century

The grey partridge (Perdix perdix Linnaeus, 1758), is a polytypic species with seven recognized subspecies, including P. p. italica (Hartert, 1917), which is endemic to Italy. Until World War II, the species was widespread across Europe but severely declined due to anthropogenic causes, jeopardizing the Italian subspecies gene pool. Genetic characterization and haplotype identification were performed by analyzing the 5'-end of the mitochondrial control region (CR). A total of 15 haplotypes were detected, seven of which were present in the population before 1915. Among them, three haplotypes were never detected again in the individuals collected after 1915. Interestingly, eight of the 15 haplotypes detected in Italian museum samples belonged exclusively to individuals collected after 1915. The obtained data highlight a high presence of specimens originating from other European populations and, despite all the conservation efforts, suggest an uncertain situation of the subspecies in Italy. This research was strongly backed up by extensive bibliographic research on historical documents, allowing the identification of hundreds of restocking events all over Italy. This is an integral part of this research and has laid the foundations for identifying and circumscribing historical periods in which introductions from the rest of Europe had different pressures, aiming to define a baseline.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Survived the Glaciations, Will They Survive the Fish? Allochthonous Ichthyofauna and Alpine Endemic Newts: A Road Map for a Conservation Strategy

The Calabrian Alpine newt (Ichthyosaura alpestris inexpectata) is a glacial relict with small and extremely localised populations in the Catena Costiera (Calabria, Southern Italy) and is considered to be "Endangered" by the Italian IUCN assessment. Climate-induced habitat loss and recent fish introductions in three lakes of the Special Area of Conservation (SAC) Laghi di Fagnano threaten the subspecies' survival in the core of its restricted range. Considering these challenges, understanding the distribution and abundance of this newt is crucial. We surveyed the spatially clustered wetlands in the SAC and neighbouring areas. First, we provide the updated distribution of this subspecies, highlighting fish-invaded and fishless sites historically known to host Calabrian Alpine newt populations and two new breeding sites that have been recently colonised. Then, we provide a rough estimate of the abundance, body size and body condition of breeding adults and habitat characteristics in fish-invaded and fishless ponds. We did not detect Calabrian Alpine newts at two historically known sites now invaded by fish. Our results indicate a reduction in occupied sites and small-size populations. These observations highlight the need for future strategies, such as fish removal, the creation of alternative breeding habitats and captive breeding, to preserve this endemic taxon.

Environmental and topographic drivers of amphibian phylogenetic diversity and endemism in the Iberian Peninsula

Understanding the ecological and evolutionary processes driving biodiversity patterns and allowing their persistence is of utmost importance. Many hypotheses have been proposed to explain spatial diversity patterns, including water-energy availability, habitat heterogeneity, and historical climatic refugia. The main goal of this study is to identify if general spatial drivers of species diversity patterns of phylogenetic diversity (PD) and phylogenetic endemism (PE) at the global scale are also predictive of PD and PE at regional scales, using Iberian amphibians as a case study. Our main hypothesis assumes that topography along with contemporary and historical climate are drivers of phylogenetic diversity and endemism, but that the strength of these predictors may be weaker at the regional scale than it tends to be at the global scale. We mapped spatial patterns of Iberian amphibians' phylogenetic diversity and endemism, using previously published phylogenetic and distribution data. Furthermore, we compiled spatial data on topographic and climatic variables related to the water-energy availability, topography, and historical climatic instability hypotheses. To test our hypotheses, we used Spatial Autoregressive Models and selected the best model to explain diversity patterns based on Akaike Information Criterion. Our results show that, out of the variables tested in our study, water-energy availability and historical climate instability are the most important drivers of amphibian diversity in Iberia. However, as predicted, the strength of these predictors in our case study is weaker than it tends to be at global scales. Thus, additional drivers should also be investigated and we suggest caution when interpreting these predictors as surrogates for different components of diversity.

Strong restructuration of skin microbiota during captivity challenges ex-situ conservation of amphibians

In response to the current worldwide amphibian extinction crisis, conservation instances have encouraged the establishment of ex-situ collections for endangered species. The resulting assurance populations are managed under strict biosecure protocols, often involving artificial cycles of temperature and humidity to induce active and overwintering phases, which likely affect the bacterial symbionts living on the amphibian skin. However, the skin microbiota is an important first line of defense against pathogens that can cause amphibian declines, such as the chytrid Batrachochytrium dendrobatidis (Bd). Determining whether current husbandry practices for assurance populations might deplete amphibians from their symbionts is therefore essential to conservation success. Here, we characterize the effect of the transitions from the wild to captivity, and between aquatic and overwintering phases, on the skin microbiota of two newt species. While our results confirm differential selectivity of skin microbiota between species, they underscore that captivity and phase-shifts similarly affect their community structure. More specifically, the translocation ex-situ is associated with rapid impoverishment, decrease in alpha diversity and strong species turnover of bacterial communities. Shifts between active and overwintering phases also cause changes in the diversity and composition of the microbiota, and on the prevalence of Bd-inhibitory phylotypes. Altogether, our results suggest that current husbandry practices strongly restructure the amphibian skin microbiota. Although it remains to be determined whether these changes are reversible or have deleterious effects on their hosts, we discuss methods to limit microbial diversity loss ex-situ and emphasize the importance of integrating bacterial communities to applied amphibian conservation.

Next-generation phylogeography of the banded newts (Ommatotriton): A phylogenetic hypothesis for three ancient species with geographically restricted interspecific gene flow and deep intraspecific genetic structure

Technological developments now make it possible to employ many markers for many individuals in a phylogeographic setting, even for taxa with large and complex genomes such as salamanders. The banded newt (genus Ommatotriton) from the Near East has been proposed to contain three species (O. nesterovi, O. ophryticus and O. vittatus) with unclear phylogenetic relationships, apparently limited interspecific gene flow and deep intraspecific geographic mtDNA structure. We use parallel tagged amplicon sequencing to obtain 177 nuclear DNA markers for 35 banded newts sampled throughout the range. We determine population structure (with Bayesian clustering and principal component analysis), interspecific gene flow (by determining the distribution of species-diagnostic alleles) and phylogenetic relationships (by maximum likelihood inference of concatenated sequence data and based on a summary-coalescent approach). We confirm that the three proposed species are genetically distinct. A sister relationship between O. nesterovi and O. ophryticus is suggested. We find evidence for introgression between O. nesterovi and O. ophryticus, but this is geographically limited. Intraspecific structuring is extensive, with the only recognized banded newt subspecies, O. vittatus cilicensis, representing the most distinct lineage below the species level. While mtDNA mostly mirrors the pattern observed in nuclear DNA, all banded newt species show mito-nuclear discordance as well.

Unraveling the Hidden Diversity of the Native White Claw Crayfish in the Iberian Peninsula

Several European freshwater crayfish species are currently included in one of the IUCN Red list categories. In the Iberian Peninsula, the native Austropotamobius pallipes species complex (the white clawed crayfish, WCC) has experienced a drastic decline since 1973. Implementing conservation management strategies for this species requires a better understanding of the patterns and structure of its genetic diversity. In this study, we assessed the levels and patterns of genetic variation in 71 populations along the whole distributional range of the WCC in the Iberian Peninsula. The two mitochondrial markers analyzed (Cytochrome Oxidase subunit I and 16S rRNA genes) indicated high levels of genetic diversity, which are significantly geographically structured in three main genetic groups, two corresponding to Northern and one to Central-Eastern and the westernmost Iberian Peninsula. The diversity found included new private haplotypes, and revealed the potential effect of paleogeographic barriers and last glaciations in the population structure observed. Current conservation and management programs for the WCC in the Iberian Peninsula should consider these three phylogeographic groups as essential management units in order to preserve the remaining genetic diversity in the species.

Batrachochytrium salamandrivorans Threat to the Iberian Urodele Hotspot

The recent introduction of the chytrid fungus Batrachochytrium salamandrivorans into northeastern Spain threatens salamander diversity on the Iberian Peninsula. We assessed the current epidemiological situation with extensive field sampling of urodele populations. We then sought to delineate priority regions and identify conservation units for the Iberian Peninsula by estimating the susceptibility of Iberian urodeles using laboratory experiments, evidence from mortality events in nature and captivity and inference from phylogeny. None of the 1395 field samples, collected between 2015 and 2021 were positive for Bsal and no Bsal-associated mortality events were recorded, in contrast to the confirmed occurrence of Bsal outbreak previously described in 2018. We classified five of eleven Iberian urodele species as highly susceptible, predicting elevated mortality and population declines following potential Bsal emergence in the wild, five species as intermediately susceptible with variable disease outcomes and one species as resistant to disease and mortality. We identified the six conservation units (i.e., species or lineages within species) at highest risk and propose priority areas for active disease surveillance and field biosecurity measures. The magnitude of the disease threat identified here emphasizes the need for region-tailored disease abatement plans that couple active disease surveillance to rapid and drastic actions.

Mitochondrial DNA diversity of the alpine newt (Ichthyosaura alpestris) in the Carpathian Basin: evidence for multiple cryptic lineages associated with Pleistocene refugia

The phylogeography and molecular taxonomy of the Alpine newt, Ichthyosaura alpestris, has been intensively studied in the past. However, previous studies did not include a comprehensive sampling from the Carpathian Basin, possibly a key region in the evolution of the species. We used a 1251 bp long fragment of the mitochondrial genome to infer the species’ evolutionary history in central-eastern Europe by assigning isolated Carpathian Basin populations from 6 regions to previously defined mtDNA lineages. We also revised the morphology-based intraspecific taxonomy of the species in the light of new genetic data. Alpine newt populations from the Carpathian Basin represented two different mitochondrial lineages. The Mátra, Bükk and Zemplén Mts populations can be assigned to the Western lineage of the nominotypical subspecies. Bakony and Őrség populations showed high haplotype diversity and formed a separate clade within the Western lineage, suggesting that the Carpathian Basin might have provided cryptic refugia for Alpine newt populations in their cold-continental forest-steppe landscapes during the younger Pleistocene. Newts from Apuseni Mts were related to the Eastern lineage but formed a distinct clade within this lineage. Considering the morphological and genetic differentiation of the Bakony and Őrség populations, consistent with a long independent evolutionary history, we propose these populations be referred to as Ichthyosaura alpestris bakonyiensis (Dely, 1964). We provide a redescription of this poorly known subspecies.

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.

Slow and steady wins the race: contrasted phylogeographic signatures in two Alpine amphibians

A deeper phylogeographic structure is expected for slow-dispersing habitat specialists compared to widespread adaptable species, especially in topographically complex regions. We tested this classic assumption by comparing the genomic (RAD-sequencing) phylogeographies of two amphibians inhabiting the Swiss Alps: the mobile, cosmopolitan common frog (Rana temporaria) against the stationary, mountain endemic Alpine salamander (Salamandra atra). Our results ran opposite of predictions: the frog displayed significantly higher genetic divergences and lower within-population variation compared to the salamander. This implies a prominent role for their distinctive glacial histories in shaping intraspecific diversity and structure: diversification and recolonization from several circum-Alpine micro-refugia for the frog versus a single refugium for the salamander, potentially combined with better population connectivity and stability. These striking differences emphasize the great variability of phylogeographic responses to the Quaternary glaciations, hence the complexity to predict general patterns of genetic diversity at the regional scale, and the forces that underlie them.

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.

Speciation patterns in the Forficula auricularia species complex: cryptic and not so cryptic taxa across the western Palaearctic region

Forficula auricularia (the European earwig) is possibly a complex of cryptic species. To test this hypothesis, we performed: (1) a phylogeographic study based on fragments of the mitochondrial COI and the nuclear ITS2 markers on a wide geographic sampling, (2) morphometric analyses of lineages present in Spain and (3) niche overlap analyses. We recovered five reciprocally monophyletic ancient phylogroups with unique historical patterns of distribution, climatic niches and diversification. External morphology was conserved and not correlated with speciation events, except in one case. Phylogenetic placement of the morphologically distinct taxon renders F. auricularia paraphyletic. Based on the congruence of the phylogenetic units defined by mtDNA and nuclear sequence data, we conclude that phylogroups have their own historical and future evolutionary trajectory and represent independent taxonomic units. Forficula auricularia is a complex of at least four species: the morphologically diagnosable Forficula aeolica González-Miguéns & García-París sp. nov., and the cryptic taxa: Forficula mediterranea González-Miguéns & García-París sp. nov., Forficula dentataFabricius, 1775stat. nov. and Forficula auriculariaLinnaeus, 1758s.s. We also provide new synonymy for F. dentata.

Limited long-distance dispersal success in a Western European fairy shrimp evidenced by nuclear and mitochondrial lineage structuring

Anostraca are known by their ability for long-distance dispersal, but the existence in several species of deep, geographically structured mtDNA lineages suggests their populations are subjected to allopatric differentiation, isolation, and prevalence of local scale dispersion. Tanymastix stagnalis is one of the most widespread species of Anostraca and previous studies revealed an unclear geographical pattern of mtDNA genetic diversity. Here, we analyze populations from the Iberian and Italian Peninsulas, Central Europe, and Scandinavia, with the aim to characterize the patterns of genetic diversity in a spatio-temporal framework using mtDNA and nuclear markers to test gene flow among close populations. For these aims we built a time-calibrated phylogeny and carried out Bayesian phylogeographic analyses using a continuous diffusion model. Our results indicated that T. stagnalis presents a deeply structured genetic diversity, including 7 ancient lineages, some of them even predating the Pleistocene. The Iberian Peninsula harbors high diversity of lineages, with strong isolation and recent absence of gene flow between populations. Dispersal at local scale seems to be the prevailing dispersal mode of T. stagnalis, which exhibits a pattern of isolation-by-distance in the Iberian Peninsula. We remark the vulnerability of most of these lineages, given the limited known geographic distribution of some of them, and the high risk of losing important evolutionary potential for the species.

Extensive cytonuclear discordance in a crested newt from the Balkan Peninsula glacial refugium

Integration of multilocus data and species distribution modelling into phylogeography allows mitochondrial DNA (mtDNA)-based scenarios to be fine-tuned. We address the question of whether extensive mtDNA substructuring in the crested newt Triturus macedonicus from the Balkan Peninsula is matched in the nuclear genome. We determine the intraspecific population structure based on 52 nuclear DNA markers and project a species distribution model on climate layers for the Last Glacial Maximum. We show that T. macedonicus accumulated nuclear DNA population structure in an area predicted to have been climatically stable during the Pleistocene, with four nuclear DNA groups in the western part of the species range. The distribution of these nuclear DNA groups shows little agreement with that of mtDNA structuring, which shows three highly distinct species-specific clades and a fourth one introgressed from another crested newt species. This cytonuclear discordance conveys that historical biogeographical scenarios based on mtDNA exclusively should be interpreted with caution. Our findings further highlight the important role the Balkan Peninsula has played in the evolution and preservation of European biodiversity.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Genes on the edge: A framework to detect genetic diversity imperiled by climate change

Ongoing global warming is disrupting several ecological and evolutionary processes, spanning different levels of biological organization. Species are expected to shift their ranges as a response to climate change, with relevant implications to peripheral populations at the trailing and leading edges. Several studies have analyzed the exposure of species to climate change but few have explored exposure at the intraspecific level. We introduce a framework to forecast exposure to climate change at the intraspecific level. We build on existing methods by combining correlative species distribution models, a model of species range dynamics, and a model of phylogeographic interpolation. We demonstrate the framework by applying it to 20 Iberian amphibian and reptile species. Our aims were to: (a) identify which species and intraspecific lineages will be most exposed to future climate change; (b) test if nucleotide diversity at the edges of species ranges are significantly higher or lower than on the overall range; and (c) analyze if areas of higher species gain, loss, and turnover coincide with those predicted for lineages richness and nucleotide diversity. We found that about 80% of the studied species are predicted to contract their range. Within each species, some lineages were predicted to contract their range, while others were predicted to maintain or expand it. Therefore, estimating the impacts of climate change at the species level only can underestimate losses at the intraspecific level. Some species had significant high amount of nucleotide at the trailing or leading edge, or both, but we did not find a consistent pattern across species. Spatial patterns of species richness, gain, loss, and turnover were fairly concurrent with lineages richness and nucleotide diversity. Our results support the need for increased attention to intraspecific diversity regarding monitoring and conservation strategies under climate change.

Geographic patterns of phenotypic diversity in incipient species of North American blister beetles (Coleoptera: Meloidae) are not determined by species niches, but driven by demography along the speciation process

The Epicauta stigmata complex is a group of blister beetles composed of three parapatric or sympatric species that occur in central Mexico to southern USA: E. stigmata, E. uniforma and E. melanochroa. These species are morphologically very similar, and are mainly distinguished by body colour differences. Here we assessed whether phenotypic divergence in coloration patterns define evolutionary units within the complex. We studied the phylogenetic relationships, demographic history and concordances between morphological and ecological traits in the group. The complex apparently had a demographic history of recent population expansion during the last glaciation period 75000 to 9500 years ago. The three species show no reciprocal monophyly, and thus their allospecificity was not confirmed. The current distribution of haplotypes and the genetic divergences in these taxa can be explained by either recent mitochondrial introgression events caused by hybridisation or by incomplete lineage sorting. Colour pattern differences in the complex are not likely a product of local selection acting over a common genetic background. We suggest that phenotypic divergence in colour patterns during an incipient speciation process might be seen as an enhancing factor of cohesion within each of the three evolutionary units.

Quaternary history, population genetic structure and diversity of the cold-adapted Alpine newt Ichthyosaura alpestris in peninsular Italy

Mediterranean peninsulas are major biodiversity hotspots, and cold-adapted species are an important component of this biodiversity. However, cold-adapted species contributed surprisingly little to our knowledge of the intimate links between Quaternary environmental changes, species' responses to these changes, and current patterns of intraspecific biodiversity. Here, we investigated the genetic structure and evolutionary history of a cold-adapted amphibian, the Alpine newt Ichthyosaura alpestris, within the Italian peninsula. Nuclear and mitochondrial markers consistently identified three distinct genetic lineages, whose divergence dates to the Early Pleistocene (1.9 and 0.8 million years ago). Our results show that the Italian peninsula provided multiple Pleistocene refugia to this cold-adapted species, and suggest that allopatric fragmentation followed by secondary admixture have been key events in the formation of its current pattern of genetic diversity. Indeed, estimates of population genetic diversity clearly identified contact populations as those achieving the highest levels of diversity. Such concordance among cold-adapted and temperate species in terms of processes triggering the formation of regional patterns of genetic diversity provides strong support for the hypothesis that gene exchange between divergent lineages, rather than long-term stability of refugial populations, has been the main step toward the formation of hotspots of intraspecific biodiversity.

Newt life after fish introduction: extirpation of paedomorphosis in a mountain fish lake and newt use of satellite pools

Fish introduction is one of the main causes of amphibian decline worldwide. It affects particularly rare aquatic phenotypes such as paedomorphs, which retain gills during the adult stage. In this context, we determined whether small wetlands, such as pools surrounding fished and fishless lakes, could sustain paedomorphic and metamorphic newts. To this end, we surveyed lakes known historically to sustain Alpine newts Ichthyosaura alpestris as well as 35 nearby pools. On the basis of the published records, the only known population exhibiting paedomorphosis in the Swiss Alps was found to be extirpated by salmonid introductions. However, the metamorphs persisted in peripheral pools, paedomorphosis was discovered at a new locality, and overwintering larvae were still present in one of the lakes. These results show the importance of conserving varied aquatic habitats such as pools in mountainous environments where the main resources can become unsuitable for amphibians because of fish introductions. Pools may also function as reservoirs in maintaining newt populations until programs to remove fish from lakes can be carried out. It is not known if paedomorphs could reappear after fish removal. However, the combined resilience of amphibians after fish removal and the genetic basis for paedomorphosis highlighted in other taxa by previous studies suggest that there is the potential to maintain this intraspecific case of diversity even after its disappearance.