Integrative inference of population history in the Ibero-Maghrebian endemic Pleurodeles waltl (Salamandridae)

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Phylogeographic and Paleoclimatic Modelling Tools Improve Our Understanding of the Biogeographic History of Hierophis viridiflavus (Colubridae)

Phylogeographic and paleoclimatic modelling studies have been combined to infer the role of Pleistocene climatic oscillations as drivers of the genetic structure and distribution of Mediterranean taxa. For the European whip snake, Hierophis viridiflavus, previous studies based on paleoclimatic modelling have depicted a low reliability in the pattern of past climatic suitability across the central Mediterranean Basin, which barely fits the species' genetic structure. In this study, we combined phylogeographic and paleoclimatic modelling tools to improve our understanding of the biogeographic history of H. viridiflavus, particularly extending the sampling and phylogeographic inferences to previously under-sampled regions. Phylogeographic analyses recovered two major clades that diverged at the beginning of the Pleistocene and had diversified in different ways by the late Pleistocene: the east clade (composed of three subclades) and the west clade (with no further structure). Paleoclimatic models highlighted the temperate character of H. viridiflavus, indicating range contractions during both the last inter-glacial and last glacial maximum periods. Range expansions from southern-located climatic refugia likely occurred in the Bølling-Allerød and Middle Holocene periods, which are supported by signals of demographic growth in the west clade and South-East-North subclade. Overall, this work improves our understanding of the historical biogeography of H. viridiflavus, providing further insights into the evolutionary processes that occurred in the Mediterranean Basin hotspot.

Projected northward shifts in eastern red-backed salamanders due to changing climate

Many species' distributions are being impacted by the acceleration of climate change. Amphibians in particular serve numerous ecosystem functions and are useful indicators of environmental change. Understanding how their distributions have been impacted by climate change and will continue to be impacted is thus important to overall ecosystem health. Plethodon cinereus (Eastern Red-Backed Salamander) is a widespread species of lungless salamander (Plethodontidae) that ranges across northeastern North America. To better understand future potential lungless salamander range shifts, we quantify environmental favorability, the likelihood of membership in a set of sites where environmental conditions are favorable for a species, for P. cinereus in multiple time periods, and examine shifts in the species' distribution. First, utilizing a large data set of georeferenced records, we assessed which bioclimatic variables were associated with environmental favorability in P. cinereus. We then used species distribution modeling for two time periods (1961-1980 and 2001-2020) to determine whether there was a regional shift in environmental favorability in the past 60 years. Models were then used to project future distributions under eight climate change scenarios to quantify potential range shifts. Shifts were assessed using fuzzy logic, avoiding thresholds that oversimplify model predictions into artificial binary outputs. We found that P. cinereus presence is strongly associated with environmental stability. There has been a substantial northward shift in environmental favorability for P. cinereus between 1961-1980 and 2001-2020. This shift is predicted to continue by 2070, with larger shifts under higher greenhouse gas emission scenarios. As climate change accelerates, it is differentially impacting species but has especially strong impacts on dispersal-limited species. Our results show substantial northward shifts in climatic favorability in the last 60 years for P. cinereus, which are likely to be exacerbated by ongoing climate change. Since P. cinereus is dispersal-limited, these models may imply local extirpations along the southern modern range with limited northward dispersal. Continued monitoring of amphibians in the field will reveal microclimatic effects associated with climate change and the accuracy of the model predictions presented here.

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.

The Strait of Gibraltar is an ineffective palaeogeographic barrier for some flightless darkling beetles (Coleoptera: Tenebrionidae: Pimelia)

