From species divergence to population structure: A multimarker approach on the most basal lineage of Salamandridae, the spectacled salamanders (genus Salamandrina) from Italy

Different Traits, Different Evolutionary Pathways: Insights from Salamandrina (Amphibia, Caudata)

Species delimitation is often based on a single or very few genetic or phenetic traits, something which leads to misinterpretations and often does not provide information about evolutionary processes. Here, we investigated the diversity pattern of multiple phenetic traits of the two extant species of Salamandrina, a genus split only after molecular traits had been studied but the two species of which are phenetically very similar. The phenetic traits we studied are size, external body shape and head colour pattern, in a model comparison framework using non-linear mixed models and unsupervised and supervised clustering. Overall, we found high levels of intra-specific variability for body size and shape, depending on population belonging and habitat, while differences between species were generally lower. The habitat the salamanders dwell in also seems important for colour pattern. Basing on our findings, from the methodological point of view, we suggest (i) to take into account the variability at population level when testing for higher level variability, and (ii) a semi-supervised learning approach to high dimensional data. We also showed that different phenotypic traits of the same organism could result from different evolutionary routes. Local adaptation is likely responsible for body size and shape variability, with selective pressures more similar across species than within them. Head colour pattern also depends on habitat, differently from ventral colour pattern (not studied in this paper) which likely evolved under genetic drift.

Past, present, and future climate space of the only endemic vertebrate genus of the Italian peninsula

The two extant Salamandrina species represent a unique case of morphology, ecology, and ethology among urodeles. The range of this genus is currently limited to Italy, where it represents the only endemic vertebrate genus, but its past range extended over a much broader area of Europe, including the Iberian and Balkan peninsulas. ENM analyses using modern occurrences of Salamandrina demonstrate that the current climate of the majority of Europe, and especially areas where fossils of this genus were found, is currently not suitable for this genus, neither was it suitable during the last 3.3 million years. This result allows possible assumptions about the climatic influence on the former extirpation of this salamander from several areas of Europe. Furthermore, it shows that, during Pliocene–Pleistocene climatic oscillations, Mediterranean peninsulas, despite being generally considered together because of similar latitude, had different potential to effectively become glacial refugia for this salamander, and possibly for other species as well. Future projections using different CO₂ emission scenarios predict that climatic suitability will be even more drastically reduced during the next 50 years, underlining once more the importance of conservation strategies and emission-reducing policies.

UCE Phylogenomics, detection of a putative hybrid population, and one older mitogenomic node age of Batrachuperus salamanders

The prevalence of incomplete lineage sorting complicates the examination of hybridization and species-level paraphyly with gene trees of a small number of loci. In Asian mountain salamanders of the genus Batrachuperus, possible hybridization and species paraphyly had been identified by utilizing mitochondrial genealogy and fixed allozyme differences. Here we sampled 2909 UCEs in 44 local populations from all six Batrachuperus species, inferred gene and species trees, compared them with mitochondrial and allozyme results, and examined the potential hybridization and species paraphyly. The clustering pattern of single-locus trees, increased proportion of heterozygous SNPs, allele frequency-based migration edge estimation, and intrapopulation long branches (as expected from an increase of genetic lineage and nucleotide diversity) support that an eastern B. karlschmidti population has experienced admixture with B. tibetanus. On the 2909-UCE concatenated and species trees, lower nodal supports were observed when similar proportions of loci agreed with alternative topologies, i.e., a reciprocal monophyly between a Pengxian lineage and the remainder of B. pinchonii (0.379) or a paraphyly of the latter with respect to Pengxian (0.362). The UCE phylogenomics agreed with the relatively recent groupings in the allozyme dendrogram. Despite incomplete lineage sorting, the mitochondrial trees were similar to the UCE trees for deeper relationships of the genus. However, one significant branch-length level discordance was identified. The branch between the common ancestor of B. daochengensis and B. yenyuanensis and common ancestor of the genus was approximately three times shorter on the mitochondrial tree than on the UCE tree, suggesting that the split of the mitochondrial lineages was likely a few million years earlier than the split of species. This finding supports considering possible ancestral polymorphism when interpreting different divergence dates estimated from mitochondrial and genome-wide data.

Reconstructing evolution at the community level: A case study on Mediterranean amphibians

Reconstructing reliable timescales for species evolution is an important and indispensable goal of modern biogeography. However, many factors influence the estimation of divergence times, and uncertainty in the inferred time trees remains a major issue that is often insufficiently acknowledged. We here focus on a fundamental problem of time tree analysis: the combination of slow-evolving (nuclear DNA) and fast-evolving (mitochondrial DNA) markers in a single time tree. Both markers differ in their suitability to infer divergences at different time scales (the 'genome-timescale-dilemma'). However, strategies to infer shallow and deep divergences in a single time tree have rarely been compared empirically. Using Mediterranean amphibians as model system that is exceptional in its geographic and taxonomic completeness of available genetic information, we analyze 202 lineages of western Palearctic amphibians across the entire Mediterranean region. We compiled data of four nuclear and five mitochondrial genes and used twelve fossil calibration points widely acknowledged for amphibian evolution. We reconstruct time trees for an extensive lineage-level data set and compare the performances of the different trees: the first tree is based on primary fossil calibration and mitochondrial DNA, while the second tree is based on a combination of primary fossil and on secondary calibrations taken from a nuclear tree using mitochondrial DNA (two-step protocol). Focusing on a set of nodes that are most likely explained by vicariance, we statistically compare the reconstructed alternative time trees by applying a biogeographical plausibility test. Our two-step protocol outperformed the alternative approach in terms of spatial and temporal plausibility. It allows us to infer scenarios for Mediterranean amphibian evolution in eight geographic provinces. We identified several tectonic and climatic events explaining the majority of Mediterranean amphibian divergences, with Plio-Pleistocene climatic fluctuations being the dominant driver for intrageneric evolution. However, often more than one event could be invoked for a specific split. We give recommendations for the use of secondary calibrations in future molecular clock analyses at the community level.

