De novo transcriptome assembly and annotation for gene discovery in Salamandra salamandra at the larval stage.

Dispersal is a key process in ecology and evolutionary biology, as it shapes biodiversity patterns over space and time. Attitude to disperse is unevenly distributed among individuals within populations, and that individual personality can have pivotal roles in the shaping of this attitude. Here, we assembled and annotated the first de novo transcriptome of the head tissues of Salamandra salamandra from individuals, representative of distinct behavioral profiles. We obtained 1,153,432,918 reads, which were successfully assembled and annotated. The high-quality of the assembly was confirmed by three assembly validators. The alignment of contigs against the de novo transcriptome led to a mapping percentage higher than 94%. The homology annotation with DIAMOND led to 153,048 (blastx) and 95,942 (blastp) shared contigs, annotated on NR, Swiss-Prot and TrEMBL. The domain and site protein prediction led to 9850 GO-annotated contigs. This de novo transcriptome represents reliable reference for comparative gene expression studies between alternative behavioral types, for comparative gene expression studies within Salamandra, and for whole transcriptome and proteome studies in amphibians.

Comparative landscape genetics reveals the evolution of viviparity reduces genetic connectivity in fire salamanders.

Evolutionary changes in reproductive mode may affect co-evolving traits, such as dispersal, although this subject remains largely underexplored. The shift from aquatic oviparous or larviparous reproduction to terrestrial viviparous reproduction in some amphibians entails skipping the aquatic larval stage and, thus, greater independence from water. Accordingly, amphibians exhibiting terrestrial viviparous reproduction may potentially disperse across a wider variety of suboptimal habitats and increase population connectivity in fragmented landscapes compared to aquatic-breeding species. We investigated this hypothesis in the fire salamander (Salamandra salamandra), which exhibits both aquatic- (larviparity) and terrestrial-breeding (viviparity) strategies. We genotyped 426 larviparous and 360 viviparous adult salamanders for 13 microsatellite loci and sequenced a mitochondrial marker for 133 larviparous and 119 viviparous individuals to compare population connectivity and landscape resistance to gene flow within a landscape genetics framework. Contrary to our predictions, viviparous populations exhibited greater differentiation and reduced genetic connectivity compared to larviparous populations. Landscape genetic analyses indicate viviparity may be partially responsible for this pattern, as water courses comprised a significant barrier only in viviparous salamanders, probably due to their fully terrestrial life cycle. Agricultural areas and, to a lesser extent, topography also decreased genetic connectivity in both larviparous and viviparous populations. This study is one of very few to explicitly demonstrate the evolution of a derived reproductive mode affects patterns of genetic connectivity. Our findings open avenues for future research to better understand the eco-evolutionary implications underlying the emergence of terrestrial reproduction in amphibians.

Allopatric diversification and evolutionary melting pot in a North African Palearctic relict: The biogeographic history of Salamandra algira.

North Africa is a climatically and topographically complex region with unique biotic assemblages resulting from the combination of multiple biogeographic realms. Here, we assess the role of climate in promoting intra-specific diversification in a Palearctic relict, the North African fire salamander, Salamandra algira, using a combination of phylogenetic and population genetic analyses, paleoclimatic modelling and niche overlap tests. We used mitochondrial DNA (Cyt-b), 9838 ddRADseq loci, and 14 microsatellite loci to characterize patterns of genetic diversity and population structure. Phylogenetic analyses recover two major clades, each including several lineages with mito-nuclear discordances suggesting introgressive patterns between lineages in the Middle Atlas, associated with a melting pot of genetic diversity. Paleoclimatic modelling identified putative climatic refugia, largely matching areas of high genetic diversity, and supports the role of aridity in promoting allopatric diversification associated with ecological niche conservatism. Overall, our results highlight the role of climatic microrefugia as drivers of populations' persistence and diversification in the face of climatic oscillations in North Africa, and stress the importance of accounting for different genomic regions when reconstructing biogeographic processes from molecular markers.

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.

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.

Inferring the shallow phylogeny of true salamanders (Salamandra) by multiple phylogenomic approaches.

