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.

No abstract

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. | 1285BURGON et al.

No abstract

No abstract