Conservation phylogeography: does historical diversity contribute to regional vulnerability in European tree frogs (Hyla arborea)?

The introduction of three cryptic tree frog species in the Dutch coastal dunes challenges conservation paradigms

Invasive species are considered one of the main drivers of the sixth mass extinction. Conservation solutions depend on whether a species is also indigenous to the country it invades (i.e., beyond its native range). In the case of invasive cryptic species, genetic tools are required to establish their identity. We illustrate these issues with the human-mediated colonization of the Dutch coastal dunes by Hyla tree frogs. Although previously assumed to concern the indigenous common tree frog H. arborea, European tree frogs comprise a complex of allopatric cryptic species, meaning the taxonomic identity of introduced Dutch populations warrants investigation. We sequence mtDNA for 173 individuals from native and introduced populations across the Netherlands and compare our dataset with hundreds of Hyla haplotypes previously barcoded in the Western Palearctic. Two of the dune populations carry an mtDNA haplotype of the native species H. arborea that occurs naturally elsewhere in the Netherlands. In contrast, mtDNA assigned to the eastern tree frog H. orientalis was detected in all three other dune populations. In one of these populations mtDNA of the Italian tree frog H. intermedia was also found. Not one, but three species of tree frogs have thus been introduced to the Dutch coastal dunes, only one of them being native to the Netherlands. This situation causes a conservation conundrum as some introduced populations are lawfully protected but could pose a threat to local biodiversity. Regarding the ‘true’ exotic tree frog species, all conservation options should be considered.

Morphometry of two cryptic tree frog species at their hybrid zone reveals neither intermediate nor transgressive morphotypes

Under incomplete reproductive isolation, secondary contact of diverged allopatric lineages may lead to the formation of hybrid zones that allow to study recombinants over several generations as excellent systems of genomic interactions resulting from the evolutionary forces acting on certain genes and phenotypes. Hybrid phenotypes may either exhibit intermediacy or, alternatively, transgressive traits, which exceed the extremes of their parents due to epistasis and segregation of complementary alleles. While transgressive morphotypes have been examined in fish, reptiles, birds, and mammals, studies in amphibians are rare. Here, we associate microsatellite-based genotypes with morphometrics-based morphotypes of two tree frog species of the Hyla arborea group, sampled across a hybrid zone in Poland, to understand whether the genetically differentiated parental species also differ in morphology between each other and their hybrids and whether secondary contact leads to the evolution of intermediate or transgressive morphotypes. Using univariate approaches, explorative multivariate methods (principal component analyses) as well as techniques with prior grouping (discriminant function analyses), we find that morphotypes of both parental species and hybrids differ from each other. Importantly, hybrid morphotypes are neither intermediate nor transgressive but found to be more similar to H. orientalis than to H. arborea.

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure, genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Performance comparison of two reduced-representation based genome-wide marker-discovery strategies in a multi-taxon phylogeographic framework

Multi-locus genetic data are pivotal in phylogenetics. Today, high-throughput sequencing (HTS) allows scientists to generate an unprecedented amount of such data from any organism. However, HTS is resource intense and may not be accessible to wide parts of the scientific community. In phylogeography, the use of HTS has concentrated on a few taxonomic groups, and the amount of data used to resolve a phylogeographic pattern often seems arbitrary. We explore the performance of two genetic marker sampling strategies and the effect of marker quantity in a comparative phylogeographic framework focusing on six species (arthropods and plants). The same analyses were applied to data inferred from amplified fragment length polymorphism fingerprinting (AFLP), a cheap, non-HTS based technique that is able to straightforwardly produce several hundred markers, and from restriction site associated DNA sequencing (RADseq), a more expensive, HTS-based technique that produces thousands of single nucleotide polymorphisms. We show that in four of six study species, AFLP leads to results comparable with those of RADseq. While we do not aim to contest the advantages of HTS techniques, we also show that AFLP is a robust technique to delimit evolutionary entities in both plants and animals. The demonstrated similarity of results from the two techniques also strengthens biological conclusions that were based on AFLP data in the past, an important finding given the wide utilization of AFLP over the last decades. We emphasize that whenever the delimitation of evolutionary entities is the central goal, as it is in many fields of biodiversity research, AFLP is still an adequate technique.

