Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae).

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Diversity of trematodes from the amphibian anomaly P hotspot: Role of planorbid snails.

Trematode infection of the second intermediate hosts can lead to changes in their fitness and, as a result, a change in the invasion rate of animal communities. It is especially pronounced during the invasion of parasite species that reduce activity due to the manipulation of hosts through the changes of their morphology and physiology. One of these cases is an anomaly P syndrome hotspot found in some populations of water frogs and toads in Europe caused by the trematode Strigea robusta metacercariae. The occurrence of pathogen and their participation in ecosystems are intrigues questions in the anomaly P phenomenon, as well as the role of planorbid snails that serve as the first intermediate hosts for many trematode species. Herein, we focused on trematodes spectra from planorbid snails and amphibians from the anomaly P hosts with the aim to undetected interactions between the pathways of parasites. Emerging cercariae of 6802 planorbid snails of dominant species (Planorbarius corneus, Planorbis planorbis, and Anisus spp.) were detected by both morphological and molecular methods in seven waterbodies in Privolzhskaya Lesostep Nature Reserve (Russia). A total of 95 sequences of 18 species were received, and 48 sequences were unique and did not present in any genetic databases. The 18 species of trematodes from snails and 14 species of trematodes from amphibian hosts (Pelophylax ridibundus; Ranidae; Anura) were detected. Three species (Echinostoma nasincovae, Tylodelphys circibuteonis and Australapatemon burti) was new for the trematode fauna of the Middle Volga River region and Russia as a whole. Eleven species of parasitic flatworms have amphibians in their life cycles and nine species used amphibians as metacercariae hosts: Echinostoma nasincovae, E. miyagawai, Echinoparyphium recurvatum, Tylodelphys circibuteonis, Neodiplostomum spathula, Paralepoderma cloacicola, Macrodera longicollis, Strigea robusta, and Strigea strigis. The occurrence of trematode species from planorbid mollusks and frogs were compared.

Effects of minimum Strigea robusta (Digenea: Strigeidae) cercariae doses and localization of cysts on the anomaly P manifestation in Pelophylax lessonae (Anura: Ranidae) tadpoles.

The anomaly P is a mass morphological anomaly found in some populations of anuran amphibians (water frogs of the genus Pelophylax and toads of the genera Bufo and Bufotes) caused by the parasitic flatworm Strigea robusta. Minimum dose of cercariae for the appearance of the anomaly P remains unknown. However, it is important information for understanding of host population dynamics after invasion and the effects of the parasite on the second intermediate hosts. Herein, the invasion properties of S. robusta in Pelophylax lessonae tadpoles (Anura: Ranidae) and minimum dose for appearance of mild and severe forms of the anomaly P syndrome were described after direct experiments with certain numbers of cercariae exposure. Experimental groups of tadpoles have been exposed to eight doses of cercariae (2, 4, 6, 8, 10, 12, 14 and 16). A total of 63.8% tadpoles survived to the end of this experiment. It was revealed that a mild form of the anomaly P (polydactyly) can appear after infection by two cercariae, while the severe form traits appear after infection by four cercariae. The mean number of detected encysted metacercariae was reached to 53.5%. Differences in infection rates can be explained by the presence of an individual immune response in tadpoles or by the presence of different genetic lineages of the parasite infecting the same snail, which have different infectious potential. Low doses of infection leading to the induction of anomalies characterize S. robusta as a highly pathogenic species for amphibian species that are susceptible to infection and show an abnormal phenotype.

Strigea robusta (Digenea: Strigeidae) infection effects on the gonadal structure and limb malformation in toad early development.