Evidences for a shared history for spectacled salamanders, haplotypes and climate

The so-called glacial refugia, formed during the Pleistocene climatic oscillations, played a major role in shaping the distribution of European species, triggering migrations or isolating populations. Many of these events were recently investigated by genetic data, mainly for the European Last Glacial stage, in the Iberic, Italian and Greek-Balkan peninsulas. The amphibian genus Salamandrina, the most ancient living salamandrid lineage, was widespread in Europe until the climatic oscillations of Miocene probably forced it to shelter in the only suitable territory at that time, the Apennines. Nowadays this genus is endemic of peninsular Italy with two parapatric species, S. perspicillata and S. terdigitata, sharing an area of secondary contact formed after the Last Glacial Maximum. Climate is generally identified as the key factor for the interpretation of genetic data. In this research, we directly measure climate influences on the two Salamandrina known species through Ensemble Modelling techniques and post-modelling GIS analyses, integrating updated genetic data in this process. Our results confirm the hypotheses of southwards (and subsequent northwards) shifts, identify glacial refugia and corridors used for the post-glacial re-colonization. Finally, we map a contact zone deserving more sampling effort to disentangle the introgression and hybridization observed.

Islands within an island: Population genetic structure of the endemic Sardinian newt, Euproctus platycephalus

The identification of historic and contemporary barriers to dispersal is central to the conservation of endangered amphibians, but may be hindered by their complex life history and elusive nature. The complementary information generated by mitochondrial (mtDNA) and microsatellite markers generates a valuable tool in elucidating population structure and the impact of habitat fragmentation. We applied this approach to the study of an endangered montane newt, Euproctus platycephalus. Endemic to the Mediterranean island of Sardinia, it is threatened by anthropogenic activity, disease, and climate change. We have demonstrated a clear hierarchy of structure across genetically divergent and spatially distinct subpopulations. Divergence between three main mountain regions dominated genetic partitioning with both markers. Mitochondrial phylogeography revealed a deep division dating to ca. 1 million years ago (Mya), isolating the northern region, and further differentiation between the central and southern regions ca. 0.5 Mya, suggesting an association with Pleistocene severe glacial oscillations. Our findings are consistent with a model of southward range expansion during glacial periods, with postglacial range retraction to montane habitat and subsequent genetic isolation. Microsatellite markers revealed further strong population structure, demonstrating significant divergence within the central region, and partial differentiation within the south. The northern population showed reduced genetic diversity. Discordance between mitochondrial and microsatellite markers at this scale indicated a further complexity of population structure, in keeping with male-biased dispersal and female philopatry. Our study underscores the need to elucidate cryptic population structure in the ecology and conservation strategies for endangered island-restricted amphibians, especially in the context of disease and climate change.

No signs of inbreeding despite long-term isolation and habitat fragmentation in the critically endangered Montseny brook newt (Calotriton arnoldi)

Endemic species with restricted geographic ranges potentially suffer the highest risk of extinction. If these species are further fragmented into genetically isolated subpopulations, the risk of extinction is elevated. Habitat fragmentation is generally considered to have negative effects on species survival, despite some evidence for neutral or even positive effects. Typically, non-negative effects are ignored by conservation biology. The Montseny brook newt (Calotriton arnoldi) has one of the smallest distribution ranges of any European amphibian (8 km²) and is considered critically endangered by the International Union for Conservation of Nature. Here we apply molecular markers to analyze its population structure and find that habitat fragmentation owing to a natural barrier has resulted in strong genetic division of populations into two sectors, with no detectable migration between sites. Although effective population size estimates suggest low values for all populations, we found low levels of inbreeding and relatedness between individuals within populations. Moreover, C. arnoldi displays similar levels of genetic diversity to its sister species Calotriton asper, from which it separated around 1.5 million years ago and which has a much larger distribution range. Our extensive study shows that natural habitat fragmentation does not result in negative genetic effects, such as the loss of genetic diversity and inbreeding on an evolutionary timescale. We hypothesize that species in such conditions may evolve strategies (for example, special mating preferences) to mitigate the effects of small population sizes. However, it should be stressed that the influence of natural habitat fragmentation on an evolutionary timescale should not be conflated with anthropogenic habitat loss or degradation when considering conservation strategies.

Extensive unidirectional introgression between two salamander lineages of ancient divergence and its evolutionary implications

Hybridization and introgression, contrary to previous beliefs, are now considered to be widespread processes even among animal species. Nonetheless, the range of their possible outcomes and roles in moulding biodiversity patterns are still far from being fully appraised. Here we investigated the pattern of hybridization and introgression between Salamandrina perspicillata and S. terdigitata, two salamanders endemic to the Italian peninsula. Using a set of diagnostic or differentiated genetic markers (9 nuclear and 1 mitochondrial), we documented extensive unidirectional introgression of S. terdigitata alleles into the S. perspicillata gene pool in central Italy, indicating that barriers against hybridization were permeable when they came into secondary contact, and despite their ancient divergence. Nonetheless, purebred S. terdigitata, as well as F1, F2, and backcrosses were not found within the hybrid zone. Moreover, Bayesian analyses of population structure identified admixed populations belonging to a differentiated gene pool with respect to both parental populations. Overall, the observed genetic structure, together with their geographic pattern of distribution, suggests that Salamandrina populations in central Italy could have entered a distinct evolutionary pathway. How far they have gone along this pathway will deserve future investigation.