Capture and return of sexual genomes by hybridogenetic frogs provide clonal genome enrichment in a sexual species.

Hybridogenesis is a reproductive tool for sexual parasitism. Hybridogenetic hybrids use gametes from their sexual host for their own reproduction, but sexual species gain no benefit from such matings as their genome is later eliminated. Here, we examine the presence of sexual parasitism in water frogs through crossing experiments and genome-wide data. We specifically focus on the famous Central-European populations where Pelophylax esculentus males (hybrids of P. ridibundus and P. lessonae) live with P. ridibundus. We identified a system where the hybrids commonly produce two types of clonal gametes (hybrid amphispermy). The haploid lessonae genome is clonally inherited from generation to generation and assures the maintenance of hybrids through a process, in which lessonae sperm fertilize P. ridibundus eggs. The haploid ridibundus genome in hybrids received from P. ridibundus a generation ago, is perpetuated as clonal ridibundus sperm and used to fertilize P. ridibundus eggs, yielding female P. ridibundus progeny. These results imply animal reproduction in which hybridogenetic taxa are not only sexual parasites, but also participate in the formation of a sexual taxon in a remarkable way. This occurs through a process by which sexual gametes are being captured, converted to clones, and returned to sexual populations in one generation.

All-male hybrids of a tetrapod Pelophylax esculentus share its origin and genetics of maintenance.

BACKGROUND: Sexual parasites offer unique insights into the reproduction of unisexual and sexual populations. Because unisexuality is almost exclusively linked to the female sex, most studies addressed host-parasite dynamics in populations where sperm-dependent females dominate. Pelophylax water frogs from Central Europe include hybrids of both sexes, collectively named P. esculentus. They live syntopically with their parental species P. lessonae and/or P. ridibundus. Some hybrid lineages consist of all males providing a chance to understand the origin and perpetuation of a host-parasite (egg-dependent) system compared to sperm-dependent parthenogenesis. METHODS: We focused on P. ridibundus-P. esculentus populations where P. ridibundus of both sexes lives together with only diploid P. esculentus males. Based on 17 microsatellite markers and six allozyme loci, we analyzed (i) the variability of individual genomes, (ii) the reproductive mode(s) of all-male hybrids, and (iii) the genealogical relationships between the hybrid and parental genomes. RESULTS: Our microsatellite data revealed that P. esculentus males bear Mendelian-inherited ridibundus genomes while the lessonae genome represents a single clone. Our data indicate that this clone did not recently originate from adjacent P. lessonae populations, suggesting an older in situ or ex situ origin. CONCLUSIONS: Our results confirm that also males can perpetuate over many generations as the unisexual lineage and successfully compete with P. ridibundus males for eggs provided by P. ridibundus females. Natural persistence of such sex-specific hybrid populations allows to studying the similarities and differences between male and female reproductive parasitism in many biological settings.

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.

The complete mitochondrial genome of Telmatobufo australis (Amphibia: Anura: Calyptocephalellidae).

The mitochondrial (mt) genome of Telmatobufo australis is a circular molecule of 17,989 bp in length, comprising 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. Gene order and content are identical to those previously reported from other neobatrachian mt genomes. Two protein-coding genes (COI and ATP6) presumably used GTG as start codons while COIII possessed an incomplete stop codon.

Genetic diversity and distribution patterns of diploid and polyploid hybrid water frog populations (Pelophylax esculentus complex) across Europe.

Polyploidization is a rare yet sometimes successful way for animals to rapidly create geno- and phenotypes that may colonize new habitats and quickly adapt to environmental changes. In this study, we use water frogs of the Pelophylax esculentus complex, comprising two species (Pelophylax lessonae, genotype LL; Pelophylax ridibundus, RR) and various diploid (LR) and triploid (LLR, LRR) hybrid forms, summarized as P. esculentus, as a model for studying recent hybridization and polyploidization in the context of speciation. Specifically, we compared the geographic distribution and genetic diversity of diploid and triploid hybrids across Europe to understand their origin, maintenance and potential role in hybrid speciation. We found that different hybrid and parental genotypes are not evenly distributed across Europe. Rather, their genetic diversity is structured by latitude and longitude and the presence/absence of parental species but not of triploids. Highest genetic diversity was observed in central and eastern Europe, the lowest in the northwestern parts of Europe. This gradient can be explained by the decrease in genetic diversity during postglacial expansion from southeastern glacial refuge areas. Genealogical relationships calculated on the basis of microsatellite data clearly indicate that hybrids are of multiple origin and include a huge variety of parental genomes. Water frogs in mixed-ploidy populations without any parental species (i.e. all-hybrid populations) can be viewed as evolutionary units that may be on their way towards hybrid speciation. Maintenance of such all-hybrid populations requires a continuous exchange of genomes between diploids and triploids, but scenarios for alternative evolutionary trajectories are discussed.