The anomaly P is a mass morphological anomaly reported in some water frog populations across Europe. It was found that polydactyly is only a mild attenuated form of heavy cases of the anomaly P syndrome, which have strong deformations of the hindlimbs and, partly, forelimbs. It was shown that the anomaly P is caused by the trematode Strigea robusta and this syndrome can be considered as a special case of strigeosis in amphibians. The anomaly P for a long time considered to be specific for water frogs of the genus Pelophylax. Herein, we describe polydactyly and heavy forms of the anomaly P syndrome in toads of the genera Bufo and Bufotes, as a result of exposure to S. robusta cercariae. A total of 150 tadpoles of Bufo bufo, 60 tadpoles of Bufotes viridis, and 60 tadpoles of Bufotes baturae were divided into five experimental and four control groups (30 tadpoles in each group). All anomalies in the toads were similar to those observed in water frogs. The survival of tadpoles in the experimental groups was 76%. The anomaly P was observed in 57.9% of toad tadpoles (51.8% of mild forms and 6.1% of heavy forms). The occurrence of the anomaly P varied among groups from 19% to 78%. Heavy forms of the anomaly P were found in all experimental groups. We described rare asymmetrical cases of the anomaly P. According to severe modification of limb morphology, we supposed changes of gonadal morphology (any modifications of the germ and somatic cells). The gonadal development of infected tadpoles was however the same as in uninfected toad tadpoles, and heterochromatin distribution within gonocytes had no differences as well. It seems like the parasite doesn't have any effect on the gonadal development of the toads. The lack of heavy forms in natural populations of toads, as well as a development of gonads were discussed.

From Gondwana to the Yellow Sea, evolutionary diversifications of true toads Bufo sp. in the Eastern Palearctic and a revisit of species boundaries for Asian lineages.

Taxa with vast distribution ranges often display unresolved phylogeographic structures and unclear taxonomic boundaries resulting in hidden diversity. This hypothesis-driven study reveals the evolutionary history of Bufonidae, covering the phylogeographic patterns found in Holarctic bufonids from the West Gondwana to the phylogenetic taxonomy of Asiatic true toads in the Eastern Palearctic. We used an integrative approach relying on fossilized birth-death calibrations, population dynamics, gene-flow, species distribution, and species delimitation modeling to resolve the biogeography of the clade and highlight cryptic lineages. We verified the near-simultaneous Miocene radiations within Western and Eastern Palearctic Bufo, c. 14.49-10.00 Mya, temporally matching with the maximum dust outflows in Central Asian deserts. Contrary to earlier studies, we demonstrated that the combined impacts of long dispersal and ice-age refugia equally contributed to the current genetic structure of Bufo in East Asia. Our findings reveal a climate-driven adaptation in septentrional Eastern Asian Bufo, explaining its range shifts toward northern latitudes. We resolve species boundaries within the Eastern Palearctic Bufo, and redefine the taxonomic and conservation units of the northeastern species: B. sachalinensis and its subspecies.

The east Asian Asiatic toad (also known by the latin name Bufo gargarizans) lives in a wide range of habitats across East Asia including forests, meadows and cultivated land. However, it remains unclear how these toads evolved and became so widespread ­ partly because it has proved difficult for researchers to clearly define the species and what distinguishes it from other closely-related species of toads (collectively known as Bufo toads). Othman et al. combined several bioinformatics techniques to study Asiatic toads and 38 other species of bufonid toads from across the globe. This approach found that Bufo toads first emerged in eastern Asia between 14 and 10 million years ago. This coincides with a point in time when large swathes of land in central Asia turned from adequate to sustain toad populations into desert, suggesting the change in climate prompted the toads to migrate eastwards from central Asia. The Bufo toads then divided into two groups of species: one in mainland East Asia and the other in Japan. Furthermore, the study revealed there is more genetic diversity ­ that is, more variety in the DNA of individuals ­ in Asiatic toads than previously thought. The findings also help to define several other species of Bufo toads more clearly and describe a new toad species restricted to the Korean Peninsula, northeastern China and the Russian Primorye region: the Sakhalin toad (Bufo sachalinensis). This work demonstrates that a large-scale study of many species across the globe can be used to understand how the species evolved and more clearly distinguish one species from another. The findings of Othman et al. will be of interest to both professional and citizen scientists interested in the natural history of Asia. Furthermore, as several species of Bufo toads are in decline in the wild, they provide evidence that may aid future efforts to conserve them.

