Mitochondrial sequence analysis of Salamandra taxa suggests old splits of major lineages and postglacial recolonizations of central Europe from distinct source populations of Salamandra salamandra

Monitoring spatiotemporal impacts of changes in land surface temperature on near eastern fire salamander (Salamandra infraimmaculata) in the Middle East

Persistence and coexistence of many pond-breeding amphibians depend on seasonality. Temperature, as a seasonal climate component, affects numerous physical and biological processes of pond-breeding amphibians. Satellite-derived land surface temperature (LST) is the radiative skin temperature of the land surface, which has received less attention in spatiotemporal seasonal habitat monitoring. The present study aims to evaluate the increasing and decreasing effects of LST trends at two levels: (1) habitat suitability and connectivity; (2) individual population sites and their longitudinal distribution (with increasing longitude). Habitat suitability modeling was conducted based on an ensemble species distribution model (eSDM). Using electrical circuit theory, the connectivity of interior and intact habitat cores was investigated. An average seasonal LST was prepared separately for each season from 2003 to 2021 and entered into Mann-Kendall (MK) analysis to determine the spatiotemporal effects of LST changes using the Z-Score (ZMK) at two confidence levels of 95 and 99%. Based on the results, in winter, 28.12% and 70.70% of the suitable habitat were affected by an increasing trend of LST at 95% and 99% confidence levels, respectively. The highest spatial overlap of the decreasing trend of LST with the suitable habitat occurred in summer and was 6.4% at the 95% confidence level and 4.2% at the 99% confidence level. Considering population site at 95% confidence interval, the increasing trend of LST was calculated to be 20.2%, 9.5%, 4.2%, and 6.3% of localities in winter, spring, summer, and autumn, respectively. At the 99% confidence level, these percentages reduced to 8.5%, 3.1%, 1%, and 1%, respectively. During winter and summer, based on the results of the longitudinal trend, an increasing trend of LST was observed in sites. Localities of Hatay and Iica village in Turkey experienced seasonally asynchronous climate change regimes. The approach used in this study allowed us to create a link between the life cycle and seasonal changes on a micro-scale (breeding sites) and macro-scale (distribution and connectivity). Findings of this paper can be effectively used by conservation managers to preserve S. infraimmaculata's metapopulation.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

Measuring athletic performance in post-metamorphic fire salamanders

OBJECTIVE

Athletic performances are dynamic movements that are physically challenging and often predict individual success in ecological contexts. They stem from a complex integration of multiple phenotypic traits-e.g., morphological, physiological and behavioural-that dictate animal survival and individual fitness. However, directly quantifying athletic performances can be particularly challenging in cryptic, slow-moving species or not very reactive in attitude. Here we present and describe a rapid, simple, and low-cost method to measure athletic performance in post-metamorphic individuals of the fire salamander Salamandra salamandra. While extremely reactive during the larval stage, adult salamanders are, in fact, cryptic and relatively slow-moving.

RESULTS

Forcing terrestrial juveniles to swim under standard, albeit ecologically plausible, laboratory conditions, and using an automatic point-mass tracking tool, we were able to measure maximal and average performance indicators of post-metamorphic individuals. This method avoids inter-individual variation in motivation, as it forces individuals to perform at their best. Moreover, with this method, measures of athletic performance will be directly comparable between larval and terrestrial stages, allowing to study the contribution of carryover effects to the wide range of processes implicated in the eco-evo-devo of athletic performance in salamanders.

Historic range dynamics in Kaiser's mountain newt (Neurergus kaiseri): Insights from phylogeographic analyses and species distribution modeling

Vulnerable Kaiser's mountain newt, Neurergus kaiseri, is endemic to highland streams, springs, and pools of the southwestern Zagros mountain, Iran. The present study aimed to use an integration of phylogeographical and species distribution modeling (SDM) approaches to provide new insights into the evolutionary history of the species throughout Quaternary climate oscillations. The phylogeographical analysis was followed by analyzing two mitochondrial DNA (mt-DNA) markers including 127 control region (D-loop) and 72 NADH dehydrogenase 2 (ND2) sequences from 15 populations in the entire species range that were obtained from GenBank. Potential recent and past distribution (the Last Glacial Maximum, LGM, 21 Kya and the Mid-Holocene, 6 Kya) reconstructed by ensemble SDM using nine algorithms with CCSM4, MIROC-ESM, and MPI-ESM-P models. N. kaiseri displayed two distinct lineages in the northern and southern regions that diverged in the Early-Pleistocene. The demographics analysis showed signs of a slight increase in effective population size for both northern and southern populations in the Mid-Pleistocene. Biogeography analysis showed that both vicariance and dispersal events played an important role in the formation of recent species distribution of N. kaiseri. Based on SDM projection onto paleoclimatic data, N. kaiseri displayed a scenario of past range expansion that followed by postglacial contraction. The models showed that the distribution range of the species may have shifted to a lower altitude during LGM while with amelioration of climatic during Mid-Holocene to recent conditions caused the species to shift to the higher altitude. The findings of the current study support the hypothesis that the Zagros mountains​ may be acting as climatic refugia and play an important role in the protection of isolated populations during climate oscillations.