Conservation genetics of yellow-bellied toads (Bombina variegata): a matter of geographical scale and isolation

Amphibian populations world-wide are threatened by declines and extinctions mainly due to habitat loss and fragmentation. Habitat fragmentation threatens the yellow-bellied toad Bombina variegata in the northern and western regions of its distribution where it is strictly protected. We studied the genetic structure and diversity of populations at three geographical scales using microsatellite loci to detect potential threats for population persistence. At the local scale, we sampled four neighbouring localities at 1–2.6 km distance to detect effects of short-term (decades) fragmentation on connectivity. At the regional scale, five additional localities in the mountains of the Westerwald (Rhineland-Palatinate, Germany) were studied at up to 50.1 km distance to analyse genetic diversity and population structure. At the continental scale, we included data from regions in the northern distribution with fragmented populations (Hesse and Lower Saxony, Germany) and more continuous populations in the South (Alsace, France; Geneva, Switzerland; Trentino, Italy) to evaluate variation of genetic diversity. At the local scale, short-term fragmentation caused significant genetic differentiation between breeding assemblages only 1.4 km apart from each other. At the regional scale, we found notable genetic distance among localities. At the continental scale, we identified Alsace, Trentino and Geneva in the South as regions with low genetic structuring and high allelic richness, and the northern remaining regions in Germany as deeply structured with reduced allelic richness. We suggest that reduced genetic diversity and habitat fragmentation in northern regions makes these populations particularly vulnerable to decline. In conclusion, informed conservation management of B. variegata should focus on measures maintaining or improving connectivity among neighbouring populations.

Mitochondrial genomics reveals the evolutionary history of the porpoises (Phocoenidae) across the speciation continuum

Historical variation in food resources is expected to be a major driver of cetacean evolution, especially for the smallest species like porpoises. Despite major conservation issues among porpoise species (e.g., vaquita and finless), their evolutionary history remains understudied. Here, we reconstructed their evolutionary history across the speciation continuum. Phylogenetic analyses of 63 mitochondrial genomes suggest that porpoises radiated during the deep environmental changes of the Pliocene. However, all intra-specific subdivisions were shaped during the Quaternary glaciations. We observed analogous evolutionary patterns in both hemispheres associated with convergent evolution to coastal versus oceanic environments. This suggests that similar mechanisms are driving species diversification in northern (harbor and Dall's) and southern species (spectacled and Burmeister's). In contrast to previous studies, spectacled and Burmeister's porpoises shared a more recent common ancestor than with the vaquita that diverged from southern species during the Pliocene. The low genetic diversity observed in the vaquita carried signatures of a very low population size since the last 5,000 years. Cryptic lineages within Dall's, spectacled and Pacific harbor porpoises suggest a richer evolutionary history than previously suspected. These results provide a new perspective on the mechanisms driving diversification in porpoises and an evolutionary framework for their conservation.

An odyssey out of Africa: an integrative review of past and present invasions by the Mediterranean tree frog (Hyla meridionalis)

An illustration of the human footprint on biodiversity are the faunal movements that have accompanied commercial and cultural exchanges between civilizations throughout history. In this article, we provide an integrative review of biogeographical and archaeological knowledge to understand these processes for the Mediterranean tree frog (Hyla meridionalis), an African species that has reached several regions of south-western Europe and associated islands. By re-analysing molecular (multilocus mitochondrial and genome-wide markers) and bioclimatic (ecological niche modelling) data, we give a comprehensive picture of the genetic diversity, structure and environmental suitability of H. meridionalis across its natural and exotic ranges, which in turn offer specific clues to the putative routes of colonization and associated events. Long-term monitoring efforts suggest northwards shifts of local range margins, potentially due to global warming, and we further demonstrate that this species is the most frequent amphibian to travel via the food supply chain in Western Europe. High dispersal ability, ecological tolerance and proximity to human settlements have made H. meridionalis a recurrent witness to the complexity and diversity of the civilizations that ruled the Western Mediterranean.