Next-generation phylogeography of the banded newts (Ommatotriton): A phylogenetic hypothesis for three ancient species with geographically restricted interspecific gene flow and deep intraspecific genetic structure.

Technological developments now make it possible to employ many markers for many individuals in a phylogeographic setting, even for taxa with large and complex genomes such as salamanders. The banded newt (genus Ommatotriton) from the Near East has been proposed to contain three species (O. nesterovi, O. ophryticus and O. vittatus) with unclear phylogenetic relationships, apparently limited interspecific gene flow and deep intraspecific geographic mtDNA structure. We use parallel tagged amplicon sequencing to obtain 177 nuclear DNA markers for 35 banded newts sampled throughout the range. We determine population structure (with Bayesian clustering and principal component analysis), interspecific gene flow (by determining the distribution of species-diagnostic alleles) and phylogenetic relationships (by maximum likelihood inference of concatenated sequence data and based on a summary-coalescent approach). We confirm that the three proposed species are genetically distinct. A sister relationship between O. nesterovi and O. ophryticus is suggested. We find evidence for introgression between O. nesterovi and O. ophryticus, but this is geographically limited. Intraspecific structuring is extensive, with the only recognized banded newt subspecies, O. vittatus cilicensis, representing the most distinct lineage below the species level. While mtDNA mostly mirrors the pattern observed in nuclear DNA, all banded newt species show mito-nuclear discordance as well.

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.

Update on Distribution and Conservation Status of Amphibians in the Democratic People's Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties.

Determining the range, status, ecology and behaviour of species from areas where surveys and samplings are uncommon or difficult to conduct is a challenge, such as in the Democratic People's Republic of Korea (DPR Korea). Here, we used genetic samples, field surveys, call recordings, photographic identification and a literature review to estimate the presence, range and status of amphibians in the DPR Korea. From our combined results and based on the IUCN Red List categories and criteria, we were able to estimate the national threat levels for most species. Our results demonstrated the presence of 18 native species and the suspected presence of Karsenia koreana and two Onychodactylus species. We reported the first record for Rana uenoi in the vicinity of Pyongyang using molecular tools and similarly confirmed the presence of Dryophytes japonicus at the same location. Based on distribution and modelling, we can expect the contact zone between species within the Rana and Onychodactylus genera to be located along the Changbai Massif, a mountain range that marks a shift in ecoregions and acts as a barrier to dispersion. The species richness was higher in the lowlands and at lower latitudes, with such areas populated by up to 11 species, while more northern regions were characterised by species richness of about half of that value. The combination of ecological models and known threats resulted in the recommendation of ten species as threatened at the national level following the IUCN Red List categories and criteria. This high number of threatened species was anticipated based on the high threat level to amphibians in bordering nations and globally. While the ecology of species in the DPR Korea is still understudied, we argue that species relying on agricultural wetlands such as rice paddies are not under imminent threat due to the enduring presence of extensive agricultural landscapes with low rates of chemical use and mechanisation. The maintenance of such landscapes is a clear benefit to amphibian species, in contrast to more industrialised agricultural landscapes in neighbouring nations. In comparison, the status of species dependent on forested habitats is unclear and threat levels are likely to be higher because of deforestation, as in neighbouring nations.

Relict groups of spiny frogs indicate Late Paleogene-Early Neogene trans-Tibet dispersal of thermophile faunal elements.