Phylogenetics of mud snakes (Squamata: Serpentes: Homalopsidae): A paradox of both undescribed diversity and taxonomic inflation

Mud snakes (Serpentes: Homalopsidae) are a family of 55 described, mainly aquatic, species primarily distributed throughout mainland Southeast Asia and the Indo-Australian Archipelago. Although they have been the focus of prior research, the basic relationships amongst genera and species remain poorly known. We used a combined mitochondrial and nuclear gene dataset to infer their phylogenetic relationships, using the highest levels of taxon and geographic sampling for any homalopsid phylogeny to date (62% generic and 62% species coverage; 140 individuals). Our results recover two reciprocally monophyletic groups: the fangless Brachyorrhos and its sister clade comprised of all rear-fanged homalopsids. Most genera and interspecific relationships were monophyletic and strongly supported, but intergeneric relationships and intraspecific population structure lack support. We find evidence of both undescribed diversity as well as cases of taxonomic inflation within several species. Tree-based species delimitation approaches (mPTP) support potential new candidate species as distinct from their conspecifics and also suggest that many named taxa may not be distinct species. Divergence date estimation and lineage-through-time analyses indicate lower levels of speciation in the Eocene, with a subsequent burst in diversification in the Miocene. Homalopsids may have diversified most rapidly during the Pliocene and Pleistocene, possibly in relation to tectonic shifts and sea-level fluctuations that took place in Sundaland and the Sahul Shelf. Our analyses provide new insights on homalopsid taxonomy, a baseline phylogeny for the family, and further biogeographic implications demonstrating how dynamic tectonics and Quaternary sea level changes may have shaped a widespread, diverse family of snakes.

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.

A case of color aberration in a fire salamander (Salamandra salamandra) larva

Color aberrations in salamanders are caused by defective chromatophores or by disrupted synthesis of the pigments, and have been previously reported in Europe. In this paper we report the first case of leucism in a fire salamander larva from a cave-dwelling population in Romania. The leucistic larva was larger and heavier compared to the other larvae from the same breeding habitat.

The evolution of pueriparity maintains multiple paternity in a polymorphic viviparous salamander

The reduction in fecundity associated with the evolution of viviparity may have far-reaching implications for the ecology, demography, and evolution of populations. The evolution of a polygamous behaviour (e.g. polyandry) may counteract some of the effects underlying a lower fecundity, such as the reduction in genetic diversity. Comparing patterns of multiple paternity between reproductive modes allows us to understand how viviparity accounts for the trade-off between offspring quality and quantity. We analysed genetic patterns of paternity and offspring genetic diversity across 42 families from two modes of viviparity in a reproductive polymorphic species, Salamandra salamandra. This species shows an ancestral (larviparity: large clutches of free aquatic larvae), and a derived reproductive mode (pueriparity: smaller clutches of larger terrestrial juveniles). Our results confirm the existence of multiple paternity in pueriparous salamanders. Furthermore, we show the evolution of pueriparity maintains, and even increases, the occurrence of multiple paternity and the number of sires compared to larviparity, though we did not find a clear effect on genetic diversity. High incidence of multiple paternity in pueriparous populations might arise as a mechanism to avoid fertilization failures and to ensure reproductive success, and thus has important implications in highly isolated populations with small broods.

Functional colour genes and signals of selection in colour-polymorphic salamanders