The rise of high-throughput sequencing techniques provides the unprecedented opportunity to analyse controversial phylogenetic relationships in great depth, but also introduces a risk of being misinterpreted by high node support values influenced by unevenly distributed missing data or unrealistic model assumptions. Here, we use three largely independent phylogenomic data sets to reconstruct the controversial phylogeny of true salamanders of the genus Salamandra, a group of amphibians providing an intriguing model to study the evolution of aposematism and viviparity. For all six species of the genus Salamandra, and two outgroup species from its sister genus Lyciasalamandra, we used RNA sequencing (RNAseq) and restriction site associated DNA sequencing (RADseq) to obtain data for: (1) 3070 nuclear protein-coding genes from RNAseq; (2) 7440 loci obtained by RADseq; and (3) full mitochondrial genomes. The RNAseq and RADseq data sets retrieved fully congruent topologies when each of them was analyzed in a concatenation approach, with high support for: (1) S. infraimmaculata being sister group to all other Salamandra species; (2) S. algira being sister to S. salamandra; (3) these two species being the sister group to a clade containing S. atra, S. corsica and S. lanzai; and (4) the alpine species S. atra and S. lanzai being sister taxa. The phylogeny inferred from the mitochondrial genome sequences differed from these results, most notably by strongly supporting a clade containing S. atra and S. corsica as sister taxa. A different placement of S. corsica was also retrieved when analysing the RNAseq and RADseq data under species tree approaches. Closer examination of gene trees derived from RNAseq revealed that only a low number of them supported each of the alternative placements of S. atra. Furthermore, gene jackknife support for the S. atra - S. lanzai node stabilized only with very large concatenated data sets. The phylogeny of true salamanders thus provides a compelling example of how classical node support metrics such as bootstrap and Bayesian posterior probability can provide high confidence values in a phylogenomic topology even if the phylogenetic signal for some nodes is spurious, highlighting the importance of complementary approaches such as gene jackknifing. Yet, the general congruence among the topologies recovered from the RNAseq and RADseq data sets increases our confidence in the results, and validates the use of phylotranscriptomic approaches for reconstructing shallow relationships among closely related taxa. We hypothesize that the evolution of Salamandra has been characterized by episodes of introgressive hybridization, which would explain the difficulties of fully reconstructing their evolutionary relationships.

Salamandra.

Beukema et al. Introduce the genus salamandra comprising the well-known fire salamanders.

Hybridization during altitudinal range shifts: nuclear introgression leads to extensive cyto-nuclear discordance in the fire salamander.

Ecological models predict that, in the face of climate change, taxa occupying steep altitudinal gradients will shift their distributions, leading to the contraction or extinction of the high-elevation (cold-adapted) taxa. However, hybridization between ecomorphologically divergent taxa commonly occurs in nature and may lead to alternative evolutionary outcomes, such as genetic merger or gene flow at specific genes. We evaluate this hypothesis by studying patterns of divergence and gene flow across three replicate contact zones between high- and low-elevation ecomorphs of the fire salamander (Salamandra salamandra) that have experienced altitudinal range shifts over the current postglacial period. Strong population structure with high genetic divergence in mitochondrial DNA suggests that vicariant evolution has occurred over several glacial-interglacial cycles and that it has led to cryptic differentiation within ecomorphs. In current parapatric boundaries, we do not find evidence for local extinction and replacement upon postglacial expansion. Instead, parapatric taxa recurrently show discordance between mitochondrial and nuclear markers, suggesting nuclear-mediated gene flow across contact zones. Isolation with migration models support this hypothesis by showing significant gene flow across all five parapatric boundaries. Together, our results suggest that, while some genomic regions, such as the mitochondria, may follow morphologic species traits and retreat to isolated mountain tops, other genomic regions, such as nuclear markers, may flow across parapatric boundaries, sometimes leading to a complete genetic merger. We show that despite high ecologic and morphologic divergence over prolonged periods of time, hybridization allows for evolutionary outcomes alternative to extinction and replacement of taxa in response to climate change.

Proteomic analysis of the skin of Chinese giant salamander (Andrias davidianus).

