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.

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.

Intraspecific variation in lizard heat tolerance alters estimates of climate impact.

Research addressing the effects of global warming on the distribution and persistence of species generally assumes that population variation in thermal tolerance is spatially constant or overridden by interspecific variation. Typically, this rationale is implicit in sourcing one critical thermal maximum (CTmax ) population estimate per species to model spatiotemporal cross-taxa variation in heat tolerance. Theory suggests that such an approach could result in biased or imprecise estimates and forecasts of impact from climate warming, but limited empirical evidence in support of those expectations exists. We experimentally quantify the magnitude of intraspecific variation in CTmax among lizard populations, and the extent to which incorporating such variability can alter estimates of climate impact through a biophysical model. To do so, we measured CTmax from 59 populations of 15 Iberian lizard species (304 individuals). The overall median CTmax across all individuals from all species was 42.8°C and ranged from 40.5 to 48.3°C, with species medians decreasing through xeric, climate-generalist and mesic taxa. We found strong statistical support for intraspecific differentiation in CTmax by up to a median of 3°C among populations. We show that annual restricted activity (operative temperature > CTmax ) over the Iberian distribution of our study species differs by a median of >80 hr per 25-km2 grid cell based on different population-level CTmax estimates. This discrepancy leads to predictions of spatial variation in annual restricted activity to change by more than 20 days for six of the study species. Considering that during restriction periods, reptiles should be unable to feed and reproduce, current projections of climate-change impacts on the fitness of ectotherm fauna could be under- or over-estimated depending on which population is chosen to represent the physiological spectra of the species in question. Mapping heat tolerance over the full geographical ranges of single species is thus critical to address cross-taxa patterns and drivers of heat tolerance in a biologically comprehensive way.

Individual and Geographic Variation of Skin Alkaloids in Three Swamp-Forest Species of Madagascan Poison Frogs (Mantella).

Seventy skins of three mantellid frog species from Madagascan swamp-forest habitats, Mantella aurantiaca, M. crocea, and M. milotympanum, were individually examined for skin alkaloids using GC/MS. These poison frogs were found to differ significantly in their alkaloid composition from species of Mantella originating from non-flooded rainforest in eastern Madagascar, which were examined in earlier work. Only 16 of the previously detected 106 alkaloids were represented among the 60 alkaloids from the swamp-forest frogs of the present study. We hypothesize this difference is related mainly to habitat but cannot exclude a phylogenetic component as the three swamp-forest species are a closely related monophyletic group. The paucity of alkaloids with unbranched-carbon skeletons (ant-derived) and the commonness of alkaloids with branched-carbon skeletons (mite-derived) indicate that oribatid mites are a major source of alkaloids in these species of mantellids. Furthermore, most of the alkaloids have an oxygen atom in their formulae. Differences in alkaloids were observed among species, populations of the same species, and habitats. In M. aurantiaca, small geographic distances among populations were associated with differences in alkaloid profiles, with a remote third site illustrating even greater differences. The present study and an earlier study of three other mantellid species suggest that oribatid mites, and not ants, are the major source of alkaloids in the species of mantellids examined thus far.

Vertebrate time-tree elucidates the biogeographic pattern of a major biotic change around the K-T boundary in Madagascar.

The geographic and temporal origins of Madagascar's biota have long been in the center of debate. We reconstructed a time-tree including nearly all native nonflying and nonmarine vertebrate clades present on the island, from DNA sequences of two single-copy protein-coding nuclear genes (BDNF and RAG1) and a set of congruent time constraints. Reconstructions calculated with autocorrelated or independent substitution rates over clades agreed in placing the origins of the 31 included clades in Cretaceous to Cenozoic times. The two clades with sister groups in South America were the oldest, followed by those of a putative Asian ancestry that were significantly older than the prevalent clades of African ancestry. No colonizations from Asia occurred after the Eocene, suggesting that dispersal and vicariance of Asian/Indian groups were favored over a comparatively short period during, and shortly after, the separation of India and Madagascar. Species richness of clades correlates with their age but those clades that have a large proportion of species diversity in rainforests are significantly more species-rich. This finding suggests an underlying pattern of continuous speciation through time in Madagascar's vertebrates, with accelerated episodes of adaptive diversification in those clades that succeeded radiating into the rainforests.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana.

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

