The relative homogeneity of pelagic environments has been regarded as the reason for the absence of hybrid zones for hybridizing planktonic Daphnia (Crustacea: Cladocera); occasional dominance of interspecific hybrids over parental species was explained by their temporal superiority in fluctuating environments. However, water bodies with spatially varying environmental conditions might facilitate the formation of hybrid zones in plankton. We studied the distribution of species and hybrids of the Daphnia longispina complex in 11 canyon-shaped reservoirs, localities characterized by horizontal environmental gradients (particularly of food supply and size-selective predation); we also analysed patterns of carapace size and fecundity among coexisting taxa. Spatial distribution of taxa agreed with their ecological characteristics; those showing different affinities along longitudinal reservoir profiles differed in size according to the presumed fish predation gradient. Only hybrids of Daphnia galeata with Daphnia cucullata and D. longispina (=hyalina) were recorded. The latter two species preferred opposite ends of gradients, such spatial segregation probably explaining the absence of their hybrids. Distributional patterns were relatively stable in two consecutive summers, apart from a substantial decline of D. galeata×cucullata in the second year. The observed pattern of a hybrid-dominated zone in intermediate conditions suggests that local Daphnia hybrid zones may indeed form within reservoirs.
The reproductive mode of cladocerans, cyclical parthenogenesis, makes this crustacean group an interesting model for ecological and evolutionary studies. The ability to switch between parthenogenetic reproduction during favourable conditions and sexual reproduction associated with the formation of diapausing stages (‘resting eggs’), and the possibility to grow clonal lineages in laboratory conditions, and therefore to disentangle the effects of environmental factors, genetic background and phenotypic variation, makes especially members of the genus Daphnia (Crustacea: Cladocera: Anomopoda) favourite study organisms.
Hybridization between different Daphnia species has been documented in several species groups (Schwenk & Spaak 1995), with the resulting reproductive modes of hybrids ranging from parthenogenetically reproducing hybrids unable to form diapausing eggs (Hebert & Finston 1996), obligate parthenogens, often polyploid, producing asexual resting eggs (e.g. Hebert et al. 1989a; Dufresne & Hebert 1994; Hebert & Finston 1996) to hybrids capable of sexual reproduction and backcrossing with the parental taxa, as observed within the Daphnia longispina complex. Hybridizing species belonging to this complex (Daphnia galeata Sars, Daphnia cucullata Sars and the pelagic form of D. longispina (O. F. Müller)=D. hyalina Leydig; for taxonomy and nomenclature, see Petrusek et al. 2008) are ecologically particularly important, as they are among the most common Daphnia species inhabiting large European permanent water bodies. The occurrence of their hybrids has been reported from many natural as well as man-made lakes (e.g. Hebert et al. 1989b; Seda et al. 2007b; Wolinska et al. 2007).
Despite potential introgression and occasional hybrid dominance, parental species within the D. longispina complex remain genetically distinct, as interspecific hybrids exhibit various aspects of hybrid breakdown (Keller et al. 2007), such as lower hatching and survival rate (Schwenk et al. 2001) and lower sexual reproductive success (Keller & Spaak 2004). It has been therefore proposed that reproductive isolation effectively exists among hybridizing species within the complex (Keller et al. 2007). Hybrids are nevertheless fully competitive with parental species when reproducing parthenogenetically, and due to asexual reproduction may at least temporarily avoid hybrid breakdown. Some Daphnia interspecific hybrids may exhibit superior population rates of increase under certain conditions, which has been shown in both the laboratory and field conditions. This phenomenon has been attributed to heterosis (Repka et al. 1999) or a combination of advantageous traits of both parental species (Spaak & Hoekstra 1995; Declerck & De Meester 2003). As a result, hybrids may occasionally be more abundant than the parental species (e.g. Spaak & Hoekstra 1993; Seda et al. 2007b).
Unlike in terrestrial habitats, which often exhibit various gradients in factors affecting the distribution of taxa and often facilitating the presence of interspecific hybrids in intermediate conditions (i.e. in spatially restricted hybrid zones; Arnold 1997), lakes and ponds tend to offer much more homogeneous conditions for plankton, and ecologically different pelagic environments are usually isolated from each other in island-like nature (Schwenk & Spaak 1995). As the space used by parental taxa and hybrids typically overlaps in the whole lake area, the coexistence of hybrids and parental species of Daphnia in the pelagic environment within any particular lake may not be appropriately modelled by hybrid zone scenarios.
The temporal hybrid superiority hypothesis (Spaak & Hoekstra 1995) has been proposed to explain the often observed hybrid dominance within localities by the fluctuation of environmental conditions that may temporarily favour hybrid genotypes. The most important factors affecting the presence of species and hybrids of the D. longispina complex seem to be size-selective fish predation pressure and food level, as it is the case for pelagic Daphnia in general (Gliwicz 2003). However, other factors such as parasites (Wolinska et al. 2006, 2007) or food quality (Seidendorf et al. 2007) may also strongly affect the patterns of species coexistence.
To study the impact of variation in environmental factors affecting the success of the parental taxa and hybrids, three potential approaches are available. First, zooplankton communities from various isolated water bodies with differing ecological conditions may be compared (e.g. Keller 2007). Second, it might be possible to track long-term changes in community structure in localities with substantial temporal variation or gradual changes in key environmental factors (such as fish stock or trophic level), e.g. by combining historical samples and data with a palaeogenetic analysis of resting eggs from the sediment (Jankowski & Straile 2003; N. Brede, C. Sandrock, D. Straile, P. Spaak, T. Jankowski, B. Streit & K. Schwenk 2007, unpublished data).
Third, we may examine the variation in the spatial distribution of the parental and hybrid taxa within water bodies in which substantial environmental gradients exist. Vertical gradients are characteristic of thermally stratified water bodies such as temperate lakes (figure 1). Hybridizing Daphnia taxa may indeed show different spatial distributions within these gradients (Weider & Stich 1992; Seda et al. 2007b), and intraspecific genetic differentiation on a vertical scale has been documented as well (Seda et al. 2007a). The vertical segregation of taxa may be caused by variable selection pressures in different depths or by individual spatial preferences. However, pelagic daphnids often exhibit diel vertical migrations across the thermocline (e.g. King & Miracle 1995; Spaak et al. 2004), showing that vertical temperature and oxygen gradients on a relatively small spatial scale do not necessarily form substantial barriers for individual movement.
Although water bodies with substantial internal horizontal gradients are apparently much rarer, they may be particularly interesting, as the dispersal of pelagic species is less limited by physical barriers and the uneven distribution of taxa in areas connected by water may suggest ‘barriers’ of ecological nature, i.e. differential selection forces. Our model systems, deep canyon-shaped reservoirs (figure 1), provide good conditions for studying this phenomenon. The well-defined main reservoir tributary as a point source of nutrient input, the unidirectional flow of water down the reservoir, occasional turbidity stress from upstream sites and significantly higher affinity of fishes to upstream locations are all significant environmental factors for the maintenance of strong heterogeneity of the selective forces responsible for the overall longitudinal heterogeneity of reservoir zooplankton (Urabe 1990).
In our recent study (Seda et al. 2007b), we demonstrated significant differences in the horizontal distribution of species and hybrids of the D. longispina complex within canyon-shaped reservoirs in Central Europe (Czech Republic). The presence of two or even three hybridizing species in a particular reservoir was common, but each taxon showed different patterns in relation to horizontal gradients. The spatial distribution of taxa was in agreement with their ecological characteristics (Gliwicz 1990; Flößner 2000): D. galeata was a ubiquitous species; D. cucullata was generally more common in the high-food and high-predation environment of upstream reservoir regions; and D. longispina had the opposite tendency. The distribution of hybrids of D. galeata with the remaining two species partly overlapped with their non-galeata parental species; the preference of D. galeata×cucullata hybrids for upstream regions was especially characteristic.
The substantial environmental gradient, together with the non-random distribution of parental taxa and zones with a preferential occurrence of hybrids, suggests that unlike in lakes studied so far, canyon-shaped reservoirs could offer favourable conditions for the emergence of local within-lake hybrid zones of zooplankters. This hypothesis may be supported by the partial overlap in the distribution of parental taxa with differing ecological requirements and the ongoing formation of hybrid genotypes, which usually only attain ecologically relevant abundances in certain parts of any given reservoir.
In this paper, we summarize our current knowledge of Daphnia hybridization in reservoirs, based on the data from two consecutive years. We discuss the patterns of taxon coexistence within environmental gradients and their temporal stability, and evaluate whether the distribution patterns re-form after the winter season, when Daphnia populations either disappear completely from the water column or at least experience severe bottlenecks. We also describe in detail a characteristic canyon-shaped reservoir with three hybridizing Daphnia species. Finally, we focus on the variation in the taxon body size and fecundity as population characteristics that probably contribute to the hybrid success, and are strongly influenced by the selection forces in reservoirs.
2. Material and methods
(a) Locality selection
To study the spatial distribution of taxa and genotypes within environmental gradients in reservoirs, it is important to ensure that local processes are crucial for shaping the observed patterns, rather than immigration or import of individuals. The primary selection criterion for our study localities was their canyon-like morphology (i.e. elongated shape along the longitudinal axis and increasing depth from the inflow regions towards the dam), which may result in longitudinal environmental gradients; however, we excluded those reservoirs in which import of zooplankton from upstream water bodies could strongly affect the Daphnia taxon composition in the inflow regions. Altogether, we sampled 11 reservoirs located in the Czech Republic, covering a substantial range of reservoir size, depth and retention time (table 1). Schematic outlines and locations of the studied reservoirs are shown in Seda et al. (2007b), where also additional details about localities are provided. We analysed the longitudinal and vertical distributions of Daphnia species and hybrids in samples from two consecutive summer seasons, collected once in July 2004 and once in July 2005.
(b) Sample collection and processing
The methods of sample collection and processing, identical for both years, were described in detail in Seda et al. (2007b). At each sampling date, we collected zooplankton by vertical hauls of plankton nets (mesh size 170 μm) from three sampling stations along the reservoir main axis, which covered the longitudinal environmental gradients: upstream near the river inflow; in the centre of the reservoir; and at the deepest part of the reservoir near the dam (figures 1a and 3). Sampling at the downstream site also focused on the vertical spatial distribution in the stratified water body, therefore we collected separately zooplankton from the epi-, meta- and hypolimnion, using a closing net for deeper layers. Prior to sampling, we measured vertical profiles of temperature and dissolved oxygen to evaluate the extent of stratification and determine the borders of the epi-, meta- and hypolimnion (figure 1b). From each site/layer, we collected a quantitative sample preserved in 4% formaldehyde solution for the analysis of zooplankton abundance (to estimate the density of all Daphnia in the respective samples, regardless of their taxon composition) and an additional sample deep-frozen on site in liquid nitrogen for genetic analyses (to identify proportions of parental species and hybrid genotypes within the D. longispina complex).
We genetically analysed approximately 40–50 randomly selected adult Daphnia females (if available) from each deep-frozen sample. Before processing, each individual was digitally photographed under the microscope for subsequent measuring of its carapace size, and the number of eggs/embryos in its brood chamber was counted.
We compared the results of two different molecular marker systems available for the discrimination of species and hybrids within the complex. Primarily, we used allozyme electrophoresis on cellulose acetate gels (Hebert & Beaton 1989), which has been successfully applied in a number of studies. Two allozyme loci were scored: sAAT, supernatant amino aspartate transferase (EC 126.96.36.199); AO, aldehyde oxidase (EC 188.8.131.52). It was shown that sAAT and AO can be used as diagnostic markers to discriminate among D. galeata, D. longispina (named D. hyalina in the respective papers) and D. cucullata, and to identify their hybrids (Wolf & Mort 1986; Gießler 1997a). Daphnia galeata is virtually fixed for the F (fast) and D. longispina for S (slow) alleles at both loci; D. cucullata is fixed for the S− (very slow) allele at sAAT. The alleged potential for the discrimination of all species using the locus AO (Gießler 1997a), however, could not be fully used. Apparently, at least in some cases, D. cucullata AO alleles may not be reliably differentiated from those of other species. In samples from the Czech Republic, bands of D. cucullata often overlapped with those of D. galeata; similarly, in a North German lake, D. cucullata could not be well differentiated at AO from D. longispina (Spaak et al. 2004).
Using the allozyme markers, we classified Daphnia individuals as belonging to parental species (carrying alleles typical of one species only) or as genotypes of hybrid origin (sharing alleles of two different species). Some individuals sharing alleles of D. galeata and D. longispina showed patterns clearly suggesting backcrossing or the formation of later-generation hybrids. However, as the relative frequencies of such genotypes were usually low, and two markers do not allow for reliable assignment to a particular hybrid class, we pooled them with apparent F1 hybrids for the purpose of this study.
Additionally, we used a DNA-based method—the restriction fragment length polymorphism (RFLP) of the nuclear ribosomal internal transcribed spacer (ITS)—to validate allozyme results on a subset of samples from 2004 (1265 individuals from 9 out of 11 studied reservoirs; table 1). This method, originally developed by Billiones et al. (2004), uses the restriction of a short part of ITS1, the 5.8S ribosomal RNA gene, and a large part of ITS2, to obtain species-specific fragment patterns that are additive in hybrid genotypes. The original ITS-RFLP protocol commonly suffered from a point mutation in D. galeata alleles, which caused their misidentification as belonging to D. cucullata; we therefore applied a newly developed alternative double-digest protocol (Skage et al. 2007) to circumvent this problem.
To allow direct comparison of allozyme and ITS-RFLP patterns, DNA was prepared from Daphnia homogenates used for the allozyme electrophoresis. A 2.5 μl aliquot of the homogenate was transferred to 30 μl of a solution containing H3 buffer and proteinase K (Schwenk et al. 1998) and incubated at 55°C for 6–10 hours; proteinase was subsequently inactivated by heating to 95°C for 10 min. The subsequent DNA amplification and restriction by overnight incubation with endonucleases MbiI and Eco52I (Fermentas) mostly followed the protocol by Skage et al. (2007). However, we used an alternative forward primer, ITS-F-New (5′-GGT AAC CGC TGA ACC TCC TTC-3′; Skage et al. 2007), which provides longer amplified fragments to reliably differentiate between D. galeata pattern and potentially uncut PCR products. The banding patterns were interpreted according to Skage et al. (2007); individuals combining fragments from two species were identified as hybrids even if bands of one species were more intense than those of another species. Occasional very weak bands were nevertheless not considered. Based on the fit between the two marker systems (see §3) and the morphology of studied individuals, allozyme data were further used for evaluating the distribution of species and hybrids in our study.
(c) Statistical analyses
Multivariate analyses of Daphnia spatial distribution were performed with the software package Canoco for Windows v. 4.5 (ter Braak & Šmilauer 2002). Original counts (number of individuals identified as respective species or hybrids) were log transformed, and standardization by sample norm was used to focus the analyses on the differences in the relative proportion of individual taxa. To summarize occurrence patterns of Daphnia taxa visually, principal component analysis (PCA) was used. To test for differences in the community composition along the horizontal profiles, partial redundancy analysis (RDA), using reservoir identity as a covariate, was performed. We pooled the data from the epi-, meta- and hypolimnion at the downstream sampling station for the analysis of horizontal distribution, weighing taxon abundances from different layers by the respective layer volume. The significance of the relationship between species composition and spatial gradient was tested using a model-based type of a Monte Carlo permutation test (ter Braak & Šmilauer 2002).
The contribution of individual factors (unique effects of locality, of longitudinal or vertical position and their interaction) was quantified using three different partial RDA models (Lepš & Šmilauer 2003). In two analyses, one of the two factors was used as the explanatory variable and the other as a covariate; in the third analysis, the interaction term was tested. The partial analyses were needed also for the two main effects due to the unbalanced nature of the sampling design.
Differences among individual taxa in the carapace size and fecundity of adult females were tested using a nested-design ANOVA model. This model included, besides the effect of taxon, also the fixed effect of habitat (epi-, hypo-, metalimnion, inflow and middle) and the random effects of reservoir identity and sampling year (nested within reservoir identity). The size of carapace and fecundity of adult females were log transformed to achieve homogeneity of variances.
Owing to the partially hierarchical nature of the ANOVA model, there was no straightforward way to perform multiple comparisons allowing the testing of differences among individual taxa. In addition, only some differences were of interest, as not all taxon pairs co-occurred in the same samples, nor did they all hybridize. Therefore, we tested the differences among D. cucullata, D. galeata and their hybrids and among D. longispina, D. galeata and their hybrids separately, by performing six pairwise comparisons followed by Holm's correction for simultaneous tests (Holm 1979).
(a) Allozymes and ITS-RFLP: marker selection
Results of the two different nuclear marker sets used to identify taxa (ITS rDNA region and allozymes) were largely in agreement (85% of the 1265 individuals analysed by both methods), although not corresponding to each other completely. In most reservoirs, the choice of marker would not substantially affect the patterns of spatial distribution of parental species and hybrids. Disagreements in identification between allozyme and ITS-RFLP patterns were almost always observed for individuals from sites where species in question co-occurred and hybridization was common. A notable exception was the Stanovice Reservoir, in which only D. galeata was detected by allozymes in 2004 but 19% of the individuals showed ITS-RFLP patterns characteristic for D. galeata×longispina hybrids. Morphological characteristics of most of these individuals agreed with the allozyme identification, which was subsequently confirmed by the detailed analysis of 12 microsatellite loci classifying all analysed individuals as D. galeata (Š. Ruthová, A. Petrusek, J. Seda, A. Thielsch, N. Brede & K. Schwenk 2008, unpublished data). This prompted us to use allozyme markers for primary taxon identification, although occasional differences between the two marker systems indicated that a certain proportion of introgressed individuals remained undetected and the extent of hybridization may have been underestimated.
(b) Patterns of taxon coexistence and their temporal stability
All the three parental Daphnia species (D. galeata, D. longispina and D. cucullata) occurred commonly in the investigated set of reservoirs (table 1), and all three actually co-occurred in three of them. Out of the three potential interspecific hybrids, however, we recorded only two; hybrids D. cucullata×longispina were never observed. Daphnia galeata×longispina strongly dominated over D. galeata×cucullata in number of individuals as well as localities occupied (656 over 51 individuals out of 4464 analysed; seven versus four localities). In three reservoirs with the presence of hybrid genotypes, we did not record one of their parental species (table 1). In one of these, the Horka Reservoir, we recorded only D. longispina and hybrid genotypes with D. galeata; the allozyme patterns (sAAT heterozygous and AO homozygous for the D. longispina marker), however, suggested that these individuals did not represent the F1 hybrid generation but later-generation hybrids or backcrosses with D. longispina.
The taxon composition of the D. longispina complex significantly differed among localities (RDA, p=0.001, pseudo-F-statistic 10.832). The results of the PCA of the same data (figure 2) show that the 11 sampled reservoirs can be divided into three groups based on the presence of particular Daphnia species and hybrids. Four reservoirs in the upper half of the ordination diagram are characterized by high proportions of D. cucullata in at least some samples (typically upstream). The isolated position of the Horka Reservoir at the right edge reflects the fact that D. galeata was completely absent from this locality (the only case among the investigated reservoirs), and only D. longispina and hybrid or backcrossed genotypes were present. The remaining six reservoirs contained a substantial proportion of D. galeata, with D. cucullata being rare or absent. The positions of sampling sites within reservoirs in figure 2 (grey squares) indicates the general tendency of D. cucullata and its hybrids to occur more at upstream sites, and the opposite tendency of D. longispina and its hybrids. The overall Daphnia taxon composition within reservoirs was very similar in 2004 and 2005 samples (table 1; figure 2). However, we observed a sharp decline in the presence of D. galeata×cucullata hybrids, which were not recorded in 2005 in three out of four reservoirs in which they had been present a year before.
To compare the extent of variability in the Daphnia taxon composition along both the horizontal and vertical gradients with the variation among localities, we decomposed the total compositional variance in an ANOVA style, using partial RDA methods (table 2). In both the cases, differences in the species composition among localities explained about two-thirds of the overall variation. The differences in the Daphnia taxon composition along the longitudinal and vertical directions explained, respectively, 23 and 16% of the variability remaining after correction for differences among reservoirs; the effects of position along the gradient as well as the interactions ‘sampling station–locality’ were always statistically significant. The residual variance, which includes the effect of season, was substantially lower for the longitudinal gradient, showing that the horizontal differentiation in the taxon composition was more stable than the vertical structure at the dam.
Apart from summarizing variability over all taxa and localities, we tested whether individual taxa exhibit significant differences in their relative abundance along both the gradients. In the longitudinal direction, the strongest response was observed for D. galeata×longispina hybrids, occurring significantly more often in the downstream locations (F=9.859, p<0.001); other two significantly responding taxa were D. longispina (F=5.496, p=0.008) and D. cucullata (F=4.803, p=0.013), the former apparently preferring the downstream and the latter the upstream sites. In reservoirs where both of these two parental species were found, they only very rarely co-occurred at the same sampling sites; if so, either one or both of them were in minor proportion (less than 5% of all Daphnia in the sample). Out of the three taxa commonly occurring at the downstream locations (D. galeata, D. longispina and their hybrids), a significant response to vertical differentiation was found for both the parental species: D. galeata typically dominating in the epilimnion and D. longispina being more common in the deeper layers (F=12.284, p<0.001 and F=7.405, p=0.002, respectively).
(c) Hybridization in the Vranov Reservoir
The Vranov Reservoir, 18 km long, up to 45 m deep, and with a theoretical retention time of 134 days, is the second longest locality in our dataset, and is the only one in which we found all the three parental species and the two hybrids. The patterns observed at this locality may be used as a characteristic example of the effect of both longitudinal and vertical environmental gradients on the Daphnia taxon composition, and on the body size and fecundity of the respective taxa. The distribution pattern of all taxa in the two seasons, shown in figure 3, illustrates not only the differing characteristics of taxon spatial distribution but also the temporal dynamics of such patterns. In 2004, the interspecific hybrids occupied an intermediate position between their parental species. However, we observed remarkable changes in the occurrence of hybrids between the two seasons. On the one hand, D. galeata×cucullata, which was common in the upstream and especially central part of the reservoir in 2004 (being actually the most common taxon at the middle sampling station), completely disappeared in 2005; on the other hand, the relative abundance of D. galeata×longispina increased in the second year.
The patterns in Daphnia size and fecundity along the horizontal axis of the Vranov Reservoir (figure 4) reflected the environmental gradients. We observed a gradual decrease in the carapace size and increase in Daphnia brood size towards the upstream region, this pattern being consistent in both years; both the sampling station and the taxon had significant effects on size and fecundity (ANOVA, p<10−5 in all cases). The D. galeata×cucullata hybrids were intermediate and significantly differed from their parental taxa in size (adjusted p<0.002) but not in fecundity. The relationships among D. galeata×longispina hybrids and their parental species were different. The size of the three taxa did not differ significantly, and the patterns of fecundity varied depending on the sampling station and date (figure 4). However, in overall comparison, coexisting D. galeata and D. longispina did not differ in brood size, and hybrids actually had slightly but significantly higher fecundity (adjusted p<0.05) than either of the parental taxa.
(d) Taxon-specific differences in carapace size and fecundity
To check whether such patterns of body size and fecundity are consistent in other localities as well, we analysed differences in carapace and brood size between coexisting and spatially overlapping taxa within the two hybridizing species pairs (D. galeata–D. cucullata and D. galeata–D. longispina), using data from all samples (4544 individuals in total). A comparison of animals from the same sampling sites and dates demonstrated highly significant effects of taxon identity on carapace size (ANOVA, p<10−15) as well as on fecundity (p<10−13).
We observed highly significant differences in each pairwise comparison among D. cucullata, D. galeata and their hybrid: the carapace size increased from D. cucullata to hybrid and to D. galeata, with the corresponding adjusted p-values of 2.64×10−8 and 6.83×10−6, respectively. However, the differences in fecundity of these taxa were not significant. In the second hybridizing pair, the general patterns were different from those recorded in the Vranov Reservoir: the carapace size of D. galeata was significantly larger than in coexisting D. longispina or their hybrids (adjusted p<10−4); the difference between D. longispina and D. galeata×longispina hybrids was not significant (adjusted p=0.12). The fecundity of the three taxa followed the same trend: D. galeata had substantially more eggs/embryos than the two other taxa (adjusted p<2.4×10−4); the difference between D. longispina and D. galeata×longispina hybrids was non-significant (adjusted p=0.10).
The taxon distribution patterns on longitudinal profiles observed in the studied reservoirs in 2004 (Seda et al. 2007b) remained relatively stable over two consecutive seasons, showing that spatially differentiated Daphnia communities re-formed after the winter bottleneck. Hybrids mostly coexisted with the parental taxa, suggesting that they are formed locally within reservoirs, as has been shown for lakes (Spaak 1997). This is also in agreement with the common presence of hybrid eggs in resting egg banks of the studied localities (I. Vaníčková, J. Seda & A. Petrusek 2007, unpublished data). In reservoirs where all the three parental species coexisted, differences in selection forces along horizontal gradients resulted in the almost complete spatial separation of two potentially hybridizing parental species (D. cucullata and D. longispina). This segregation, acting as an ecological reproductive barrier, was the most likely cause of the absence of D. cucullata×longispina hybrids in our samples. In general, the coexistence of these two species, which prefer opposite ends of predation and trophic gradients, is less likely than the coexistence of either of them with D. galeata. This may explain the relative scarcity or absence of D. cucullata×longispina hybrids in many lakes (e.g. Wolf & Mort 1986; Hebert et al. 1989b; Spaak 1997), although locally they may be more common (Gießler 1997b).
The numerical dominance of D. galeata×longispina over D. galeata×cucullata in our samples, however, is unlikely to be explained by the different degrees of coexistence of the two taxa. Both the parental species pairs co-occur, and while species of the former show different preferences on the vertical profile (figure 3; see also Seda et al. 2007b), this certainly does not prevent successful hybridization; on the contrary, D. galeata×longispina hybrids were common and apparently successful in reservoirs. Hybridization between D. galeata and D. cucullata is common as well; in a preliminary analysis of resting egg banks from those reservoirs where all the three parental species coexisted (I. Vaníčková, J. Seda & A. Petrusek 2007, unpublished data), hybrids with D. cucullata were actually found more frequently than those with D. longispina, and D. galeata×cucullata hybrid eggs were common in the resting egg bank of the Kníničky Reservoir, where such hybrids were apparently absent from the active population.
A reduced hatching success of hybrids (Schwenk et al. 2001; Keller et al. 2007) combined with the competition with parental species may partly explain this absence of hybrids in the water column. However, our long-term data on the clonal composition from the Římov Reservoir (J. Seda, J. Macháček, A. Petrusek & K. Schwenk 2007, unpublished data), as well as the temporal pattern of taxon dominance in the Vranov Reservoir (figure 3), suggest another factor contributing to the scarcity of D. galeata×cucullata hybrids in our samples: apparently, they are less likely to overwinter than D. galeata×longispina, and new hybrid genotypes have to be recruited from resting egg banks. The higher sensitivity of D. galeata×cucullata to winter conditions may be due to their smaller size and therefore lower efficiency under low food supply in winter, or because shallower parts of reservoirs, where they occur more often, provide a less suitable environment in winter and the bottlenecks in Daphnia populations are much more severe there than in the deeper parts (J. Seda, J. Macháček, A. Petrusek & K. Schwenk 2007, unpublished data). By contrast, at least some successful D. galeata×longispina clones probably survive in the deep lacustrine parts of reservoirs and contribute to the next season's community; the long-term survival of such hybrid genotypes has been documented in various European lakes (e.g. Spaak & Hoekstra 1993; Jankowski & Straile 2004). However, the apparently strong decline of D. galeata×cucullata hybrids in 2005 may have also been due to within-season fluctuations of taxon abundance (e.g. Keller & Spaak 2004), which would not be detected by our single summer sampling, or due to their distribution being restricted to the non-sampled regions of the studied localities.
It is important to note that the presence of substantial gradients in the Daphnia taxon composition did not depend on the reservoir flushing rate, although we would not expect them to develop in localities with very fast water flow. Gradients corresponding to species' ecological requirements were observed in reservoirs with a short retention time (average value below one month, e.g. Kníničky), in those with intermediate values (four to six months, e.g. Vranov and Vír) as well as in one with a retention time well over 1 year (Želivka; see table 1). Of all factors affecting the spatial heterogeneity of the Daphnia taxonomic composition, we presume the most crucial one to be the intensity of size-selective fish predation and different susceptibilities of different taxa to it (Spaak & Hoekstra 1997; Spaak & Boersma 2006). Fish predation pressure physically limits the size of Daphnia occurring in different parts of the water body, and only when the predator pressure is relaxed under a certain threshold does interspecific competition among larger taxa come into play. Fishes are preferentially present in the upstream reservoir regions (Vašek et al. 2004); this fact is reflected in the small size of Daphnia present there (figure 4). The tendencies in longitudinal arrangement of the taxa from the hybridizing pair D. galeata–cucullata agree well with significant differences in carapace size of the respective taxa, favouring the smallest Daphnia upstream.
The dominance of larger Daphnia further downstream corresponds with their higher filtering efficiency (Gliwicz 2003), small D. cucullata or its hybrids being probably outcompeted by larger taxa. Predation pressure may contribute not only to the absence of D. longispina and its hybrids in the upstream regions of eutrophic reservoirs but also to the vertical distribution of these taxa at the dam. Although the differences in size of coexisting taxa are rather low, hardly giving one taxon a selective advantage against predators, species and hybrids differ in predator susceptibility due to their differing vertical distribution in the stratified water column. Additionally, differences in life-history traits (Macháček & Seda 2007, 2008) certainly contribute to changes in the community and clonal composition; the responses to food quality and quantity nevertheless highly depend on individual genotypes even within taxa (Seidendorf et al. 2007). Fecundity of at least some hybrid clones may certainly exceed their coexisting parental species (figure 4); nevertheless, the brood size variation seems to be locality dependent—in Vranov, D. galeata×longispina hybrids exhibited higher fecundity than parental taxa, but, overall, such a pattern was not observed. This might be due to different genetic backgrounds of hybrids in different reservoirs, or due to their response to temporally or spatially varying environmental conditions (Spaak & Hoekstra 1995).
The existence of longitudinal environmental gradients in reservoirs may improve conditions for the presence of hybrids. Characteristics of interspecific hybrids, such as the intermediate size of D. galeata×cucullata, may promote their intermediate spatial distribution in comparison with the parental species in the way observed in the Vranov Reservoir. The presence of a clear hybrid-dominated zone on the horizontal axis of this reservoir suggests that hybrid dominance may not only be a temporal phenomenon, but also have a specific spatial aspect. We hypothesize that somewhere along reservoir axes, an environment favouring intermediate phenotypes may exist over most of the growing season. Areas of hybrid dominance would certainly be limited to only a certain section of the reservoir longitudinal profile, spatially more restricted than regions where parental taxa possibly coexist. Additionally, we presume that such zones may shift along the reservoir axis, depending on the stability of environmental gradients and the flow rate in the reservoir. Three sampling stations in our project only very crudely covered the longitudinal heterogeneity in reservoirs, so it is possible that denser sampling would reveal more areas where the relative abundance of hybrid genotypes exceed that of the parental species.
Hybrid zones in plankton, if developed, are definitely less stable and a much rarer phenomenon than those in the terrestrial environment. The three-dimensional nature of the pelagic environment allowing relatively free movement of zooplankters and passive dispersal due to wind-induced mixing of water masses and water currents tend to disrupt heterogeneities in the community composition necessary for forming such zones. For taxa entering diapause during their life cycle, such as cladocerans, the process is further complicated by the dispersal of the dormant stages (ephippia in Daphnia), which may be deposited in other parts of the water body than where they are formed; hybrid genotypes may therefore be hatching in conditions substantially different from those where the distribution of their parental species overlap. However, our observations suggest that despite all mixing, strong environmental gradients in reservoirs may exert selection pressures sufficient for hybrid zone formation.
The long-term hybrid presence in reservoirs may have important consequences on the gene flow among parental species during the induction of sexual reproduction, as zones of hybrid dominance may increase the likelihood of the production of later-generation hybrid or backcrossed genotypes. It is not unlikely that the reservoir environment, which offers suitable conditions to species adapted to strikingly different environments, is also suitable for introgression among species to a higher degree than in other, more homogeneous lakes. Although our data are not detailed enough to analyse the extent of introgression, it is worth noting that the disagreement between ITS-RFLP and the allozyme pattern was more often found for individuals from localities with more coexisting taxa, and that the fit between the identification obtained from ITS-RFLP and from other markers was lower in our studied Czech reservoirs than in other localities such as mountain or boreal lakes (Skage et al. 2007).
Patterns of introgression and maintenance of polymorphism of the ITS rDNA region are nevertheless complicated by the multi-copy character of this marker, and its concerted evolution through gene conversion (Arnheim 1983). At least in the Stanovice Reservoir, for which detailed microsatellite data were available (Š. Ruthová, A. Petrusek, J. Seda, A. Thielsch, N. Brede & K. Schwenk 2008, unpublished data), the apparently common ‘hybrid’ patterns detected by ITS-RFLP could not be due to contemporary local hybridization. In this case, the alleged longispina-specific ITS allele must had been introgressed into D. galeata genotypes earlier, and probably persisted in a certain proportion of ITS copies in the genome. Apparently, a detailed comparison of different marker systems, and the use of a wide array of polymorphic markers, preferably microsatellites (Brede et al. 2006), would be needed to reveal the fine-scale patterns of local introgression processes.
This study conformed to all relevant legislation of the Czech Republic.
Our study of reservoir Daphnia was funded by the Czech Science Foundation (project 206/04/0190), Grant Agencies of the Charles University (GAUK 114807) and of the Academy of Sciences of the Czech Republic (AVOZ 60170517), and the Czech Ministry of Education (MSM0021620828 and MSM6007665801). Kateřina Kolářová and Barbora Horová helped during allozyme analyses and Ivana Vaníčková provided results from resting egg banks. Comments from Derek Taylor and one anonymous reviewer helped to improve the manuscript, and David Hardekopf made language corrections. We appreciate also the support of State Water Authority managers during sampling.
One contribution of 16 to a Theme Issue ‘Hybridization in animals: extent, processes and evolutionary impact’.
- © 2008 The Royal Society