Coloration has been associated with multiple biologically relevant traits that drive adaptation and diversification in many taxa. However, despite the great diversity of colour patterns present in amphibians the underlying molecular basis is largely unknown. Here, we use insight from a highly colour-variable lineage of the European fire salamander (Salamandra salamandra bernardezi) to identify functional associations with striking variation in colour morph and pattern. The three focal colour morphs-ancestral black-yellow striped, fully yellow and fully brown-differed in pattern, visible coloration and cellular composition. From population genomic analyses of up to 4,702 loci, we found no correlations of neutral population genetic structure with colour morph. However, we identified 21 loci with genotype-phenotype associations, several of which relate to known colour genes. Furthermore, we inferred response to selection at up to 142 loci between the colour morphs, again including several that relate to coloration genes. By transcriptomic analysis across all different combinations, we found 196 differentially expressed genes between yellow, brown and black skin, 63 of which are candidate genes involved in animal coloration. The concordance across different statistical approaches and 'omic data sets provide several lines of evidence for loci linked to functional differences between colour morphs, including TYR, CAMK1 and PMEL. We found little association between colour morph and the metabolomic profile of its toxic compounds from the skin secretions. Our research suggests that current ecological and evolutionary hypotheses for the origins and maintenance of these striking colour morphs may need to be revisited.

Hybridizing salamanders experience accelerated diversification

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

A reference DNA barcode library for Austrian amphibians and reptiles

In the last few years, DNA barcoding became an established method for species identification in biodiversity inventories and monitoring studies. Such studies depend on the access to a comprehensive reference data base, covering all relevant taxa. Here we present a comprehensive DNA barcode inventory of all amphibian and reptile species native to Austria, except for the putatively extinct Vipera ursinii rakosiensis and Lissotriton helveticus, which has been only recently reported for the very western edge of Austria. A total of 194 DNA barcodes were generated in the framework of the Austrian Barcode of Life (ABOL) initiative. Species identification via DNA barcodes was successful for most species, except for the hybridogenetic species complex of water frogs (Pelophylax spp.) and the crested newts (Triturus spp.), in areas of sympatry. However, DNA barcoding also proved powerful in detecting deep conspecific lineages, e.g. within Natrix natrix or the wall lizard (Podarcis muralis), resulting in more than one Barcode Index Number (BIN) per species. Moreover, DNA barcodes revealed the presence of Natrix helvetica, which has been elevated to species level only recently, and genetic signatures of the Italian water frog Pelophylax bergeri in Western Austria for the first time. Comparison to previously published DNA barcoding data of European amphibians and reptiles corroborated the results of the Austrian data but also revealed certain peculiarities, underlining the particular strengths and in the case of the genus Pelophylax also the limitations of DNA barcoding. Consequently, DNA barcoding is not only powerful for species identification of all life stages of most Austrian amphibian and reptile species, but also for the detection of new species, the monitoring of gene flow or the presence of alien populations and/or species. Thus, DNA barcoding and the data generated in this study may serve both scientific and national or even transnational conservation purposes.

Phylogeographical structure of the pygmy shrew: revisiting the roles of southern and northern refugia in Europe

Southern and northern glacial refugia are considered paradigms that explain the complex phylogeographical patterns and processes of European biota. Here, we provide a revisited statistical phylogeographical analysis of the pygmy shrew Sorex minutus Linnaeus, 1766 (Eulipotyphla, Soricidae), examining its genetic diversity, genetic differentiation and demographic history in the Mediterranean peninsulas and in Western and Central Europe. The results showed support for genetically distinct and diverse phylogeographical groups consistent with southern and northern glacial refugia, as expected from previous studies. We also identified geographical barriers concordant with glaciated mountain ranges during the Last Glacial Maximum (LGM), early diversification events dated between the Late Pleistocene and Early Holocene for the main phylogeographical groups, and recent (post-LGM) patterns of demographic expansions. This study is the most comprehensive investigation of this species to date, and the results have implications for the conservation of intraspecific diversity and the preservation of the evolutionary potential of S. minutus.

'Mainland-island' population structure of a terrestrial salamander in a forest-bocage landscape with little evidence for in situ ecological speciation