The complete mitochondrial genome of the Critically Endangered Saba Green Iguana, Iguana iguana (Squamata: Iguanidae).

The populations of native iguanas in the Caribbean Lesser Antilles are threatened by the wide occurrence and spread of non-native iguanas. Until recently, competitive hybridization was not believed to threaten the Saba Green Iguana, a subpopulation of Iguana iguana (Linnaeus, 1758) from the island of Saba. However, the arrival of non-native iguanas has put the native population at risk, leading to a change in the conservation status of the Saba Green Iguana to Critically Endangered, according to guidelines from the International Union for the Conservation of Nature. Here, we generated the complete mitogenome of the Saba Green Iguana using Oxford Nanopore long-read technology. The mitogenome is 16,626 bp long and has 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a control region (1194 bp). Noteworthy, this is only the second published mitogenome for the Iguana iguana species complex, despite the known high intraspecific genetic variation.

Molecular phylogeny and biogeography of Malagasy frogs of the genus Gephyromantis.

The genus Gephyromantis is a clade within the Malagasy-Comoroan family Mantellidae composed of rainforest frogs that live and breed to varying degrees independently from water. Based on DNA sequences of five mitochondrial and five nuclear genes we inferred the phylogeny of these frogs with full taxon coverage at the species level. Our preferred consensus tree from a partitioned Bayesian analysis of 5843 base pairs of 51 nominal and candidate species supports various major clades within the genus although the basal relationships among these remain unresolved. The data provide strong evidence for the monophyly of the subgenera Gephyromantis (after exclusion of Gephyromantis klemmeri), Laurentomantis, Vatomantis, and Phylacomantis. Species assigned to the subgenus Duboimantis belong to two strongly supported clades of uncertain relationships. G. klemmeri, previously in the subgenus Gephyromantis, was placed with high support sister to the Laurentomantis clade, and the Laurentomantis+G. klemmeri clade was sister to Vatomantis. A reconstruction of ancestral distribution areas indicates a diversification of several subgenera in the northern biogeographic regions of Madagascar and the dispersal out of northern Madagascar for several clades.

Discordant patterns of nuclear and mitochondrial introgression in Iberian populations of the European common frog (Rana temporaria).

Amphibians often show complex histories of intraspecific and interspecific genetic introgression, which might differ in mitochondrial and nuclear genes. In our study of the genetic differentiation of the European common frog, Rana temporaria (159 specimens from 23 populations were analyzed for 24 presumptive allozyme loci; 82 specimens were sequenced for a 540-bp fragment of the mitochondrial 16S rRNA gene), multilocus correspondence analysis (CA) and Bayesian assignment tests of the nuclear data were concordant in identifying 2 population groups corresponding to 1) the Pyrenees in the east and 2) the Galicia and Asturias regions in the west, the latter corresponding to the subspecies R. temporaria parvipalmata. Geographically intermediate populations were genetically intermediate in the allozyme CA and, less clearly in the Bayesian assignment, with mitochondrial haplotypes exclusively belonging to the parvipalmata group. This indicates different degrees of introgression in the mitochondrial and nuclear genomes. Although Pyrenean high-altitude populations are morphologically distinct from low-altitude populations, these 2 groups were not separate clusters in any analysis. This suggests that the morphological differences may be due to fast adaptation to elevational gradients, likely under maintenance of gene flow, and that the underlying genetic changes are not detectable by the analyzed markers. We argue that a parsimonious explanation for the observed pattern along the east-west axis in northern Spain may be competition between invading and resident populations, with no need to invoke selection. However, in order to conclusively rule out selective processes, additional and finer scale data are required to test for asymmetric mating preference/behaviour, sex-biased gene flow, or sex-biased survival of potential hybrids.

Reconstruction of the climate envelopes of salamanders and their evolution through time.