LTR retroelements are intrinsic components of transcriptional networks in frogs.

BACKGROUND: LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology. RESULTS: Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active. CONCLUSIONS: The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple "junk DNA" or "harmful genomic parasites" but essential components of the transcriptomic machinery in vertebrates.

Batrachochytrium dendrobatidis in Germany: distribution, prevalences, and prediction of high risk areas.

In Germany, the pathogenic fungus Batrachochytrium dendrobatidis (Bd) was detected in 11 indigenous frog species, 4 newt species, and 1 salamander species in 64 out of the 181 locations (35%) investigated. Among the 3450 samples collected between 2003 and 2011, 284 (8.2%) were positive for Bd infections. The highest prevalences were observed in Alytes obstetricans (17.8% of individuals, 20% of populations), followed by Ichthyosaura alpestris (14.7%, 22.2%), Bombina variegata (13.9%, 38.5%), and water frogs comprising 2 species, Pelophylax lessonae and P. ridibundus, and their hybrid form P. esculentus (13.5%, 29.0%). Bd is widespread; areas of higher prevalence were detected in eastern, southeastern, western, and southwestern Germany. Our data indicate that drift fencing of amphibians is not a risk factor for the anthropogenic spread of Bd. Although chytridiomycosis outbreaks have never been observed in Germany, it cannot be excluded that Bd infections affect the dynamics of local amphibian populations. Among the questions still to be answered is whether juveniles are more susceptible to Bd infections than adults. Further work, especially long-term observations including capture-mark-recapture studies, is required to clarify the impact Bd has on amphibians in Germany and Central Europe.

West Africa - a safe haven for frogs? A sub-continental assessment of the chytrid fungus (Batrachochytrium dendrobatidis).

A putative driver of global amphibian decline is the panzootic chytrid fungus Batrachochytrium dendrobatidis (Bd). While Bd has been documented across continental Africa, its distribution in West Africa remains ambiguous. We tested 793 West African amphibians (one caecilian and 61 anuran species) for the presence of Bd. The samples originated from seven West African countries - Bénin, Burkina Faso, Côte d'Ivoire, Ghana, Guinea, Liberia, Sierra Leone - and were collected from a variety of habitats, ranging from lowland rainforests to montane forests, montane grasslands to humid and dry lowland savannahs. The species investigated comprised various life-history strategies, but we focused particularly on aquatic and riparian species. We used diagnostic PCR to screen 656 specimen swabs and histology to analyse 137 specimen toe tips. All samples tested negative for Bd, including a widespread habitat generalist Hoplobatrachus occipitalis which is intensively traded on the West African food market and thus could be a potential dispersal agent for Bd. Continental fine-grained (30 arc seconds) environmental niche models suggest that Bd should have a broad distribution across West Africa that includes most of the regions and habitats that we surveyed. The surprising apparent absence of Bd in West Africa indicates that the Dahomey Gap may have acted as a natural barrier. Herein we highlight the importance of this Bd-free region of the African continent - especially for the long-term conservation of several threatened species depending on fast flowing forest streams (Conraua alleni ("Vulnerable") and Petropedetes natator ("Near Threatened")) as well as the "Critically Endangered" viviparous toad endemic to the montane grasslands of Mount Nimba (Nimbaphrynoides occidentalis).

Balancing a cline by influx of migrants: a genetic transition in water frogs of eastern Greece.