The Chinese giant salamander (Andrias davidianus), renowned as a living fossil, is the largest and longest-lived amphibian species in the world. Its skin has developed mucous gland which could secrete a large amount of mucus under the scraping and electric stimulation, and the molting is the degraded skin stratum corneum. Although several proteomic studies have focused on functional proteomes of mammalian and frog skin, the skin proteome of Chinese giant salamander has not yet been carefully studied. To establish the functional skin proteome of Chinese giant salamander, two-dimensional gel electrophoresis (2DE) and mass spectrometry (MS) were applied to detect the composition and relative abundance of the proteins in the skin, mucus and molting. Our findings indicated that 249 proteins were identified in the skin, 155 proteins in the mucus, and 97 proteins in the molting. Furthermore, Gene Ontology (GO) analysis showed that these proteins participated in various physiological activities, including extracellular matrix organization, defense, immune response, wound healing, respiration, etc. In conclusion, the proteomic results provide new insight in the aspects of the proteomes in the skin, mucus and the molting of Chinese giant salamander. BIOLOGICAL SIGNIFICANCE: This was the first study to examine the protein expression abundance in the skin, mucus and molting of Chinese giant salamander by a proteomics approach. Meantime, the identification of a more global proteome in normal skin may provide a basis for characterizing and comparing the skin proteomes from other amphibian species.

Landscape influences on dispersal behaviour: a theoretical model and empirical test using the fire salamander, Salamandra infraimmaculata.

When populations reside within a heterogeneous landscape, isolation by distance may not be a good predictor of genetic divergence if dispersal behaviour and therefore gene flow depend on landscape features. Commonly used approaches linking landscape features to gene flow include the least cost path (LCP), random walk (RW), and isolation by resistance (IBR) models. However, none of these models is likely to be the most appropriate for all species and in all environments. We compared the performance of LCP, RW and IBR models of dispersal with the aid of simulations conducted on artificially generated landscapes. We also applied each model to empirical data on the landscape genetics of the endangered fire salamander, Salamandra infraimmaculata, in northern Israel, where conservation planning requires an understanding of the dispersal corridors. Our simulations demonstrate that wide dispersal corridors of the low-cost environment facilitate dispersal in the IBR model, but inhibit dispersal in the RW model. In our empirical study, IBR explained the genetic divergence better than the LCP and RW models (partial Mantel correlation 0.413 for IBR, compared to 0.212 for LCP, and 0.340 for RW). Overall dispersal cost in salamanders was also well predicted by landscape feature slope steepness (76%), and elevation (24%). We conclude that fire salamander dispersal is well characterised by IBR predictions. Together with our simulation findings, these results indicate that wide dispersal corridors facilitate, rather than hinder, salamander dispersal. Comparison of genetic data to dispersal model outputs can be a useful technique in inferring dispersal behaviour from population genetic data.

Nuclear and mitochondrial multilocus phylogeny and survey of alkaloid content in true salamanders of the genus Salamandra (Salamandridae).

The genus Salamandra represents a clade of six species of Palearctic salamanders of either contrasted black-yellow, or uniformly black coloration, known to contain steroidal alkaloid toxins in high concentrations in their skin secretions. This study reconstructs the phylogeny of the genus Salamandra based on DNA sequences of segments of 10 mitochondrial and 13 nuclear genes from 31 individual samples representing all Salamandra species and most of the commonly recognized subspecies. The concatenated analysis of the complete dataset produced a fully resolved tree with most nodes strongly supported, suggesting that a clade composed of the Alpine salamander (S. atra) and the Corsican fire salamander (S. corsica) is the sister taxon to a clade containing the remaining species, among which S. algira and S. salamandra are sister species. Separate analyses of mitochondrial and nuclear data partitions disagreed regarding basal nodes and in the position of the root but concordantly recovered the S. atra/S. corsica as well as the S. salamandra/S. algira relationship. A species-tree analysis suggested almost simultaneous temporal splits between these pairs of species, which we hypothesize was caused by vicariance events after the Messinian salinity crisis (from late Miocene to early Pliocene). A survey of toxins with combined gas chromatography/mass spectroscopy confirmed the presence of samandarine and/or samandarone steroidal alkaloids in all species of Salamandra as well as in representatives of their sister group, Lyciasalamandra. Samandarone was also detected in lower concentrations in other salamandrids including Calotriton, Euproctus, Lissotriton, and Triturus, suggesting that the presence and possible biosynthesis of this alkaloid is plesiomorphic within the Salamandridae.

The more the better - polyandry and genetic similarity are positively linked to reproductive success in a natural population of terrestrial salamanders (Salamandra salamandra).