Ecological niches evolve through time, but at different rates and to different degrees. An integrated approach using diverse databases, methods, and analytical tools is used to estimate climate envelopes for species of salamanders (family Salamandridae). These species, which range widely across the Holarctic and have a rich and long fossil record, are used to probe the evolutionary dynamics of niches studied in a phylogenetic context through time and across space. Climate data and statistical methods are used to estimate niche dimensions related to precipitation and temperature, in both a phylogenetic and ecogeographic context. Using phylogenetic methods, climate envelopes are estimated for segments of a time-calibrated phylogenetic tree of salamandrids, to explore how far back in time it is reasonable to make such estimates. Our research strategy illuminates some of the challenges and limitations of the available data and methods and identifies potential solutions, such as the need for physiological and behavioral data that may help to better define species' niches or the development of novel evolutionary models that account for paleoclimatic data. We explore and analyze limits to the application of currently available methodologies. Organisms known to have evolved slowly and conservatively, and which are ectotherms and likely to be profoundly affected by climatic variables, may be the most useful for studies of this kind.

Vast underestimation of Madagascar's biodiversity evidenced by an integrative amphibian inventory.

Amphibians are in decline worldwide. However, their patterns of diversity, especially in the tropics, are not well understood, mainly because of incomplete information on taxonomy and distribution. We assess morphological, bioacoustic, and genetic variation of Madagascar's amphibians, one of the first near-complete taxon samplings from a biodiversity hotspot. Based on DNA sequences of 2,850 specimens sampled from over 170 localities, our analyses reveal an extreme proportion of amphibian diversity, projecting an almost 2-fold increase in species numbers from the currently described 244 species to a minimum of 373 and up to 465. This diversity is widespread geographically and across most major phylogenetic lineages except in a few previously well-studied genera, and is not restricted to morphologically cryptic clades. We classify the genealogical lineages in confirmed and unconfirmed candidate species or deeply divergent conspecific lineages based on concordance of genetic divergences with other characters. This integrative approach may be widely applicable to improve estimates of organismal diversity. Our results suggest that in Madagascar the spatial pattern of amphibian richness and endemism must be revisited, and current habitat destruction may be affecting more species than previously thought, in amphibians as well as in other animal groups. This case study suggests that worldwide tropical amphibian diversity is probably underestimated at an unprecedented level and stresses the need for integrated taxonomic surveys as a basis for prioritizing conservation efforts within biodiversity hotspots.

Speciation in little: the role of range and body size in the diversification of Malagasy mantellid frogs.

BACKGROUND: The rate and mode of lineage diversification might be shaped by clade-specific traits. In Madagascar, many groups of organisms are characterized by tiny distribution ranges and small body sizes, and this high degree of microendemism and miniaturization parallels a high species diversity in some of these groups. We here investigate the geographic patterns characterizing the radiation of the frog family Mantellidae that is virtually endemic to Madagascar. We integrate a newly reconstructed near-complete species-level timetree of the Mantellidae with georeferenced distribution records and maximum male body size data to infer the influence of these life-history traits on each other and on mantellid diversification. RESULTS: We reconstructed a molecular phylogeny based on nuclear and mitochondrial DNA for 257 species and candidate species of the mantellid frog radiation. Based on this phylogeny we identified 53 well-supported pairs of sister species that we used for phylogenetic comparative analyses, along with whole tree-based phylogenetic comparative methods. Sister species within the Mantellidae diverged at 0.2-14.4 million years ago and more recently diverged sister species had geographical range centroids more proximate to each other, independently of their current sympatric or allopatric occurrence. The largest number of sister species pairs had non-overlapping ranges, but several examples of young microendemic sister species occurring in full sympatry suggest the possibility of non-allopatric speciation. Range sizes of species included in the sister species comparisons increased with evolutionary age, as did range size differences between sister species, which rejects peripatric speciation. For the majority of mantellid sister species and the whole mantellid radiation, range and body sizes were associated with each other and small body sizes were linked to higher mitochondrial nucleotide substitution rates and higher clade diversity. In contrast, small range sizes were unexpectedly associated with a slow-down of mitochondrial substitution rates. CONCLUSIONS: Based on these results we define a testable hypothesis under which small body sizes result in limited dispersal capabilities and low physiological tolerances, causing smaller and more strongly fragmented ranges. This can be thought to facilitate reproductive isolation and thus favor speciation. Contrary to the expectation of the faster speciation of such microendemic phenotype species, we only found small body sizes of mantellid frogs to be linked to higher diversification and substitution rates, but not small range sizes. A joint analysis of various species-rich regional anuran radiations might provide enough species with all combinations of range and body sizes for a more conclusive test of this hypothesis.