The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Because it is indicative of reproductive isolation, the amount of genetic introgression across secondary contact zones is increasingly considered in species delimitation. However, patterns of admixture at range margins can be skewed by the regional dynamics of hybrid zones. In this context, we posit an important role for phylogeographic history: hybrid zones located within glacial refugia (putatively formed during the Late-Pleistocene) should be better defined than those located in post-glacial or introduced ranges (putatively formed during the Holocene and the Anthropocene). We test this hypothesis in a speciation continuum of tree frogs from the Western Palearctic (Hyla), featuring ten identified contacts between species spanning Plio-Pleistocene to Miocene divergences. We review the rich phylogeographic literature of this group and examine the overlooked transition between H. arborea and H. molleri in Western France using a multilocus dataset. Our comparative analysis supports a trend that contacts zones resulting from post-glacial expansions and human translocations feature more extensive introgression than those established within refugial areas. Integrating the biogeographic history of incipient species, i.e. their age since first contact together with their genetic divergence, thus appears timely to draw sound evolutionary and taxonomic inferences from patterns of introgression across hybrid zones.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

Phylogeography of Aegean green toads (Bufo viridis subgroup): continental hybrid swarm vs. insular diversification with discovery of a new island endemic

Background

Debated aspects in speciation research concern the amount of gene flow between incipient species under secondary contact and the modes by which post-zygotic isolation accumulates. Secondary contact zones of allopatric lineages, involving varying levels of divergence, provide natural settings for comparative studies, for which the Aegean (Eastern Mediterranean) geography offers unique scenarios. In Palearctic green toads (Bufo viridis subgroup or Bufotes), Plio-Pleistocene (~ 2.6 Mya) diverged species show a sharp transition without contemporary gene flow, while younger lineages, diverged in the Lower-Pleistocene (~ 1.9 Mya), admix over tens of kilometers. Here, we conducted a fine-scale multilocus phylogeographic analysis of continental and insular green toads from the Aegean, where a third pair of taxa, involving Mid-Pleistocene diverged (~ 1.5 Mya) mitochondrial lineages, earlier tentatively named viridis and variabilis, (co-)occurs.

Results

We discovered a new lineage, endemic to Naxos (Central Cyclades), while coastal islands and Crete feature weak genetic differentiation from the continent. In continental Greece, both lineages, viridis and variabilis, form a hybrid swarm, involving massive mitochondrial and nuclear admixture over hundreds of kilometers, without obvious selection against hybrids.

Conclusions

The genetic signatures of insular Aegean toads appear governed by bathymetry and Quaternary sea level changes, resulting in long-term isolation (Central Cyclades: Naxos) and recent land-bridges (coastal islands). Conversely, Crete has been isolated since the end of the Messinian salinity crisis (5.3 My) and Cretan populations thus likely result from human-mediated colonization, at least since Antiquity, from Peloponnese and Anatolia. Comparisons of green toad hybrid zones support the idea that post-zygotic hybrid incompatibilities accumulate gradually over the genome. In this radiation, only one million years of divergence separate a scenario of complete reproductive isolation, from a secondary contact resulting in near panmixia.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1179-0) contains supplementary material, which is available to authorized users.

Using conventional F-statistics to study unconventional sex-chromosome differentiation

Species with undifferentiated sex chromosomes emerge as key organisms to understand the astonishing diversity of sex-determination systems. Whereas new genomic methods are widening opportunities to study these systems, the difficulty to separately characterize their X and Y homologous chromosomes poses limitations. Here we demonstrate that two simple F-statistics calculated from sex-linked genotypes, namely the genetic distance (F_st) between sexes and the inbreeding coefficient (F_is) in the heterogametic sex, can be used as reliable proxies to compare sex-chromosome differentiation between populations. We correlated these metrics using published microsatellite data from two frog species (Hyla arboreaand Rana temporaria), and show that they intimately relate to the overall amount of X–Y differentiation in populations. However, the fits for individual loci appear highly variable, suggesting that a dense genetic coverage will be needed for inferring fine-scale patterns of differentiation along sex-chromosomes. The applications of these F-statistics, which implies little sampling requirement, significantly facilitate population analyses of sex-chromosomes.

Hidden diversity and phylogeographic history provide conservation insights for the edible seaweed Sargassum fusiforme in the Northwest Pacific

