Hybridizing Daphnia communities from ten neighbouring lakes: spatio-temporal dynamics, local processes, gene flow and invasiveness
© Yin et al.; licensee BioMed Central Ltd. 2014
Received: 10 December 2013
Accepted: 4 April 2014
Published: 12 April 2014
In natural communities of cyclical parthenogens, rapid response to environmental change is enabled by switching between two reproduction modes. While long periods of asexual reproduction allow some clones to outcompete others, and may result in “clonal erosion”, sexual reproduction restores genetic variation in such systems. Moreover, sexual reproduction may result in the formation of interspecific hybrids. These hybrids can then reach high abundances, through asexual clonal reproduction. In the present study, we explored genetic variation in water fleas of the genus Daphnia. The focus was on the short-term dynamics within several clonal assemblages from the hybridizing Daphnia longispina complex and the impact of gene flow at small spatial scales.
Daphnia individuals belonged either to the parental species D. galeata and D. longispina, or to different hybrid classes, as identified by 15 microsatellite markers. The distribution and genotypic structure of parental species, but not hybrids, corresponded well with the geographical positions of the lakes. Within parental species, the genetic distance among populations of D. galeata was lower than among populations of D. longispina. Moreover, D. galeata dominance was associated with higher phosphorous load. Finally, there was no evidence for clonal erosion.
Our results suggest that the contemporary structure of hybridizing Daphnia communities from ten nearby lakes is influenced by colonization events from neighbouring habitats as well as by environmental factors. Unlike the parental species, however, there was little evidence for successful dispersal of hybrids, which seem to be produced locally. Finally, in contrast to temporary Daphnia populations, in which a decrease in clonal diversity was sometimes detectable over a single growing season, the high clonal diversity and lack of clonal erosion observed here might result from repeated hatching of sexually produced offspring. Overall, our study provides insights into spatio-temporal dynamics in a hybridizing Daphnia species complex in a recently established lake system, and relates genetic similarities of populations to a scenario of secondary invasion enhanced by environmental factors.
KeywordsCyclical parthenogenesis Daphnia longispina complex Clonal richness Microsatellite Population structure
In natural communities, cyclically parthenogenetic reproduction is a complex and adaptable reproductive strategy with strong implications for ecological and evolutionary processes, such as genetic diversity, clonal erosion, interspecific hybridization and gene flow. In this mode of reproduction, one or more generations of parthenogenetically produced females alternate with a sexual generation formed by males producing sperm, and females producing haploid eggs. In cladoceran Daphnia, for example, females usually clone themselves by producing parthenogenetic daughters. Only when unfavourable conditions arise, such as food shortage, overcrowding, or change in temperature (e.g. ), do Daphnia individuals switch to the production of males and sexual haploid eggs that require fertilization, and then diapause. The general advantages of this reproductive mode are, on the one hand, enhanced genotypic diversity after sexual recombination and, on the other, high population growth rates during asexual periods (e.g. ).
The sexual phase of Daphnia’s reproductive cycle results not only in the production of new genotypes, but also in the formation of interspecific hybrids if closely related species co-occur (e.g. ). Via asexual propagation, these hybrids can then sometimes reach high abundances (e.g. [4, 5]). Such hybridization is especially common among the members of the Daphnia longispina complex and widespread among the species D. galeata and D. longispina (taxonomy revised in ). Most previous work on spatial patterns in the D. longispina complex explored Daphnia communities across large geographical regions (e.g. [5, 7]). For example, Keller et al.  found that although D. galeata and D. longispina are present on both sides of the Alps, D. galeata is more common in lakes south of the Alps, whereas D. longispina dominates in the north. In addition to this apparent geographical gradient, differences in environmental conditions also play a role . However, little is known about spatial variation in the hybridizing D. longispina complex over small scales, where dispersal among neighbouring habitats might be of great importance. The main stage for Daphnia’s passive dispersal involves sexual diapause eggs being carried by water birds or transported by wind . Since parental species produce more conspecific- than interspecific-crosses (therefore, more diapause eggs represent parental than hybrid genotypes), and since the “hatching success” of hybrid diapause eggs is inherently lower than of parental offspring [9, 10], parental species are expected to have higher gene flow capacities than hybrids .
In cyclical parthenogens, long periods of strict asexual reproduction might result in declines of clonal diversity. This process has been described in life history models as “clonal erosion” [11, 12]. In the field, such a decrease in clonal diversity over the growing season has been reported for populations of Daphnia magna (e.g. ) and of cyclically parthenogenetic rotifers (e.g. ). With respect to Daphnia, evidence for clonal erosion comes from populations inhabiting temporary ponds (e.g. ). In these habitats, sexual and asexual periods are strictly separated, with sexually produced offspring hatching from diapause eggs mainly in spring, giving rise to a large number of clonal lineages. Due to subsequent clonal competition, genetic diversity decreases gradually. In contrast, previous studies of Daphnia populations from permanent lakes did not detect a consistent decrease in clonal diversity over a growing season (e.g. [4, 9]). In these and many other studies, however, the “clones” were defined by few allozyme loci and might actually have represented clonal groups instead of single lineages . This methodological constraint could bias the results in such a way that a potential decrease in genotypic diversity might have been undetectable. Indeed, clonal erosion might also shape the genetic structure in overwintering Daphnia populations, as argued in a study applying higher resolution microsatellite markers . However, this conclusion remains speculative, because it was derived from an analysis of genotypic structure of Daphnia populations only at the end of the growing season. In a more detailed summer-autumn field survey of the D. longispina complex in two large lakes, clonal diversity was traced over time by genotyping Daphnia at 10 polymorphic microsatellite loci . This set of markers is able to detect true clonal lineages , and the clonal diversity was found to be roughly constant . In that study , the exclusion of the spring period, during which thousands of new genotypes are expected to hatch from diapause eggs, may have prevented the detection of an initial drop in diversity (e.g. [13, 19]). Consequently, to rigorously test for clonal erosion in permanent habitats, a survey must start early in the year when genetic diversity is expected to be highest, and changes must be followed over the whole growing season until late autumn, otherwise a decline might remain undetected.
In the present study, we explored spatial patterns and temporal dynamics in the genotypic structure of Daphnia communities from several lakes inhabited by taxa of the D. longispina complex. This was undertaken over a small spatial scale, with study sites separated by 30 km at most. We applied high-resolution genotyping (15 microsatellite loci). Specifically, we tested if the relative frequencies of parental species and their hybrids correspond with the geographical proximity and/or environmental descriptors of the lakes. Then, we looked for changes in taxon composition between spring and autumn seasons. In addition, we estimated gene flow and compared genotypic similarities among populations belonging to either parental species or hybrids, in order to find evidence for invasiveness and/or local production of hybrid taxa. Finally, some lakes were tracked throughout one growing season to determine if clonal diversity decreases over time (i.e., clonal erosion).
Spatial patterns in taxon composition
Environmental preferences of parental taxa
Temporal changes in taxon composition
Over time, Daphnia taxon composition changed significantly between spring and autumn in two of the seven lakes tested. In FELD, the community changed from complete dominance by F1 hybrids to complete dominance by D. galeata; in LUSS, D. longispina increased its abundance in comparison to the F1 hybrids (Figure 2).
Genetic distances among spatially and temporally isolated populations of parental species and hybrids
Gene flow and local production of hybrids
The mean number of migrants among lakes (all pairwise comparisons per taxon) differed between the parental species and hybrids (df = 2, F = 4.73, P = 0.023). Specifically, there were significantly more migrants among D. galeata populations (Nm = 4.26, SD = 4.63, N = 6) than among D. longispina (Nm = 0.84, SD = 0.28, N = 9) or hybrids (Nm = 0.37, SD = 0.10, N = 6), meaning that gene flow was largest for D. galeata (Post-Hoc test: P < 0.05). Then, evidence for local production of hybrids was found. Specifically, the hybrids from LUSS (i.e. a lake where hybrids coexisted with parental species) were more likely assigned to the respective local parental D. longispina population than to populations of D. longispina from other lakes (df = 17, t = -14.2, P < 0.001).
Test for clonality
Among 1856 individuals with complete multilocus genotypes (MLG) profiles (or data missing solely at the SwiD2 locus) 867 unique MLGs were detected. The hypothesis that individuals sharing identical MLGs were of sexual origin was rejected (P < 0.001 in all cases), indicating that individuals with identical MLGs can be treated as the same clone.
Temporal changes in clonal composition
In the present study, we surveyed spatio-temporal dynamics in hybridizing communities of Daphnia originating from ten neighbouring lakes. This is a first insight into genetic diversity, clonal erosion, hybridization and gene flow on such a small spatial scale which is tightly linked to the adaptable reproductive strategy in species complexes formed by cyclical parthenogens.
Spatial pattern and importance of environmental factors on taxon composition
A geographical separation was observed between communities dominated by D. galeata or D. longispina despite presumably ongoing dispersal between the neighbouring lakes. In a first step, lakes have been possibly colonized by extant neighbouring populations in the past. For example, Lußsee (LUSS) could have been colonized by D. longispina from Langwieder See (LANG), a lake created six decades earlier, and only 100 m away (see Figure 2). In Daphnia, and other zooplankton species, the advantage of first arrivals is relatively strong, due to their rapid population growth rate (e.g. [22, 23]). Indeed, from studies of the colonization of artificial ponds, it has been reported that Daphnia populations were founded by only a few genotypes [24, 25]. Moreover, after initial colonization, no new alleles were observed in these populations, indicating limited establishment success for genotypes arriving later . However, it has also been observed that changes in environmental conditions might result in species shifts, especially in the D. longispina communities (e.g. ). Thus, independent from geographical proximity, differences in lake environmental conditions may explain part of the observed spatial variation in taxon composition. Indeed, PCA and corresponding statistical tests showed that D. galeata occurred in more eutrophic (and older) habitats, and D. longispina in less eutrophic (and younger) ones. This pattern is consistent with a postulate of Flößner and Kraus  that D. galeata is promoted by eutrophication. Moreover, by comparing the physical parameters and water quality of several lakes, Keller et al.  concluded that D. longispina prefers habitats with low phosphorus load. Species-specific contrasting preferences were confirmed by the detection of a gradual replacement of D. longispina with D. galeata during periods of eutrophication , and by follow-up experimental tests showing differential fitness of these two parental species under various phosphorus loads .
Surprisingly, although D. galeata occurred exclusively in the oldest lakes, the genetic distance among populations of D. galeata was ~1/5th of that among populations of D. longispina, implying lower divergence time for D. galeata. This could possibly be explained by the fact that D. galeata spread later and, being more invasive and a better colonizer as suggested by the highest number of migrants, replaced former residents. Assuming that newly established artificial lake habitats are oligotrophic in the beginning and therefore first colonized by D. longispina, eutrophication might have allowed D. galeata to invade (see ). Overall, in our system the geographical separation between D. galeata and D. longispina dominated communities is consistent with the environmental conditions of the occupied habitats. Thus, geographical and environmental forces are not exclusive, and both most likely play a role in shaping Daphnia hybridizing communities.
Temporal pattern and lack of clonal erosion
The relative proportion of Daphnia taxa changed significantly between spring and autumn in two of the seven tested lakes. In addition, there was significant replacement of clones over time in two out of the three lakes for which a complete set of time series data was available. The observed changes may be due to the impact of environmental variation on the relative fitness of parental species and hybrids, as well as the fitness of specific clones within species, as has been shown in previous experimental surveys (e.g. [29–31]). However, although extinction of clones due to selection and random events should lead to an overall decrease in clonal diversity towards the end of the growing season , no signs of clonal erosion were detected in our study. Consistent with previous reports for permanent lakes (e.g. [4, 9]), in our system the indices for clonal diversity remained roughly constant over the growing season, both for D. galeata and D. longispina populations (hybrids were not tested because of their small sample size). The presence of clonal erosion in Daphnia populations may thus depend on the habitat type [11, 32]. Specifically, in contrast to temporary pond populations with clear separation of sexual and asexual phases in the Daphnia reproductive cycle (e.g. ), and for which the importance of clonal erosion has been documented (e.g. [13, 34]), the two modes of reproduction often overlap in the larger permanent lakes (e.g. [9, 10, 35]). This can contribute to the maintenance of high clonal diversity throughout the year. Moreover, given that fitness of genotypes in the D. longispina complex varies strongly with environment (e.g. [30, 31]), seasonal environmental changes may favour the maintenance of a high number of genotypes through fluctuating selection pressures.
Local production of hybrids
In contrast to parental species, neither the frequency of hybrids in the community, nor the clonal composition of hybrid populations was associated with the geographical neighbourhood of the lakes. Also, the UPGMA analysis showed that the highest genetic distance was detected among populations belonging to hybrids. Gene flow among hybrids was negligible and hybrids were more related to parents from the same lake than to parents from other lakes. These results complement evidence from previous studies showing that D. galeata × D. longispina hybrids have impaired hatching success from diapause eggs [9, 10], possibly restricting their gene flow. This all suggests that these hybrids are rather locally produced, by multiple hybridization events .
In summary, by applying high-resolution microsatellite markers, we addressed spatial and temporal patterns in the hybridizing D. longispina complex on a small scale, in ten recently established lakes. We found that the distribution and genotypic structure of parental Daphnia species was consistent with the geographical position and the differences in environmental attributes of the lakes; in particular, the presence of D. galeata was related to higher phosphorous load. Furthermore, the genetic distance among populations of D. galeata was lower than among populations of D. longispina, implying lower divergence time for D. galeata. The explanation for this may be that D. galeata spread later, being more invasive and a better colonizer, and then replaced former residents in lakes with high trophy level. Finally, no clonal erosion was detected across the studied lakes, suggesting that this phenomenon may depend on the habitat type and the extent of hatching of sexual offspring. Overall, the Daphnia system is a useful model to gain a better understanding of dynamics in hybridizing communities formed by cyclical parthenogens.
Local administration has been contacted at the beginning of the study. Collection of zooplankton (Daphnia) did not require specific permissions as these samples were obtained from unprotected lakes that are open for public activities. Our study did not involve the use or collection of endangered or protected species.
Study sites and field sampling
Daphnia individuals were sampled from ten flooded gravel pits in and around Munich (Germany; Additional file 1: Table S1). The lakes were created between 1930 and 2000, most are small and shallow, and the maximum distance between them is about 30 km. Geographical coordinates, lake age, lake size, maximum depth, total phosphorus concentration (Ptot), nitrogen concentration (N) and Daphnia densities (averaged across monthly densities from April to November 2011) are given in the supplementary material (Additional file 2: Table S2). The phosphorus and nitrogen data were provided by Andreas Scholz from the Wasserwirtschaftsamt München (http://www.wwa-m.bayern.de/). These measurements were taken in spring, from the surface to the bottom of the lake at one-meter intervals (average Ptot and N values were then calculated). Zooplankton samples were collected monthly from April to November 2011 (eight samples per lake; sampling interval was about 30 days). A 95-μm plankton net was hauled through the whole water column at two different sites per lake within the deep basin (depth was measured using a portable depth sounder); the two samples were then pooled and preserved in 95% ethanol. A random subsample was used for density counts (only adult females were taken into consideration). Then, ca. 50 adult females from the D. longispina complex were chosen randomly per sample for genotyping. Since Daphnia densities were sometimes very low, a complete time series of eight samples could be genotyped for only three lakes (FASA, LERC and WALD; for lake abbreviations see Additional file 1: Table S1). Two to seven samples were genotyped for the seven remaining lakes, resulting in 1934 Daphnia individuals being genotyped in total.
The DNA of each individual Daphnia was extracted and genotyped at 15 microsatellite markers , in two multiplex polymerase chain reactions, following a protocol described elsewhere . PCR products were analysed on an ABI PRISM 3700 capillary sequencer using a LIZ 500 labelled size standard. Genotypes were scored using GeneMapper version 3.7 (Applied Biosystems). Alleles at each locus were defined by their fragment length (in base pairs). Across different runs of genotyping, the consistency of alleles was checked, with locus-specific patterns of one reference genotype used in each run. This enabled us to adjust alleles with small differences in fragment length among different runs. In addition, there was no evidence of scoring errors due to stuttering, large allele dropout, or the presence of null alleles , as indicated by respective tests using MICRO-CHECKER 2.2.3 (104 permutations, ).
To display the genetic relatedness of multilocus genotypes (MLGs), factorial correspondence analysis (FCA) was applied in GENETIX 4.05 , on all unique MLGs together with 49 well-defined reference genotypes (same set of genotypes as in ), representing three parental species (D. cucullata, D. galeata, D. longispina) and their interspecific hybrids. It was confirmed that only D. galeata and D. longispina parental species were present in our dataset. Bayesian statistics were then used (in NewHybrids 1.1, ) to assign individuals to one of six possible classes (i.e., two parental species, F1 and F2 hybrids and both backcrosses). The probability threshold for assignments was set to 95% (106 iterations after a burn-in of length 106). All the following calculations of taxon-specific parameters are based on the taxon assignment in NewHybrids.
Spatial patterns in taxon composition
To test if the taxonomic similarity of Daphnia communities corresponds with the geographical proximity of the lakes, the correlation between pairwise Euclidean distances based on the taxon composition and respective pairwise geographical distances between sampling sites was evaluated using a Mantel test (5000 permutations, Past v2.07, ). Two estimates of Euclidean distances were derived in Past v2.07 for the spring samples (i.e., April or May), based on either the frequency, or the presence/absence, of five Daphnia groups. These groups included all classes as identified by NewHybrids (i.e., D. galeata, D. longispina, F1 hybrids and F2 hybrids; backcrosses were not detected) and an additional group of individuals that could not be identified by this method. All ten lakes were included in these tests.
Environmental preferences of parental taxa
A principal components analysis (PCA) of five environmental variables (i.e., lake age, lake size, maximum depth, total phosphorus concentration (Ptot), nitrogen concentration (N)) and mean Daphnia density (see Additional file 2: Table S2) was used to relate Daphnia species occurrence to combinations of these parameters. Variables that did not conform to normality were log-transformed. To test if Daphnia species were nonrandomly distributed along each component of the PCA, a one-way ANOVA was used with the presence/absence of each Daphnia species (i.e., D. galeata or D. longispina) as a main factor. The species presence/absence category was assigned based on the whole-year dataset (Additional file 1: Table S1). All ten lakes were included in these tests.
Temporal changes in taxon composition
The relative frequencies of five Daphnia groups (i.e., four classes: D. galeata, D. longispina, F1 hybrids and F2 hybrids identified by NewHybrids, and one group of individuals not identified by this method) were compared between spring (April or May) and autumn (October or November or, in the case of FELD, September), by applying a Monte Carlo approach with 105 simulation runs . Since separate tests were run for each of the seven lakes tested here (BOHM, FASA, FELD, HEIM, LERC, LUSS and WALD), sequential Bonferroni corrections  were applied when interpreting the results.
Genetic distances among spatially and temporally isolated populations of parental species and hybrids
Pairwise Nei’s genetic distances  were calculated in GENALEX 6 , among populations belonging to D. galeata, D. longispina or hybrids (F1 and F2 NewHybrids classes were pooled into one group of “hybrids”, because of the low sample size of the F2 class). Population samples originating from different months were considered separately, and only populations with sample size N ≥ 6 were included. Then, genetic similarities among population samples were displayed using the Unweighted Pair-Group Method with Arithmetic Mean (UPGMA), as calculated in MEGA 4 .
Gene flow and local production of hybrids
The number of migrants between lakes was calculated per generation (Nm) to estimate gene flow among populations of the two parental species (i.e., D. galeata and D. longispina) and the hybrids (F1 and F2 NewHybrids classes were pooled into one group of “hybrids”, because of the low sample size of the F2 class), respectively (in GENETIX 4.05). The whole-year dataset was pooled per taxon and lake. For each of the three taxa, only lake samples with N > 10 were considered to represent a local population and included in this analysis. Then, one-way ANOVA and Post-Hoc Tests were used to test for taxon-specific differences in the number of migrants (in SPSS 20.0).
In addition, it has been tested if hybrids are more likely formed by the parental species inhabiting their home lake (i.e., local production of hybrids) than by the parental species from other lakes. Thus, for each individual hybrid genotype, the log-likelihood with which it was assigned to parental groups (specific parental species from the home lake versus pooled populations of the same parental species from the other lakes) was calculated in Arlequin 3.0 . To distinguish statistical differences in respective log-likelihoods of hybrid assignments, a paired t-test (103 bootstrap replications) was run in SPSS 20.0. Although hybrids and respective parental species co-existed in FELD and LUSS, the test above could only be performed for LUSS hybrids, because only one hybrid genotype was detected in FELD.
Test for clonality
To test the null hypothesis that an MLG encountered more than once was the result of sexual recombination, instead of clonal propagation, the Psex value was calculated in GENCLONE 2.0 . The calculations were performed individually for each abundant taxon and lake. Only individuals characterized at all 15 microsatellite loci (or with missing data at a single locus, SwiD2; many D. longispina individuals could not be amplified at this locus) were included in these and the follow-up analyses of clonal structure (i.e., 1856 of 1934 genotyped Daphnia, see Additional file 1: Table S1).
Temporal changes in clonal composition
To test if differences in clonal composition increase with the time interval between samples, the clonal turnover rate was calculated between sampling dates (for all possible pairwise comparisons), per taxon and lake, by computing complements of the Morisita-Horn index (1-MH ; see also ) using the program SPADE  and 104 bootstrap replications. This turnover rate takes into account not only changes in the presence of individual MLGs but also changes in MLG frequencies. Turnover rate varies from 0 (complete similarity) to 1 (no similarity). Then, the correlation between the pairwise turnover rate in clonal composition and the corresponding pairwise time differences between sampling dates (calculated in monthly intervals) was evaluated by a Mantel test (5000 permutations, Past v2.07). Finally, evidence was sought for the existence of an overall decrease in clonal diversity during the course of the growing season (i.e., clonal erosion). Clonal diversity was measured per taxon and sampling date as: (i) clonal richness R, and (ii) the inverse of the Simpson index. Clonal richness R was calculated as R = (G-1) / (N-1), where G is the number of genotypes and N represents sample size . R varies from 0 (all individuals belong to one clone) to 1 (all individuals are different). The inverse of the Simpson index was calculated as 1/(∑Pi2), where Pi is the proportion of the ith clone. This index thus includes additional information about the distribution of clones within samples. The minimum value of this index is 1 (all individuals belong to one clone) whereas the maximum value is the number of unique clones in the sample. Linear regression analyses were applied to test for the effect of time on clonal diversity (i.e., either clonal richness R or the inverse of the Simpson index) using R-software . The tests were run only on those populations for which samples spanning the entire growing season were available: D. galeata from FASA, D. galeata from LERC and D. longispina from WALD (a minimum sample size was set of 30 individuals).
We would like to thank Rita Jaenichen who organized and carried out the zooplankton sampling. Further, Rita Jaenichen and Sabine Radetzki helped with the genotyping whereas Monika Poxleitner counted the zooplankton samples for density estimates. Finally, we thank Andreas Scholz from Wasserwirtschaftsamt München who provided environmental data, and Thomas Hovestadt and two anonymous reviewers for comments on the manuscript and Mark Phillipo for the linguistic help. This research was funded by the German Science Foundation grant to JW (WO 1587⁄4–1) and the National Natural Science Foundation of China grant to MY (31200284).
- Pijanowska J, Stolpe G: Summer diapause in Daphnia as a reaction to the presence of fish. J Plankton Res. 1996, 18 (8): 1407-1412. 10.1093/plankt/18.8.1407.View ArticleGoogle Scholar
- Barton NH, Charlesworth B: Why sex and recombination?. Science. 1998, 281 (5385): 1986-1990.PubMedView ArticleGoogle Scholar
- Hebert PDN: Interspecific hybridization between cyclic parthenogens. Evolution. 1985, 39: 216-220. 10.2307/2408534.View ArticleGoogle Scholar
- Spaak P: Temporal changes in the genetic structure of the Daphnia species complex in Tjeukemeer, with evidence for backcrossing. Heredity. 1996, 76: 539-548. 10.1038/hdy.1996.77.View ArticleGoogle Scholar
- Keller B, Wolinska J, Manca M, Spaak P: Spatial, environmental and anthropogenic effects on the taxon composition of hybridizing Daphnia. Phil Trans R Soc B. 2008, 363 (1505): 2943-2952. 10.1098/rstb.2008.0044.PubMedPubMed CentralView ArticleGoogle Scholar
- Petrusek A, Hobaek A, Nilssen JP, Skage M, Cerny M, Brede N, Schwenk K: A taxonomic reappraisal of the European Daphnia longispina complex (Crustacea, Cladocera, Anomopoda). Zool Scr. 2008, 37 (5): 507-519. 10.1111/j.1463-6409.2008.00336.x.View ArticleGoogle Scholar
- Petrusek A, Seda J, Macháček J, Ruthová Š, Šmilauer P: Daphnia hybridization along ecological gradients in pelagic environments: the potential for the presence of hybrid zones in plankton. Phil Trans R Soc B. 2008, 363 (1505): 2931-2941. 10.1098/rstb.2008.0026.PubMedPubMed CentralView ArticleGoogle Scholar
- Havel JE, Shurin JB: Mechanisms, effects, and scales of dispersal in freshwater zooplankton. Limnol Oceanogr. 2004, 49 (4): 1229-1238.View ArticleGoogle Scholar
- Keller B, Spaak P: Nonrandom sexual reproduction and diapausing egg production in a Daphnia hybrid species complex. Limnol Oceanogr. 2004, 49 (4): 1393-1400.View ArticleGoogle Scholar
- Keller B, Wolinska J, Tellenbach C, Spaak P: Reproductive isolation keeps hybridizing Daphnia species distinct. Limnol Oceanogr. 2007, 52 (3): 984-991. 10.4319/lo.2007.52.3.0984.View ArticleGoogle Scholar
- de Meester L, Vanoverbeke J, de Gelas K, Ortells R, Spaak P: Genetic structure of cyclic parthenogenetic zooplankton populations - a conceptual framework. Arch Hydrobiol. 2006, 167 (1–4): 217-244.View ArticleGoogle Scholar
- Vanoverbeke J, de Meester L: Clonal erosion and genetic drift in cyclical parthenogens - the interplay between neutral and selective processes. J Evol Biol. 2010, 23 (5): 997-1012. 10.1111/j.1420-9101.2010.01970.x.PubMedView ArticleGoogle Scholar
- Carvalho GR, Crisp DJ: The clonal ecology of Daphnia magna. I. Temporal changes in the clonal structure of a natural population. J Anim Ecol. 1987, 56: 453-468. 10.2307/5060.View ArticleGoogle Scholar
- Ortells R, Gomez A, Serra M: Effects of duration of the planktonic phase on rotifer genetic diversity. Arch Hydrobiol. 2006, 167: 203-216. 10.1127/0003-9136/2006/0167-0203.View ArticleGoogle Scholar
- Thielsch A, Brede N, Petrusek A, de Meester L, Schwenk K: Contribution of cyclic parthenogenesis and colonization history to population structure in Daphnia. Mol Ecol. 2009, 18 (8): 1616-1628. 10.1111/j.1365-294X.2009.04130.x.PubMedView ArticleGoogle Scholar
- Hamrova E, Mergeay J, Petrusek A: Strong differences in the clonal variation of two Daphnia species from mountain lakes affected by overwintering strategy. BMC Evol Biol. 2011, 11: 231-10.1186/1471-2148-11-231.PubMedPubMed CentralView ArticleGoogle Scholar
- Yin M, Petrusek A, Seda J, Wolinska J: Fine-scale temporal and spatial variation of taxon and clonal structure in the Daphnia longispina hybrid complex in heterogeneous environments. BMC Evol Biol. 2012, 12: 12-10.1186/1471-2148-12-12.PubMedPubMed CentralView ArticleGoogle Scholar
- Yin M, Wolinska J, Gießler S: Clonal diversity, clonal persistence and rapid taxon replacement in natural populations of species and hybrids of the Daphnia longispina complex. Mol Ecol. 2010, 19 (19): 4168-4178. 10.1111/j.1365-294X.2010.04807.x.PubMedView ArticleGoogle Scholar
- Hairston NG, Hansen AM, Schaffner WR: The effect of diapause emergence on the seasonal dynamics of a zooplankton assemblage. Freshwat Biol. 2000, 45 (2): 133-145. 10.1046/j.1365-2427.2000.00386.x.View ArticleGoogle Scholar
- Anderson EC, Thompson EA: A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002, 160 (3): 1217-1229.PubMedPubMed CentralGoogle Scholar
- Nei M: Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978, 89 (3): 583-590.PubMedPubMed CentralGoogle Scholar
- de Meester L, Gomez A, Okamura B, Schwenk K: The monopolization hypothesis and the dispersal-gene flow paradox in aquatic organisms. Acta Oecol. 2002, 23 (3): 121-135. 10.1016/S1146-609X(02)01145-1.View ArticleGoogle Scholar
- Mergeay J, de Meester L, Eggermont H, Verschuren D: Priority effects and species sorting in a long paleoecological record of repeated community assembly through time. Ecology. 2011, 92 (12): 2267-2275. 10.1890/10-1645.1.PubMedView ArticleGoogle Scholar
- Ortells R, Vanoverbeke J, Louette G, de Meester L: Colonization of Daphnia magna in a newly created pond: founder effects and secondary immigrants. Hydrobiologia. 2014, 723 (1): 167-179. 10.1007/s10750-013-1593-7.View ArticleGoogle Scholar
- Louette G, Vanoverbeke J, Ortells R, de Meester L: The founding mothers: the genetic structure of newly established Daphnia populations. Oikos. 2007, 116 (5): 728-741. 10.1111/j.0030-1299.2007.15664.x.View ArticleGoogle Scholar
- Brede N, Sandrock C, Straile D, Spaak P, Jankowski T, Streit B, Schwenk K: The impact of human-made ecological changes on the genetic architecture of Daphnia species. Proc Natl Acad Sci U S A. 2009, 106 (12): 4758-4763. 10.1073/pnas.0807187106.PubMedPubMed CentralView ArticleGoogle Scholar
- Flößner D, Kraus K: On the taxonomy of the Daphnia hyalina-galeata complex (Crustacea: Cladocera). Hydrobiologia. 1986, 137: 97-115. 10.1007/BF00004206.View ArticleGoogle Scholar
- Spaak P, Fox J, Hairston NG: Modes and mechanisms of a Daphnia invasion. P Roy Soc B-Biol Sci. 2012, 279 (1740): 2936-2944. 10.1098/rspb.2012.0280.View ArticleGoogle Scholar
- Spaak P, Vanoverbeke J, Boersma M: Predator induced life history changes and the coexistence of five taxa in a Daphnia species complex. Oikos. 2000, 89 (1): 164-174. 10.1034/j.1600-0706.2000.890118.x.View ArticleGoogle Scholar
- Wolinska J, Bittner K, Ebert D, Spaak P: The coexistence of hybrid and parental Daphnia: the role of parasites. P Roy Soc B-Biol Sci. 2006, 273: 1977-1983. 10.1098/rspb.2006.3523.View ArticleGoogle Scholar
- Seidendorf B, Boersma M, Schwenk K: Evolutionary stoichiometry: the role of food quality for clonal differentiation and hybrid maintenance in a Daphnia species complex. Limnol Oceanogr. 2007, 52: 385-394. 10.4319/lo.2007.52.1.0385.View ArticleGoogle Scholar
- de Meester L: Local genetic differentiation and adaptation in freshwater zooplankton populations: patterns and processes. Ecoscience. 1996, 3 (4): 385-399.Google Scholar
- Hebert PDN: Enzyme variability in natural populations of Daphnia magna. III. Genotypic frequencies in intermittent populations. Genetics. 1974, 77: 335-341.PubMedPubMed CentralGoogle Scholar
- Ruvinsky AO, Perelygin AA, Lobkov YI, Belyaev DK: Factors organizing and maintaining polymorphism in a cyclic parthenogenetic species - Daphnia pulex. Heredity. 1986, 57: 15-22. 10.1038/hdy.1986.81.View ArticleGoogle Scholar
- Hebert PDN: Genotypic characteristics of cyclic parthenogens and their obligately asexual derivates. The Evolution of sex and its Consequences. Edited by: Stearns S. 1987, Basel: Birkhäuser Verlag, 175-195.View ArticleGoogle Scholar
- Spaak P: Hybridization in the Daphnia galeata complex: are hybrids locally produced?. Hydrobiologia. 1997, 360: 127-133. 10.1023/A:1003157117667.View ArticleGoogle Scholar
- Brede N, Thielsch A, Sandrock C, Spaak P, Keller B, Streit B, Schwenk K: Microsatellite markers for European Daphnia. Mol Ecol Notes. 2006, 6 (2): 536-539. 10.1111/j.1471-8286.2005.01218.x.View ArticleGoogle Scholar
- Selkoe KA, Toonen RJ: Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers - appendix. Ecol Lett. 2006, 9 (5): 615-629. 10.1111/j.1461-0248.2006.00889.x.PubMedView ArticleGoogle Scholar
- van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4 (3): 535-538. 10.1111/j.1471-8286.2004.00684.x.View ArticleGoogle Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR. 5000, Montpellier (France): 1996–2004, Université de Montpellier II, Available from URL: http://www.univ-montp2.fr/Google Scholar
- Hammer Ø, Harper DAT, Ryan PD: PAST paleontological statistics software package for education and data analysis. Palaeontol Electroni. 2001, 4: 1-9.Google Scholar
- Sham PC, Curtis D: Monte-Carlo tests for associations between disease and alleles at highly polymorphic loci. Ann Hum Genet. 1995, 59: 97-105. 10.1111/j.1469-1809.1995.tb01608.x.PubMedView ArticleGoogle Scholar
- Rice WR: Analyzing tables of statistical tests. Evolution. 1989, 43: 223-225. 10.2307/2409177.View ArticleGoogle Scholar
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6 (1): 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform. 2005, 1: 47-50.Google Scholar
- Haond A, Belkhir K: GENCLONE: a computer program to analyse genotypic data, test for clonality and describe spatial clonal organization. Mol Ecol Notes. 2007, 7 (1): 15-17.View ArticleGoogle Scholar
- Magurran AE: Ecological Diversity and its Measurement. 1988, Cambridge: University PressView ArticleGoogle Scholar
- Seda J, Petrusek A, Machacek J, Smilauer P: Spatial distribution of the Daphnia longispina species complex and other planktonic crustaceans in the heterogeneous environment of canyon-shaped reservoirs. J Plankton Res. 2007, 29 (7): 619-628. 10.1093/plankt/fbm044.View ArticleGoogle Scholar
- Chao A, Chazdon RL, Colwell RK, Shen TJ: A new statistical approach for assessing similarity of species composition with incidence and abundance data. Ecol Lett. 2005, 8 (2): 148-159.View ArticleGoogle Scholar
- Dorken ME, Eckert CG: Severely reduced sexual reproduction in northern populations of a clonal plant, Decodon verticillatus (Lythraceae). J Ecol. 2001, 89 (3): 339-350. 10.1046/j.1365-2745.2001.00558.x.View ArticleGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2009, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.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.