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.

Electronic supplementary material

The online version of this article (10.1186/s13293-018-0172-z) contains supplementary material, which is available to authorized users.

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.

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.

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.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

Photo-physical properties of 2-(1-ethynylpyrene)-adenosine: influence of hydrogen bonding on excited state properties

The photo-physical properties of 2-(1-ethynylpyrene)-adenosine (PyA), a fluorescent probe for RNA dynamics, were examined by solvation studies. The excited-state dynamics display the influence of the vicinity on the spectral features. Combining improved transient absorption and streak camera measurements along with a new analysis method provide a detailed molecular picture of the photophysics. After intramolecular vibrational energy redistribution (IVR), two distinct states are observed. Solvent class (protic/aprotic) and permittivity strongly affect the properties of these states and their population ratio. As a result their emission spectrum is altered, while the fluorescence quantum yield and the overall lifetime remain nearly unchanged. Consequently, the hitherto existing model of the photophysics is herein refined and extended. The findings can serve as basis for improving the information content of measurements with PyA as a label in RNA.

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.

Batrachochytrium dendrobatidis in amphibians of Cameroon, including first records for caecilians

Amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) has been hypothesised to be an indigenous parasite of African amphibians. In Cameroon, however, previous surveys in one region (in the northwest) failed to detect this pathogen, despite the earliest African Bd having been recorded from a frog in eastern Cameroon, plus one recent record in the far southeast. To reconcile these contrasting results, we present survey data from 12 localities across 6 regions of Cameroon from anurans (n = 1052) and caecilians (n = 85) of ca. 108 species. Bd was detected in 124 amphibian hosts at 7 localities, including Mt. Oku, Mt. Cameroon, Mt. Manengouba and lowland localities in the centre and west of the country. None of the hosts were observed dead or dying. Infected amphibian hosts were not detected in other localities in the south and eastern rainforest belt. Infection occurred in both anurans and caecilians, making this the first reported case of infection in the latter order (Gymnophiona) of amphibians. There was no significant difference between prevalence and infection intensity in frogs and caecilians. We highlight the importance of taking into account the inhibition of diagnostic qPCR in studies on Bd, based on all Bd-positive hosts being undetected when screened without bovine serum albumin in the qPCR mix. The status of Bd as an indigenous, cosmopolitan amphibian parasite in Africa, including Cameroon, is supported by this work. Isolating and sequencing strains of Bd from Cameroon should now be a priority. Longitudinal host population monitoring will be required to determine the effects, if any, of the infection on amphibians in Cameroon.

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.

New records of Batrachochytrium dendrobatidis in Chilean frogs

We used molecular techniques to examine 11 species of frogs in 6 localities in southern Chile to ascertain the incidence of the chytrid fungus Batrachochytrium dendrobatidis (Bd). We detected the fungus in 2 localities (Coñaripe and Raúl Marín Balmaceda) in 3 species: Batrachyla leptopus, Pleurodema thaul and Rhinoderma darwinii. Our findings expand the list of Bd hosts to include B. leptopus and P. thaul and extend the spatial distribution in Chile to include the southernmost Bd record at Raúl Marín Balmaceda.

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.

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.

Widespread unidirectional transfer of mitochondrial DNA: a case in western Palaearctic water frogs

Interspecies transfer of mitochondrial (mt) DNA is a common phenomenon in plants, invertebrates and vertebrates, normally linked with hybridization of closely related species in zones of sympatry or parapatry. In central Europe, in an area north of 48 degrees N latitude and between 8 degrees and 22 degrees E longitude, western Palaearctic water frogs show massive unidirectional introgression of mtDNA: 33.7% of 407 Rana ridibunda possessed mtDNA specific for Rana lessonae. By contrast, no R. lessonae with R. ridibunda mtDNA was observed. That R. ridibunda with introgressed mitochondrial genomes were found exclusively within the range of the hybrid Rana esculenta and that most hybrids had lessonae mtDNA (90.4% of 335 individuals investigated) is evidence that R. esculenta serves as a vehicle for transfer of lessonae mtDNA into R. ridibunda. Such introgression has occurred several times independently. The abundance and wide distribution of individuals with introgressed mitochondrial genomes show that R. lessonae mt genomes work successfully in a R. ridibunda chromosomal background despite their high sequence divergence from R. ridibunda mtDNAs (14.2-15.2% in the ND2/ND3 genes). Greater effectiveness of enzymes encoded by R. lessonae mtDNA may be advantageous to individuals of R. ridibunda and probably R. esculenta in the northern parts of their ranges.

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.

Divergence in mitochondrial DNA of Near Eastern water frogs with special reference to the systematic status of Cypriote and Anatolian populations (Anura, Ranidae)

Abstract

Water frogs from Anatolia, Syria, Jordan, and central Asia (Turkmenistan, Uzbekistan, Kazakhstan) were compared on the basis of the complete mitochondrial (mt) ND3 gene (340 bp), two flanking mt t-RNA gene fragments (26 bp), and a 374 bp fragment of the mt 12S rRNA gene. A total of 27 haplotypes were found among the investigated individuals. Anatolian water frogs differed from Syrian and Jordanian Rana bedriagae by 2.2-3.4% of the analysed sites. The observed divergence (2.8-4.1%) between the Cypriote water frogs and frogs from the surrounding mainland (southern Turkey, west Syria) was in the same range as between R. bedriagae and European R. ridibunda (3.1-3.9%). These results suggest that neither the Cypriote nor the Anatolian water frogs represent R. bedriagae. Furthermore, maximum parsimony and neighbour-joining trees showed a clear subdivision of Asian water frogs into three Anatolian lineages, two central Asian lineages, a Cypriote lineage and the bedriagae lineage. In all trees the Cypriote lineage branches off first and a clade formed by two Anatolian lineages is placed as the sister group of water frogs from Turkmenistan and Uzbekistan, whereas the phylogenetic positions of R. bedriagae, the Ceyhan lineage, the Kazakhstan lineage and R. ridibunda remain unclear.