Combining phylogeography and landscape genetics to infer the evolutionary history of a short-range Mediterranean relict, Salamandra salamandra longirostris

Temporal Range Dynamics of the Lataste's Viper (Vipera latastei Boscá, 1878) in Doñana (Spain): Insights into Anthropogenically Driven Factors

Doñana (southern Spain), a region of notable biodiversity richness, is highly threatened by ongoing landscape transformation and climate change. We investigated the local effect of these anthropogenic factors on the temporal range dynamics of Lataste's viper (Vipera latastei), an Iberian endemic Mediterranean reptile that has apparently become rare over the years in Doñana. Using ecological niche-based models, based on climatic and remote sensing variables, we analyzed historical (1959-1999) and contemporary (2000-2022) records of the species to assess range shifts and identify environmental factors that may influence them. Our results show that V. latastei is mostly restricted to the coastal region of Doñana and that one temperature variable is the most important factor explaining this distribution pattern in both periods. Additional climatic and vegetation variables play a role in its historical distribution, but they become less important in contemporary times, suggesting a niche simplification over time. We found 30.5% of reduction in the species suitable area from historical to contemporary conditions, a reduction that would be even greater (83.37%) in the absence of niche shift. These findings underscore the species' heightened vulnerability to ongoing environmental changes and highlight the urgent need for targeted conservation strategies.

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.

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.

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.

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.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Niche models at inter- and intraspecific levels reveal hierarchical niche differentiation in midwife toads

Variation and population structure play key roles in the speciation process, but adaptive intraspecific genetic variation is commonly ignored when forecasting species niches. Amphibians serve as excellent models for testing how climate and local adaptations shape species distributions due to physiological and dispersal constraints and long generational times. In this study, we analysed the climatic factors driving the evolution of the genus Alytes at inter- and intraspecific levels that may limit realized niches. We tested for both differences among the five recognized species and among intraspecific clades for three of the species (Alytes obstetricans, A. cisternasii, and A. dickhilleni). We employed ecological niche models with an ordination approach to perform niche overlap analyses and test hypotheses of niche conservatism or divergence. Our results showed strong differences in the environmental variables affecting species climatic requirements. At the interspecific level, tests of equivalence and similarity revealed that sister species were non-identical in their environmental niches, although they neither were entirely dissimilar. This pattern was also consistent at the intraspecific level, with the exception of A. cisternasii, whose clades appeared to have experienced a lower degree of niche divergence than clades of the other species. In conclusion, our results support that Alytes toads, examined at both the intra- and interspecific levels, tend to occupy similar, if not identical, climatic environments.