BACKGROUND: The Himalaya-Tibet orogen (HTO) presents an outstanding geologically active formation that contributed to, and fostered, modern Asian biodiversity. However, our concepts of the historical biogeography of its biota are far from conclusive, as are uplift scenarios for the different parts of the HTO. Here, we revisited our previously published data set of the tribe Paini extending it with sequence data from the most western Himalayan spiny frogs Allopaa and Chrysopaa and using them as an indirect indicator for the potential paleoecological development of Tibet. METHODS: We obtained sequence data of two mitochondrial loci (16S rRNA, COI) and one nuclear marker (Rag1) from Allopaa samples from Kashmir Himalaya as well as Chrysopaa sequence data from the Hindu Kush available from GenBank to complement our previous data set. A Maximum likelihood and dated Bayesian gene tree were generated based on the concatenated data set. To resolve the inconsistent placement of Allopaa, we performed different topology tests. RESULTS: Consistent with previous results, the Southeast Asian genus Quasipaa is sister to all other spiny frogs. The results further reveal a basal placement of Chrysopaa relative to Allopaa and Nanorana with an estimated age of ca. 26 Mya. Based on the topology tests, the phylogenetic position of Allopaa as a sister clade to Chaparana seems to be most likely, resulting in a paraphyletic genus Nanorana and a separation from the latter clade around 20 Mya, although a basal position of Allopaa to the genus Nanorana cannot be entirely excluded. Both, the placements of Chrysopaa and Allopaa support the presence of basal Paini lineages in the far northwestern part of the HTO, which is diametrically opposite end of the HTO with respect to the ancestral area of spiny frogs in Southeast Asia. These striking distributional patterns can be most parsimoniously explained by trans-Tibet dispersal during the late Oligocene (subtropical Chrysopaa) respectively early Miocene (warm temperate Allopaa). Within spiny frogs, only members of the monophyletic Nanorana+Paa clade are adapted to the colder temperate climates, indicating that high-altitude environments did not dominate in the HTO before ca. 15 Mya. Our results are consistent with fossil records suggesting that large parts of Tibet were characterized by subtropical to warm temperate climates at least until the early Miocene. They contradict prevalent geological models of a highly uplifted late Paleogene proto-Plateau.

Hybridization and introgression between toads with different sex chromosome systems.

The growing interest in the lability of sex determination in non-model vertebrates such as amphibians and fishes has revealed high rates of sex chromosome turnovers among closely related species of the same clade. Can such lineages hybridize and admix with different sex-determining systems, or could the changes have precipitated their speciation? We addressed these questions in incipient species of toads (Bufonidae), where we identified a heterogametic transition and characterized their hybrid zone with genome-wide markers (RADseq). Adult and sibship data confirmed that the common toad B. bufo is female heterogametic (ZW), while its sister species the spined toad B. spinosus is male heterogametic (XY). Analysis of a fine scale transect across their parapatric ranges in southeastern France unveiled a narrow tension zone (∼10 km), with asymmetric mitochondrial and nuclear admixture over hundreds of kilometers southward and northward, respectively. The geographic extent of introgression is consistent with an expansion of B. spinosus across B. bufo's former ranges in Mediterranean France, as also suggested by species distribution models. However, widespread cyto-nuclear discordances (B. spinosus backrosses carrying B. bufo mtDNA) run against predictions from the dominance effects of Haldane's rule, perhaps because Y and W heterogametologs are not degenerated. Common and spined toads can thus successfully cross-breed despite fundamental differences in their sex determination mechanisms, but remain partially separated by reproductive barriers. Whether and how the interactions of their XY and ZW genes contribute to these barriers shall provide novel insights on the debated role of labile sex chromosomes in speciation.

Strigea robusta causes polydactyly and severe forms of Rostand's anomaly P in water frogs.