Although classically thought to be rare, female polyandry is widespread and may entail significant fitness benefits. If females store sperm over extended periods of time, the consequences of polyandry will depend on the pattern of sperm storage, and some of the potential benefits of polyandry can only be realized if sperm from different males is mixed. Our study aimed to determine patterns and consequences of polyandry in an amphibian species, the fire salamander, under fully natural conditions. Fire salamanders are ideal study objects, because mating, fertilization and larval deposition are temporally decoupled, females store sperm for several months, and larvae are deposited in the order of fertilization. Based on 18 microsatellite loci, we conducted paternity analysis of 24 female-offspring arrays with, in total, over 600 larvae fertilized under complete natural conditions. More than one-third of females were polyandrous and up to four males were found as sires. Our data clearly show that sperm from multiple males is mixed in the female's spermatheca. Nevertheless, paternity is biased, and the most successful male sires on average 70% of the larvae, suggesting a 'topping off' mechanism with first-male precedence. Female reproductive success increased with the number of sires, most probably because multiple mating ensured high fertilization success. In contrast, offspring number was unaffected by female condition and genetic characteristics, but surprisingly, it increased with the degree of genetic relatedness between females and their sires. Sires of polyandrous females tended to be genetically similar to each other, indicating a role for active female choice.

Sex-specific estimates of dispersal show female philopatry and male dispersal in a promiscuous amphibian, the alpine salamander (Salamandra atra).

Amphibians display wide variations in life-history traits and life cycles that should prove useful to explore the evolution of sex-biased dispersal, but quantitative data on sex-specific dispersal patterns are scarce. Here, we focused on Salamandra atra, an endemic alpine species showing peculiar life-history traits. Strictly terrestrial and viviparous, the species has a promiscuous mating system, and females reproduce only every 3 to 4 years. In the present study, we provide quantitative estimates of asymmetries in male vs. female dispersal using both field-based (mark-recapture) and genetic approaches (detection of sex-biased dispersal and estimates of migration rates based on the contrast in genetic structure across sexes and age classes). Our results revealed a high level of gene flow among populations, which stems exclusively from male dispersal. We hypothesize that philopatric females benefit from being familiar with their natal area for the acquisition and defence of an appropriate shelter, while male dispersal has been secondarily favoured by inbreeding avoidance. Together with other studies on amphibians, our results indicate that a species' mating system alone is a poor predictor of sex-linked differences in dispersal, in particular for promiscuous species. Further studies should focus more directly on the proximate forces that favour or limit dispersal to refine our understanding of the evolution of sex-biased dispersal in animals.

Phylogeny, life history, and ecology contribute to differences in amphibian susceptibility to ranaviruses.

Research that identifies the potential host range of generalist pathogens as well as variation in host susceptibility is critical for understanding and predicting the dynamics of infectious diseases within ecological communities. Ranaviruses have been linked to amphibian die-off events worldwide with the greatest number of reported mortality events occurring in the United States. While reports of ranavirus-associated mortality events continue to accumulate, few data exist comparing the relative susceptibility of different species. Using a series of laboratory exposure experiments and comparative phylogenetics, we compared the susceptibilities of 19 amphibian species from two salamander families and five anurans families for two ranavirus isolates: frog virus 3 (FV3) and an FV3-like isolate from an American bullfrog culture facility. We discovered that ranaviruses were capable of infecting 17 of the 19 larval amphibian species tested with mortality ranging from 0 to 100%. Phylogenetic comparative methods demonstrated that species within the anuran family Ranidae were generally more susceptible to ranavirus infection compared to species from the other five families. We also found that susceptibility to infection was associated with species that breed in semi-permanent ponds, develop rapidly as larvae, and have limited range sizes. Collectively, these results suggest that phylogeny, life history characteristics, and habitat associations of amphibians have the potential to impact susceptibility to ranaviruses.

Salamandra salamandra (Amphibia: Caudata: Salamandridae) in Portugal: not all black and yellow.

The fire salamander complex is quite diverse in the Iberian Peninsula where nine subspecies of Salamandra salamandra are currently recognized. Here, we analysed the geographical distribution of the subspecies S. s. gallaica and S. s. crespoi using partial sequences of the mitochondrial cytochrome b gene of 168 individuals from 12 locations in Portugal. Our results support the existence of a deep lineage divergence between the two subspecies, with non-overlapping geographical distributions except in two contact zones: one in Sesimbra on the western coast, and another in Alcoutim on the southeastern border with Spain. Moreover, S. s. crespoi displays signs of gene flow among the sampled locations whereas S. s. gallaica shows evidence of some restriction to gene flow. Present-day genetic make-up of S. s. gallaica and S. s. crespoi is a result of past historical events, fine-tuned by contemporary Iberian geoclimate. Humid mountain areas were found to harbour increased genetic diversity possibly acting as past refugia during drier interglacial periods. To analyse wider geographical patterns and lineage splitting events within S. salamandra we performed a Bayesian dating analysis completing our data set with previously published sequences. The observed divergences were associated to successive biogeographic scenarios, and to other Iberian species showing similar trends.