Adaptation to different ecological environments can, through divergent selection, generate phenotypic and genetic differences between populations, and eventually give rise to new species. The fire salamander (Salamandra salamandra) has been proposed to represent an early stage of ecological speciation, driven by differential habitat adaptation through the deposition and development of larvae in streams versus ponds in the Kottenforst near Bonn (Germany). We set out to test this hypothesis of ecological speciation in an area different from the one where it was raised and we took the opportunity to explore for drivers of genetic differentiation at a landscape scale. A survey over 640 localities demonstrated the species' presence in ponds and streams across forests, hilly terrain and areas with hedgerows ('bocage'). Genetic variation at 14 microsatellite loci across 41 localities in and around two small deciduous forests showed that salamander effective population sizes were higher in forests than in the bocage, with panmixia in the forests (Fst < 0.010) versus genetic drift or founder effects in several of the small and more or less isolated bocage populations (Fst > 0.025). The system fits the 'mainland-island' metapopulation model rather than indicating adaptive genetic divergence in pond versus stream larval habitats. A reanalysis of the Kottenforst data indicated that microsatellite genetic variation fitted a geographical rather than an environmental axis, with a sharp transition from a western pond-breeding to an eastern, more frequently stream-breeding group of populations. A parallel changeover in mitochondrial DNA exists but remains to be well documented. The data support the existence of a hybrid zone following secondary contact of differentiated lineages, more so than speciation in situ.

A new subspecies of African fire salamander Salamandra algira (Urodela, Salamandridae) from the Middle Atlas Mountains, Morocco

A new subspecies within the Salamandra algira species complex from north-western Africa is described. Previous molecular analysis showed that the populations of S. algira splendens from north-western Morocco consisted of two well supported clades: clade 1 distributed in the Rif Mountains, from Chefchaouen (type locality) to Al Hoceima, and clade 2, located southern from clade 1 being isolated in the northern and central Middle Atlas Mountains. Clade 2 is herein described as a distinct subspecies: Salamandra algira atlantica ssp. nov. based on morphological data, allopatric range and molecular divergence. This new subspecies shows an uncorrected pairwise distance of 0.0265 from clade 1 based on cytochrome b DNA sequences. Salamandra algira atlantica ssp. nov. is a slender and large sized salamander with a highly variable colouration pattern. It can be distinguished from S. algira splendens by the greater proportion of coral red in the background colouration, being the only known subspecies of S. algira in which coral red can exceed the proportion of black. Variable number (0-5) of yellow to golden yellow dorsal blotches, but usually in lower numbers than the nominotypical subspecies. Salamandra algira atlantica ssp. nov. inhabits subhumid to humid forests and karstic systems at mid to high elevations. We briefly discuss the phylogenetic and taxonomic issues among the genus Salamandra which encompasses more valid species than currently recognised.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Detection of elusive fire salamander larvae (Salamandra salamandra) in streams via environmental DNA

In the face of the global biodiversity crisis, the monitoring of species richness and diversity is experiencing an increased demand entailing a raise in cost and time investment. The analysis of species-specific DNA fragments in environmental samples (eDNA) such as from water or soil, facilitate the molecular detection of species without the specific sampling of individuals. The invasive chytrid fungus Batrachochytrium salamandrivorans (Bsal) is infecting natural fire salamander populations (Salamandra salamandra) and causes chytridiomycosis resulting in infrequent regional extinctions of populations across Central Europe. With regard to the expanding distribution of Bsal over the last years, cost-effective monitoring of fire salamanders is important for the conservation of this species. Based on a real-time quantitative PCR (qPCR) assay, we developed a new protocol to detect S. salamandra larvae in streams via eDNA, using species-specific primers of the mitochondrial control region (D-loop). We tested the efficiency of qPCR primer sets for six combinations of DNA extraction kits coupled with subsequent PCR inhibitor removal kits for obtaining qPCR-detectable S. salamandra eDNA from water filters, that were taken both from natural streams and artificial water tanks in the laboratory as positive controls. We found that the DNeasy Blood & Tissue Kit in combination with the DNeasy PowerClean CleanUp Kit performed best for detecting salamander larvae from natural streams. Our experimental protocol paves the way for resource-saving approaches to monitor S. salamandra larvae, but also confirms the limits to this eDNA approach in that it requires optimized laboratory protocols.

Hybridization and extensive mitochondrial introgression among fire salamanders in peninsular Italy

