Evolution of serum albumin intron-1 is shaped by a 5' truncated non-long terminal repeat retrotransposon in western Palearctic water frogs (Neobatrachia)

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.

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.

Hybridogenesis in the Water Frogs from Western Russian Territory: Intrapopulation Variation in Genome Elimination

Hybridogenesis in an interspecific hybrid frog is a coupling mechanism in the gametogenic cell line that eliminates the genome of one parental species with endoduplication of the remaining genome of the other parental species. It has been intensively investigated in the edible frog Pelophylax kl. esculentus (RL), a natural hybrid between the marsh frog P. ridibundus (RR) and the pool frog P. lessonae (LL). However, the genetic mechanisms involved remain unclear. Here, we investigated the water frogs in the western Russian territory. In three of the four populations, we genetically identified 16 RL frogs living sympatrically with the parental LL species, or with both parental species. In addition, two populations contained genome introgression with another species, P. bedriagae (BB) (a close relative of RR). In the gonads of 13 RL frogs, the L genome was eliminated, producing gametes of R (or R combined with the B genome). In sharp contrast, one RL male eliminated the L or R genome, producing both R and L sperm. We detected a variation in genome elimination within a population. Based on the genetic backgrounds of RL frogs, we hypothesize that the introgression of the B genome resulted in the change in choosing a genome to be eliminated.

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.

Preliminary genetic data suggest the occurrence of the Balkan water frog, Pelophylax kurtmuelleri, in southwestern Poland

Recent molecular studies have detected the occurrence of exotic water frog species (Pelophylaxsp.) in central and western European populations. Here, we report genetic evidence for the occurrence of the Balkan water frog,Pelophylax kurtmuelleri, in southwestern Poland. We found a high frequency of an allele of serum albumin intron-1 and a mitochondrial cytochromebhaplotype specific for this southern taxon in frogs from the Barycz river drainage system. We interpret this finding as evidence of admixture betweenP. kurtmuelleriand the localridibundus-esculentuswater frog population. The origin of the exoticP. kurtmuellerimitochondrial and nuclear alleles in southwestern Poland could be due to (i) hybridization after a human-mediated introduction ofP. kurtmuelleri, (ii) the persistence of ancestral polymorphism in central EuropeanP. ridibundus, or (iii) hybridization betweenP. kurtmuelleriandP. ridibundusin the Balkans followed by the northward expansion of admixedP. ridibundus. Identical mtDNA haplotypes found in southwestern Poland and localities on the borders between Greece, Albania and Macedonia suggest that this region harboured the source population ofP. kurtmuelleriat the studied site.

Is premeiotic genome elimination an exclusive mechanism for hemiclonal reproduction in hybrid males of the genus Pelophylax?

Background

The ability to eliminate a parental genome from a eukaryotic germ cell is a phenomenon observed mostly in hybrid organisms displaying an alternative propagation to sexual reproduction. For most taxa, the underlying cellular pathways and timing of the elimination process is only poorly understood. In the water frog hybrid Pelophylax esculentus (parental taxa are P. ridibundus and P. lessonae) the only described mechanism assumes that one parental genome is excluded from the germline during metamorphosis and prior to meiosis, while only second genome enters meiosis after endoreduplication. Our study of hybrids from a P. ridibundus—P. esculentus-male populations known for its production of more types of gametes shows that hybridogenetic mechanism of genome elimination is not uniform.

Results

Using comparative genomic hybridization (CGH) on mitotic and meiotic cell stages, we identified at least two pathways of meiotic mechanisms. One type of Pelophylax esculentus males provides supporting evidence of a premeiotic elimination of one parental genome. In several other males we record the presence of both parental genomes in the late phases of meiotic prophase I (diplotene) and metaphase I.

Conclusion

Some P. esculentus males have no genome elimination from the germ line prior to meiosis. Considering previous cytological and experimental evidence for a formation of both ridibundus and lessonae sperm within a single P. esculentus individual, we propose a hypothesis that genome elimination from the germline can either be postponed to the meiotic stages or absent altogether in these hybrids.

Electronic supplementary material

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

Modelling interspecific hybridization with genome exclusion to identify conservation actions: the case of native and invasive Pelophylax waterfrogs

Interspecific hybridization occurs in nature but can also be caused by human actions. It often leads to infertile or fertile hybrids that exclude one parental genome during gametogenesis, escaping genetic recombination and introgression. The threat that genome-exclusion hybridization might represent on parental species is poorly understood, especially when invasive species are involved. Here, we show how to assess the effects of genome-exclusion hybridization and how to elaborate conservation actions by simulating scenarios using a model of nonintrogressive hybridization. We examine the case of the frog Pelophylax ridibundus, introduced in Western Europe, which can hybridize with the native Pelophylax lessonae and the pre-existing hybrid Pelophylax esculentus, maintained by hybridogenesis. If translocated from Southern Europe, P. ridibundus produces new sterile hybrids and we show that it mainly threatens P. esculentus. Translocation from Central Europe leads to new fertile hybrids, threatening all native waterfrogs. Local extinction is demographically mediated via wasted reproductive potential or via demographic flow through generations towards P. ridibundus. We reveal that enlarging the habitat size of the native P. lessonae relative to that of the invader is a promising conservation strategy, avoiding the difficulties of fighting the invader. We finally stress that nonintrogressive hybridization is to be considered in conservation programmes.