Variation patterns of allozymes and of ND3 haplotypes of mitochondrial DNA reveal a zone of genetic transition among western Palearctic water frogs extending across northeastern Greece and European Turkey. At the western end of the zone, allozymes characteristic of Central European frogs known as Pelophylax ridibundus predominate, whereas at the eastern end, alleles characteristic of western Anatolian water frogs (P. cf. bedriagae) prevail. The ND3 haplotypes reveal 2 major clades, 1 characteristic of Anatolian frogs, the other of European; the European clade itself has distinct eastern and western subclades. Both the 2 major clades and the 2 subclades overlap within the transition zone. Using Bayesian model selection methods, allozyme data suggest considerable immigration into the Nestos River area from eastern and western populations. In contrast, the ND3 data suggest that migration rates are so high among all locations that they form a single panmictic unit; the best model for allozymes is second best for mitochondrial DNA (mtDNA). Nuclear markers (allozymes), which have roughly 4 times as deep a coalescent history as mtDNA data and thus may reflect patterns over a longer time, indicate that eastern and western refugial populations have expanded since deglaciation (in the last 10,000 years) and have met near the Nestos River, whereas the mtDNA with its smaller effective population size has already lost the signal of partitioning into refugia.

Delimitation of the Thoracosphaeraceae (Dinophyceae), including the calcareous dinoflagellates, based on large amounts of ribosomal RNA sequence data.

The phylogenetic relationships of the Dinophyceae (Alveolata) are not sufficiently resolved at present. The Thoracosphaeraceae (Peridiniales) are the only group of the Alveolata that include members with calcareous coccoid stages; this trait is considered apomorphic. Although the coccoid stage apparently is not calcareous, Bysmatrum has been assigned to the Thoracosphaeraceae based on thecal morphology. We tested the monophyly of the Thoracosphaeraceae using large sets of ribosomal RNA sequence data of the Alveolata including the Dinophyceae. Phylogenetic analyses were performed using Maximum Likelihood and Bayesian approaches. The Thoracosphaeraceae were monophyletic, but included also a number of non-calcareous dinophytes (such as Pentapharsodinium and Pfiesteria) and even parasites (such as Duboscquodinium and Tintinnophagus). Bysmatrum had an isolated and uncertain phylogenetic position outside the Thoracosphaeraceae. The phylogenetic relationships among calcareous dinophytes appear complex, and the assumption of the single origin of the potential to produce calcareous structures is challenged. The application of concatenated ribosomal RNA sequence data may prove promising for phylogenetic reconstructions of the Dinophyceae in future.

A new PCR-RFLP-based method for an easier systematic affiliation of European water frogs.

We describe a non-invasive, PCR-RFLP-based method that allows reliable determination of the European water frog species Pelophylax lessonae and Pelophylax ridibundus and the hybrid form Pelophylax esculentus. Maximum-likelihood analysis of ITS2 sequences revealed two robust monophyletic clades corresponding to water frogs of the P. lessonae and P. ridibundus groups. Three restriction enzymes (KpnI, HaeII, and SmaI) were used to digest three conserved ITS2 domains. Taxonomic identification was unambiguous; the three restriction enzymes gave the same results. A French reference sample was identified using allozyme electrophoresis. Our PCR-RFLP method confirmed circa 83% of identification of the allozyme method. We conclude that the difference between identifications was caused by introgression.

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.

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