Discordance between mitochondrial and nuclear patterns of population genetic structure is providing key insights into the eco-evolutionary dynamics between and within species, and their assessment is highly relevant to biodiversity monitoring practices based on DNA barcoding approaches. Here, we investigate the population genetic structure of the fire salamander Salamandra salamandra in peninsular Italy. Both mitochondrial and nuclear markers clearly identified two main population groups. However, nuclear and mitochondrial zones of geographic transition between groups were located 600 km from one another. Recent population declines in central Italy partially erased the genetic imprints of past hybridization dynamics. However, the overall pattern of genetic variation, together with morphological and fossil data, suggest that a rampant mitochondrial introgression triggered the observed mitonuclear discordance, following a post-glacial secondary contact between lineages. Our results clearly show the major role played by reticulate evolution in shaping the structure of Salamandra salamandra populations and, together with similar findings in other regions of the species' range, contribute to identify the fire salamander as a particularly intriguing case to investigate the complexity of mechanisms triggering patterns of mitonuclear discordance in animals.

Plasticity and evolutionary divergence in gene expression associated with alternative habitat use in larvae of the European Fire Salamander

Transcriptomes of organisms reveal differentiation associated with the use of different habitats. However, this leaves open how much of the observed differentiation can be attributed to genetic differences or to transcriptional plasticity. In this study, we disentangle causes of differential gene expression in larvae of the European fire salamander from the Kottenforst forest in Germany. Larvae inhabit permanent streams and ephemeral ponds and represent an example of a young evolutionary split associated with contrasting ecological conditions. We hypothesized that adaptation towards differences in water temperature plays a role because the thermal regime between stream and pond habitats differs notably. Tissue samples from tail fins of larvae were collected to study gene expression using microarrays. We found ample evidence for differentiation among larvae occupying different habitats in nature with 2,800 of 11,797 genes being differentially expressed. We then quantified transcriptional plasticity towards temperature and genetic differentiation based on controlled temperature laboratory experiments. Gene-by-environment interactions modelling revealed that 28% of the gene expression divergence observed among samples in nature could be attributed to plasticity related to water temperature. Expression patterns of only a small number of 101 genes were affected by the genotype. Our analysis demonstrates that effects of environmental factors must be taken into account to explain variation of gene expression in salamanders in nature. Notwithstanding, it provides first evidence that genetic factors determined gene expression divergence between pond and stream ecotypes and could be involved in adaptive evolution.

Parallel habitat acclimatization is realized by the expression of different genes in two closely related salamander species (genus Salamandra)

The utilization of similar habitats by different species provides an ideal opportunity to identify genes underlying adaptation and acclimatization. Here, we analysed the gene expression of two closely related salamander species: Salamandra salamandra in Central Europe and Salamandra infraimmaculata in the Near East. These species inhabit similar habitat types: 'temporary ponds' and 'permanent streams' during larval development. We developed two species-specific gene expression microarrays, each targeting over 12 000 transcripts, including an overlapping subset of 8331 orthologues. Gene expression was examined for systematic differences between temporary ponds and permanent streams in larvae from both salamander species to establish gene sets and functions associated with these two habitat types. Only 20 orthologues were associated with a habitat in both species, but these orthologues did not show parallel expression patterns across species more than expected by chance. Functional annotation of a set of 106 genes with the highest effect size for a habitat suggested four putative gene function categories associated with a habitat in both species: cell proliferation, neural development, oxygen responses and muscle capacity. Among these high effect size genes was a single orthologue (14-3-3 protein zeta/YWHAZ) that was downregulated in temporary ponds in both species. The emergence of four gene function categories combined with a lack of parallel expression of orthologues (except 14-3-3 protein zeta) suggests that parallel habitat adaptation or acclimatization by larvae from S. salamandra and S. infraimmaculata to temporary ponds and permanent streams is mainly realized by different genes with a converging functionality.

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.

Phylogeny, age structure, growth dynamics and colour pattern of the Salamandra algira algira population in the Edough Massif, northeastern Algeria

The distribution of the Fire Salamander in North Africa is discontinuous and the Edough Peninsula, Algeria, is considered as the eastern edge of the distribution area. In the current study, we establish a description of the Salamandra algira algira population in its type locality. In this context, an analysis of the mitochondrial DNA D-loop of 47 sequences comes to confirm the phylogenetic status of our population with regard to the other Algerian and Moroccan populations. Also, we used the skeletochronological method for establishing the age structure of the population. Maximum longevity reached 18 years, with a high frequency of young adults, which suggests a good survival of the juveniles. The growth of males is faster than that of the females, although the maximum size of the males is 180 mm, while that of the females is 210 mm. The Edough’s salamander’s phenotype is characterized by multiple small white spots dispersed in different parts of the body (belly, sides, legs and throat) and a high number of large red spots. These red spots are surrounded by a ring of small white spots on the lower part of the body and sometimes on the legs, thus creating specific eyespots that are often aligned along the lower sides.