Understanding the evolutionary processes that have created diversity and the genetic potential of species to adapt to environmental change is an important premise for biodiversity conservation. Herein, we used mitochondrial trnW‐L and cox3 and plastid rbcL‐S data sets to analyze population genetic variation and phylogeographic history of the brown alga Sargassum fusiforme, whose natural resource has been largely exterminated in the Asia–Northwest Pacific in the past decades. Phylogenetic trees and network analysis consistently revealed three major haplotype groups (A, B, and C) in S. fusiforme, with A and B distributed in the Japan‐Pacific coast. Group C consisted of three subgroups (C1, C2, and C3) which were distributed in the Sea of Japan, the Yellow–Bohai Sea, and East China Sea, respectively. Isolation‐with‐migration (IMa) analysis revealed that the three groups diverged approximately during the mid‐Pleistocene (c. 756–1,224 ka). Extended Bayesian skyline plots (EBSP) showed that groups A and B underwent relatively long‐term stable population size despite a subsequent rapid demographic expansion, while subgroups C2 and C3 underwent a sudden expansion at c. 260 ka. F_ST and AMOVA detected low population‐level genetic variation and high degrees of divergence between groups. The cryptic diversity and phylogeographic patterns found in S. fusiforme not only are essential to understand how environmental shifts and evolutionary processes shaped diversity and distribution of coastal seaweeds but also provide additional insights for conserving and managing seaweed resources and facilitate predictions of their responses to future climate change and habitat loss.

Phylogeography and postglacial colonization of Central Europe by Anguis fragilis and Anguis colchica

Five slow-worm species are distributed in the Palearctic region. Two species, Anguis fragilis and A. colchica, are native to Central Europe, where only limited information on the phylogeography of the genus exists. Here, we examined the situation in Poland and surrounding regions, where a mitochondrial contact zone between the species is expected. We used new mitochondrial DNA sequences and available published data from Central Europe and the northern Balkans. Haplotypes of both species were recorded in the study area. Anguis fragilis is represented by a single haplogroup, while A. colchica by three haplogroups. This suggests four independent sources/refugia for postglacial colonization of northern Central Europe. The mitochondrial contact zone seems to mirror the borders between lowlands of the North European Plain and East European Plain, and the south-eastern Poland uplands, while the Vistula River does not constitute a barrier. The presence of both species, A. fragilis and A. colchica, in Poland should also be considered for protection by the Polish conservation legislation.

Phylogeography reveals an ancient cryptic radiation in East-Asian tree frogs (Hyla japonica group) and complex relationships between continental and island lineages

Background

In contrast to the Western Palearctic and Nearctic biogeographic regions, the phylogeography of Eastern-Palearctic terrestrial vertebrates has received relatively little attention. In East Asia, tectonic events, along with Pleistocene climatic conditions, likely affected species distribution and diversity, especially through their impact on sea levels and the consequent opening and closing of land-bridges between Eurasia and the Japanese Archipelago. To better understand these effects, we sequenced mitochondrial and nuclear markers to determine phylogeographic patterns in East-Asian tree frogs, with a particular focus on the widespread H. japonica.

Results

We document several cryptic lineages within the currently recognized H. japonica populations, including two main clades of Late Miocene divergence (~5 Mya). One occurs on the northeastern Japanese Archipelago (Honshu and Hokkaido) and the Russian Far-East islands (Kunashir and Sakhalin), and the second one inhabits the remaining range, comprising southwestern Japan, the Korean Peninsula, Transiberian China, Russia and Mongolia. Each clade further features strong allopatric Plio-Pleistocene subdivisions (~2-3 Mya), especially among continental and southwestern Japanese tree frog populations. Combined with paleo-climate-based distribution models, the molecular data allowed the identification of Pleistocene glacial refugia and continental routes of postglacial recolonization. Phylogenetic reconstructions further supported genetic homogeneity between the Korean H. suweonensis and Chinese H. immaculata, suggesting the former to be a relic population of the latter that arose when the Yellow Sea formed, at the end of the last glaciation.

Conclusions

