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.

Hagfish intestinal antimicrobial peptides are ancient cathelicidins.

Three potent broad-spectrum antimicrobial peptides (HFIAP-1, -2, and -3) isolated from intestinal tissues of Myxine glutinosa (Atlantic hagfish) are identified as ancient members of the cathelicidin family of antimicrobial peptides, hitherto known only from mammals. In situ hybridization reveals that HFIAPs are produced in nests of myeloid cells within the loose connective tissue of the gut wall, a tissue reminiscent of both gut-associated lymphoid tissue (GALT) and vertebrate spleen. We suggest that this tissue organization provides local defense of the hagfish gastrointestinal tract via innate immunity and possibly served as the architectural plan upon which the adaptive immune system evolved.

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.

Comparative analysis of mitochondrial genomes in Bombina (Anura; Bombinatoridae).

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.

Phylogeography of the fire-bellied toads Bombina: independent Pleistocene histories inferred from mitochondrial genomes.

The fire-bellied toads Bombina bombina and Bombina variegata, interbreed in a long, narrow zone maintained by a balance between selection and dispersal. Hybridization takes place between local, genetically differentiated groups. To quantify divergence between these groups and reconstruct their history and demography, we analysed nucleotide variation at the mitochondrial cytochrome b gene (1096 bp) in 364 individuals from 156 sites representing the entire range of both species. Three distinct clades with high sequence divergence (K2P = 8-11%) were distinguished. One clade grouped B. bombina haplotypes; the two other clades grouped B. variegata haplotypes. One B. variegata clade included only Carpathian individuals; the other represented B. variegata from the southwestern parts of its distribution: Southern and Western Europe (Balkano-Western lineage), Apennines, and the Rhodope Mountains. Differentiation between the Carpathian and Balkano-Western lineages, K2P approximately 8%, approached interspecific divergence. Deep divergence among European Bombina lineages suggests their preglacial origin, and implies long and largely independent evolutionary histories of the species. Multiple glacial refugia were identified in the lowlands adjoining the Black Sea, in the Carpathians, in the Balkans, and in the Apennines. The results of the nested clade and demographic analyses suggest drastic reductions of population sizes during the last glacial period, and significant demographic growth related to postglacial colonization. Inferred history, supported by fossil evidence, demonstrates that Bombina ranges underwent repeated contractions and expansions. Geographical concordance between morphology, allozymes, and mtDNA shows that previous episodes of interspecific hybridization have left no detectable mtDNA introgression. Either the admixed populations went extinct, or selection against hybrids hindered mtDNA gene flow in ancient hybrid zones.

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.

GYNOGENETIC REPRODUCTION IN HYBRID MOLE SALAMANDERS (GENUS AMBYSTOMA).

Ambystoma platineum, a unisexual clonal triploid taxon of mole salamander, originated by hybridization between the Mendelian species A. jeffersonianum and A. laterale. Studies of lampbrush chromosomes indicated that A. platineum reproduces gynogenetically, that is, sperm from a sexual host species is required to activate egg development but makes no genetic contribution to the developing embryo. Nevertheless, electrophoretic diversity in populations of some hybrid Ambystoma suggested continual in situ recreation of unisexual hybrids and bidirectional gene exchange between the parental species and the hybrids. A. platineum usually lives with, and is sexually dependent on, one of its parental species, A. jeffersonianum. In central Indiana, however, A. platineum populations have shifted their host dependency to A. texanum. Such A. texanum-dependent populations of A. platineum provide an almost ideal system for studying reproductive mode in A. platineum, because both replacement of a jeffersonianum or laterale genome of A. platineum by a texanum genome, and movement of genes from A. platineum to the host species, A. texanum, would be readily detected by electrophoretic markers. Our samples of A. texanum provided no evidence for the transfer of jeffersonianum or laterale genes into A. texanum. Similarly, among 32 A. platineum sampled from six localities in east-central Illinois and central Indiana, we find no texanum alleles, and thus no evidence for genome replacement. The one diploid hybrid individual contained only a jeffersonianum and a laterale genome; because of the absence of either parental species from these populations, this hybrid could only have come from a diploid ovum produced by A. platineum. Both morphometric and electrophoretic results for the two tetraploid individuals indicate that they resulted from fertilization of triploid oocytes of A. platineum by sperm of A. texanum. Because genome replacement in A. texanum-dependent populations of A. platineum is irreversible, the persistence of A. platineum in A. texanum-dependent populations demonstrates conclusively that the major mode of reproduction in A. platineum populations is clonal: A. platineum produces mainly triploid eggs that develop gynogenetically.