A simplified molecular method for distinguishing among species and ploidy levels in European water frogs (Pelophylax)

Western Palearctic water frogs in the genus Pelophylax are a set of morphologically similar anuran species that form hybridogenetic complexes. Fully reliable identification of species and especially of hybrid ploidy depends on karyological and molecular methods. In central Europe, native water frog populations consist of the Pelophylax esculentus complex, that is, P. lessonae (LL), P. ridibundus (RR) and the hybrid form P. esculentus that can have different karyotypes (RL, LLR and RRL). We developed existing molecular methods further and propose a simple PCR method based on size-differences in the length of the serum albumin intron-1 and the RanaCR1, a non-LTR retrotransposon of the chicken repeat (CR) family. This PCR yields taxon-specific banding patterns that can easily be screened by standard agarose gel electrophoresis and correctly identify species in all of the 160 samples that had been identified to karyotype with other methods. To distinguish ploidy levels in LR, LLR and RRL specimens, we used the ratio of the peak heights of the larger (ridibundus specific) to the smaller (lessonae specific) bands of fluorescently labelled PCR products resolved on a capillary DNA sequencer and obtained a correct assignment of the karyotype in 93% of cases. Our new method will cut down time and expenses drastically for a reliable identification of water frogs of the P. esculentus complex and potentially for identification of other hybridogenetic complexes and/or taxa, and it even serves as a good indicator of the ploidy status of hybrid individuals.

RrS1-like sequences of water frogs from Central Europe and around the Aegean Sea: chromosomal organization, evolution, possible function

RrS1-like sequences of water frogs (genus Pelophylax) display varied genomic organization, whereas the centromeric hybridization pattern reveals species-specific differences. Using fluorescent in situ hybridization, Pelophylax cf. bedriagae, Pelophylax kurtmuelleri, and Pelophylax ridibundus showed a hybridization signal at centromeres of chromosomes 1-5, but in P. kurtmuelleri the medium-small chromosome labeled was 10 rather than 8. Pelophylax cretensis had almost 16 of 26 centromeres labeled, as did Pelophylax lessonae from Poland when its chromosomes are hybridized with a homologous probe. When StuI-digested genomic DNA was hybridized with RrS1 probe, hybridization ladders for P. ridibundus from Poland have evenly spaced steps (about 100 bp) of uniform intensity from about 200 bp upward. Steps in hybridization ladders from circum-Aegean taxa vary in intensity: larger, odd-numbered steps are often fainter. A strong double band (800/900 bp) in Anatolian P. cf. bedriagae, emphasized by a weak 700 bp band, distinguishes them from P. kurtmuelleri from the Peloponnisos, in which the 900 bp band is almost absent. The ladder in P. cretensis lacks odd-numbered steps. A and B repeats, observed originally within the RrS1 satellite of P. ridibundus, occur also in the circum-Aegean frogs and in P. lessonae, Pelophylax epeiroticus, Pelophylax saharicus, and Pelophylax shqipericus. It is plausible that AB dimers or ABB trimers rather than A or B monomers correspond to functional/evolutionary units. The presence of regions similar to yeast CDEs and mammalian CENP-B boxes suggests a role for RrS1 sequences in centromere organization.

Phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs have been determined by geological processes and climate change in the Late Cenozoic

AIM: Our aims were to assess the phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs and to estimate divergence times using different geological scenarios. We related divergence times to past geological events and discuss the relevance of our data for the systematics of eastern Mediterranean water frogs. LOCATION: The eastern Mediterranean region. METHODS: Genetic diversity and divergence were calculated using sequences of two protein-coding mitochondrial (mt) genes: ND2 (1038 bp, 119 sequences) and ND3 (340 bp, 612 sequences). Divergence times were estimated in a Bayesian framework under four geological scenarios representing alternative possible geological histories for the eastern Mediterranean. We then compared the different scenarios using Bayes factors and additional geological data. RESULTS: Extensive genetic diversity in mtDNA divides eastern Mediterranean water frogs into six main haplogroups (MHG). Three MHGs were identified on the Anatolian mainland; the most widespread MHG with the highest diversity is distributed from western Anatolia to the northern shore of the Caspian Sea, including the type locality of Pelophylax ridibundus. The other two Anatolian MHGs are restricted to south-eastern Turkey, occupying localities west and east of the Amanos mountain range. One of the remaining three MHGs is restricted to Cyprus; a second to the Levant; the third was found in the distribution area of European lake frogs (P. ridibundus group), including the Balkans. MAIN CONCLUSIONS: Based on geological evidence and estimates of genetic divergence we hypothesize that the water frogs of Cyprus have been isolated from the Anatolian mainland populations since the end of the Messinian salinity crisis (MSC), i.e. since c. 5.5-5.3 Ma, while our divergence time estimates indicate that the isolation of Crete from the mainland populations (Peloponnese, Anatolia) most likely pre-dates the MSC. The observed rates of divergence imply a time window of c. 1.6-1.1 million years for diversification of the largest Anatolian MHG; divergence between the two other Anatolian MHGs may have begun about 3.0 Ma, apparently as a result of uplift of the Amanos Mountains. Our mtDNA data suggest that the Anatolian water frogs and frogs from Cyprus represent several undescribed species.