Patterns of divergence and diversity were likely triggered by Miocene tectonic activities and Quaternary climatic fluctuations (including glaciations), causing the formation and disappearance of land-bridges between the Japanese islands and the continent. Overall, this resulted in a ring-like diversification of H. japonica around the Sea of Japan. Our findings urge for important taxonomic revisions in East-Asian tree frogs. First, they support the synonymy of H. suweonensis (Kuramoto, 1980) and H. immaculata (Boettger, 1888). Second, the nominal H. japonica (Günther, 1859) represents at least two species: an eastern (new taxon A) on the northern Japanese and Russian Far East islands, and a southwestern species (n. t. B) on southern Japanese islands and possibly also forming continental populations. Third, these continental tree frogs may also represent an additional entity, previously described as H. stepheni Boulenger, 1888 (senior synonym of H. ussuriensis Nikolskii, 1918). A complete revision of this group requires further taxonomic and nomenclatural analyses, especially since it remains unclear to which taxon the species-epitheton japonica corresponds to.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0814-x) contains supplementary material, which is available to authorized users.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Multilocus phylogeography of the common midwife toad, Alytes obstetricans (Anura, Alytidae): Contrasting patterns of lineage diversification and genetic structure in the Iberian refugium

Recent investigations on the evolutionary history of the common midwife toad (Alytes obstetricans) revealed high levels of geographically structured genetic diversity but also a situation where delineation of major historical lineages and resolution of their relationships are much more complex than previously thought. We studied sequence variation in one mitochondrial and four nuclear genes throughout the entire distribution range of all recognized A. obstetricans subspecies to infer the evolutionary processes that shaped current patterns of genetic diversity and population subdivision. We found six divergent, geographically structured mtDNA haplogroups diagnosing population lineages, and varying levels of admixture in nuclear markers. Given the timeframe inferred for the splits between major lineages, the climatic and environmental changes that occurred during the Pleistocene seem to have shaped the diversification history of A. obstetricans. Survival of populations in allopatric refugia through the Ice Ages supports the generality of the "refugia-within-refugia" scenario for the Iberian Peninsula. However, lineages corresponding to subspecies A. o. almogavarii, A. o. pertinax, A. o. obstetricans, and A. o. boscai responded differently to Pleistocene climatic oscillations after diverging from a common ancestor. Alytes o. obstetricans expanded northward from a northern Iberian refugium through the western Pyrenees, leaving a signal of contrasting patterns of genetic diversity, with a single mtDNA haplotype north of the Pyrenees from SW France to Germany. Both A. o. pertinax and A. o. boscai are widespread and genetically diverse in Iberia, the latter comprising two divergent lineages with a long independent history. Finally, A. o. almogavarii is mostly restricted to the north-eastern corner of Iberia north of the Ebro river, with additional populations in a small region in south-eastern France. This taxon exhibits unparalleled levels of genetic diversity and little haplotype sharing with other lineages, suggesting a process of incipient speciation.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

Background

Hybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.

Results

We found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.

Conclusions

Our findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0385-2) contains supplementary material, which is available to authorized users.

Effect of biogeographic history on population vulnerability in European amphibians

The genetic diversity of populations, which contributes greatly to their adaptive potential, is negatively affected by anthropogenic habitat fragmentation and destruction. However, continental-scale losses of genetic diversity also resulted from the population expansions that followed the end of the last glaciation, an element that is rarely considered in a conservation context. We addressed this issue in a meta-analysis in which we compared the spatial patterns of vulnerability of 18 widespread European amphibians in light of phylogeographic histories (glacial refugia and postglacial routes) and anthropogenic disturbances. Conservation statuses significantly worsened with distances from refugia, particularly in the context of industrial agriculture; human population density also had a negative effect. These findings suggest that features associated with the loss of genetic diversity in post-glacial amphibian populations (such as enhanced fixation load or depressed adaptive potential) may increase their susceptibility to current threats (e.g., habitat fragmentation and pesticide use). We propose that the phylogeographic status of populations (i.e., refugial vs. post-glacial) should be considered in conservation assessments for regional and national red lists.

Genetic consequences of population expansions and contractions in the common hippopotamus (Hippopotamus amphibius) since the Late Pleistocene