BACKGROUND: Cases of polydactyly in natural populations of amphibians have attracted great interest from biologists. At the end of the 1940s, the French biologist Jean Rostand discovered a polymorphic syndrome in some water frog (Anura: Pelophylax) populations that included polydactyly and some severe morphological anomalies (he called it 'anomaly P'). The cause of this anomaly remains unknown for 70 years. In a previous study, we obtained anomaly P in the laboratory in tadpoles of water frogs that developed together with molluscs Planorbarius corneus (Mollusca: Gastropoda) collected in the field. We thus proposed the 'trematode hypothesis', according to which the infectious agent responsible for anomaly P is a trematode species. METHODS: Metacercariae from tadpoles with anomaly P were identified using ITS2 gene sequencing as Strigea robusta (Trematoda: Strigeidae). To verify teratogenic features of the species, cercariae of S. robusta were tested for the possibility to cause anomalies. Identification of cercariae species was made using morphological and molecular methods (sequencing of ITS2 and 28S rRNA). The tadpoles were exposed to parasites at four doses of cercariae (control, low, medium and high) and divided into two groups: "early" (at 25-27 Gosner stages) and "late" (at 29-34 Gosner stages) exposure. RESULTS: A total of 58 (72.5%) tadpoles survived until metamorphosis under the dose-dependent experiment with the trematode S. robusta. Differences in the survival rates were observed between the exposed and unexposed tadpoles both in the group of "early" tadpoles and "late" tadpoles. The exposure of tadpoles to the cercariae of S. robusta induced anomaly P in 82% of surviving tadpoles. The severe forms developed only in "early" stages under all doses of cercariae exposure. Polydactyly predominantly developed in the "late" stages; under a light exposure dose, polydactyly also developed in "early" tadpoles. Laboratory-hatched tadpoles reared together with infected snails had different rates of survival and complexity of deformations associated with the period of coexistence. CONCLUSIONS: The experiments with direct cercariae exposure provide compelling evidence that S. robusta leads to anomaly P in tadpoles of water frogs. The manifestation of anomaly P turned out to be dependent on the stage of development, cercariae dose, and the location of the cysts.

Reconstruction of past distribution for the Mongolian toad, Strauchbufo raddei (Anura: Bufonidae) using environmental modeling.

The use of ecological models enables determining the current distribution of species, but also their past distribution when matching climatic conditions are available. In specific cases, they can also be used to determine the likelihood of fossils to belong to the same species-under the hypothesis that all individuals of a species have the same ecological requirements. Here, using environmental modeling, we reconstructed the distribution of the Mongolian toad, Strauchbufo raddei, since the Last Glacial Maximum and thus covering the time period between the Late Pleistocene and the Holocene. We found the range of the species to have shifted over time, with the LGM population clustered around the current southern range of the species, before expanding east and north during the Pleistocene, and reaching the current range since the mid-Holocene. Finally, we determined that the ecological conditions during the life-time of the mid-Pleistocene fossils attributed to the species in Europe were too different from the one of the extant species or fossils occurring at the same period in Asia to belong to the same species.

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.

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.

Variation in hybridogenetic hybrid emergence between populations of water frogs from the Pelophylax esculentus complex.

Many closely related species are capable of mating to produce hybrid offspring, which are usually sterile. Nevertheless, altering the gametogenesis of hybrid offspring can rescue hybrids from sterility by enabling asexual reproduction. Hybridogenesis is one of the most complicated asexual reproductive modes, and it includes drastic genome reorganization only in the germline; this is achieved through elimination of one parental genome and duplication of the remaining one to restore diploid chromosomal set and overcome blocks in meiotic progression. We investigated a model of hybridogenesis, namely, water frogs from the Pelophylax esculentus complex, for the emergence of asexual reproduction. Further, we assessed the impact of its asexual reproduction on the maintenance of interspecies hybrids from two populations on the western edge of the P. esculentus range, in which hybrids coexist with either both parental species or with only one parental species. After analysing tadpole karyotypes, we conclude that in both studied populations, the majority of diploid hybrid males produced haploid gametes with the P. ridibundus genome after elimination of the P. lessonae genome. Hybrid females exhibited problems with genome elimination and duplication; they usually produced oocytes with univalents, but there were observations of individual oocytes with 13 bivalents and even 26 bivalents. In some hybrid tadpoles, especially F1 crosses, we observed failed germ cell development, while in tadpoles from backcrosses, germ cells were normally distributed and contained micronuclei. By identifying chromosomes present in micronuclei, we estimated that the majority of tadpoles from all crosses were able to selectively eliminate the P. lessonae chromosomes. According to our results, hybridogenesis in hybrids can appear both from crosses of parental species and crosses between sexual species with hybrid individuals. The ability to eliminate a genome and perform endoreplication to ensure gamete formation differed between male and female hybrids from the studied populations. Some diploid hybrid females can rarely produce not only haploid gametes but also diploid gametes, which is a crucial step in the formation of triploid hybrids.

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.