The geographic distribution of a species is shaped by its biology and by environmental and palaeogeographic factors that interact at different spatial-temporal scales, which leads to distributions and diversification patterns observed between and within lineages. The darkling beetle genus Pimelia has been diversifying for more than 31.2 Mya showing different colonization patterns after the opening of the Gibraltar Strait 5 Mya. Three of the 14 subgenera of Pimelia have populations on both sides of the Strait. Through extensive sampling and the analysis of three molecular markers, we determine levels of intra- and interspecific genetic variation, identify evolutionary lineages in subgenera, estimate their temporal origin and distribution ranges and discuss the historical basis for the geographic and diversification patterns of Pimelia around the Strait. This single geographical feature acted both as a barrier and as a dispersal route for different Pimelia species. The Strait has represented a strong barrier for the subgenus Magrebmelia since the Middle Miocene. However, the subgenera Amblyptera and Amblypteraca share repetitive signatures of post-Messinian colonization across the Strait, possibly driven by stochastic or ‘catastrophic’ events such as tsunamis. Our demographic analyses support Wallace’s hypothesis on insect dispersal stochasticity. Some taxonomic changes, including the designation of a lectotype for Pimelia maura, are also proposed.

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.

Niche differentiation in a postglacial colonizer, the bank vole Clethrionomys glareolus

Species-level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species-wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid-Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.

Physical and ecological isolation contribute to maintain genetic differentiation between fire salamander subspecies

Landscape features shape patterns of gene flow among populations, ultimately determining where taxa lay along the continuum between panmixia to complete reproductive isolation. Gene flow can be restricted, leading to population differentiation in two non-exclusive ways: "physical isolation", in which geographic distance in combination with the landscape features restricts movement of individuals promoting genetic drift, and "ecological isolation", in which adaptive mechanisms constrain gene flow between different environments via divergent natural selection. In central Iberia, two fire salamander subspecies occur in parapatry across elevation gradients along the Iberian Central System mountains, while in the adjacent Montes de Toledo Region only one of them occurs. By integrating population and landscape genetic analyses, we show a ubiquitous role of physical isolation between and within mountain ranges, with unsuitable landscapes increasing differentiation between populations. However, across the Iberian Central System, we found strong support for a significant contribution of ecological isolation, with low genetic differentiation in environmentally homogeneous areas, but high differentiation across sharp transitions in precipitation seasonality. These patterns are consistent with a significant contribution of ecological isolation in restricting gene flow among subspecies. Overall, our results suggest that ecological divergence contributes to reduce genetic admixture, creating an opportunity for lineages to follow distinct evolutionary trajectories.

Evaluating surrogates of genetic diversity for conservation planning

Protected-area systems should conserve intraspecific genetic diversity. Because genetic data require resources to obtain, several approaches have been proposed for generating plans for protected-area systems (prioritizations) when genetic data are not available. Yet such surrogate-based approaches remain poorly tested. We evaluated the effectiveness of potential surrogate-based approaches based on microsatellite genetic data collected across the Iberian Peninsula for 7 amphibian and 3 reptilian species. Long-term environmental suitability did not effectively represent sites containing high genetic diversity (allelic richness). Prioritizations based on long-term environmental suitability had similar performance to random prioritizations. Geographic distances and resistance distances based on contemporary environmental suitability were not always effective surrogates for identification of combinations of sites that contain individuals with different genetic compositions. Our results demonstrate that population genetic data based on commonly used neutral markers can inform prioritizations, and we could not find an adequate substitute. Conservation planners need to weigh the potential benefits of genetic data against their acquisition costs.

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo)

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

Rattlesnake (Crotalus spp.) distribution and diversity in Zacatecas, Mexico

Mexico is home to a large number of reptile species and has one of the greatest diversities of venomous snakes, among which the rattlesnakes pertaining to the Crotalus genus stand out. Out of more than 40 species in the country, nine are found in Zacatecas: C.aquilus, C.atrox, C.basiliscus, C.lepidus, C.molossus, C.polystictus, C.pricei, C.scutulatus and C.willardi. Although these reptiles are important, due to their relevance in terms of ecology, cultural use and public health, their conservation is impacted by multiple factors, such as habitat fragmentation and indiscriminate killing. Thus, most species within this genus are found in some type of risk category at both the national and international level. The purpose of this study was to determine the potential distribution and diversity of rattlesnakes at the municipal level in the understudied state of Zacatecas. To do this, we analyzed and described the global distribution of nine rattlesnake species by building species distribution models, which determined their potential distribution based on a set of ecological variables and presence records. The resulting models were used to assess the diversity of rattlesnake species potentially present in each municipality within the state. Thirty-nine (67.24 %) out of fifty-eight municipalities registered at least one rattlesnake species. Fresnillo, Sombrerete and Valparaíso were some of the municipalities showing greatest diversity. Moreover, C.atrox, C.lepidus, C.molossus and C.scutulatus were the most widely found species in the state. On the other hand, C.basiliscus, C.polystictus, C.pricei and C.willardi were rarely spotted and so, information on their distribution patterns within Zacatecas is limited. Finally, the areas having the largest potential for the distribution of these species were defined. These findings should make field work much more time- and cost-effective, facilitating the collection of in situ data that are useful for management and conservation plans of these species in Zacatecas.

