Acknowledging more biodiversity without more species.

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

Phylogenetic insights into Central European Chorthippus and Pseudochorthippus (Orthoptera: Acrididae) species using ddRADseq data.

The evolution of several orthopteran groups, especially within the grasshopper family Acrididae, remains poorly understood. This is particularly true for the subfamily Gomphocerinae, which comprises cryptic sympatric and syntopic species. Previous mitochondrial studies have highlighted major discrepancies between taxonomic and phylogenetic hypotheses, thereby emphasizing the necessity of genome-wide approaches. In this study, we employ double-digest restriction site-associated DNA sequencing (ddRADseq) to reconstruct the evolution of Central European Chorthippus and Pseudochorthippus species, especially C.smardai, P.tatrae and the C.biguttulus group. Our phylogenomic analyses recovered deep discordance with mitochondrial DNA barcoding, emphasizing its unreliability in Gomphocerinae grasshoppers. Specifically, our data robustly distinguished the C.biguttulus group and confirmed the distinctiveness of C.eisentrauti, also shedding light on its presence in the Berchtesgaden Alps. Moreover, our results support the reclassification of C.smardai to the genus Pseudochorthippus and of P.tatrae to the genus Chorthippus. Our study demonstrates the efficiency of high-throughput genomic methods such as RADseq without prior optimization to elucidate the complex evolution of grasshopper radiations with direct taxonomic implications. While RADseq has predominantly been utilized for population genomics and within-genus phylogenomics, its application extends to resolve relationships between deeply-diverged clades representative of distinct genera.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture.

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Phylogenomic insights into the diversity and evolution of Palearctic vipers.

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of âˆ¼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions.

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Countrywide screening supports model-based predictions of the distribution of Batrachochytrium dendrobatidis in Ukraine.

The chytrid Batrachochytrium dendrobatidis (Bd) is a widespread fungus causing amphibian declines across the globe. Although data on Bd occurrence in Eastern Europe are scarce, a recent species distribution model (SDM) for Bd reported that western and north-western parts of Ukraine are highly suitable to the pathogen. We verified the SDM-predicted range of Bd in Ukraine by sampling amphibians across the country and screening for Bd using qPCR. A total of 446 amphibian samples (tissue and skin swabs) from 11 species were collected from 36 localities. We obtained qPCR-positive results for 33 samples including waterfrogs (Pelophylax esculentus complex) and fire- and yellow-bellied toads (Bombina spp.) from 8 localities. We found that Bd-positive localities had significantly higher predicted Bd habitat suitability than sites that were pathogen-free. Amplification and sequencing of the internal transcribed spacer (ITS) region of samples with the highest Bd load revealed matches with ITS haplotypes of the globally distributed BdGPL strain, and a single case of the BdASIA-2/BdBRAZIL haplotype. We found that Bd was non-randomly distributed across Ukraine, with infections present in the western and north-central forested peripheries of the country with a relatively cool, moist climate. On the other hand, our results suggest that Bd is absent or present in low abundance in the more continental central, southern and eastern regions of Ukraine, corroborating the model-predicted distribution of chytrid fungus. These areas could potentially serve as climatic refugia for Bd-susceptible amphibian hosts.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation.

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes).

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide a comprehensive phylogeographic framework for the A. obstetricans complex, as well as a fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

When Sex Chromosomes Recombine Only in the Heterogametic Sex: Heterochiasmy and Heterogamety in Hyla Tree Frogs.

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs.

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean.

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.

Population genomics of an exceptional hybridogenetic system of Pelophylax water frogs.