Call a spade a spade: taxonomy and distribution of Pelobates, with description of a new Balkan endemic.

The genomic era contributes to update the taxonomy of many debated terrestrial vertebrates. In an accompanying work, we provided a comprehensive molecular assessment of spadefoot toads (Pelobates) using genomic data. Our results call for taxonomic updates in this group. First, nuclear phylogenomics confirmed the species-level divergence between the Iberian P.cultripes and its Moroccan relative P.varaldii. Second, we inferred that P.fuscus and P.vespertinus, considered subspecies until recently, feature partial reproductive isolation and thus deserve a specific level. Third, we evidenced cryptic speciation and diversification among deeply diverged lineages collectively known as Pelobatessyriacus. Populations from the Near East correspond to the Eastern spadefoot toad P.syriacus sensu stricto, which is represented by two subspecies, one in the Levant (P.s.syriacus) and the other in the rest of the range (P.s.boettgeri). Populations from southeastern Europe correspond to the Balkan spadefoot toad, P.balcanicus. Based on genetic evidence, this species is also polytypic: the nominal P.b.balcanicus inhabits the Balkan Peninsula; a new subspecies P.b.chloeae ssp. nov. appears endemic to the Peloponnese. In this paper, we provide an updated overview of the taxonomy and distribution of all extant Pelobates taxa and describe P.b.chloeae ssp. nov.

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.

Phylogeny of spiny frogs Nanorana (Anura: Dicroglossidae) supports a Tibetan origin of a Himalayan species group.

Recent advances in the understanding of the evolution of the Asian continent challenge the long-held belief of a faunal immigration into the Himalaya. Spiny frogs of the genus Nanorana are a characteristic faunal group of the Himalaya-Tibet orogen (HTO). We examine the phylogeny of these frogs to explore alternative biogeographic scenarios for their origin in the Greater Himalaya, namely, immigration, South Tibetan origin, strict vicariance. We sequenced 150 Nanorana samples from 62 localities for three mitochondrial (1,524 bp) and three nuclear markers (2,043 bp) and complemented the data with sequence data available from GenBank. We reconstructed a gene tree, phylogenetic networks, and ancestral areas. Based on the nuDNA, we also generated a time-calibrated species tree. The results revealed two major clades (Nanorana and Quasipaa), which originated in the Lower Miocene from eastern China and subsequently spread into the HTO (Nanorana). Five well-supported subclades are found within Nanorana: from the East, Central, and Northwest Himalaya, the Tibetan Plateau, and the southeastern Plateau margin. The latter subclade represents the most basal group (subgenus Chaparana), the Plateau group (Nanorana) represents the sister clade to all species of the Greater Himalaya (Paa). We found no evidence for an east-west range expansion of Paa along the Himalaya, nor clear support for a strict vicariance model. Diversification in each of the three Himalayan subclades has probably occurred in distinct areas. Specimens from the NW Himalaya are placed basally relative to the highly diverse Central Himalayan group, while the lineage from the Tibetan Plateau is placed within a more terminal clade. Our data indicate a Tibetan origin of Himalayan Nanorana and support a previous hypothesis, which implies that a significant part of the Himalayan biodiversity results from primary diversification of the species groups in South Tibet before this part of the HTO was uplifted to its recent heights.