Toward reliable habitat suitability and accessibility models in an era of multiple environmental stressors

Global biodiversity declines, largely driven by climate and land-use changes, urge the development of transparent guidelines for effective conservation strategies. Species distribution modeling (SDM) is a widely used approach for predicting potential shifts in species distributions, which can in turn support ecological conservation where environmental change is expected to impact population and community dynamics. Improvements in SDM accuracy through incorporating intra- and interspecific processes have boosted the SDM field forward, but simultaneously urge harmonizing the vast array of SDM approaches into an overarching, widely adoptable, and scientifically justified SDM framework. In this review, we first discuss how climate warming and land-use change interact to govern population dynamics and species' distributions, depending on species' dispersal and evolutionary abilities. We particularly emphasize that both land-use and climate change can reduce the accessibility to suitable habitat for many species, rendering the ability of species to colonize new habitat and to exchange genetic variation a crucial yet poorly implemented component of SDM. We then unite existing methodological SDM practices that aim to increase model accuracy through accounting for multiple global change stressors, dispersal, or evolution, while shifting our focus to model feasibility. We finally propose a roadmap harmonizing model accuracy and feasibility, applicable to both common and rare species, particularly those with poor dispersal abilities. This roadmap (a) paves the way for an overarching SDM framework allowing comparison and synthesis of different SDM studies and (b) could advance SDM to a level that allows systematic integration of SDM outcomes into effective conservation plans.

An odyssey out of Africa: an integrative review of past and present invasions by the Mediterranean tree frog (Hyla meridionalis)

An illustration of the human footprint on biodiversity are the faunal movements that have accompanied commercial and cultural exchanges between civilizations throughout history. In this article, we provide an integrative review of biogeographical and archaeological knowledge to understand these processes for the Mediterranean tree frog (Hyla meridionalis), an African species that has reached several regions of south-western Europe and associated islands. By re-analysing molecular (multilocus mitochondrial and genome-wide markers) and bioclimatic (ecological niche modelling) data, we give a comprehensive picture of the genetic diversity, structure and environmental suitability of H. meridionalis across its natural and exotic ranges, which in turn offer specific clues to the putative routes of colonization and associated events. Long-term monitoring efforts suggest northwards shifts of local range margins, potentially due to global warming, and we further demonstrate that this species is the most frequent amphibian to travel via the food supply chain in Western Europe. High dispersal ability, ecological tolerance and proximity to human settlements have made H. meridionalis a recurrent witness to the complexity and diversity of the civilizations that ruled the Western Mediterranean.

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.

Phylogeny of the Eurasian Wren Nannus troglodytes (Aves: Passeriformes: Troglodytidae) reveals deep and complex diversification patterns of Ibero-Maghrebian and Cyrenaican populations

