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.

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.

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.

Deleterious alleles and differential viability in progeny of natural hemiclonal frogs.

Abstract.-Spontaneous deleterious mutations are expected to accumulate through Muller's ratchet in clonally reproducing organisms and may lead to their extinction. We study deleterious mutations and their effects in a system of European frogs. Rana esculenta (RL), natural hybrids R. ridibunda (RR) X R. lessonae (LL), reproduce hemiclonally; both sexes exclude the L genome in the germ line and produce unrecombined R gametes; hybridity is restored each generation by matings of RL with coexisting LL. Different allozyme-defined hybrid hemiclones (R genome haplotypes) are thought to have originated independently from primary hybridizations RR x LL. Natural matings between two hybrids usually lead to inviable RR tadpoles. This inviability is thought to result from unmasked deleterious alleles on the clonally transmitted R genomes. Most simply it reflects homozygosity for recessive deleterious alleles at particular loci; alternatively (consistent with absence of RR adults in multiclonal populations) it may reflect hemiclone-specific sets of incompletely recessive deleterious mutations that cumulatively cause inviability when two such genomes are combined. If inviability results from the former, progeny of two hybrids of different hemiclones, whether allopatric or coexisting, should be viable, because it is improbable that their R genomes share recessive deleterious alleles at the same set of loci; if inviability results from the latter, progeny of hybrids of different hemiclones should be inviable, especially when hybrid lineages are old. We tested these hypotheses in artificial crosses, using frogs from three regions: hemiclonal hybrids outside R. ridibunda's range from northern Switzerland (two abundant coexisting allozyme-defined hemiclones; estimated lineage age < or = 5,000 generations) and from Sicily, Italy (one hemiclone; estimated age > or = 25,000 generations) and R. ridibunda from Poland. We generated RR progeny, which we reared under benign conditions in the laboratory, by crossing (1) two hybrids from the same region (H x H local); (2) two hybrids from different regions (H X H foreign); (3) hybrids and R. ridibunda (H X R); and (4) two R. ridibunda (R X R). Survival to metamorphosis was similar and high for R x R, H X H foreign, and H X R, whereas all tadpoles of H X H local died before metamorphosis. This supports the hypothesis that homozygosity for recessive deleterious mutations at particular loci causes inviability. Crosses within and between the two coexisting hemiclones from Switzerland were, however, equally inviable. This result may reflect episodic sexual recombination in RR progeny from exceptional successful interclonal hybrid X hybrid matings, followed by matings of such RR with LL. This process would both slow down or halt Muller's ratchet and disrupt genetic independence of coexisting hemiclones, so that the same remaining deleterious R alleles could exist in different allozyme-defined hemiclones. Whereas all data are consistent with the prediction of Muller's ratchet operating on clonally transmitted R genomes of natural hybrid lineages, they are insufficient to demonstrate such operation, because deleterious recessives that mutated after clone formation and those that preexisted in the R. ridibunda source populations that formed the hemiclonal lineages are not distinguished. The possibility of episodic sexual recombination must be carefully taken into account when studying Muller's ratchet in natural populations of this Rana system.

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.

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.

GEOLOGICALLY DATED SEA BARRIERS CALIBRATE A PROTEIN CLOCK FOR AEGEAN WATER FROGS.