BACKGROUND: Hybridogenesis can represent the first stage towards hybrid speciation where the hybrid taxon eventually weans off its parental species. In hybridogenetic water frogs, the hybrid Pelophylax kl. esculentus (genomes RL) usually eliminates one genome from its germline and relies on its parental species P. lessonae (genomes LL) or P. ridibundus (genomes RR) to perpetuate in so-called L-E and R-E systems. But not exclusively: some all-hybrid populations (E-E system) bypass the need for their parental species and fulfill their sexual cycle via triploid hybrid frogs. Genetic surveys are essential to understand the great diversity of these hybridogenetic dynamics and their evolution. Here we conducted such study using RAD-sequencing on Pelophylax from southern Switzerland (Ticino), a geographically-isolated region featuring different assemblages of parental P. lessonae and hybrid P. kl. esculentus. RESULTS: We found two types of hybridogenetic systems in Ticino: an L-E system in northern populations and a presumably all-hybrid E-E system in the closely-related southern populations, where P. lessonae was not detected. In the latter, we did not find evidence for triploid individuals from the population genomic data, but identified a few P. ridibundus (RR) as offspring from interhybrid crosses (LR × LR). CONCLUSIONS: Assuming P. lessonae is truly absent from southern Ticino, the putative maintenance of all-hybrid populations without triploid individuals would require an unusual lability of genome elimination, namely that P. kl. esculentus from both sexes are capable of producing gametes with either L or R genomes. This could be achieved by the co-existence of L- and R- eliminating lineages or by "hybrid amphigamy", i. e. males and females producing sperm and eggs among which both genomes are represented. These hypotheses imply that polyploidy is not the exclusive evolutionary pathway for hybrids to become reproductively independent, and challenge the classical view that hybridogenetic taxa are necessarily sexual parasites.

Sex-Chromosome Recombination in Common Frogs Brings Water to the Fountain-of-Youth.

According to the canonical model of sex-chromosome evolution, the degeneration of Y or W chromosomes (as observed in mammals and birds, respectively) results from an arrest of recombination in the heterogametic sex, driven by the fixation of sexually antagonistic mutations. However, sex chromosomes have remained homomorphic in many lineages of fishes, amphibians, and nonavian reptiles. According to the "fountain-of-youth" model, this homomorphy results from occasional events of sex reversal. If recombination arrest in males is controlled by maleness per se (and not by genotype), then Y chromosomes are expected to recombine in XY females, preventing their long-term degeneration. Here, we provide field support for the fountain-of-youth, by showing that sex-chromosome recombination in Rana temporaria only depends on phenotypic sex: naturally occurring XX males show the same restriction of recombination as XY males (average map length ∼2 cM), while XY females recombine as much as XX females (average map length ∼150 cM). Our results challenge several common assumptions regarding the evolution of sex chromosomes, including the role of sexually antagonistic genes as drivers of recombination arrest, and that of chromosomal inversions as underlying mechanisms.

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.

Diversification and speciation in tree frogs from the Maghreb (Hyla meridionalis sensu lato), with description of a new African endemic.

Comparative molecular studies emphasized a new biogeographic paradigm for the terrestrial fauna of North Africa, one of the last uncharted ecoregions of the Western Palearctic: two independent east-west divisions across the Maghreb. Through a comprehensive phylogeography, we assessed how this model suits the genetic diversification documented for the tree frog Hyla meridionalis sensu lato. Analyses of mtDNA variation and thousands of nuclear loci confirmed the old split (low-Pliocene) between Tunisian and Algerian populations. These lineages meet but barely admix in the eastern Maghreb (Algerian-Tunisian border), a sign of putatively advanced reproductive isolation. In the western Maghreb, we report a Pleistocene divergence between Moroccan and Algerian populations. Tree frogs thus follow both predictions: a double east-west break that gave rise to two suture zones characteristic of North-African phylogeography. Moreover, some intraspecific mtDNA variation is not mirrored by the nuclear data, emphasizing that evolutionary units should always be designated by multilocus approaches. Last but not least, we describe the Tunisian lineage as a new species endemic to Africa.

Fifteen shades of green: The evolution of Bufotes toads revisited.