A multigenic perspective on phylogenetic relationships in the largest family of salamanders, the Plethodontidae.

Despite several recent studies, the phylogeny of plethodontid salamanders is not yet fully resolved and the phylogenetic positions of several key genera, especially Aneides, Hemidactylium, Hydromantes and Karsenia, are contentious. Here we present a combined dataset of complete mitochondrial genomes and three nuclear loci for 20 species (16 genera) of plethodontids, representing all major clades in the family. The combined dataset without mitochondrial third codon positions provides a fully resolved, statistically well-supported tree. In this topology two major clades are recovered. A northern clade includes Aneides, Desmognathus, Ensatina, Hydromantes, Karsenia, Phaeognathus and Plethodon, with Plethodon being the sister taxon to the rest of the clade. Hydromantes and Karsenia are sister taxa, and Aneides is recovered as the sister taxon to Ensatina. Desmognathus+Phaeognathus form the sister taxon to Aneides+Ensatina. An eastern/southern clade comprises two subclades. One subclade, the spelerpines (Eurycea, Gyrinophilus, Pseudotriton, Stereochilus, Urspelerpes) is the sister taxon to a subclade comprising Hemidactylium, Batrachoseps and the tropical plethodontids (represented by Bolitoglossa, Nototriton and Thorius). In this topology Hemidactylium is well-supported as the sister taxon to Batrachoseps. Only when mitochondrial third codon positions are included using maximum likelihood analysis is Hemidactylium recovered as the sister taxon to Batrachoseps+tropical genera. Hypothesis testing of alternative topologies supports these conclusions. On the basis of these results we propose a conservative taxonomy for Plethodontidae.

Deep genealogical lineages in the widely distributed African helmeted terrapin: evidence from mitochondrial and nuclear DNA (Testudines: Pelomedusidae: Pelomedusa subrufa).

We investigated the phylogeographic differentiation of the widely distributed African helmeted terrapin Pelomedusa subrufa based on 1503 base pairs of mitochondrial DNA (partial cyt b and ND4 genes with adjacent tRNAs) and 1937 bp of nuclear DNA (partial Rag1, Rag2, R35 genes). Congruent among different analyses, nine strongly divergent mitochondrial clades were found, representing three major geographical groupings: (1) A northern group which includes clades I from Cameroon, II from Ghana and Ivory Coast, III from Benin, Burkina Faso and Niger, IV from the Central African Republic, and V from Kenya, (2) a northeastern group consisting of clades VI from Somalia, and VII from Saudi Arabia and Yemen, and (3) a southern group comprising clade VIII from Botswana, the Democratic Republic of Congo, Madagascar and Malawi, and clade IX from South Africa. Malagasy and continental African populations were not clearly differentiated, indicating very recent arrival or introduction of Pelomedusa in Madagascar. The southern group was in some phylogenetic analyses sister to Pelusios, rendering Pelomedusa paraphyletic with respect to that genus. However, using partitioned Bayesian analyses and sequence data of the three nuclear genes, Pelomedusa was monophyletic, suggesting that its mitochondrial paraphyly is due to either ancient introgressive hybridization or phylogenetic noise. Otherwise, nuclear sequence data recovered a lower level of divergence, but corroborated the general differentiation pattern of Pelomedusa as revealed by mtDNA. This, and the depth of the divergences between clades, indicates ancient differentiation. The divergences observed fall within, and in part exceed considerably, the differentiation typically occurring among chelonian species. To test whether Pelomedusa is best considered a single species composed of deep genealogical lineages, or a complex of up to nine distinct species, we suggest a future taxonomic revision that should (1) extend the geographical sampling of molecular data, specifically focusing on contact zones and the possible sympatric occurrence of lineages without admixture, and (2) evaluate the morphology of the various genealogical lineages using the type specimens or topotypical material of the numerous junior synonyms of P. subrufa.

