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.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

The incidence of the anomaly P syndrome in water frogs (Anura, Ranidae, Pelophylax) from the Middle Volga River (Russia)

The anomaly P is a widespread morphological anomaly, which occurs in some groups of amphibians, caused by the trematode parasite Strigea robusta (Digenea: Strigeidae). This anomaly has been previously recorded in water frogs of the genus Pelophylax and toads of the genera Bufo and Bufotes. The anomaly P includes symmetrical polydactyly cases as a mild attenuated form of the complex syndrome, which in severe cases includes strong deformations of hindlimbs and forelimbs. Strigea robusta has a complex 3-host life cycle using planorbid mollusks as the first intermediate hosts, amphibian larvae as the second intermediate hosts, and anatid birds as the definitive hosts. Herein, we described new records of the anomaly P syndrome in water frogs of the genus Pelophylax from the northeastern parts of their ranges. Symmetrical polydactyly (as a mild form of the anomaly P syndrome) was found in 30 individuals of three species of water frogs from seven localities: in 25 individuals of P. lessonae from four waterbodies, in four individuals of P. ridibundus from three waterbodies, and one individual of P. esculentus. In Gusevo pond, three individuals of P. lessonae with severe cases of the syndrome were found. This is the first record of the anomaly P in reliably identified hybridogenetic edible frogs (P. esculentus) that have been identified in nature. Additionally, we provided new data about the occurrence of the anomaly P and the prevalence of the trematode S. robusta in mollusks taken from two water bodies where anomalous water frogs were found.

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.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1482-4) contains supplementary material, which is available to authorized users.

An extinct vertebrate preserved by its living hybridogenetic descendant

Hybridogenesis is a special mode of hybrid reproduction where one parental genome is eliminated and the other is transmitted clonally. We propose that this mechanism can perpetuate the genome of extinct species, based on new genetic data from Pelophylax water frogs. We characterized the genetic makeup of Italian hybridogenetic hybrids (P. kl. hispanicus and esculentus) and identified a new endemic lineage of Eastern-Mediterranean origin as one parental ancestor of P. kl. hispanicus. This taxon is nowadays extinct in the wild but its germline subsists through its hybridogenetic descendant, which can thus be considered as a "semi living fossil". Such rare situation calls for realistic efforts of de-extinction through selective breeding without genetic engineering, and fuels the topical controversy of reviving long extinct species. "Ghost" species hidden by taxa of hybrid origin may be more frequent than suspected in vertebrate groups that experienced a strong history of hybridization and semi-sexual reproduction.

The first taxonomic revaluation of the Iranian water frogs of the genus Pelophylax (Anura: Ranidae) using sequences of the mitochondrial genome

The Eurasian water frog species and their geographic ranges have undergone considerable changes in the last four decades, but the Iranian populations have largely remained unknown. All the Iranian populations of water frogs, despite their vast distribution range have attributed to a single species: Rana ridibunda. In order to understand the phylogenetic relationships and taxonomic status of water frogs of Iran, we collected samples from many populations across the country and used the mitochondrial DNA sequence variation. A data set with a final sequence length of 616 nucleotides was generated for Cyt b from 70 individuals of Pelophylax in which there are 422 invariable sites, 174 variable sites of which 123 were parsimony informative. In total, 43 haplotypes were found (Hd: 0.9752). The result demonstrated that, two major clades with strong support can be identified within the Iranian water frogs. One of these clades that include north western and southwestern populations forms a monophyletic group along with P. bedriagae samples from Turkey. The second clade consists of water frog populations of north and northeastern parts of Iran which in turn is subdivided into two subclades. Inclusion of water frog samples from adjacent areas showed that the second clade of our study is, most likely, a distinct taxonomic entity at species rank with its two subclades indicating two diagnosable subspecies for the clade. In conclusion, we suggest that two distinct species, P. bedriagae and Pelophylax sp., with its two subspecies, should be identified as water frogs of Iran. In Addition, another traditionally reported water frog of Iran, P.ridibundus, most likely should be excluded from the Iranian water frog's checklist.

Extensive mitochondrial heteroplasmy in hybrid water frog (Pelophylax spp.) populations from Southeast Europe

Water frogs of the genus Pelophylax (previous Rana) species have been much studied in Europe for their outstanding reproductive mechanism in which sympatric hybridization between genetically distinct parental species produces diverse genetic forms of viable hybrid animals. The most common hybrid is P. esculentus that carries the genomes of both parental species, P. ridibundus and P. lessonae, but usually transfers the whole genome of only one parent to its offsprings (hybridogenesis). The evolutionary cost of transfer of the intact genome and hence the hemiclonal reproduction is the depletion of heterozygosity in the hybrid populations. Pelophylax esculentus presents an excellent example of the long‐term sustained hybridization and hemiclonal reproduction in which the effects of the low genetic diversity are balanced through the novel mutations and periodic recombinations. In this study, we analyzed the mitochondrial (mt) and microsatellites DNA variations in hybrid Pelophylax populations from southern parts of the Pannonian Basin and a north–south transect of the Balkan Peninsula, which are home for a variety of Pelophylax genetic lineages. The mtDNA haplotypes found in this study corresponded to P. ridibundus and P. epeiroticus of the Balkan – Anatolian lineage (ridibundus–bedriagae) and to P. lessonae and a divergent lessonae haplotype of the lessonae lineage. The mtDNA genomes showed considerable intraspecific variation and geographic differentiation. The Balkan wide distributed P. ridibundus was found in all studied populations and its nuclear genome, along with either the lessonae or the endemic epeiroticus genome, in all hybrids. An unexpected finding was that the hybrid populations were invariably heteroplasmic, that is, they contained the mtDNA of both parental species. We discussed the possibility that such extensive heteroplasmy is a result of hybridization and it comes from regular leakage of the paternal mtDNA from a sperm of one species that fertilizes eggs of another. In this case, the mechanisms that protect the egg from heterospecific fertilization and further from the presence of sperm mtDNA could become compromised due to their differences and divergence at both, mitochondrial and nuclear DNA. The heteroplasmy once retained in the fertilized egg could be transmitted by hybrid backcrossing to the progeny and maintained in a population over generations. The role of interspecies and heteroplasmic hybrid animals due to their genomic diversity and better fitness compare to the parental species might be of the special importance in adaptations to miscellaneous and isolated environments at the Balkan Peninsula.