How little do we know about the reproductive mode in the north African salamander, Salamandra algira? Pueriparity in divergent mitochondrial lineages of S. a. tingitana

Salamandra algirais one of the few species of vertebrates that displays intraspecific variation in reproductive strategies, making it a valuable study model in reproductive mode evolution. How pueriparity (viviparity) inS. algiraarose remains unclear mainly due to insufficient information on the distribution of both reproductive modes in the species. This note summarizes what is known about pueriparous reproduction inS. algiraand adds new data on the distribution and geographic limits of both reproductive strategies and mitochondrial lineages withinS. a. tingitanaacross the Western Rif in Morocco, with possible implications for our understanding of the ecological context of pueriparity. The occurrence of pueriparity in two mitochondrial sublineages suggests multiple events of pueriparity evolution inS. algira, or a model of mitochondrial capture through hybridization and introgression across the contact zone.

The complete mitochondrial genome of Salamandra salamandra (Amphibia: Urodela: Salamandridae)

The complete mitochondrial genome of Salamandra salamandra was reconstructed using bycatch sequences from an RNAseq library generated from muscle tissue. This study shows the potential of using data originally produced for transcriptome assembly to additionally generate complete mitochondrial genomes. The resulting mitogenome was circular, consisted of 16,331 bp and followed the standard vertebrate gene order. Subsequent analysis of the current mitochondrial genome sequence for the species as listed on GenBank highlighted a species misidentification.

Population Structure and Evolution after Speciation of the Hokkaido Salamander (Hynobius retardatus)

The Hokkaido salamander (Hynobius retardatus) is endemic to Hokkaido Island, Japan, and shows intriguing flexible phenotypic plasticity and regional morphological diversity. However, to date, allozymes and partial mitochondria DNA sequences have provided only an outline of its demographic histories and the pattern of its genetic diversification. To understand the finer details of the population structure of this species and its evolution since speciation, we genotyped five regional populations by using 12 recently developed microsatellite polymorphic markers. We found a clear population structure with low gene flow among the five populations, but a close genetic relationship between the Teshio and Kitami populations. Our demographic analysis suggested that Teshio and Erimo had the largest effective population sizes among the five populations. These findings regarding the population structure and demography of H. retardatus improve our understanding of the faunal phylogeography on Hokkaido Island and also provide fundamental genetic information that will be useful for future studies.

Morphological variation in salamanders and their potential response to climate change

Despite the recognition that some species might quickly adapt to new conditions under climate change, demonstrating and predicting such a fundamental response is challenging. Morphological variations in response to climate may be caused by evolutionary changes or phenotypic plasticity, or both, but teasing apart these processes is difficult. Here, we built on the number of thoracic vertebrae (NTV) in ectothermic vertebrates, a known genetically based feature, to establish a link with body size and evaluate how climate change might affect the future morphological response of this group of species. First, we show that in old-world salamanders, NTV variation is strongly related to changes in body size. Secondly, using 22 salamander species as a case study, we found support for relationships between the spatial variation in selected bioclimatic variables and NTV for most of species. For 44% of species, precipitation and aridity were the predominant drivers of geographical variation of the NTV. Temperature features were dominant for 31% of species, while for 19% temperature and precipitation played a comparable role. This two-step analysis demonstrates that ectothermic vertebrates may evolve in response to climate change by modifying the number of thoracic vertebrae. These findings allow to develop scenarios for potential morphological evolution under future climate change and to identify areas and species in which the most marked evolutionary responses are expected. Resistance to climate change estimated from species distribution models was positively related to present-day species morphological response, suggesting that the ability of morphological evolution may play a role for species' persistence under climate change. The possibility that present-day capacity for local adaptation might help the resistance response to climate change can be integrated into analyses of the impact of global changes and should also be considered when planning management actions favouring species persistence.

Phylogeny, phylogeography, and evolution in the Mediterranean region: News from a freshwater mussel (Potomida, Unionida)