The Mediterranean Basin represents a Global Biodiversity Hotspot where many organisms show high inter- and intraspecific differentiation. Extant phylogeographic patterns of terrestrial circum-Mediterranean faunas were mainly shaped through Pleistocene range shifts and range fragmentations due to retreat into different glacial refugia. Thus, several extant Mediterranean bird species have diversified by surviving glaciations in different hospitable refugia and subsequently expanded their distribution ranges during the Holocene. Such a scenario was also suggested for the Eurasian Wren (Nannus troglodytes) despite the lack of genetic data for most Mediterranean subspecies. Our phylogenetic multi-locus analysis comprised 18 out of 28 currently accepted subspecies of N. troglodytes, including all but one subspecies which are present in the Mediterranean Basin. The resulting phylogenetic reconstruction dated the onset of the entire Holarctic radiation of three Nannus species to the early Pleistocene. In the Eurasian Wren, two North African subspecies represented separate basal lineages from the Maghreb (N. t. kabylorum) and from the Libyan Cyrenaica (N. t. juniperi), being only distantly related to other Mediterranean populations. Although N. troglodytes appeared to be paraphyletic with respect to the Nearctic Winter Wren (N. hiemalis), respective nodes did not receive strong statistical support. In contrast, paraphyly of the Ibero-Maghrebian taxon N. t. kabylorum was strongly supported. Southern Iberian populations of N. t. kabylorum did not clade with Maghrebian populations of the same subspecies but formed a sister clade to a highly diverse European clade (including nominate N. t. troglodytes and eight further taxa). In accordance with a pattern also found in other birds, Eurasian populations were split into a western clade (Europe, Caucasus) and an eastern clade (Central Asia, Sino-Himalayas, East Asia). This complex phylogeographic pattern revealed cryptic diversification in N. troglodytes, especially in the Iberio-Maghrebian region.

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.

Patterns of morphological diversification in giant Berberomeloe blister beetles (Coleoptera: Meloidae) reveal an unexpected taxonomic diversity concordant with mtDNA phylogenetic structure

Delimiting species boundaries is a complex challenge usually hindered by overlooked morphological diversification or misinterpretation of geographically structured phenotypic variability. Independent molecular data are extremely useful to characterize and understand such morphological diversity. Morphological and molecular variability of the non-phoretic and apterous, widely distributed, giant blister beetles of the genus Berberomeloe, were investigated within and between lineages across most of the distributional range of the genus. We used two mtDNA gene fragments to characterize genetic variability and to produce a time-calibrated phylogeny of the genus. Our results reveal several mitochondrial lineages, allopatrically, parapatrically and sympatrically distributed. Most clades are not distinguishable between each other based on morphometrics. However, no morphometric overlap is observed between two closely related clades, one of them occurring in sympatry with a distantly congeneric species (B. insignis), suggesting that sympatry could trigger morphological diversification. Although most species share a morphometric space, they can be morphologically identified by a combination of easily observed characteristic qualitative features. Based on the concordance between mtDNA clades and morphological units, we describe six new species of Berberomeloe (B. castuo sp. nov., B. comunero sp. nov., B. indalo sp. nov, B. yebli sp. nov., B. payoyo sp. nov. and B. tenebrosus sp. nov.), revalidate two taxa (B. maculifrons comb. nov. and B. laevigatus comb. nov.) and redefine B. majalis.

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.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean F_ST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

Comparative landscape genetics of pond-breeding amphibians in Mediterranean temporal wetlands: The positive role of structural heterogeneity in promoting gene flow

Comparative landscape genetics studies can provide key information to implement cost-effective conservation measures favouring a broad set of taxa. These studies are scarce, particularly in Mediterranean areas, which include diverse but threatened biological communities. Here, we focus on Mediterranean wetlands in central Iberia and perform a multi-level, comparative study of two endemic pond-breeding amphibians, a salamander (Pleurodeles waltl) and a toad (Pelobates cultripes). We genotyped 411 salamanders from 20 populations and 306 toads from 16 populations at 18 and 16 microsatellite loci, respectively, and identified major factors associated with population connectivity through the analysis of three sets of variables potentially affecting gene flow at increasingly finer levels of spatial resolution. Topographic, land use/cover, and remotely sensed vegetation/moisture indices were used to derive optimized resistance surfaces for the two species. We found contrasting patterns of genetic structure, with stronger, finer scale genetic differentiation in Pleurodeles waltl, and notable differences in the role of fine-scale patterns of heterogeneity in vegetation cover and water content in shaping patterns of regional genetic structure in the two species. Overall, our results suggest a positive role of structural heterogeneity in population connectivity in pond-breeding amphibians, with habitat patches of Mediterranean scrubland and open oak woodlands ("dehesas") facilitating gene flow. Our study highlights the usefulness of remotely sensed continuous variables of land cover, vegetation and water content (e.g., NDVI, NDMI) in conservation-oriented studies aimed at identifying major drivers of population connectivity.