A 5' truncated non-LTR CR1-like retrotransposon, named RanaCR1, was identified in the serum albumin intron-1 (SAI-1) of at least seven species of western Palearctic water frogs (WPWF). Based on sequence similarity of the carboxy-terminal region (CTR) of ORF2 and/or the highly conserved 3' untranslated region (3' UTR), RanaCR1-like elements occur also in the genome of Xenopus tropicalis and Rana temporaria. Unlike other CR1 elements, RanaCR1 contains a CA microsatellite in its 3' UTR. The low nucleotide diversity of the 3' UTR compared to the CTR and to SAI-1 suggests that this region still plays a role in WPWF, either as a structure-stabilizing element, or within a species-specific transcriptional network. Length variation of water frog SAI-1 sequences is caused by deletions that extend in some cases beyond the 5' or 3' ends of RanaCR1, probably a result of selection for structural and functional stability of the primary transcript. The impact of RanaCR1 on SAI-1 evolution is also indicated by the significant negative correlation between the length of both SAI-1 and RanaCR1 and the percentage GC content of RanaCR1. Both SAI-1 and RanaCR1 sequences support the sister group relationship of R. perezi and R. saharica, which are placed in the phylogenetic tree at a basal position, the sister clade to other water frog taxa. It also supports the monophyly of the R. lessonae group; of Anatolian water frogs (R. cf. bedriagae), which are not conspecific with R. bedriagae, and of the European ridibunda group. Within the ridibunda clade, Greek frogs are clearly separated, supporting the hypothesis that Balkan water frogs represent a distinct species. Frogs from Atyrau (Kazakhstan), the type locality of R. ridibunda, were heterozygous for a ridibunda and a cf. bedriagae specific allele.

Gametogenesis of intergroup hybrids of hemiclonal frogs.

European water frog hybrids Rana esculenta (R. ridibundaxR. lessonae) reproduce hemiclonally, by hybridogenesis: in the germ line they exclude the genome of one parental species and produce haploid gametes with an unrecombined genome of the other parental species. In the widespread L-E population system, both sexes of hybrids (E) coexist with R. lessonae (L). They exclude the lessonae genome and produce ridibunda gametes. In the R-E system, hybrid males coexist with R. ridibunda (R); they exclude either their ridibunda or their lessonae genome and produce sperm with a lessonae or with a ridibunda genome or a mixture of both kinds of sperm. We examined 13 male offspring, 12 of which were from crosses between L-E system and R-E system frogs. All were somatically hybrid. With one exception, they excluded the lessonae genome in the germ line and subsequently endoreduplicated the ridibunda genome. Spermatogonial metaphases contained a haploid or a diploid number of ridibunda chromosomes, identified through in situ hybridization to a satellite DNA marker, and by spermatocyte I metaphases containing a haploid number of ridibunda bivalents. The exception, an F1 hybrid between L-E system R. lessonae and R-E system R. ridibunda, was not hybridogenetic, showed no genome exclusion, and evidenced a disturbed gametogenesis resulting from the combination of two heterospecific genomes. None of the hybridogenetic hybrids showed any cell lines excluding the ridibunda genome, the pattern most frequent in hybrids of the R-E system, unique to that system, and essential for its persistence. A particular combination of R-E system lessonae and R-E system ridibunda genomes seems necessary to induce the R-E system type of hemiclonal gametogenesis.

Phylogeny of calcareous dinoflagellates as inferred from ITS and ribosomal sequence data.

The phylogenetic relationships of calcareous dinoflagellates (i.e., Calciodinellaceae and Thoracosphaera) are investigated. Molecular data from the ribosomal 5.8S rRNA and highly conserved motifs of the ITS1 show Calciodinellaceae s.l. to be monophyletic when few non-calcareous taxa are included. They segregate into three monophyletic assemblages in a molecular analysis that considers the 5.8S rRNA and both the Internal Transcribed Spacer regions ITS1 and ITS2: a clade comprising species of Ensiculifera and Pentapharsodinium (E/P-clade), Scrippsiella s.l. (including fossil-based taxa such as Calciodinellum and Calcigonellum), and a heterogeneous group (T/P-clade) of calcareous (e.g., Thoracosphaera) and non-calcareous taxa (e.g., the highly toxic Pfiesteria). The potential to produce calcareous structures is considered as apomorphic within alveolates, and non-calcareous taxa nesting with calcareous dinoflagellates may have reduced calcification secondarily. Molecular results do not contradict general evolutionary scenarios provided by previous morphological (mainly paleontological) investigations.

Secondary structure models of the nuclear internal transcribed spacer regions and 5.8S rRNA in Calciodinelloideae (Peridiniaceae) and other dinoflagellates.

Secondary structure models of the 5.8S rRNA and both internal transcribed spacers (ITS1 and ITS2) are proposed for Calciodinelloideae (Peridiniaceae) and are also plausible for other dinoflagellates. The secondary structure of the 5.8S rRNA corresponds to previously developed models, with two internal paired regions and at least one 5.8S rRNA-28S rRNA interaction. A general secondary structure model of ITS1 for Calciodinelloideae (and other dinoflagellates), consisting of an open multibranch loop with three major helices, is proposed. The homology of these paired regions with those found in other taxa, published in previous studies (e.g. yeast, green algae and Platyhelmithes) remains to be determined. Finally, a general secondary structure model of ITS2 for Calciodinelloideae (and other dinoflagellates) is reconstructed. Based on the 5.8S rRNA-28S rRNA interaction, it consists of a closed multibranch loop, with four major helices. At least helix III and IV have homology with paired regions found in other eukaryotic taxa (e.g. yeast, green algae and vertebrates). Since the secondary structures of both ITS regions are more conserved than the nucleotide sequences, their analysis helps in understanding molecular evolution and increases the number of structural characters. Thus, the structure models developed in this study may be generally useful for future phylogenetic analyses.