Testing species-level diversification hypotheses in Madagascar: the case of microendemic Brookesia leaf chameleons.

Madagascar's flora and fauna are remarkable both for their diversity and supraspecific endemism. Moreover, many taxa contain large numbers of species with limited distributions. Several hypotheses have been proposed to explain this high level of microendemism, including 1) riverine barrier, 2) mountain refuge, and 3) watershed contraction hypotheses, the latter 2 of which center on fragmentation due to climatic shifts associated with Pliocene/Pleistocene glaciations. The Malagasy leaf chameleon genus Brookesia is a speciose group with a high proportion of microendemic taxa, thus making it an excellent candidate to test these vicariance scenarios. We used mitochondrial and nuclear sequence data to construct a Brookesia phylogeny, and temporal concordance with Pliocene/Pleistocene speciation scenarios was tested by estimating divergence dates using a relaxed-clock Bayesian method. We strongly reject a role for Pliocene/Pleistocene climatic fluctuations in species-level diversification of Brookesia. We also used simulations to test the spatial predictions of the watershed contraction model in a phylogenetic context, independent of its temporal component, and found no statistical support for this model. The riverine barrier model is likewise a qualitatively poor fit to our data, but some relationships support a more ancient mountain refuge effect. We assessed support for the 3 hypotheses in a nonphylogenetic context by examining altitude and species richness and found a significant positive correlation between these variables. This is consistent with a mountain refuge effect but does not support the watershed contraction or riverine barrier models. Finally, we find repeated higher level east-west divergence patterns 1) between the 2 sister clades comprising the Brookesia minima group and 2) within the clade of larger leaf chameleons, which shows a basal divergence between western and eastern/northern sister clades. Our results highlight the central role of phylogeny in any meaningful tests of species-level diversification theories.

Post-mating clutch piracy in an amphibian.

Female multiple mating and alternative mating systems can decrease the opportunity for sexual selection. Sperm competition is often the outcome of females mating with multiple males and has been observed in many animals, and alternative reproductive systems are widespread among species with external fertilization and parental care. Multiple paternity without associated complex behaviour related to mating or parental care is also seen in simultaneously spawning amphibians and fishes that release gametes into water. Here we report 'clutch piracy' in a montane population of the common frog Rana temporaria, a reproductive behaviour previously unknown in vertebrates with external fertilization. Males of this species clasp the females and the pair deposits one spherical clutch of eggs. No parental care is provided. 'Pirate' males search for freshly laid clutches, clasp them as they would do a female and fertilize the eggs that were left unfertilized by the 'parental' male. This behaviour does not seem to be size-dependent, and some males mate with a female and perform clutch piracy in the same season. Piracy affected 84% of the clutches and in some cases increased the proportion of eggs fertilized, providing direct fitness benefits both for the pirate males and the females. Sexual selection--probably caused by a strong male-biased sex ratio--occurs in this population, as indicated by size-assortative mating; however, clutch piracy may reduce its impact. This provides a good model to explore how alternative mating strategies can affect the intensity of sexual selection.

Rapid diversification and dispersal during periods of global warming by plethodontid salamanders.