Over the past two decades, an increasing amount of phylogeographic work has substantially improved our understanding of African biogeography, in particular the role played by Pleistocene pluvial-drought cycles on terrestrial vertebrates. However, still little is known on the evolutionary history of semi-aquatic animals, which faced tremendous challenges imposed by unpredictable availability of water resources. In this study, we investigate the Late Pleistocene history of the common hippopotamus (Hippopotamus amphibius), using mitochondrial and nuclear DNA sequence variation and range-wide sampling. We documented a global demographic and spatial expansion approximately 0.1-0.3 Myr ago, most likely associated with an episode of massive drainage overflow. These events presumably enabled a historical continent-wide gene flow among hippopotamus populations, and hence, no clear continental-scale genetic structuring remains. Nevertheless, present-day hippopotamus populations are genetically disconnected, probably as a result of the mid-Holocene aridification and contemporary anthropogenic pressures. This unique pattern contrasts with the biogeographic paradigms established for savannah-adapted ungulate mammals and should be further investigated in other water-associated taxa. Our study has important consequences for the conservation of the hippo, an emblematic but threatened species that requires specific protection to curtail its long-term decline.

First-generation linkage map for the European tree frog (Hyla arborea) with utility in congeneric species

Background

Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications.

Findings

We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement.

Conclusions

The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-850) contains supplementary material, which is available to authorized users.

Sex-chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea)

Occasional XY recombination is a proposed explanation for the sex-chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW-European Hyla arborea populations identified male-specific alleles at sex-linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex-linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex-chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY-recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).

Range-wide sex-chromosome sequence similarity supports occasional XY recombination in European tree frogs (Hyla arborea)

In contrast with mammals and birds, most poikilothermic vertebrates feature structurally undifferentiated sex chromosomes, which may result either from frequent turnovers, or from occasional events of XY recombination. The latter mechanism was recently suggested to be responsible for sex-chromosome homomorphy in European tree frogs (Hyla arborea). However, no single case of male recombination has been identified in large-scale laboratory crosses, and populations from NW Europe consistently display sex-specific allelic frequencies with male-diagnostic alleles, suggesting the absence of recombination in their recent history. To address this apparent paradox, we extended the phylogeographic scope of investigations, by analyzing the sequences of three sex-linked markers throughout the whole species distribution. Refugial populations (southern Balkans and Adriatic coast) show a mix of X and Y alleles in haplotypic networks, and no more within-individual pairwise nucleotide differences in males than in females, testifying to recurrent XY recombination. In contrast, populations of NW Europe, which originated from a recent postglacial expansion, show a clear pattern of XY differentiation; the X and Y gametologs of the sex-linked gene Med15 present different alleles, likely fixed by drift on the front wave of expansions, and kept differentiated since. Our results support the view that sex-chromosome homomorphy in H. arborea is maintained by occasional or historical events of recombination; whether the frequency of these events indeed differs between populations remains to be clarified.

Stronger transferability but lower variability in transcriptomic- than in anonymous microsatellites: evidence from Hylid frogs

A simple way to quickly optimize microsatellites in nonmodel organisms is to reuse loci available in closely related taxa; however, this approach can be limited by the stochastic and low cross-amplification success experienced in some groups (e.g. amphibians). An efficient alternative is to develop loci from transcriptome sequences. Transcriptomic microsatellites have been found to vary in their levels of cross-species amplification and variability, but this has to date never been tested in amphibians. Here, we compare the patterns of cross-amplification and levels of polymorphism of 18 published anonymous microsatellites isolated from genomic DNA vs. 17 loci derived from a transcriptome, across nine species of tree frogs (Hyla arborea and Hyla cinerea group). We established a clear negative relationship between divergence time and amplification success, which was much steeper for anonymous than transcriptomic markers, with half-lives (time at which 50% of the markers still amplify) of 1.1 and 37 My, respectively. Transcriptomic markers are significantly less polymorphic than anonymous loci, but remain variable across diverged taxa. We conclude that the exploitation of amphibian transcriptomes for developing microsatellites seems an optimal approach for multispecies surveys (e.g. analyses of hybrid zones, comparative linkage mapping), whereas anonymous microsatellites may be more informative for fine-scale analyses of intraspecific variation. Moreover, our results confirm the pattern that microsatellite cross-amplification is greatly variable among amphibians and should be assessed independently within target lineages. Finally, we provide a bank of microsatellites for Palaearctic tree frogs (so far only available for H. arborea), which will be useful for conservation and evolutionary studies in this radiation.