Reliable estimates of phylogenetic relationships and divergence times are a crucial requirement for many evolutionary studies, but are usually difficult because fossils are scarce and their interpretation is often uncertain. Frogs are fresh water animals that generally are unable to cross salt water barriers (their skin is readily permeable to both salt and water). The geologically determined ages of salt water barriers that isolate related frog populations thus provide an independent measure of the minimum date of genetic divergence between pairs of such populations. For the genetically well-studied western Palearctic water frogs (Rana esculenta group), the Aegean region provides an ideal area for determining the relationship between genetic divergence and time of spatial isolation, using a nested set of geologically determined isolation times (12,000 yr, 200,000 yr, 1.8 Myr, 2-3 Myr, and 5.2 Myr). Using 31 electrophoretic loci for 33 pairs of neighboring frog populations, a linear relationship between geologically determined isolation time and Hillis' modified Nei genetic distance was found: D*Nei = (0.04 ± 0.01) + (0.10 ± 0.01) isolation time [Myr] corresponding to an average divergence rate ("molecular clock" pace) of 0.10 D*Nei /Myr (0.10 DNei /Myr). This rate is in the range of previous estimates reported for protein electrophoretic data; the value is conservative because relatively few of the loci used are "fast evolvers" (13%; sAAT, ALB, EST-5, MPI). Removing these fast evolvers from the analysis results in 0.08 D*Nei /Myr (0.08 DNei /Myr). The confidence limits for estimation of the divergence time given the genetic distance are large, but unusually narrow for this kind of study; they permit us to estimate divergence times during the Pliocene and Miocene. Few previous studies, including sequence analyses, have provided reasonable estimates of divergence time for the Pliocene. A test using the outgroup taxa Rana perezi and Rana saharica (also isolated for 5.2 Myr by the Strait of Gibraltar) fits the calibration well: observed genetic Nei distance D*Nei = 0.55, expected D*Nei = 0.56. The calculated divergence times, based on this absolute molecular clock, suggest a series of speciation events after the Messinian (5.2 Myr), possibly triggered by the rapid ecological changes accompanying the desiccation and refilling of the Mediterranean Basin.

COMPETITION AMONG TADPOLES OF COEXISTING HEMICLONES OF HYBRIDOGENETIC RANA ESCULENTA: SUPPORT FOR THE FROZEN NICHE VARIATION MODEL.

Vertebrate animals reproducing without genetic recombination typically are hybrids, which have large ranges, are locally abundant, and live in disturbed or harsh habitats. This holds for the hemiclonal hybridogenetic frog Rana esculenta: it is widespread in Europe and commonly is found in disturbed habitats such as gravel pits. We hypothesize that its widespread occurrence may either be the result of natural selection for a single hemiclone acting as a broadly adapted "general-purpose" genotype, or of interclonal selection, which maintains multiple hemiclones that each are relatively narrowly adapted and perform differently across environments, that is, the Frozen Niche Variation model. We tested these competing hypotheses using 1000-L outdoor artificial ponds to rear tadpoles of the parental species (Rana lessonae [LL] and Rana ridibunda [RR]) alone, and each of three hemiclones of Rana esculenta (GUT1, GUT2, GUT3) alone, and in mixed hemiclonal populations from hatching to metamorphosis. Tadpoles of three coexisting hemiclones from a single natural population (near Gütighausen, Switzerland) were reared in both two- and three-way mixtures in equal total numbers at high and low density. For each species and hemiclone, the proportion of tadpoles metamorphosing decreased as the density of tadpoles increased, with the three hemiclones spanning the range of values exhibited by the two parental species. LL and GUT1 tadpoles produced the highest proportion of metamorphs, whereas tadpoles of RR produced the fewest metamorphs at both densities. GUT1 tadpoles also produced the largest metamorphs at low density, GUT2 and GUT3 tadpoles produced smaller metamorphs than did GUT1 tadpoles at the low density, but the three hemiclones did not differ from each other at high density. The parental species (LL and RR) were intermediate in metamorphic size to the hemiclones at low density, but all genotypes converged on a similar size at high density. Length of the larval period also was affected by density, but its effect was dependent on genotype. GUT1 tadpoles had the shortest larval period at the low density, but larval period was longer and not different between GUT1, GUT3, and LL at high density. RR tadpoles had the longest larval period at both densities. The most dramatic results were that three genotypes (GUT1, GUT2, and RR) maintained rank order and increased days to metamorphosis from low to high density, whereas two genotypes (GUT3 and LL) changed rank order and decreased days to metamorphosis from low to high density. Mixtures of hemiclones in two- and three-way combinations facilitated the proportion of tadpoles metamorphosing for GUT1 and GUT2 at both densities, but only at the low density for GUT3 tadpoles. Results from this experiment are incompatible with the General-Purpose Genotype model as a global explanation of hybrid abundance in these frogs. Alternatively, the Frozen Niche Variation prediction of general performance superiority of clonal mixtures relative to single clone populations is strongly supported. The data confirm that fitness advantages of hemiclones change, depending on the environment, such that in temporally and spatially heterogeneous habitats like ponds, frequency-dependent selection among hemiclones may promote coexistence in hemiclonal assemblages. Yet, differential dispersal or colonization ability and historical factors affecting hemiclone distribution may also be important in shaping patterns of clonal coexistence.