A phylogeny and timescale derived from analyses of multilocus nuclear DNA sequences for Holarctic genera of plethodontid salamanders reveal them to be an old radiation whose common ancestor diverged from sister taxa in the late Jurassic and underwent rapid diversification during the late Cretaceous. A North American origin of plethodontids was followed by a continental-wide diversification, not necessarily centered only in the Appalachian region. The colonization of Eurasia by plethodontids most likely occurred once, by dispersal during the late Cretaceous. Subsequent diversification in Asia led to the origin of Hydromantes and Karsenia, with the former then dispersing both to Europe and back to North America. Salamanders underwent rapid episodes of diversification and dispersal that coincided with major global warming events during the late Cretaceous and again during the Paleocene-Eocene thermal optimum. The major clades of plethodontids were established during these episodes, contemporaneously with similar phenomena in angiosperms, arthropods, birds, and mammals. Periods of global warming may have promoted diversification and both inter- and transcontinental dispersal in northern hemisphere salamanders by making available terrain that shortened dispersal routes and offered new opportunities for adaptive and vicariant evolution.

Challenges in estimating ancestral state reconstructions: the evolution of migration in Sylvia warblers as a study case.

Our current understanding of how species have evolved is mainly based on comparative phylogenetic methods, which use phylogenies to infer the evolution of traits. The development of ancestral state reconstruction (ASR) methods has provided the tools to reconstruct trait evolution, which are widely used in fields like evolutionary biology, macroecology and paleontology. As there are different elements involved in those analyses, with different levels of uncertainty (i.e. relating to branch length estimation, trait coding, statistical framework, taxon sampling or software), the various combinations of these elements likely have a strong impact on the reconstruction of the evolution of traits, potentially leading to opposite conclusions. To assess the impact of these different elements in ASR, we performed a set of analyses, including all possible combinations of such elements and using the evolution of migratory behavior in Sylvia warblers as a case study, which was coded as a continuous or as a discrete character. Our results show that taxon sampling, character coding, tree shape, statistical framework and software all significantly affect ASR, both individually and in combination. Not all reconstructed tree nodes show discrepancies, but in the critical ones most pairwise comparisons between analyses lead to conflicting and unexpectedly antagonistic results (zero migration vs fully migratory), thus challenging interpretations of trait evolution. We propose some possible solutions to partly inform decisions, involving the method selection and the incorporation of biological or fossil evidence regarding how traits evolve, but our results demand serious rethinking about how the research community currently uses ASR.

ACDC, a global database of amphibian cytochrome-b sequences using reproducible curation for GenBank records.

Genetic data are a crucial and exponentially growing resource across all biological sciences, yet curated databases are scarce. The widespread occurrence of sequence and (meta)data errors in public repositories calls for comprehensive improvements of curation protocols leading to robust research and downstream analyses. We collated and curated all available GenBank cytochrome-b sequences for amphibians, a benchmark marker in this globally declining vertebrate clade. The Amphibia's Curated Database of Cytochrome-b (ACDC) consists of 36,514 sequences representing 2,309 species from 398 genera (median = 2 with 50% interquartile ranges of 1-7 species/genus). We updated the taxonomic identity of >4,800 sequences (ca. 13%) and found 2,359 (6%) conflicting sequences with 84% of the errors originating from taxonomic misidentifications. The database (accessible at https://doi.org/10.6084/m9.figshare.9944759 ) also includes an R script to replicate our study for other loci and taxonomic groups. We provide recommendations to improve genetic-data quality in public repositories and flag species for which there is a need for taxonomic refinement in the face of increased rate of amphibian extinctions in the Anthropocene.

Hidden in plain sight: a new frog species of the genus Blommersia from the oceanic island of Mayotte, Comoros archipelago.

The amphibian fauna of the western Indian ocean volcanic island of Mayotte is currently constituted by two species belonging to two genera of the anuran family Mantellidae: Blommersia transmarina and Boophis nauticus. These were recently described after intense fieldwork on the herpetofauna of the island. We here describe a third new species of frog from Mayotte, based on morphological and molecular data, that occurs in sympatry with the others and was utterly unnoticed until now. Genetic analyses of the16S rRNA gene, including all described and several undescribed species of the genus Blommersia from Madagascar and Mayotte, confirms that the new species is the sister species of Blommersia transmarina. Both species show apparent morphological differences as well as different life histories, ecology and genetics that confirm Blommersia nataliae sp. nov. as a new species. We propose an IUCN Red List status of Critically Endangered for B. nataliae sp. nov.