The Potomida genus (Bivalvia, Unionida) has a Circum-Mediterranean distribution and like other freshwater mussel species, its populations have suffered dramatic declines. Although this genus is currently considered as monotypic, it has a long history of taxonomic revisions and presently many aspects of its systematics and evolutionary history are unclear. We sampled a total of 323 individuals from 39 different sites across the Potomida genus distribution, and sequenced two mitochondrial (16S rDNA and Cytochrome c Oxidase Subunit I) and one nuclear (28S rDNA) genes to clarify its phylogeny and phylogeographic history. Our results show that the genus includes two well-supported clades, one comprising solely the western Mediterranean species Potomida littoralis, and the other including two eastern Mediterranean species, the Greek endemic P. acarnanica and the Anatolian and Middle Eastern P. semirugata. We suggest that Potomida started radiating during the upper Miocene, and that both vicariance and dispersal events shaped the diversification and distribution of the genus along the Mediterranean region. P. littoralis is further divided in two mitochondrial lineages, one restricted to Europe and the other occurring mostly in North Africa. Moreover, some European basins present both lineages in sympatry. The conservation status of the three recognized species should be reevaluated, particularly P. acarnanica, since it is restricted to two Greek river basins presenting a high risk of extinction. Overall, our results clarify some important gaps in knowledge concerning the phylogeny, phylogeography and evolution of the Potomida genus in the Mediterranean region with important taxonomical, ecological and conservational implications.

Population Genetic Structure of the Endangered Kaiser's Mountain Newt, Neurergus kaiseri (Amphibia: Salamandridae)

Species often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in an endemic and critically endangered stream breeding mountain newt, Neurergus kaiseri, within its entire range in southwestern Iran. We identified two geographic regions based on phylogenetic relationships using Bayesian inference and maximum likelihood of 779 bp mtDNA (D-loop) in 111 individuals from ten of twelve known breeding populations. This analysis revealed a clear divergence between northern populations, located in more humid habitats at higher elevation, and southern populations, from drier habitats at lower elevations regions. From seven haplotypes found in these populations none was shared between the two regions. Analysis of molecular variance (AMOVA) of N. kaiseri indicates that 94.03% of sequence variation is distributed among newt populations and 5.97% within them. Moreover, a high degree of genetic subdivision, mainly attributable to the existence of significant variance among the two regions is shown (θ_CT = 0.94, P = 0.002). The positive and significant correlation between geographic and genetic distances (r = 0.61, P = 0.002) following controlling for environmental distance suggests an important influence of geographic divergence of the sites in shaping the genetic variation and may provide tools for a possible conservation based prioritization policy for the endangered species.

Phylogeographic relationships and shallow mitochondrial divergence of Algerian populations of Salamandra algira

The North African fire salamander, Salamandra algira, is distributed in Algeria, Morocco and Ceuta (Spanish territory located on the north coast of Africa), but until now rather limited information has been available on the populations across the Algerian part of its range. We here provide a first analysis of the phylogeography of this species in Algeria, based on 44 samples from populations distributed across 15 localities in Central Algeria. We sequenced three segments of mitochondrial DNA, covering 12S rRNA, cytochrome b (Cytb) and the D-loop. The mtDNA sequences of the Algerian samples were strongly different from the Moroccan populations occurring west of the Moulouya River (corresponding to the subspecies S. a. tingitana and S. a. splendens) but sister to the genetically rather similar population from the Beni Snassen Massif in eastern Morocco (subspecies S. algira spelaea). Among the Algerian specimens studied, those from the westernmost site, Chrea Massif, were the sister clade to the remaining populations, and the overall genetic divergence was low, with a maximum of five mutational steps in a 295 bp fragment of cytochrome b.

Genetic structure of the fire salamander Salamandra salamandra in the Polish Sudetes

We analysed genetic variation within and differentiation between nineteen populations of the fire salamander Salamandra salamandra inhabiting the north-eastern margin of the species range in the Sudetes Mountains (south-western Poland). The results were compared with those obtained recently for the Polish part of the Carpathians. Variation of 10 nuclear microsatellite loci was analysed in 744 individuals to estimate genetic structure, gene flow, isolation and to test for a geographic gradient of genetic variation. Mitochondrial DNA control region (D-loop) of 252 specimens from all localities was used to identify the origin of populations currently inhabiting its north-eastern range. We found little genetic differentiation among populations in the Sudetes indicating substantial recent or ongoing gene flow. The exceptions were one isolated peripheral population located outside the continuous distribution range which displayed extremely reduced genetic variation probably due to a combination of long term isolation and low population size, and one population located at the eastern margin of the Polish Sudetes. Populations inhabiting the Sudetes and the Carpathians formed two separate clusters based on microsatellite loci. In accordance with available phylogeographic information, single mitochondrial haplotype (type IIb) fixed in all populations indicates their origin from a single refugium and may suggest colonization from the Balkan Peninsula. The analysis of geographic gradient in variation showed its decline in the westerly direction suggesting colonization of Poland from the east, however, alternative scenarios of postglacial colonization could not be rejected with the available data.