Tracing the first step to speciation: ecological and genetic differentiation of a salamander population in a small forest.

Mechanisms and processes of ecologically driven adaptive speciation are best studied in natural situations where the splitting process is still occurring, i.e. before complete reproductive isolation is achieved. Here, we present a case of an early stage of adaptive differentiation under sympatric conditions in the fire salamander, Salamandra salamandra, that allows inferring the underlying processes for the split. Larvae of S. salamandra normally mature in small streams until metamorphosis, but in an old, continuous forest area near Bonn (the Kottenforst), we found salamander larvae not only in small streams but also in shallow ponds, which are ecologically very different from small streams. Common-environment experiments with larvae from both habitat types reveal specific adaptations to these different ecological conditions. Mitochondrial and microsatellite analyses show that the two ecologically differentiated groups also show signs of genetic differentiation. A parallel analysis of animals from a neighbouring much larger forest area (the Eifel), in which larvae mature only in streams, shows no signs of genetic differentiation, indicating that gene flow between ecologically similar types can occur over large distances. Hence, geographical factors cannot explain the differential larval habitat adaptations in the Kottenforst, in particular since adult life and mating of S. salamandra is strictly terrestrial and not associated with larval habitats. We propose therefore that the evolution of these adaptations was coupled with the evolution of cues for assortative mating which would be in line with models of sympatric speciation that suggest a co-evolution of habitat adaptations and associated mating signals.

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.

Representatives of the genus Salamandra occur in Europe, Northern Africa and the Near East. Many local variants are known but species and subspecies status of these is still a matter of dispute. We have analysed samples from locations covering the whole expansion range of Salamandra by sequence analysis of mitochondrial D-loop regions. In addition, we have calibrated the rate of divergence of the D-loop on the basis of geologically dated splits of the closely related genus Euproctus. Phylogenetic analysis of the sequences suggests that six major monophyletic groups exist (S. salamandra, S. algira, S. infraimmaculata, S. corsica, S. atra and S. lanzai) which have split between 5 and 13 million years ago (Ma). We find that each of the Salamandra species occupies a distinct geographical area, with the exception of S. salamandra. This species occurs all over Europe from Spain to Greece, suggesting that it was the only species that has recolonized Central Europe after the last glaciation. The occurrence of specific east and west European haplotypes, as well as allozyme alleles in the S. salamandra populations suggests that this recolonization has started from at least two source populations, possibly originating in the Iberian peninsula and the Balkans. Two subpopulations of S. salamandra were found that are genetically very distinct from the other populations. One lives in northern Spain (S. s. bernardezi) and one in southern Italy (S. s. gigliolii). Surprisingly, the mitochondrial lineages of these subpopulations group closer together than the remainder S. salamandra lineages. We suggest that these populations are remnants of a large homogeneous population that had colonized Central Europe in a previous interglacial period, approximately 500 000 years ago. Animals from these populations were apparently not successful in later recolonizations. Still, they have maintained their separate genetic identity in their areas, although they are not separated by geographical barriers from very closely related neighbouring populations.

Gypsy/Ty3-like elements in the genome of the terrestrial Salamander hydromantes (Amphibia, Urodela).

We have studied a family of long repetitive DNA sequences (Hsr1) interspersed in the large genome of the European plethodontid salamander Hydromantes. The sequence analysis of a 5-kb fragment (Hsr1A) of one member has revealed significant similarities with amino acidic domains of retroviruses and retrotransposons. The similarity of the reverse transcriptase domain and the gene organization identifies Hsr1A as a member of the gypsy/Ty3 class of retrotransposons. We hypothesize that Hsr1 sequences are vestiges of an invasion of the Hydromantes genome that occurred early in the evolutionary history of these European plethodontids. About 10(6) Hsr1 sequences are present in the large Hydromantes genome. This is the highest number of copies so far discovered for retrotransposon-like elements in eukaryote organisms.