The radiation of Palearctic green toads (Bufotes) holds great potential to evaluate the role of hybridization in phylogeography at multiple stages along the speciation continuum. With fifteen species representing three ploidy levels, this model system is particularly attractive to examine the causes and consequences of allopolyploidization, a prevalent yet enigmatic pathway towards hybrid speciation. Despite substantial efforts, the evolutionary history of this species complex remains largely blurred by the lack of consistency among the corresponding literature. To get a fresh, comprehensive view on Bufotes phylogeography, here we combined genome-wide multilocus analyses (RAD-seq) with an extensive compilation of mitochondrial, genome size, niche modelling, distribution and phenotypic (bioacoustics, morphometrics, toxin composition) datasets, representing hundreds of populations throughout Eurasia. We provide a fully resolved nuclear phylogeny for Bufotes and highlight exceptional cyto-nuclear discordances characteristic of complete mtDNA replacement (in 20% of species), mitochondrial surfing during post-glacial expansions, and the formation of homoploid hybrid populations. Moreover, we traced the origin of several allopolyploids down to species level, showing that all were exclusively fathered by the West Himalayan B. latastii but mothered by several diploid forms inhabiting Central Asian lowlands, an asymmetry consistent with hypotheses on mate choice and Dobzhansky-Muller incompatibilities. Their intermediate call phenotypes potentially allowed for rapid reproductive isolation, while toxin compositions converged towards the ecologically-closest parent. Across the radiation, we pinpoint a stepwise progression of reproductive isolation through time, with a threshold below which hybridizability is irrespective of divergence (<6My), above which species barely admix and eventually evolve different mating calls (6-10My), or can successfully cross-breed through allopolyploidization (>15My). Finally, we clarified the taxonomy of Bufotes (including genetic analyses of type series) and formally described two new species, B. cypriensis sp. nov. (endemic to Cyprus) and B. perrini sp. nov. (endemic to Central Asia). Embracing the genomic age, our framework marks the advent of a new exciting era for evolutionary research in these iconic amphibians.

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.

Introgression across Hybrid Zones Is Not Mediated by Large X-Effects in Green Toads with Undifferentiated Sex Chromosomes.

Divergence between incipient species remains an incompletely understood process. Hybrid zones provide great research potential, reflecting natural organismal genomic interactions and gene evolution in a variety of recombinants over generations. While sex chromosomes are known evolutionary drivers of reproductive isolation, empirical population genetics has mostly examined species with heteromorphic sex chromosomes. We recently reported restricted introgression at sex-linked markers in an amphibian system with homomorphic sex chromosomes (Hyla), consistent with a large X-effect, designating a greater role of sex chromosomes in driving hybrid incompatibilities. Here, using a similar approach, we examined two hybrid zones of Palearctic green toads (Bufo viridis subgroup), involving several lineages that arose at different times and form secondary contacts. We find no evidence for differential introgression of sex-linked versus autosomal markers across both zones. This absence of large X-effects in Bufo indicates that, unlike in Hyla, hybrid incompatibilities may not result from the faster-heterogametic sex and faster-male aspects of Haldane's rule. The recent suppression of XY recombination in Hyla but not in Bufo may have driven greater divergence between Hyla sex chromosomes, causing stronger reproductive isolation. Alternatively, stronger linkage among Hyla's sex-linked markers could restrict introgression. We hypothesize that the degree of sex-specific recombination may condition the importance of homomorphic sex chromosomes in speciation.

Howling from the past: historical phylogeography and diversity losses in European grey wolves.

Genetic bottlenecks resulting from human-induced population declines make alarming symbols for the irreversible loss of our natural legacy worldwide. The grey wolf (Canis lupus) is an iconic example of extreme declines driven by anthropogenic factors. Here, we assessed the genetic signatures of 150 years of wolf persecution throughout the Western Palaearctic by high-throughput mitochondrial DNA sequencing of historical specimens in an unprecedented spatio-temporal framework. Despite Late Pleistocene bottlenecks, we show that historical genetic variation had remained high throughout Europe until the last several hundred years. In Western Europe, where wolves nearly got fully exterminated, diversity dramatically collapsed at the turn of the twentieth century and recolonization from few homogeneous relict populations induced drastic shifts of genetic composition. By contrast, little genetic displacement and steady levels of diversity were maintained in Eastern European regions, where human persecution had lesser effects on wolf demography. By comparing prehistoric, historic and modern patterns of genetic diversity, our study hence traces the timeframe and the active human role in the decline of the grey wolf, an emblematic yet controversial animal which symbolizes the complex relationship between human societies and nature conservation.