Intercohort size structure dynamics of fire salamander larvae in ephemeral habitats: a mesocosm experiment

The size structure of a larval population facilitates interaction asymmetries that, in turn, influence the dynamics of size-structure. In species that exhibit conspicuous aggressive interactions, the competitive effects of the smaller individuals may be overlooked. We manipulated initial size differences between two larval cohorts and young-cohort density of Salamandra infraimmaculata in mesocosms to determine: (1) whether young individuals function primarily as prey or as competitors of older and larger individuals; (2) the resulting dynamics of size variation; and (3) recruitment to the postmetamorph population. Intercohort size differences generally remained constant over time at low young-cohort densities, but reduced over time at high densities due to retardation of the old-cohort growth rate. This suggests a competitive advantage to the young cohort that outweighs the interference advantage of older cohorts previously documented in this species. The increase in mortality from desiccation due to high young-cohort density was an order of magnitude greater in the old cohort than in the young-cohort, further indicating size-dependent vulnerability to competition. However, the conditions least favorable to most of the old-cohort larvae (large size difference and high young-cohort density) promoted cannibalism. Among cannibals, mortality and time to metamorphosis decreased and sizes at metamorphosis increased substantially. Thus, a balance between the competitive advantage to young cohorts, and the interference and cannibalism advantage to old cohorts shapes larval size-structure dynamics. Larval densities and individual expression of cannibalism can shift this balance in opposite directions and alter relative recruitment rates from different cohorts.

Divergence preceding island formation among Aegean insular populations of the freshwater snail genus Pseudorientalia (Caenogastropoda: Truncatelloidea)

Freshwater snails that inhabit islands are excellent model organisms for testing relationships between geological events and phylogeography, especially in the Aegean region. Although many Aegean islands were searched in the present study, species of the genus Pseudorientalia were only found on Lesvos, Samos, and Chios. Phylogenetic relationships between specimens living on these three islands were analysed using COI and 16S rRNA molecular markers and morphological data. A high level of diversity was found between islands. Genetic distances between clades showed differences high enough for the samples from different islands to be considered distinct species (p-distance: 0.105-0.133). These results are also supported by obvious morphological differences in shell morphology between islands. The mean divergence time between the Lesvos clade and Samos/Chios clade was 24.13 ± 3.30 Mya; between the Samos and Chios clades the divergence time was 14.80 ± 1.11 Mya. Our data suggest that high divergence may have occurred between Pseudorientalia populations during the Upper and Middle Miocene, when the Aegean region was part of a united landmass. It is possible that the observed highly divergent Pseudorientalia clades are relicts of high regional diversity that existed in the past.

Genetic structure and differentiation of the fire salamander Salamandra salamandra at the northern margin of its range in the Carpathians

Amphibian populations occurring at the margin of the species range exhibit lower genetic variation due to strong genetic drift and long-term isolation. Limited mobility and site fidelity together with habitat changes may accelerate genetic processes leading to local extinction. Here, we analyze genetic variation of the fire salamander subspecies Salamandra s. salamandra inhabiting the Outer Carpathian region in Poland, at the northern border of its distribution. Nuclear DNA polymorphism based on 10 microsatellite loci of 380 individuals sampled in 11 populations were analysed to measure gene flow between subpopulations and possible long-term isolation. Mitochondrial DNA control region analysis among 17 individuals representing 13 localities was used to detect the origin of populations which colonized Northern Europe after the last glaciation. Overall, pairwise FST’s and AMOVA test of ‘among group’ variation showed little differences in the allele frequencies and relatively high local gene flow. However, Bayesian clustering results revealed subtle structuring between eastern and western part of the studied region. Two extreme marginal populations from the Carpathian Piedmont revealed reduced genetic variation which may be attributed to strong influence of genetic drift. Only one mitochondrial DNA haplotype (type IIb) was found in all individuals and suggest that after the Last Glacial Maximum Salamandra salamandra migrated to the North-Western Europe from the single glacial refugium placed in the Balkan Peninsula.