Post-zygotic selection against parental genotypes during larval development maintains all-hybrid populations of the frog Pelophylax esculentus
© Reyer et al. 2015
Received: 6 March 2015
Accepted: 2 June 2015
Published: 4 July 2015
Hybridization between two species usually leads to inviable or infertile offspring, due to endogenous or exogenous selection pressures. Yet, hybrid taxa are found in several plant and animal genera, and some of these hybrid taxa are ecologically and evolutionarily very successful. One example of such a successful hybrid is the water frog, Pelophylax esculentus which originated from matings between the two species P. ridibundus (genotype RR) and P. lessonae (LL). At the northern border of the distribution all-hybrid populations consisting of diploid (LR) and one or two triploid (LLR, LRR) frog types have been established. Here, the hybrid has achieved reproductive independence from its sexual ancestors and forms a self-sustaining evolutionary unit. Based on the gamete production of these hybrids, certain mating combinations should lead to LL and RR offspring, but these parental forms are absent among the adults.
In order to investigate the mechanisms that maintain such an all-hybrid system, we performed a field study and a crossing experiment. In the field we sampled several ponds for water frog larvae at different developmental stages. Genotype compositions were then analysed and life-history differences between the genotypes examined. In the experiment we crossed diploid and triploid males and females from different ponds and determined fertilization success as well as development speed and survival rates of the offspring under high, medium and low food availability. In both parts of the study, we found numerous LL and RR offspring during the egg and early larval stages; but the frequency of these parental genotypes decreased drastically during later stages. In natural ponds almost all of them had disappeared already before metamorphosis; under the more benign experimental conditions the last ones died as juveniles during the following year.
From the combined results we conclude that the absence of parental genotypes in all-hybrid populations is due to post-zygotic selection against them, rather than to pre-zygotic mechanisms that might prevent their formation in the first place. For this post-zygotic selection, genetic mechanisms resulting from low genetic diversity and fixation of deleterious mutations seem to be a more likely explanation than ecological factors.
Hybridization, the production of offspring between individuals of different species, is quite common in nature (reviewed by [1, 2]), and the number of recognized taxa that are of hybrid origin is increasing, mostly due to the improvement of molecular techniques . Yet, the evolutionary importance of hybridization has been strongly debated among biologists for many decades. Whereas botanists have long accepted hybridization as an important evolutionary force in speciation , zoologists traditionally neglected such an influence and considered it an evolutionary dead end . Animal hybrids are usually thought to have lower fitness than the parental species due to endogenous (genetically based) or exogenous (environmentally based) selection acting upon them. Endogenous selection is based on the concurrence of two different genomes, which often results in developmental instability, sterility or even death of the newly formed hybrid [6, 7]. Exogenous selection against hybrids can result from the fact that their morphological, physiological and/or behavioural traits are often intermediate between those of their two parental species [8–11]. This will leave them at a disadvantage in the parental habitats; successful establishment may only be possible where intermediate habitat conditions exist .
However, hybrid fitness in relation to fitness of the parental forms can vary greatly, and cases of hybrid superiority to at least one parental form have been well documented (reviewed in ). Such hybrid superiority is often caused by heterosis effects and therefore restricted to (or at least highest in) the F1 generation . As a result, the hybrid taxon will only persist under conditions of repeated primary hybridization. To be of evolutionary importance the hybrid either has to backcross with at least one of the parental species and introduce new genes into the parental gene pool; or the hybrid has to become an evolutionary significant unit (ESU) of its own by either outcompeting the parental species and overtaking their niches or by developing adaptations to a new specific niche and, as a result, restrict gene flow from other lineages within the higher organizational level of the species . Thanks to increased research in this field and improved genetic tools, meanwhile several hybrid taxa have been identified that are ecologically very successful and seem to have reached an evolutionarily old age [16–18].
The palaearctic water frog complex (Pelophylax)
One of the systems where the hybrid taxon has existed over a long time period, i.e., 10,000-50,000 years , and is ecologically successful (in the sense that it is geographically widespread) is the European water frog Pelophylax esculentus, which was formed through repeated hybridization between two parental species P. lessonae (genotype LL) and P. ridibundus (genotype RR). The hybrid nature of P. esculentus was first shown by Berger [20, 21] through biometric analyses and breeding experiments; further investigations revealed that its usual reproductive mode is hybridogenetic . Hybridogenesis is characterized by premeiotic exclusion of one parental genome during gametogenesis, followed by endomitosis of the remaining genome (i.e., its duplication prior to cell division) and its clonal transmission to gametes that all carry identical chromosomes [23, 24]. Offspring from matings between hybrids usually do not survive, due to the accumulation of mutations on the clonally inherited genomes, which can impair development through homozigosity of recessive deleterious alleles at particular loci and/or general deterioration from high mutational load, independent of homozigosity [19, 25, 26]. Therefore, production of viable hybrid offspring requires that P. esculentus acts as a sexual parasite that restores hybridity by mating with the respective parental species (the sexual host) whose genome was excluded. Thus, P. esculentus represent ‘hemiclones’  in which selection against deleterious mutations is weak because the clonal genome transmitted by the parasite is constantly sheltered by the sexual genome of the host .
Whether or not genome exclusion is induced and, if so, which parental genome is excluded and to what extent may vary with population composition, ploidy and geographic area [28, 29]. Where the L-genome is eliminated and the R-genome transmitted, P. esculentus occur in sympatry and mate with the parental species P. lessonae (LE-system); where the reverse pattern exists (elimination of R, clonal transmission of L), hybrids co-occur and mate with P. ridibundus (RE-system) . In both these systems, the vast majority of the hybrids are diploid (genotype LR); and if triploids are present at all, they occur in proportions of less than 15 % (for details see Appendix 1 ). Models have shown that stability of these mixed hybrid-parental systems is very sensitive to several factors, such as mating preference, female fecundity and larval performance of the involved taxa [32–34].
Offspring with parental genotypes are not produced from the very beginning, because of assortative mating and/or unsuccessful fertilization of eggs (pre-zygotic selection).
All genotypes are initially produced, but those of the parental species do not survive to the adult stage (post-zygotic selection).
To distinguish between these two hypotheses we combined results from a field study, where we sampled 12 ponds for water frog offspring at different developmental stages, with results from an experiment, where we raised tadpoles emerging from artificial crosses between adult frogs of known genotypes. In both parts of the study, we analyzed the genotype composition at different developmental stages to test whether the relative frequencies of the genotypes changed from eggs through tadpoles and metamorphs to 1-year old juveniles. In the experimental part, we further identified the produced gametes by comparing the genotypes of parents and offspring and recorded genotype- and sex-specific fertilization rates. We also raised the tadpoles under three feeding regimes with high, medium and low food availability, respectively, and tested for a genotype x environment effect to see whether genotype-specific development varies with environmental conditions.
Ours is not the first study to document changes in larval genotype proportions in all-hybrid populations (e.g., ), but it is the first one to systematically monitor such changes from the egg to the adult stage under natural and experimental conditions for all five genotypes that theoretically could be produced in such populations (LL, LLR, LR, LRR and RR). Moreover, with the exception of one study that used microsatellite markers , earlier studies on larval genotypes were based on measures that cannot reliably distinguish between genotypes, either because the distributions of the respective parameters overlap as for egg size, morphometric indices and erythrocyte size or the resolution is not fine enough as in allozyme analyses [44–46]. By combining microsatellite and flow cytometry analysis we were able to overcome these limitations and achieve unambiguous genotype identification (even of rare forms) at various developmental stages.
Survival and development in natural ponds
Results from general linear models relating genotype proportions to pond of origin and four developmental stages (eggs, tadpoles, metamorphs and juveniles)
Genotype proportions did not differ among ponds, except for LR (Table 1). This difference was mainly due to one pond 102 which had very few LR individuals in all samples. LR proportions in this pond were significantly lower than those in three other ponds (111, 134 and 138; all P ≤ 0.023), whereas all other ponds did not differ in pairwise comparisons. However, sample size in pond 102 was very low, so that the difference may be due to stochasticity.
Results from general linear models for a) tadpole, b) metamorph and c) juvenile development in natural ponds
Developmental stage 
Later in the development, there were hardly any genotypes present, except the three hybrid types (cf. Fig. 2). Among metamorphs, we found overall significant differences regarding snout-vent length and weight (Table 2b). Pairwise comparisons showed that LRR metamorphs were significantly heavier than LR individuals (P = 0.002) and larger than both LR and LLR (both P ≤ 0.044), whereas other pairwise differences were not significant (all P ≥ 0.275). Among juveniles, the genotype effect was not significant (Table 2c), but at both stages size and body mass of individuals differed among ponds.
Fertilization success and hatching rate
Results from a general linear model for fertilization success and hatching rate in relation to female and male genotype, their interaction and the effect of crossing type, i.e., within versus between pond crossings
Source of variation
Genotype F (pond F)
Genotype M (pond M)
Genotype F x genotype M
Survival and development of tadpoles and metamorphs
Results from general linear models testing for differences in a) tadpole survival and development (according to ) until metamorphosis and b) froglet survival and body mass until July of the following year
Offspring Type x Treatment
Our results from the 12 natural ponds show that initially all offspring genotypes that are possible (based on the existing gamete types) are actually produced. At the fertilized egg stage about 50 % of the sample was comprised of genotypes that existed among adults only very rarely (LLRR hybrids and uncertain genotypes) or not at all, such as the parental forms (LL, LLL, RR and RRR) (cf. Fig. 1); but these genotypes gradually disappeared during subsequent stages. From the metamorph stage on, the genotype composition in all ponds was that of typical EE-systems, i.e., with almost exclusively LR, LLR and LRR hybrids among the adults in ratios that are very similar to those of juveniles and adults (Fig. 2). Therefore, we conclude that the existence and maintenance of all-hybrid populations is due to post-zygotic selection against parental forms (hypothesis 2) during the tadpole stage, rather than to pre-zygotic selection like assortative mating and/or unsuccessful fertilization (hypothesis 1).
Pre-zygotic selection: assortative mating
The absence of assortative mating in all-hybrid populations of the EE-system is supported by two other studies. In two all-hybrid populations the relative frequencies of the various male x female mating combinations did not differ from those calculated from the genotype proportions under the assumption of random mating . Also, mate choice experiments revealed no preference of diploid and triploid females for males of any of the three hybrid types . These results contrast with those from similar experiments in LE populations: when given a choice between LL and LR males or their calls, LR and LL females both preferred LL males [49–52]. The difference in female choice between the two systems can be explained in both proximate and ultimate ways. In terms of proximate mechanisms, the three hybrid types (LR, LLR, LRR) are more similar to each other in size, coloration and male advertisement calls than diploid hybrids (LR) and the parental species (LL) are, thus making discrimination and choice more difficult [53, 54]. In terms of ultimate reasons, genetic fixation of mate preferences in a particular genotype should be impossible in diploid-triploid all-hybrid populations, because suitable mating partners (i.e., the ones guaranteeing viable progeny types) alternate each generation: diploid LR females producing diploid eggs should choose triploid LLR or LRR males; the resulting triploid daughters should choose diploid LR males etc. . In contrast, in LE-systems, mating with LL males is strongly selected for in LR females, because it is always the only way to produce viable offspring [32, 33].
Pre-zygotic selection: gamete types and fertilization success
Pre-zygotic selection against certain gamete types and their fusion as a reason for the lack of P. lessonae and P. ridibundus among adults of all-hybrid populations can also be dismissed. Although, for technical reasons, we had to deduce the produced gamete types from genotyping tadpoles that were 18 days old (see Methods) and, thus, cannot rule out that genotype proportions at this stage differed somewhat from those among the originally produced gametes and zygotes, this potential bias cannot explain the observed population composition. Haploid L and R eggs and sperm, as well as offspring with parental genotypes were produced in large numbers, both in natural ponds (Figs. 2 and 3) and in the crossing experiment (Fig. 5) and with no effect of crossing type; i.e., fertilization success and hatching rate were independent of whether the crossed males and females came from the same or from different ponds. They were also independent of whether diploid LR-females exclusively produced the typical diploid LR- or in addition varying amounts of haploid R-eggs (Table 3). The simultaneous production of LR- and R-eggs by the same diploid LR-female with apparently little fertility impairment has also been described in a number of other P. esculentus studies (reviewed by [35, 36]).
This differs from the situation in diploid males where average fertilization and hatching rates were lower than for all other parent types, with values decreasing from males in pond 011 that produced only R sperm through those from pond 089 producing predominantly R- and a small proportion of LR sperm to males from pond 001 with all three gamete types (L, R and LR) in fairly high proportions (Fig. 4). These results strongly suggest that the poor reproductive success in these hybrid males is due to disturbed gametogenesis, as has been already suggested previously [56–59]. In extreme cases, this can lead to the development of abnormal gonads and germ cells  and complete male sterility [35, 61].
The link between reproductive success and the extent of gametogenesis problems is further supported by our finding that the high rates of fertilization and hatching in triploid hybrids of both types (LLR and LRR) are paralleled by no aberrant gamete forms in males and only very few in females (Fig. 3). The reason why gametogenesis seems to be more consistent and less disturbed in triploid than in diploid frogs lies in the simpler cytogenetic mechanism of hemiclonal inheritance through so-called “meiotic hybridogenesis” [62, 63]: after the elimination of the unmatched chromosome (“homogenizing elimination”), in triploids “no compensatory duplication of the remaining genetic material is necessary, as it is in diploids” , and the two homologous chromosomes can go through a normal meiosis.
Although gamete variability of triploids is low among individuals from the same pond, it can differ between populations from different geographical areas, even to the extent that in some populations LLR males typically produce diploid LL instead of haploid L sperm [41, 59, 65, 66]. This and the fact that varying proportions of diploid eggs and sperm have also been found in other studies (e.g., [56, 67]) probably reflects different abilities of the R-genome, the L-genome or both to induce and/or resist genome exclusion (e.g., [68–70]). These differences may result from different primary hybridizations, ecological adaptations and/or reflect variable outcomes of a genomic conflict between segregation distorters (or meiotic driver genes) in the excluding genome and modifier genes in the host genome that counteract meiotic driver genes [71, 72]. Production of diploid gametes may also represent an adaptation of the host genome to prevent exclusion : if, for instance, in an LE-system P. lessonae males produced LL rather than L sperm that fertilise R eggs of diploid P. esculentus, then genome exclusion would (probably) reverse from L in LR offspring to R in LLR offspring. Although we know of no such LL sperm production in P. lessonae, it just may have escaped our attention because of its rareness. This, however, does not exclude the possibility that it can be successful. Production of LL sperm is the usual pattern in LLR males from some populations where it most likely arose from a single event of L genome doubling .
Production of polyploid gametes that upon fusion with haploid ones will lead to unusual offspring types such as the very rare LLL, RRR and the about 7 % uncertain genotypes found in this study or the pentaploid LLLRR reported by .
Progeny with deviations from the parental chromosome proportions expected for diploid (13 + 13) and triploids (13 + 26) (cf. ).
It seems plausible to assume that such aberrant offspring encounter developmental problems and are not viable or fertile. But why do the normal parental genotypes (LL and RR) gradually disappear towards sexual maturity? Why is there post-zygotic selection against the parental species in EE-systems, rather than against hybrids as in many other systems? Below, we discuss two potential explanations: environmental factors and genetic mechanisms.
Post-zygotic selection: environmental factors
Survival and development of Pelophylax tadpoles is definitely affected by environmental conditions. This is indicated by various lines of evidence. First, survival rate and developmental stage of tadpoles and the size and body mass of metamophs differed among ponds (Table 2). Second, development improved with improving environmental conditions: it was worst in natural ponds where LL and RR offspring did not even survive to the metamorph stage, whereas under the more benign conditions of the experiments development improved with increasing food availability, and several parental offspring types reached the froglet stage (also found by ) and a few even the stage of 1-year old juveniles. Besides food, water temperature has been identified as a factor contributing to developmental differences. Under cold conditions, as they exist at the northern range of the water frog distribution, hybrids not only survive better, they also metamorph faster than the parental species [78, 79] which can have a beneficial influence on survival at later life ages .
Based on our study, however, neither food nor temperature can explain the results. The parental genotypes disappeared from all investigated natural ponds, independent of their specific combination of biotic and abiotic conditions (for details see ). And although tadpoles in the experiments were raised under different water temperatures (colder in Sweden than in Zürich) and at three different food levels, parental genotypes performed worse than hybrid types under all conditions of our study and a previous one by . This lack of a genotype x treatment interaction (cf. Table 4) indicates that at least the environmental conditions tested here cannot be made responsible for the differential offspring survival that leads to the gradual disappearance of LL and RR animals from the population.
Post-zygotic selection: genetic mechanisms
In terms of genetic post-zygotic selection mechanisms, inviability, reduced fertility and other developmental abnormalities in offspring resulting from hybrid-hybrid matings are generally attributed to an accumulation of deleterious mutations on clonally transmitted genomes through Muller’s ratchet . Such mutational load can act in two not mutually exclusive ways: 1) through homozigosity of recessive deleterious alleles at particular loci or 2) by the cumulative load, and hence general deterioration, of the clonal genomes, independent of homozigosity [26, 36]. Crossing experiments that explicitly addressed these two explanations provided more support for the first mechanism [19, 25, 26]. At a first glance neither mechanism seems to offer an explanation for the gradual disappearance of LL and RR offspring from hybrid x hybrid matings, because reproduction is basically sexual: both the L and the R genome are regularly recombined in triploid individuals, the L when in LLR and the R when in LRR [39, 41, 83]. Therefore, deleterious mutations can regularly be purged and should not accumulate and get fixed. Nevertheless, homozygosity for deleterious mutations may contribute to the mortality of LL and RR progeny, since fixation and low genetic diversity does exist in our Scandinavian all-hybrid populations, at least at microsatellite loci [81, 83–85]. At present, we do not know what mechanism is responsible for that low genetic diversity. Maybe it is due to repeated population bottlenecks and founder effects after glacial periods when water frogs expanded their range from southern refuge areas northwards [86, 87].
The formation of diploid LR-sperms by hybrid males observed in our study has been described earlier. Some researchers have dispatched them as unimportant due to their larger size and supposed lower swimming speed and fertilisation success compared to haploid sperm [56, 67, 88]. Christiansen  on the other hand highlighted the importance of diploid LR sperm for the structure of all-hybrid populations after comparing predictions from a theoretical model with empirical data. She found that LR sperm production in LRR-rich populations was sufficiently high (22 %) to explain the observed proportion of LRR males. Their formation previously had not been understood, because the usual genesis of LRR from LR eggs and R sperm yields only females. Also, the fertilization of LR-eggs with LR-sperm results in symmetrical viable tetraploid offspring (LLRR), which were found in both, natural ponds and our crossing experiment; and in nature at least a few of them survive to adults of both sexes . Since symmetrical polyploids do not encounter the meiotic problems that triploids and other asymmetrical polyploids are facing during gametogeneis, LLRR hybrids can, at least theoretically, have an important evolutionary perspective by opening the possibility for hybrid speciation. Speciation via polyploid hybrids is common in plants [12, 89, 90], but examples from the animal kingdom are scarce, especially when it comes to vertebrates [1, 2, 18, 91, 92]. Unfortunately, we do not know yet whether and how regularly male and female tetraploids produce the diploid LR gametes that would be required for perpetuating a tetraploid line. The only LLRR male for which gamete production is known, produced haploid R sperm and a few diploid cells of unknown genotypic composition . But even if LLRR individuals of both sexes regularly produced LR gametes, Christiansen’s  model predicts that a stable tetraploid population would only be established under a more than a twofold increase in survival or reproductive output of both male and female LLRR. This is an unlikely scenario and may explain why, so far, no all-hybrid Pelophylax population with a high proportion of tetraploids has been found, despite extensive sampling all over Europe . Yet, theoretically, conditions could be fulfilled under certain ecological conditions, as exemplified by the Iberian minnow Squalius alburnoides where the proportion of tetraploids varies markedly between populations living in different habitats [18, 93].
It has been shown, however, that triploid forms producing diploid gametes in one sex and haploid ones in the other sex can act as a stepping stone towards tetraploidization (triploid bridge; [93–95]). They even may establish stable bisexually reproducing all-triploid populations, as  have found in Batura toads (Bufo pseudoraddei baturae). Thus, under certain genetic and ecological conditions, hybrids can become evolutionary significant units (ESUs) in the sense that some lineages demonstrate “highly restricted gene flow from other such lineages within the higher organizational level (lineage) of the species” . Such is also illustrated by the fact that triploid individuals in all-hybrid Pelophylax esculentus populations originate from different gamete combinations: in most populations haploid sperm fuses with diploid eggs, whereas in some others diploid sperm fertilizes haploid eggs .
For two reasons, the existence of all-hybrid populations of the hemiclonal frog Pelophylax esculentus has long remained a puzzle. First, it was believed that its hybridogenetic mode of reproduction forces P. esculentus (the sexual parasite) to always live in sympatry and mate with one of its parental species (the sexual host). This puzzle was solved when it was detected that triploid hybrids can serve as sexual hosts for diploid hybrids and vice versa. Second, analysis of gamete production patterns and mating behavior in diploid and triploid P. esculentus revealed, that offspring of both parental species should regularly be produced in all-hybrid populations; yet, they do not exist among adults. In this study, we solved this second puzzle through a field study and a crossing experiment. From the combined results we conclude that the absence of parental genotypes among adults of all-hybrid populations is due to post-zygotic selection against them, rather than to pre-zygotic mechanisms that might prevent their formation in the first place. For this post-zygotic selection, genetic mechanisms seem to be a more likely explanation than ecological factors. At present, we do not know the nature of these mechanisms, but fixation of deleterious mutations and low genetic diversity resulting from repeated population bottlenecks are likely candidates. Given, that all-hybrid populations with diploid and triploid frogs can be perpetuated in the absence of the parental species, P. esculentus can be viewed as an evolutionary significant unit that may be on its way towards hybrid speciation via different possible evolutionary trajectories .
Our analysis of larval and juvenile development under natural conditions is based on repeated sampling in 12 ponds in Skåne (Southern Sweden) at four different stages: eggs, tadpoles, metamorphs and one-year old juveniles. Eggs were collected between June 1 and June 10, 2003; tadpoles were obtained about 6 weeks later (July 16–21) by catching a random sample with a dip net. Metamorphs were caught by hand between August 5 and August 12 and juveniles in spring and summer of the following year. Originally, we had aimed at collecting at least five egg clutches, 25 tadpoles and 20 metamorphs per pond. Unfortunately, these sample sizes could not be achieved in all ponds, partly because of low numbers (perhaps due to high fish abundance) and partly due to problems in finding the existing eggs, tadpoles and metamorphs when the water was muddy and/or vegetation was dense.
For investigating larval development under experimental conditions, we crossed adult frogs from three ponds in Skåne that were chosen based on different population composition, assessed in 2002 and 2003. Pond 001 was LLR-dominated, 089 LRR-dominated and in pond 011 the two triploid genotypes occurred in equal proportions. A genotype was considered dominant when adult frogs of this type constituted more than 50 % of the adult population in both years. While LLR frogs are slightly biased towards males, LRR frogs are heavily sex-biased towards females, making it difficult to encounter LRR males. All frogs were caught between May 15 and May 19, 2004, at night by hand and kept at the field station of the University of Lund at 10 °C prior to the crossing.
Since both parts of the study hinge on correct genotype determination, we combined several techniques to obtain accurate results: microsatellite analysis on tissue, flow cytometry on blood, and morphometric indices. For analyzing the genotypes in the egg stage, we raised a subsample of 15 individuals per clutch at Stensoffa, the field station of the University of Lund, under ad libitum food conditions until July 22. This was necessary because analyzable amounts of blood and tissue can only be collected once the tadpoles have reached a certain size (~50 days old). Both sets of tadpoles (those raised from the eggs and those sampled from ponds) were staged for their development according to  and then killed with a solution of 3-aminobenzoic acid ethyl ester methanesulfonate (MS-222, 5 g/l, Sigma A5040), because it was not possible to obtain enough blood from living animals. Tissue was collected by cutting off part of the tadpole tail. From metamorphs and juveniles we took a toe clip for microsatellite analysis and a blood sample for flow cytometry analysis. Blood was obtained by cutting the web of a hind foot and collecting the emerging drop with a heparinized capillary. All blood samples were stored in a flow rate calibration (FRC) solution , and all tissue samples were kept in 70 % ETOH until lab analysis for the genotype determination was done. Additionally, we measured snout-vent length (SVL) and body mass of the metamorphs and juveniles and in adults also the tibia length (TL) and the callus length (CL).
Genome composition of the to-be crossed adults was initially determined based on morphological indices calculated from SVL, TL and CL (for details see [20, 35, 36]) and flow cytometry and later (i.e., after the crossing) confirmed by microsatellite analysis on tissue samples. For this analysis, we used seven polymorphic microsatellite loci: Ca1b5, Ca5, Ca18 , Res16 , Ca1b6, Re1CAGA10, Ga1a19 . Alleles at these loci are species-specific, i.e., they can unambiguously be assigned to either the L- or the R-genome. Additionally, the loci Ca1b5, Ca1b6, Ga1a19 and Res16 show gene dosage  which - in addition to flow cytometry - provided further information about the exact genotype. Details of the protocol for microsatellite analysis are given by . Flow cytometry on blood samples allowed us to distinguish between LR, LLR, LRR and other ploidy levels (e.g., LLRR, LLL and RRR), because L- and R-genomes have different amounts of DNA ; the protocol we used is described in . Because the triploid parental types LLL and RRR were rare (0.4 % of the whole sample), we included them for the analysis in the “normal” parental genotypes (LL and RR, respectively).
If flow cytometry and all microsatellite loci showed the same result, the individual was clearly assigned to one genotype. In some cases, however, the results were ambiguous or contradictory, even after reanalysis, e.g., when some microsatellite loci showed a LLR genotype pattern but others a LRR, or when flow cytometry and most microsatellite loci indicated an LLR genotype, but one locus showed only LR, so there was one allele missing. Cases showing repeatedly such contradictory results were categorized as “uncertain” genotypes. They may have been aneuploid (i.e., have missing or additional chromosome fractions) as suggested by  or have had the same number of chromosomes as normal diploid and triploid individuals (2n = 26 or 3n = 39), but in an unusual composition of L or R chromosomes, e.g., LR: L = 12 and R = 14, instead of 13 each . Because of these uncertainties such offspring were not considered in most analyses.
On May 23 (1 day after the crossings), the fertilized eggs were transferred from the petri dishes into 1–liter tubs containing aged tap water and then stored in a room at approximately 20 °C. Unfertilized eggs or embryos that stop development can cause degradation of water quality and were therefore removed during regular checking of the tubs. Water was changed twice a week for the same reasons. The larvae were kept indoors until 18 days after fertilization (June 9), when all had reached at least stage 25 which is characterized by the disappearance of external gills . The surviving larvae were weighed, staged and counted, and hatching rate was calculated as the number of hatched larvae relative to the number of fertilized eggs.
On the same day (June 9), we created from each cross four sets of randomly selected tadpoles without any obviously lethal morphological abnormalities and/or feeding and swimming impairments. Set one was used for later genotype determination through microsatellites and flow cytometry as described above. For this analysis, we used 25–56 offspring from each parent individual. The only exception was an LR male from pond 001, for which only two larvae could be analyzed, due to extremely low fertilization success of this male.
Tadpoles of set two were transferred to outdoor tubs at the Stensoffa Fieldstation in Sweden, while those of sets three and four were transported on June 10 by car to the University of Zürich (Switzerland), where on June 11 they were also transferred to outdoor tubs. In both localities, we randomly assigned 15 tadpoles from each cross to tubs of 60 l (Sweden) and 90 l (Zürich), respectively. These tubs had been filled 5–6 weeks earlier with water, provided with a handful of dried leaves and 1–3 snails (Lymnaea sp.), and inoculated with phyto- and zooplankton to create a self-sustaining aquatic community . Tadpoles of set two (Sweden) were fed every other day ad libitum with rabbit chow (high food treatment); those of set three in Zürich were fed twice a week by adding to each tub 0.5 g rabbit chow (medium food), while those of set four in Zürich could only feed on algae growing naturally in the tub (low food). All tubs were arranged on meadows in a random design and covered with lids to prevent colonization by invertebrate predators.
Eleven percent of the crosses produced not enough viable tadpoles to stock every tub with 15 individuals; in these cases we transferred the available number, ranging from 5 to 14 tadpoles. Five crosses could not be included in further analyses because none of the tadpoles had survived. Larval survival rates under the three food treatments were determined as the proportion of individuals surviving until nine (Sweden) and ten (Zürich) weeks, respectively, after transferring tadpoles to the outdoor tubs. (The 1 week shorter period in Sweden was necessary to avoid time conflicts between raising tadpoles and sampling in natural ponds.) Until then they were counted once a week in Zürich and, on average, every 9 days in Sweden. At the end of this period, we also determined the developmental stage of the surviving animals according to . At this point, the high and low food treatments were terminated and tadpoles preserved for subsequent genotype determination.
Tadpoles from the medium food treatment were returned to their respective tubs and pieces of floating wood added to each tub to allow animals to leave the water after metamorphosis. From then on tubs were checked daily and new metamorphs collected. They were transferred to 1 liter plastic indoor boxes (still separated by cross) with a bottom of wet soil, grass and a few pieces of wood for hiding. In their new environment, animals were kept at room temperature and fed with live crickets twice a week until December 15 when the temperature was gradually lowered to 5 °C over a period of 5 days. After hibernation until early April 2005, temperature was raised to 15 °C within 5 days. All froglets that had survived the winter were transferred to cross-specific 5 liter outdoor boxes, with soil, grass and wooden pieces in one part and water in the other. Here they were kept until July 2005, receiving live crickets twice a week. On June 2, 2005, survival rate and average body mass was determined for each cross, and survival was determined again in July. Thereafter, the experiment was ended by euthanizing the remaining individuals, because due to their origin from artificial crossing between frogs from a foreign country they could not be released into nature.
For our samples from natural populations we calculated the pond-specific proportions of all genotypes for each of the four offspring stages (eggs, tadpoles, metamorphs and juveniles) and tested for differences between them by means of general linear models (PROC GLM), with offspring stage and pond as independent factors. Differences in developmental rates of tadpoles (expressed by Gosner stage, ) and in snout-vent length and body mass of metamorphs and juveniles were also analyzed with general linear models.
In our crossing experiment, the produced gamete types were identified by comparing the genotypes of the tadpoles with those of the males and females that were crossed. If an individual produced more than one gamete type, the proportions were calculated as the number of offspring from each gamete type relative to the total number of offspring from this individual.
The early developmental variables (fertilization success and hatching rate) were analyzed with general linear models to test for the effects of female genotype, male genotype, their interaction and crossing type. Genotypes were nested within pond to account for different origins of the same genotypes. Because of missing genotypes, the experiment was not completely balanced and we therefore used Type III sums of squares. The category “crossing type” indicates whether the parents that were mated originated from the same pond (“within ponds”) or from different ponds (“between ponds”).
General linear models were also used to relate larval survival and developmental stage in outdoor tubs, as well as metamorph survival and body mass, to offspring type, food treatment, their interaction and crossing type. For these analyses we used means per cross, including only those crossings where all offspring had the same genotype. For all GLMs we applied post-hoc pairwise Bonferroni tests to investigate which of the genotypes differed. In case of small sample sizes (e.g., few surviving LLRR or LL tadpoles) we pooled offspring types into two categories: hybrids and parental animals. Proportions like genotype ratios in natural ponds and fertilization success, hatchling rate and tadpole and metamorph survival in the experiment were arcsine-square root transformed prior to analyses . All statistical analyses were performed with SYSTAT Version 11 (Systat Software, Inc. 2008).
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad repository: https://datadryad.org/resource/doi:10.5061/dryad.6rg87.
The authors thank S. Röthlisberger for help with microsatellite analyses. The group of R. Zinkernagel at the University Hospital Zürich kindly provided us access to their FACSCalibur. B. Baldetorp and M. Johansson of the Lund University Hospital as well as E. Niederer of IBT at the ETH Zürich provided advice and help for flow cytometry analysis. In counting and staging hundreds of tadpoles we were assisted by C. Haeberli, E. Embrechts and B. Schnüriger. This work was supported by SNSF grant No. 31–64004.00 to H.-U. Reyer. Catching and crossing of frogs was granted by the Swedish authority Länsstyrelsen i Skåne Län through permits no. 552-18591-02, 522-9286-03 and 552-6571-04.
- Arnold ML. Natural hybridization and evolution. Oxford: Oxford University Press; 1997Google Scholar
- Mallet J. Hybrid speciation. Nature. 2007;446:279–83.PubMedView ArticleGoogle Scholar
- Seehausen O. Hybridization and adaptive radiation. Trends Ecol Evol. 2004;19(4):198–207.PubMedView ArticleGoogle Scholar
- Grant V. Plant speciation. 2nd ed. New York: Columbia University Press; 1981.Google Scholar
- Mayr E. Animal species and evolution. Cambridge: Belknap Press of Harvard University Press; 1963.View ArticleGoogle Scholar
- Bronson CL, Grubb TC, Braun MJ. A test of the endogenous and exogenous selection hypotheses for the maintenance of a narrow avian hybrid zone. Evolution. 2003;57(3):630–7.PubMedView ArticleGoogle Scholar
- Peterson MA, Monsen KJ, Pedersen H, McFarland T, Bearden J. Direct and indirect analysis of the fitness of Chrysochus (Coleoptera : Chrysomelidae) hybrids. Biol J Linnean Soc. 2005;84(2):273–86.View ArticleGoogle Scholar
- Grant PR, Grant BR. Hybridization of bird species. Science. 1992;256:193–7.PubMedView ArticleGoogle Scholar
- Mavarez J, Salazar CA, Bermingham E, Salcedo C, Jiggins CD, Linares M. Speciation by hybridization in Heliconius butterflies. Nature. 2006;441(7095):868–71.PubMedView ArticleGoogle Scholar
- Mendelson TC, Shaw KL. Close-range acoustic signaling and mate choice in Hawaiian crickets (Gryllidae : Laupala). Behav Ecol Sociobiol. 2006;59(6):770–6.View ArticleGoogle Scholar
- Meyer A, Salzburger W, Schartl M. Hybrid origin of a swordtail species (Teleostei : Xiphophorus clemenciae) driven by sexual selection. Mol Ecol. 2006;15(3):721–30.PubMedView ArticleGoogle Scholar
- Anderson E, Stebbins GL. Hybridzation as an evolutionary stimulus. Evolution. 1954;8(4):378–88.View ArticleGoogle Scholar
- Arnold ML, Hodges SA. Are natural hybrids fit or unfit relative to their parents? Trends Ecol Evol. 1995;10(2):67–71. doi:10.1016/S0169-5347(00)88979-X.PubMedView ArticleGoogle Scholar
- Falconer DS, Mackay TFC. Introduction to Quantitative Genetics. 4th ed. Harlow, England: Pearson Education Ltd.; 1996.Google Scholar
- Fraser DJ, Bernatchez L. Adaptive evolutionary conservation: towards a unified concept for defining conservation units. Mol Ecol. 2001;10(12):2741–52.PubMedView ArticleGoogle Scholar
- Hedges SB, Bogart JP, Maxson LR. Ancestry of unisexual salamanders. Nature. 1992;356:708–10.PubMedView ArticleGoogle Scholar
- Quattro JM, Avise JC, Vrijenhoek RC. Mode of origin and sources of genotypic diversity in triploid gynogenetic fish clones (Poeciliopsis, Poeciliidae). Genetics. 1992;130:621–8.PubMed CentralPubMedGoogle Scholar
- Alves MJ, Coelho MM, Collares-Pereira MJ. Evolution in action through hybridisation and polyploidy in an Iberian freshwater fish: a genetic review. Genetica. 2001;111(1–3):375–85.PubMedView ArticleGoogle Scholar
- Guex GD, Hotz H, Semlitsch RD. Deleterious alleles and differential viability in progeny of natural hemiclonal frogs. Evolution. 2002;56(5):1036–44.PubMedView ArticleGoogle Scholar
- Berger L. Some characteristics of the crosses within Rana esculenta complex in post larval development. Ann Zool Warszawa. 1970;27:373–416.Google Scholar
- Berger L. Embryonla and larval development of F-1 generation of green frogs: Different combinations of Rana lessonae, Rana esculenta and Rana ridibunda characters. Acta Zool Cracov. 1967;12(7):123–60.Google Scholar
- Tunner HG. Das Albumin und andere Bluteiweisse bei Rana ridibunda Pallas, Rana lessonae Camerano, Rana esculenta Linné und deren Hybriden. Z Zool Syst Evol. 1973;11:219–33.View ArticleGoogle Scholar
- Schultz RJ. Hybridization, unisexuality and polyploidy in the teleost Poeciliopsis (Poeciliidae) and other vertebrates. Am Nat. 1969;103:605–19.View ArticleGoogle Scholar
- Dawley RM. An introduction to unisexual vertebrates. In: Dawley RM, Bogart JP, editors. Evolution and ecology of unisexual vertebrates. Albany New York, USA: New York State Museum; 1989. p. 1–18.Google Scholar
- Vorburger C. Fixation of deleterious mutations in clonal lineages: evidence from hybridogenetic frogs. Evolution. 2001;55(11):2319–32.PubMedView ArticleGoogle Scholar
- Vorburger C, Schmeller DS, Hotz H, et al. Maske Damage: Mutational Load in Hemiclonal Water Frogs. In: Schön I, Martens K, van Dijk P, editors. Lost Sex – The Evolutionary Biology of parthenogenesis. Heidelberg: Springer; 2009. p. 433–46.View ArticleGoogle Scholar
- Vrijenhoek RC, Angus RA, Schultz RJ. Variation and Heterozygosity in Sexually Vs Clonally Reproducing Populations of Poeciliopsis. Evolution. 1977;31(4):767–81.View ArticleGoogle Scholar
- Berger L. Systematic and genetic system in water frogs of Europe. Przegl Zool. 1983;27:47–62.Google Scholar
- Vinogradov AE, Borkin LJ, Gunther R, Rosanov JM. Two germ cell lineages with genomes of different species in one and the same animal. Hereditas. 1991;114:245–51.PubMedView ArticleGoogle Scholar
- Uzzell T, Berger L. Electrophoretic phenotypes of Rana ridibunda, Rana lessonae, and their hybridogenetic associate Rana esculenta. Proc Acad Nat Sci Philadelphia. 1975;127:13–24.Google Scholar
- Hoffmann A, Plötner J, Pruvost NBM et al. Genetic diversity and distribution of diploid and polyploid hybrid water frog populations (Pelophylax esculentus complex) across Europe. BMC Mol Ecol. 2015 (in press).Google Scholar
- Som C, Anholt BR, Reyer HU. The effect of assortative mating on the coexistence of a hybridogenetic waterfrog and its sexual host. Am Nat. 2000;156(1):34–46.PubMedView ArticleGoogle Scholar
- Hellriegel B, Reyer HU. Factors influencing the composition of mixed populations of a hemiclonal hybrid and its sexual host. J Evol Biol. 2000;13(6):906–18.View ArticleGoogle Scholar
- Reyer HU, Walti MO, Battig I, Altwegg R, Hellriegel B. Low proportions of reproducing hemiclonal females increase the stability of a sexual parasite-host system (Rana esculenta, R-lessonae). J Anim Ecol. 2004;73(6):1089–101.View ArticleGoogle Scholar
- Die GR, Europas W. Neue Brehm-Bücherei. A. Ziemsen Verlag: Wittenberg Lutherstadt; 1990.Google Scholar
- Plötner J. Die westpaläarktischen Wasserfrösche. Laurenti-Verlag: Bielefeld; 2005.Google Scholar
- Jakob C, Arioli M, Reyer H-U. Ploidy composition in all-hybrid frog populations in relation to ecological conditions. Evol Ecol Res. 2010;12(5):633–52.Google Scholar
- Uzzell T, Berger L, Günther R. Diploid and triploid progeny from a diploid female of Rana esculenta (Amphibia Salientia). Proceedings of the Academy of Natural Sciences of Philadelphia. 1975;127:81–91.Google Scholar
- Günther R, Uzzell T, Berger L. Inheritance patterns in triploid Rana “esculenta” (Amphibia, Salientia). Mitteilungen des Zoologischen Museums Berlin. 1979;55:35–57.Google Scholar
- Christiansen DG. Gamete types, sex determination and stable equilibria of all-hybrid populations of diploid and triploid edible frogs (Pelophylax esculentus). BMC Evol Biol. 2009;9(135):1–16.Google Scholar
- Pruvost NBM, Hoffmann A, Reyer H-U. Gamete production patterns, ploidy, and population genetics reveal evolutionary significant units in hybrid water frogs (Pelophylax esculentus). Ecol Evol. 2013;3(9):2933–46.PubMed CentralPubMedView ArticleGoogle Scholar
- Eikhorst R. Die Verteilung von diploiden und triploiden Larven des Teichfrosches Rana 59esculenta Linnaeus, 1758 in einer reinen Bastardpopulation (Anura: Ranidae). Salamandra. 1988;24:59–68.Google Scholar
- Christiansen DG. A microsatellite-based method for genotyping diploid and triploid water frogs of the Rana esculenta hybrid complex. Mol Ecol Notes. 2005;5(1):190–3.View ArticleGoogle Scholar
- Der ER. Laich des Teichfrosches Rana esculenta Linnaeus, 1758 in einer reinen Bastardpopulation (Anura: Ranidae). Salamandra. 1987;23:122–31.Google Scholar
- Pagano A, Joly P. Limits of the morphometric method for field identification of water frogs. Alytes (Paris). 1999;16(3–4):130–8.Google Scholar
- Schmeller D, Crivelli A, Veith M. Is triploidy indisputably determinable in hybridogenetic hybrids by planimetric analyses of eryhtrocytes? Mitt zool Mus Berlin. 2001;77(1):71–7.Google Scholar
- Günther R, Plötner J. Mating pattern in pure hybrid populations of water frogs, Rana kl. esculenta (Anura, Ranidae). Alytes. 1989;8:90–8.Google Scholar
- Rondinelli B. Female choice in all-hybrid populations of Rana esculenta [Master thesis]: Zurich: University of Zurich; 2006.Google Scholar
- Abt G, Reyer H-U. Mate choice and fitness in a hybrid frog: Rana esculenta females prefer Rana lessonae males over their own. Behav Ecol Sociobiol. 1993;32:221–8.View ArticleGoogle Scholar
- Roesli M, Reyer HU. Male vocalization and female choice in the hybridogenetic Rana lessonae/Rana esculenta complex. Anim Behav. 2000;60:745–55.PubMedView ArticleGoogle Scholar
- Engeler B, Reyer HU. Choosy females and indiscriminate males: mate choice in mixed populations of sexual and hybridogenetic water frogs (Rana lessonae, Rana esculenta). Behav Ecol. 2001;12(5):600–6.View ArticleGoogle Scholar
- Reyer HU, Frei G, Som C. Cryptic female choice: frogs reduce clutch size when amplexed by undesired males. P Roy Soc Lond B Bio. 1999;266(1433):2101–7.View ArticleGoogle Scholar
- Schmeller DS, O’hara R, Kokko H. Male adaptive stupidity: male mating pattern in hybridogenetic frogs. Evol Ecol Res. 2005;7(7):1039–50.Google Scholar
- Hoffmann A, Reyer HU. Genomic effects on advertisement call structure in diploid and triploid hybrid waterfrogs (Anura, Pelophylax esculentus). BMC Ecol. 2013;13:47. doi:10.1186/1472-6785-13-47.PubMed CentralPubMedView ArticleGoogle Scholar
- Som C, Reyer HU. Demography and evolution of pure hybridogenetic frog (Rana esculenta) populations. Evol Ecol Res. 2006;8(7):1235–48.Google Scholar
- Uzzell T, Günther R, Berger L. Rana ridibunda and Rana esculenta: A leaky hybridogenetic system (Amphibia, Salientia). P Acad Nat Sci Phila. 1977;128(9):147–71.Google Scholar
- Tunner HG, Heppich-Tunner S. Genome exclusion and two strategies of chromosome duplication in oogenesis of a hybrid frog. Naturwissenschaften. 1991;78(1):32–4.View ArticleGoogle Scholar
- Ogielska M. Rana esculenta developmental syndrom: Fates of abnormal embryos from the first cleavage until spontaneous death. Zool Polon. 1994;39:447–59.Google Scholar
- Günther R. Über die verwandtschaftlichen Beziehungen zwischen den europäischen Grünfröschen und den Bastardcharakter von Rana esculenta L. (Anura). Zool Anz. 1973;190:250–85.Google Scholar
- Ogielska M, Bartmanska J. Development of testes and differentiation of germ cells in water frogs of the Rana esculenta-complex (Amphibia, Anura). Amphibia-Reptilia. 1999;20(3):251–63.View ArticleGoogle Scholar
- Günther R. Untersuchungen der Meiose bei Männchen von Rana ridibunda Pall., Rana lessonae Cam. und der Bastardform “Rana esculenta” L. Biol Zentralbibl. 1975;94:277–94.Google Scholar
- Alves MJ, Coelho MM, Collares-Pereira MJ. Diversity in the reproductive modes of females of the Rutilus alburnoides complex (Teleostei, Cyprinidae): a way to avoid the genetic constraints of uniparentalism. Mol Biol Evol. 1998;15:1233–42.View ArticleGoogle Scholar
- Morishima K, Yoshikawa H, Arai K. Meiotic hybridogenesis in triploid Misgurnus loach derived from a clonal lineage. Heredity. 2008;100(6):581–6.PubMedView ArticleGoogle Scholar
- Vinogradov AE, Borkin LJ, Guenther R, Rosanov JM. Genome elimination in diploid and triploid Rana esculenta males: cytological evidence from DNA flow cytometry. Genome. 1990;33(5):619–27.PubMedView ArticleGoogle Scholar
- Graf JD, Polls PM. Evolutionary genetics of the Rana esculenta complex. In: Dawley RM, Bogart JP, editors. Ecology and Evolution of Unisexual Vertebrates. Albany, New York: New York State Museum; 1989. p. 289–301.Google Scholar
- Mikulíček, Kotlík P. Two water frog populations from western Slovakia consisting of diploid females and diploid and triploid males of the hybridogenetic hybrid Rana esculenta (Anura, Ranidae). Mitt Mus Nat kd Berlin, Zool Reihe. 2001;77:59–64.Google Scholar
- Rybacki M. Diploid males of Rana esculenta from natural populations in Poland producing diploid spermatozoa. Zool Polon. 1994;39:517–8.Google Scholar
- Hotz H, Mancino G, Bucciinnocenti S, Ragghianti M, Berger L, Uzzell T. Rana ridibunda varies geographically in inducing clonal gametogenesis in interspecies hybrids. J Exp Zool. 1985;236(2):199–210.View ArticleGoogle Scholar
- Guerrini F, Bucci S, Ragghianti M, Mancino G, Hotz H, Uzzell T, et al. Genomes of two water frog species resist germ line exclusion in interspecies hybrids. J Exp Zool. 1997;279(2):163–76.PubMedView ArticleGoogle Scholar
- Ragghianti M, Bucci S, Marracci S, Casola C, Mancino G, Hotz H, et al. Gametogenesis of intergroup hybrids of hemiclonal frogs. Genet Res. 2007;89(1):39–45.PubMedView ArticleGoogle Scholar
- Joly P. The future of the selfish hemiclone: a Neodarwinian approach to water frog evolution. Mitt Mus Nat kd Berlin, Zool Reihe. 2001;77:31–8.Google Scholar
- Hurst LD, Atlan A, Bengtsson BO. Genetic conflicts. Q Rev Biol. 1996;71(3):317–64.PubMedView ArticleGoogle Scholar
- Pruvost NBM, Mikulicek P, Choleva L, Reyer H-U. Contrasting reproductive strategies of triploid hybrid males in vertebrate mating systems. J Evol Biol. 2014;28:189–204.View ArticleGoogle Scholar
- Schmeller DS. Tying ecology and genetics of hemiclonally reproducing waterfrogs (Rana, Anura). Ann Zool Fennici. 2004;41:681–7.Google Scholar
- Hermaniuk A, Pruvost NBM, Kierzkowski P, Ogielska M. Genetic and cytogenetic characteristics of pentaploidy in water frogs. Herpetologica. 2013;69:36–45.View ArticleGoogle Scholar
- Ogielska M, Kierzkowski P, Rybacki M. DNA content and genome composition of diploid and triploid water frogs belonging to the Rana esculenta complex (Amphibia, Anura). Can J Zool. 2004;82(12):1894–901.View ArticleGoogle Scholar
- Berger L. On the origin of genetic systems in European water frog hybrids. Zool Polon. 1990;35(1–4):5–32.Google Scholar
- Negovetic S, Anholt BR, Semlitsch RD, Reyer HU. Specific responses of sexual and hybridogenetic European waterfrog tadpoles to temperature. Ecology. 2001;82(3):766–74.View ArticleGoogle Scholar
- Pruvost NBM, Hollinger D, Reyer H-U. Genotype-temperature interactions on larval performance shape population structure in hybridogenetic water frogs (Pelophylax esculentus complex). Funct Ecol. 2013;27:459–71.View ArticleGoogle Scholar
- Altwegg R, Reyer H-U. Patterns of natural selection on size at metamorphosis in water frogs. Evolution. 2003;57(4):872–82.PubMedView ArticleGoogle Scholar
- Christiansen DG, Fog K, Pedersen BV, Boomsma JJ. Reproduction and hybrid load in all-hybrid populations of Rana esculenta water frogs in Denmark. Evolution. 2005;59(6):1348–61.PubMedView ArticleGoogle Scholar
- Muller HJ. The relation of recombination to mutational advance. Mutat Res. 1964;1:2–9.View ArticleGoogle Scholar
- Christiansen DG, Reyer H-U. From clonal to sexual hybrids: genetic recombination via triploids in all-hybrid populations of water frogs. Evolution. 2009;63(7):1754–68.PubMedView ArticleGoogle Scholar
- Arioli M. Reproductive patterns and population genetics in pure hybridogenetic water frog populations of Rana esculenta. [PhD thesis]: Zurich: University of Zurich; 2007.Google Scholar
- Arioli M, Jakob C, Reyer HU. Genetic diversity in water frog hybrids (Pelophylax esculentus) varies with population structure and geographic location. Mol Ecol. 2010;19(9):1814–28.PubMedView ArticleGoogle Scholar
- Hewitt GM. Post-glacial re-colonization of European biota. Biol J Linnean Soc. 1999;68:87–112.View ArticleGoogle Scholar
- Zeisset I, Beebee TJ. Determination of biogeographical range: an application of molecular phylogeography to the European pool frog Rana lessonae. Proc Biol Sci. 2001;268(1470):933–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Tunner HG. Evidence for genomic imprinting in unisexual triploid hybrid frogs. Amphibia-Reptilia. 2000;21(2):135–41.View ArticleGoogle Scholar
- Rieseberg LH. Hybrid origins of plant species. Ann Revi Ecol Syst. 1997;28:359–89.View ArticleGoogle Scholar
- Rieseberg LH, Willis JH. Plant speciation. Science. 2007;317:910–4.PubMed CentralPubMedView ArticleGoogle Scholar
- Schwenk K, Brede N, Streit B. Introduction. Extent, processes and evolutionary impact of interspecific hybridization in animals. Philos Trans Royal Soc London [Biol]. 2008;363:2805–11.View ArticleGoogle Scholar
- Vrijenhoek RC. Polyploid Hybrids: Multiple origins of a treefrog species. Curr Biol. 2006;16(7):R245–7.PubMedView ArticleGoogle Scholar
- Cunha C, Doadrio I, Coelho MM. Speciation towards tetraploidization after intermediate processes of non-sexual reproduction. Philos Trans Royal Soc London [Biol]. 2008;363:2921–9. doi:10.1098/rstb.2008.0048.View ArticleGoogle Scholar
- Ramsey J, Schemske DW. Pathways, mechanisms, and rates of polyploid formation in flowering plants. Annu Rev Ecol Syst. 1998;29:467–501.View ArticleGoogle Scholar
- Mable BK. Why polyploidy is rarer in animals than in plants:/’ myths and mechanisms. Biol J Linnean Soc. 2004;82:453–66.View ArticleGoogle Scholar
- Stöck M, Lamatsch DK, Steinlein C, Epplen JT, Grosse WR, Hock R, et al. A bisexually reproducing all-triploid vertebrate. Nat Genet. 2002;30(3):325–8.PubMedView ArticleGoogle Scholar
- Gosner KL. A simplified table for staging anuran embryos and larvae with note on identification. Herpetologica. 1960;16:183–90.Google Scholar
- Storie I, Sawle A, Goodfellow K, et al. Flow rate calibration I: A novel approach for performing absolute cell counts. Cytometry. 2003;55B:1–7.View ArticleGoogle Scholar
- Garner TWJ, Gautschi B, Röthlisberger S, Reyer HU. A set of CA repeat microsatellite markers derived from the pool frog. Rana lessonae Mol Ecol. 2000;9(12):2173–5.View ArticleGoogle Scholar
- Zeisset I, Rowe G, Beebee TJC. Polymerase chain reaction primers for microsatellite loci in the north European water frogs Rana ridibunda and R. lessonae. Mol Ecol. 2000;9(8):1173–4.PubMedView ArticleGoogle Scholar
- Berger L, Rybacki M, Hotz H. Artificial fertilization of water frogs. Amphibia-Reptilia. 1994;15:408–13.View ArticleGoogle Scholar
- Reyer H-U, Bättig I. Identification of reproductive status in female frogs – A quantitative comparison of nine methods. Herpetologica. 2004;60:349–57.View ArticleGoogle Scholar
- Reyer HU, Niederer B, Hettyey A. Variation in fertilisation abilities between hemiclonal hybrid and sexual parental males of sympatric water frogs (Rana lessonae, R. esculenta, R. ridibunda). Behav Ecol Sociobiol. 2003;54(3):274–84.View ArticleGoogle Scholar
- Semlitsch RD. Asymmetric competition in mixed populations of tadpoles of the hybridogenetic Rana esculenta complex. Evolution. 1993;47:510–9.View ArticleGoogle Scholar
- Sokal RR, Rohlf FJ. Biometry. San Francisco: W. H., Freeman & Co.; 1969.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.