Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Who Eats Whom in a Pool? A Comparative Study of Prey Selectivity by Predatory Aquatic Insects

  • Jan Klecka ,

    kleckj01@prf.jcu.cz

    Affiliations Department of Ecosystems Biology, Faculty of Science, University of South Bohemia, České Budějovice, Czech Republic, Laboratory of Theoretical Ecology, Biology Centre of the Academy of Sciences of the Czech Republic, v.v.i, Institute of Entomology, České Budějovice, Czech Republic

  • David S. Boukal

    Affiliations Department of Ecosystems Biology, Faculty of Science, University of South Bohemia, České Budějovice, Czech Republic, Laboratory of Theoretical Ecology, Biology Centre of the Academy of Sciences of the Czech Republic, v.v.i, Institute of Entomology, České Budějovice, Czech Republic

Abstract

Predatory aquatic insects are a diverse group comprising top predators in small fishless water bodies. Knowledge of their diet composition is fragmentary, which hinders the understanding of mechanisms maintaining their high local diversity and of their impacts on local food web structure and dynamics. We conducted multiple-choice predation experiments using nine common species of predatory aquatic insects, including adult and larval Coleoptera, adult Heteroptera and larval Odonata, and complemented them with literature survey of similar experiments. All predators in our experiments fed selectively on the seven prey species offered, and vulnerability to predation varied strongly between the prey. The predators most often preferred dipteran larvae; previous studies further reported preferences for cladocerans. Diet overlaps between all predator pairs and predator overlaps between all prey pairs were non-zero. Modularity analysis separated all primarily nectonic predator and prey species from two groups of large and small benthic predators and their prey. These results, together with limited evidence from the literature, suggest a highly interconnected food web with several modules, in which similarly sized predators from the same microhabitat are likely to compete strongly for resources in the field (observed Pianka’s diet overlap indices >0.85). Our experiments further imply that ontogenetic diet shifts are common in predatory aquatic insects, although we observed higher diet overlaps than previously reported. Hence, individuals may or may not shift between food web modules during ontogeny.

Introduction

Who eats whom and how much? Answering this seemingly simple question is vital for the understanding of processes structuring animal communities. Data on prey selectivity are crucial for mapping the topology of food webs and predicting the effects of species invasions and extinctions on food web structure and stability [1][4]. Data on the diets of different predators are also required to quantify resource partitioning, which can underlie their coexistence [5], [6]. Yet for many food webs, publication of detailed data on the trophic links is sacrificed to achieve more compact description of the often complex food web topology [3], and data coverage varies across habitat types. In freshwater, food webs in standing fishless water bodies have been much less studied than those in streams and lakes (see [7]). Different physical factors and biotic interactions shape the communities in these habitat types, and many species are present in only one of them [8]. Conclusions drawn from the studies of food webs in streams and lakes may thus have only limited applicability to small standing waters without fish. For example, predator-prey body mass ratios differ across habitat types and taxonomic groups of consumers, which may have important implications for food web stability because predator-prey body mass ratios affect interaction strengths [9], [10].

Higher trophic levels in small standing waters are occupied by anurans and aquatic insects (e.g., [11]). The top predators include, at least in the temperate zone, mainly dragonfly and damselfly larvae (Odonata), diving beetles (Coleoptera: Dytiscidae) and bugs (Heteroptera: Nepomorpha). All three groups are speciose and diversified [12] and tens of species can coexist locally (e.g., [11], [13]). They have been traditionally considered as generalist predators [14][16], most likely because of the paucity of experimental data. However, many empirical studies suggest that these predators frequently prefer some prey over others (e.g., [17][20]). Their prey selectivity may lead to cascading effects in the food web [21] and contribute to the maintenance of high levels of biodiversity in standing waters.

Although predatory aquatic insects have been studied for decades, their feeding relationships are surprisingly little known apart from a few model taxa. A synthesis of their prey selectivity is missing and available data need to be described in detail. Hence, it cannot be assessed to what extent the mechanisms of selective predation and resource/habitat niche partitioning promote the diversity of communities in small water bodies, e.g. through food web compartmentalization [22], [23]. Neither do we know which predators have the largest impact on food web structure and which prey are keystone species supporting a disproportionate number of predators in these habitats.

To help answer these questions, we carried out a series of multiple-choice predation experiments with common predator and prey species that often coexist in pools and other small standing water bodies in central Europe. Experiments are the only viable option to compare the diets and prey selectivity across all these predators, as bugs and diving beetle larvae are suctorial. Gut content analyses based on morphological identification of the remains in the gut of the predator are hence applicable only to dragonfly larvae [24], [25] and adult diving beetles. Even when gut contents can be analysed, the estimates of prey selectivity and consumption rates may be severely biased by the fact that different types of food may take very different times to pass through the gut (e.g., [26]). Moreover, neither stable isotope analysis [27][29] nor gut contents can reveal predator selectivity in the absence of detailed data on available prey [18].

In this paper, we summarize our experimental results on selective predation by diving beetles, bugs and odonate larvae together with previously published experiments. We subsequently discuss the importance of diet overlaps, varying vulnerability of prey and ontogenetic diet shifts for the structuring of food webs in small fishless water bodies. The influence of body size and other trophic traits on the strength of predation links will be thoroughly analysed elsewhere (Klecka & Boukal, in prep.).

Methods

Ethics Statement

No specific permits were required for fieldwork as the sampled localities are not protected or privately owned. The use of tadpoles in the experiment was permitted by the regional authority (permit no. KUJCK 12524/2010 OZZL/2/Do) and the Ministry of Education of the Czech Republic (permit no. 7947/2010 30). No permit was needed for the use of invertebrates in the experiment because none of the species is protected.

Laboratory Experiment

We performed multiple choice predation experiments in an experimentally assembled, semi-natural food web with nine regionally common species (13 different stages) of predatory aquatic insects (Table 1) and seven prey species (Table 2). We also used different stages of three predators to study ontogenetic diet shifts. We were not able to cover more species or stages due to limited time available for the experiment, constrained chiefly by the availability of small tadpoles. Acilius and Libellula were chosen because they were among the most abundant species in the field and multiple stages were available simultaneously during the experiment. Dytiscus was used because both larvae and adults are voracious predators [30][32] that may even cause trophic cascades [21]; understanding the differences in their diets could help assess their potentially contrasting impacts on prey populations.

Experiments were carried out in May and June 2010 in a climate room with a regular temperature cycle (day: max. 22°C, night: min. 18°C) and 18 L:6D photoperiod. All animals were collected at various sites in South Bohemia (Czech Republic) and acclimated for 2–5 days prior to experiments. Predators were kept individually in small containers (0.25–0.7 l) and fed daily ad libitum with prey different from those used in the experiments (mainly larvae of Trichoptera). Each predator was starved for 24 hours prior to the experiment to standardize its hunger level. Prey were kept in larger containers (2–20 l) and supplied with abundant natural food (decaying plant material, detritus, algae etc.). Prey individuals which were unused or survived the experiment were released to their natural habitat.

Experiments were performed in translucent whitish plastic boxes (bottom dimensions 24×16 cm) filled with 2.5 litres of tap water (depth ca. 8 cm) aged for one day. The vessels had no substrate on the bottom; four narrow stripes of white plastic mesh suspended vertically in the water column provided simple perching sites. In each replicate, all prey individuals were released first (six Rana tadpoles, six Lymnaea, 10 Chironomus, 10 Cloeon, 10 Culex, 10 Asellus and 30 Daphnia; the densities were within the range of natural densities observed in small pools in the field); the predator was added after several minutes. Each experiment was left undisturbed to run for 24 hours, after which we counted all surviving prey; hence, dead prey were not replaced during the experiment and we did not collect data on the predation sequence.

All individual predators and prey were used only once. Natural mortality of prey, evaluated in four control trials run in the experimental vessels without a predator using the same prey combination and density as in the experiments with predators, was low (Table 2). To account for its potential impact on the results, mean number of dead prey in control trials was subtracted from prey missing at the end of each experiment. No dead uneaten prey was found in the predation experiments suggesting that wasteful killing [33] did not occur.

All predators and 20 randomly chosen individuals of each prey species were preserved in 80% ethanol and their body length excluding appendages was measured to nearest 0.1 mm (Table 1 and 2). We also classified their microhabitat use during the experiments. Almost no individuals of any species used the perching sites except Coenagrion larvae. Only two microhabitats were thus recognized: water column (including perching sites) and bottom. We refer to the second- and third-instar beetle larvae as L2 and L3, respectively. Last-instar larvae of dragonflies and damselflies are referred to as F-0 and larvae of the second instar before the last as F-2. Instar numbers are omitted throughout the text if only one instar was investigated and names are abbreviated to the genus except if multiple species from the same genus are discussed.

Data Analyses

Analyses were carried out and figures made in R 2.11.0 [34] unless stated otherwise. Selectivity of individual predators was evaluated using Manly’s selectivity index α [35], [36]:(1)where ni0 is the initial number of prey items of type i, ri is the number of prey items of type i consumed by the predator and m is the number of prey types used in the experiment. Occasionally the predator consumed all individuals of the most preferred prey. To calculate Manly’s α in these cases, eq. (1) was modified by adding one individual of this prey to the corresponding ni0 and nj0. This assumes that the added individual would have survived, and the corresponding estimate of αi is slightly conservative. Values of αi for individual prey species were compared with values expected for no selectivity using separate t-tests as recommended by Manly [37]. For presentation, values of αi were converted into electivity indices [36]. The indices for individual prey types range from −1 (prey absent in diet) to +1 (prey representing 100% of diet), with a value of 0 corresponding to unselective feeding. Diet breadth of a predator was defined as the number of prey types with electivity index larger than −1 (i.e. it only excluded prey that was never consumed). Numbers of prey consumed by individual predators and corresponding values of Manly’s alpha are shown in Table S1 (Supporting Information).

Pairwise diet overlaps of predators were calculated using Pianka’s overlap index [6] in Ecosim 7.0 [38]:(2)where Oij = Oji is the diet overlap between predator species i and j and symbol αpk denotes value of Manly’s α for prey type k consumed by predator p (p = i or j). Value of Pianka’s overlap index Oij = Oji = 1 means that the diet of the two predators is identical; the lower the value, the less similar their diets. We also modified eq. (2) to calculate overlaps Pij = Pji in predation pressure between prey species i and j by replacing Manly’s α with prey mortality:(3)where the mortality mph of prey type p (p = i or j) consumed by predator k is calculated as the ratio of consumed prey individuals, rpk, over the initial number of prey, np0. We used diet dissimilarity Dij = 1– Oij and a predator dissimilarity index Dij = 1– Pij as input data in Ward’s hierarchical clustering (stats package for R; [34]) and in nonmetric multidimensional scaling (NMDS, vegan 1.17-2 package for R; [39]) to classify and visualize the overlaps.

Finally, we analyzed the modularity of our experimental food web with two trophic levels to complement the analysis of similarities of predators based on their diet overlaps and similarities of prey based on shared predators. The aim of the analysis was to test whether our experimental food web consists of distinct modules characterized by numerous (strong) interactions within modules and few (weak) interactions among modules [22], [23], [40]. Our experiments provided us with quantitative data on the strength of predator-prey interactions, representing a weighted bipartite network. Hence, we used an algorithm for weighted networks [41] implemented in bipartite 1.17 package for R [42] to detect modules. Our experiments yielded two alternative measures of interaction strength, predator selectivity (Manly’s α) and prey mortality. We used both measures to assess the robustness of the results; the original data on the scale between 0 and 1 were multiplied by 100 and rounded to the nearest integer before each analysis.

Review of Published Experiments

To complement our results, we reviewed previous laboratory experiments on prey selectivity with the same three groups of predators and various aquatic invertebrates used as live prey. We first searched Web of Science and Zoological Record using search phrases ‘predation AND taxon’, ‘foraging AND taxon’, ‘predator AND taxon’ and ‘prey AND taxon’, where ‘taxon’ stands for appropriate names of the predators and prey at various taxonomical levels. The results were complemented by a thorough search of the references in the relevant papers and of publications citing these papers. Only studies using more than one prey species for a predator were included, i.e. we omitted studies of stage or size selectivity. We also excluded studies dealing only with vertebrate prey (tadpoles and fish fry) and studies on cannibalism and intraguild predation. For each experiment, we noted the predator and prey taxa and developmental stages, experimental setup, method of data analysis and the main results (preferred prey or lack of selective feeding). We further classified the predators and prey as occupants of water column or the benthic microhabitat as in our experiment; we mostly used known information on their microhabitat use because most experiments did not specify this behaviour. We also noted the use of any habitat structure such as bottom substrate and natural or artificial vegetation to assess the impact of refuges and perching sites on the results.

To test which prey (at the level of order) are more preferred, we pooled all published experiments except those with prey from a single order and ranked each prey from the most to least preferred. The matrix of the prey ranks in each individual experiments based on the entire dataset had 85% empty cells as most experiments used only 2–3 prey types, thereby precluding the use of the method of analysis of incomplete ranking data as described in [43]. A necessary condition for a meaningful analysis is at least ∼50% non-empty cells [43], which we could achieve only by restricting the dataset to studies involving only Cladocera, Diptera and/or Ephemeroptera. However, these studies clearly showed that Diptera were more preferred than Ephemeroptera and slightly less preferred than Cladocera, making the analysis redundant. We therefore simply scored the preference for each taxon in each experiment on binary scale (1 = most preferred prey and 0 = all other prey in a given experiment) and compared the probability of being the most preferred prey taxon, using a generalized linear model with quasi-binomial distribution. This allowed us to compare the preferences across all prey taxa. We further used multiple comparisons of means for generalized linear models in multcomp package for R [44] to perform post-hoc pairwise comparisons of preference between different prey taxa.

To test for microhabitat association between predators and their preferred prey, we counted experiments that identified one or multiple prey from a single microhabitat (benthic/water column) as the most preferred and had at least one non-preferred prey from the other habitat. We used the resulting 2×2 contingency table to test the microhabitat association with a one-tailed Fisher’s exact test.

Results

Laboratory Experiment

We focus on the following six aspects of our experimentally assembled food web: selectivity of predators, diet overlaps of different predators, ontogenetic diet shifts, prey vulnerability, predator overlaps of different prey and food web modularity.

All predator species fed selectively but differed in their level of specialization (Figure 1). Adult diving beetles (Acilius and Hydaticus) and adult Notonecta were most selective, having only one preferred prey and at most one neutrally selected prey type (i.e., consumed proportionally to its abundance). Two more predators, Libellula F-2 and Acilius L2 larvae, were fairly specialized with one preferred prey and two prey with neutral preference. Most other predators (Acilius L3 larvae, Dytiscus and Ilyocoris adults and Coenagrion, Sympetrum and Libellula F-0 larvae) significantly preferred two prey and had neutral preference to one more prey type. Dytiscus and Anax larvae were least selective. Dytiscus larvae strongly preferred and nearly depleted three prey types (Asellus, Chironomus and Rana), while Anax larvae consumed five out of seven prey species at least proportionally to their abundance and significantly preferred two of them (Chironomus and Culex). Diet breadth was related to but not identical with the preference patterns. It ranged from all seven prey in Anax larvae to four prey in adult Ilyocoris and in Coenagrion larvae. Diet breadth of the other predators was five or six prey types.

thumbnail
Figure 1. Prey selectivity of predatory aquatic insects.

Mean values ± SE of electivity index are plotted. Positive values indicate preferred prey. Prey species with electivity values significantly different from zero (P<0.05, Holm’s correction of P-levels within each predator species was used) are marked by asterisk. Predator stage as in Table 1. Panels are sorted taxonomically: A–F = Coleoptera, G–H = Heteroptera and I–M = Odonata.

https://doi.org/10.1371/journal.pone.0037741.g001

Consequently, diet overlaps of predator pairs varied greatly between 0.09 (Acilius L2 and adult Hydaticus) and 0.99 (Sympetrum and Libellula F-0; Table 3). Cluster analysis suggested four predator groups (Figure 2) with strong pairwise overlaps (≥0.80 except between Coenagrion and Acilius L3) within each group (Table 3). The first group comprises medium-sized benthic predators with strong preference for Chironomus and avoidance of Asellus (adult Acilius, adult Hydaticus and Libellula F-2 larvae; diet overlap 0.89–0.98). The second group consists of larger benthic odonate larvae which fed mainly on both species of dipteran larvae and neutrally selected Asellus (Anax, Sympetrum and Libellula F-0 larvae; diet overlap 0.89–0.99). Another group of large-bodied predators foraging on the bottom contains larvae and adults of Dytiscus and adult Ilyocoris (diet overlap 0.82–0.95). They all consumed large numbers of Asellus and Chironomus, although Ilyocoris did not feed on tadpoles. The fourth group (Coenagrion larvae, L2 and L3 larvae of Acilius and adult Notonecta; diet overlap 0.66–0.97) foraged mainly in the water column and near the surface and fed mainly on Culex or, as in Coenagrion larvae, on Daphnia and Culex. We call these predators nektonic hereafter.

thumbnail
Figure 2. Similarity of diets of the predators used in the experiment.

A. Ward’s hierarchical clustering of the diet dissimilarities Dij; height = value of clustering criterion for the particular cluster. B. Diet dissimilarities Dij visualised by nonmetric multidimensional scaling (NMDS, stress = 3.99), with groups identified by the cluster analysis highlighted; position of each species in the diagram corresponds to the centre of its label.

https://doi.org/10.1371/journal.pone.0037741.g002

We tested for an ontogenetic diet shift (ODS) in two diving beetles (L2 and L3 larvae and adults of Acilius, L3 larvae and adults Dytiscus) and one dragonfly (F-2 and F-0 larvae of Libellula). In all three species, diets differed significantly between the stages (Table 4 and Figure 1). Significant diet shifts occurred mainly in the preferred prey. We observed strong ODS in Acilius from Culex (preferred by L2 larvae) to Cloeon and Culex (preferred by L3 larvae) and subsequently to Chironomus (preferred by adults). Moreover, adult Acilius also fed on Rana tadpoles, which were never eaten by the larvae. Diet overlap was therefore much lower between the adults and larvae (0.41 and 0.55) than between the two larval instars (0.89; Table 3). Diet overlaps within the other two species were high (0.91 and 0.92) and the resulting ODS mainly quantitative (Figure 1).

thumbnail
Table 4. Tests of ontogenetic diet shifts (pairwise comparisons based on t-test).

https://doi.org/10.1371/journal.pone.0037741.t004

Vulnerability of each prey species to predation differed significantly across all predators (proportion of prey individuals consumed during experiment; GLM with quasi-binomial distribution, P<0.0001 in all cases; Figure 3). Chironomus and Culex were most vulnerable overall (average mortality 51% and 46%, respectively). Either of them was the most preferred prey for each predator (Figure 1) except Dytiscus larvae. Three other species were highly vulnerable only to a subset of predators (mortalities of Rana tadpoles: 100% from Dytiscus larvae, 54% from Anax larvae; Asellus: 92% from Dytiscus larvae, 68% from Ilyocoris; Cloeon: 73% from Acilius L3 larvae). The least consumed prey was Lymnaea, eaten only by Anax larvae and adult Dytiscus (ca. 10% mortality from either predator; Figure 3).

thumbnail
Figure 3. Mortality of individual prey species subjected to different predators.

Prey mortality is expressed as proportion of individuals eaten during the experiment (mean ± SE). Dotted vertical lines represent the overall observed mortality averaged across all combinations of prey and predator species. Predator stage as in Table 1; predators ordered taxonomically as in Table 1 and Figure 1.

https://doi.org/10.1371/journal.pone.0037741.g003

Predator overlaps among prey species varied between 0.13 (Daphnia and Lymnaea) and 0.83 (Asellus and Rana; Table 5). Cluster analysis identified three prey groups (Figure 4): the largely invulnerable Lymnaea, one group of larger benthic prey with mostly shared predators (Rana, Asellus and Chironomus; predator overlap 0.67–0.83) and another group of smaller, non-benthic prey (Cloeon, Culex and Daphnia; predator overlap 0.39–0.78). Pairwise overlaps in predators within the prey groups were thus on average lower than diet overlaps within the predator groups.

thumbnail
Figure 4. Similarity of predation pressure among prey species used in the experiment.

A. Ward’s hierarchical clustering of the dissimilarities of predation pressure Dij; height = value of clustering criterion for the particular cluster. B. Predation pressure dissimilarities Dij visualised by nonmetric multidimensional scaling (NMDS, stress = 2.29), with groups identified by the cluster analysis highlighted; position of each species in the diagram corresponds to the centre of its label.

https://doi.org/10.1371/journal.pone.0037741.g004

Finally, modularity analysis identified three modules in our experimental food web. A nektonic module containing four predators (Notonecta, L2 and L3, Acilius larvae and Coenagrion) and three prey (Daphnia, Cloeon and Culex) is identical to the combination of the respective predator and prey groups identified by cluster analysis (Figures 2 and 4). The other two modules involve benthic prey and predators (Figure 5). Both measures of interaction strength yielded the same results (Figure 5A and 5B), suggesting that the conclusions are robust. In addition, predators in the two “benthic” modules correspond well to the results of the cluster analysis, which subdivided one of the modules into two clusters and otherwise assigned only one species (Anax) differently (Figures 2 and 5). The benthic prey modules differ from the results of cluster analysis only by isolating the strongly linked Chironomus rather than excluding the weakly linked Lymnaea from the remaining three benthic prey (Figures 4 and 5).

thumbnail
Figure 5. Modularity of the experimental food web.

Modules identified by modularity analysis displayed as boxes; symbol size corresponds to predator-prey interaction strength. A. Predator preference (Manly’s α) used as measure of interaction strength. B. Prey mortality used as measure of interaction strength.

https://doi.org/10.1371/journal.pone.0037741.g005

Review of Published Experiments

Thirty-five studies reporting 59 experiments with more than 40 predator species satisfy the predefined criteria (Supporting Information Table S2). Similar numbers of experiments used diving beetles (n = 19), bugs (n = 22) and odonate larvae (n = 18) but one model taxon prevails in each group: Dytiscus in diving beetles, Notonecta in bugs and Anax in dragonflies. The prey included more than 50 species ranging from protozoans to amphibians. Surprisingly rare were studies in which the prey included Isopoda (four experiments), Oligochaeta (one) and studies comparing invertebrate and vertebrate prey (five). Prey composition varied greatly among studies and rarely was diverse enough to represent a semi-natural mixture. Up to 21 (mostly 2–6) prey types were offered; more than six prey types usually involved experiments with some prey types representing multiple size classes of the same species.

thumbnail
Table 5. Pairwise overlaps of predator assemblages associated with different prey species.

https://doi.org/10.1371/journal.pone.0037741.t005

Almost all papers reported distinct selectivity of the predator towards some of the prey (Supporting Information Table S2). Taken together, they reveal large and significant differences between preferences for different prey taxa (GLM with quasi-binomial distribution, F = 5.01, P = 0.00003; Figure 6). Cladocera were most preferred in 21 out of 27 experiments in which they were used together with alternative prey from a different order; most of these experiments used Daphnia (preferred in 19 out of 22 experiments). Various dipterans were also frequently preferred (19 out of 29 experiments); larvae of Culicidae were favoured in 12 out of 28 experiments and Chironomidae in 10 out of 13 experiments (in six of these cases, both families were tested together). Ten other prey taxa were preferred in at least one experiment, with Copepoda (two of 10 cases), Rotifera (one of four cases) and Ephemeroptera (two of 12 cases) among the least preferred. Five taxa were never preferred: Heteroptera (n = 7 experiments), Ostracoda (n = 6), Odonata (n = 4), Hydrachnida (n = 1) and Turbellaria (n = 1). Post-hoc pairwise comparisons of the simplified data on predator preferences suggested that (i) the preference for the three taxa never preferred in multiple experiments is significantly lower (t-test, P<0.05) than for the remaining taxa and (ii) the preference for Cladocera, Trichoptera and Diptera is significantly higher than that for Copepoda, Rotifera and Ephemeroptera (Figure 6). Due to the small sample sizes, we could not compare if the preferences differed among predator taxa or if a pre-existing bias in the selection of prey combinations affected the results.

thumbnail
Figure 6. Relative vulnerability of prey taxa most often used in previous experiments.

Based on data in Supporting Information Table S1. Preference = proportion of experiments in which the prey taxon was most preferred (mean ± SE). Number of experiments that included a given prey taxon is stated above each error bar. Dotted line separates prey types used in >3 experiments (to the left) and <3 experiments (to the right). Horizontal grey lines denote prey groups that do not differ significantly according to multiple comparisons of means (P>0.05).

https://doi.org/10.1371/journal.pone.0037741.g006

Among 33 experiments with benthic predators, benthic prey was preferred over prey in the water column in 10 experiments and the reverse was found in 4 experiments; the remaining experiments used only prey occupying a single microhabitat or the experiment found conditional preference for either type of prey. The proportion of clear outcomes was even lower among the 24 experiments with predators foraging in the water column, with one and three experiments respectively reporting preference for benthic and nektonic prey. This is suggestive of an overall preference for prey from the same microhabitat, but the results are inconclusive (P = 0.14, one-tailed Fisher’s exact test).

Taken together, the results of previous experiments show that larvae of dragonflies and damselflies, aquatic bugs and diving beetles are selective predators. In general, they seem to feed most heavily on cladocerans and on the larvae of Diptera (mosquitoes and chironomids). The preferences are at least partly driven by overlapping microhabitat use: predators preferring cladocerans and mosquito larvae (e.g., Notonecta) usually forage in the water column, while chironomid larvae are generally preferred by benthic predators (Supporting Information Table S2).

Discussion

Diet overlaps and Coexistence of Multiple Predators

Predators that coexist in the same habitat need to occupy different space, time and/or food niche [6], [45]. Resource partitioning and restricted diet overlap leading to distinct food niches supposedly drive predator coexistence in lizards (e.g., [6], [45]), fish [46], [47] and carnivorous mammals [48], [49], although some recent studies question the importance of food niche partitioning (e.g., [50]). Data on invertebrates are scarcer and more controversial (e.g., [51][54]) and lacking for aquatic insects in standing waters. Although the link between niche separation and coexistence should be rigorously tested [55], measuring diet overlaps among co-occurring species provides important insights into the potential contribution of resource partitioning to long-term coexistence.

Our experimental food web consisted of three modules determined by individual body size and microhabitat use. Nektonic predators that utilized mainly the water column (Acilius and Coenagrion larvae and Notonecta) and fed mainly on nektonic prey (Daphnia and Culex) had limited diet overlap (0.09–0.82, mean 0.47) with bottom-foraging “benthic” predators (adult diving beetles, Ilyocoris and Libellula and Dytiscus larvae; Figure 2). On the contrary, predators sharing the same microhabitat had a strong, typically >0.80 overlap in diets. This result likely holds across taxa. For example, diet preferences of other nektonic backswimmers (genera Anisops and Buenoa) are similar to Notonecta, while the more sedentary bugs from the families Pleidae, Belostomatidae and Naucoridae seem to prefer benthic prey (Supporting Information Table S2).

A formal test supporting the idea of microhabitat-use driven modularity by the literature data was inconclusive, mainly because many previous experiments, especially with nektonic predators, used prey from only one microhabitat. Among the four nektonic predators with a clear preference, only larvae of Acilius preferred benthic prey [56]. This could have arisen if the individual predators perceived the benthic habitat and water column as one habitat, e.g. in shallow water, and chose prey according to some other criteria. We could not assess this phenomenon due to paucity of data; the impact of water depth on diet overlaps between nektonic and benthic predators deserves further study.

Furthermore, large predators are probably less constrained by microhabitats than small-bodied predators, and their diet is driven primarily by high metabolic demands and the need for high energy intake rates. Feeding links of such predators may thus provide connections between separate food web modules. In our experiments, Anax larvae were the least selective and consumed all prey species as the only predator. Overall, larvae of diving beetles (Dytiscus) and dragonflies (Anax) are known as voracious predators of tadpoles [30][32], [57][59] and other large prey including smaller conspecifics and intraguild prey [60][63]. This effect may not be universal: some Dytiscus species have a specialized diet, such as large caddisfly larvae [30], [58].

Prey selectivity of smaller predators has been less studied (Supporting Information Table S2). It is incompletely understood apart from damselflies, which are known to feed mostly on zooplankton and are thus linked within a nektonic module ([20], [64]; references in Supporting Information Table S2). Only few studies focused on the diets and prey selectivity of medium-sized dragonflies, which are among classic taxa used in various ecological experiments (Supporting Information Table S2). The two species in our study, Libellula and Sympetrum, fed mainly on smaller benthic prey (Chironomus) and mosquito larvae, which is in line with previous results [56].

In other words, similarly sized predators foraging in the same microhabitat will often–but not always–have strongly overlapping diets. These predators are “generalists in a narrow sense”, i.e. have a broad diet conditional on their size and foraging microhabitat, and can be subject to intense indirect competitive interactions with potentially fundamental consequences for the entire food web structure. For example, food depletion in predatory aquatic insects may increase cannibalism and intraguild predation [60], which are both common among larvae of diving beetles [63] and odonates [60][62]. Intraguild predation of dytiscid larvae by odonates may cause negative correlations between odonate and dytiscid densities [65]. However, indirect competitive interactions are difficult to measure in the field and have been reported by very few studies on predatory aquatic insects. Data in [66] suggest that adult diving beetles are not food limited and hence protected from exploitative competition. That study found significant density-dependent mortality, possibly caused by competition for food or cannibalism, only in the larvae.

Diet overlaps have not been formally calculated in predatory aquatic insects before. Values found in our experiment mostly fall within the range known in other taxa, although Pianka’s indices of 0.66–0.99 (mostly≥0.87) within each of the predator clusters identified in our experiment are unusually high. Such nearly complete diet overlaps are uncommon in both terrestrial and aquatic vertebrate predators (e.g., [48], [67][70]), and overlap index as low as 0.76 has been implicated in species replacement driven by food competition [71]. Coexistence of predators with so strongly overlapping diets requires additional mechanisms such as exploitation of different size classes or stages of the shared prey [72], [73].

The observed overlaps would likely decrease with a broader range of prey species that are unfeasible to test in the laboratory, but we believe that the decrease would be limited given the broad coverage of prey sizes and functional types. Moreover, we have tested only one size class for each prey and thus cannot establish if size selectivity or other factors–such as differences in the diet concerning prey not included in the experiments, different time and/or space niches, and strong intraspecific competition or cannibalism–help these predators coexist in the same habitat. The importance of apparent competition and interference mediated by overlapping diets for population dynamics of predatory aquatic insects thus requires further study.

Coexistence on local spatial scales, e.g. through diversification of diets, might also lead to coexistence at larger scales. We speculate that this mechanism could contribute to high regional diversity of diving beetles relative to the other predators (e.g., Czech Republic: ∼120 species of diving beetles, ∼60 species of odonates and ∼40 species of aquatic bugs; [74], [75]). For example, we found that adult Hydaticus consumed mainly Chironomus larvae, while experiments with other similarly sized genera and different prey sets found preference for Daphnia [19]. Larvae of Agabus and Acilius are efficient predators of mosquito larvae [76], [77], which is consistent with our results on Acilius canaliculatus. Some diving beetles even prefer dead prey [78]. They may be opportunistic scavengers that exploit yet another food niche.

Ontogenetic Diet Shifts

The concept of ontogenetic diet shift (ODS) is rapidly becoming a central theme in studies of aquatic food webs. Most animals grow substantially during their development and body size is now recognized as a key driver of predator-prey relationships, particularly in the aquatic environment [2], [79]. Predation pressure on prey assemblages can thus change considerably as the predators grow, which can have both ecological and evolutionary consequences [80]. In addition, diet shifts may release individual predators from intraspecific competition for food in the same way as reduced diet overlaps decrease apparent competition between species. ODSs are well documented in various holometabolous taxa with complex life histories [81]. Shifting diets may also reflect changes in foraging (micro)habitat and behaviour. They occur in Notonecta bugs [82][84] and in larval odonates in both running [24] and standing waters [18]. Odonate larvae begin to feed on rotifers and even protozoans after hatching [85][87] and later switch to larger benthic prey.

We detected more or less pronounced ODSs in all three predators for which we tested more than one life stage (Acilius, Dytiscus and Libellula), indicating that ODSs are widespread in predatory aquatic insects. ODSs should be particularly common in larvae and adults of diving beetles, which differ greatly in morphology, foraging habitat and behaviour. Surprisingly, observed diet overlaps between larvae and adults (∼0.4–0.5 in Acilius and ∼0.9 in Dytiscus) were quite high. Diet overlaps between last-instar larvae and adults calculated as in [81] (0.42 in Acilius and 0.29 in Dytiscus) greatly exceed the previously reported range (0–0.08) for metamorphic species and are closer to the values typically ascribed to continually growing, gape-limited predators [81]. This suggests that limited diet overlaps may require non-overlapping habitats (as in odonates) or the presence of other mechanisms absent in the diving beetles. Diet shifts occurred also between consecutive larval instars (Acilius and Libellula), even if they were smaller and the overlaps in diet (∼0.9) fell within those observed in the four predator clusters.

We conclude that ODSs in predatory insects in small fishless pools sometimes allow species to move between food web modules during ontogeny (larvae and adults of Acilius). However, ODSs may not always be strong enough to release successive instars/stages from intraspecific competition and the predator may remain in one food web module throughout most of its aquatic life (as in Dytiscus and larvae of Libellula). Further experiments and observations are needed to quantify the overlaps across multiple instars/stages for multiple predators and resolve this issue.

Vulnerability of Prey Species to Predation

Last but not least, we showed that vulnerability to predation varies greatly among common prey found in fishless pools and other standing water bodies of the temperate zone in Europe. The most vulnerable prey were dipteran larvae (Culex and Chironomus), which were also preferred by predators in most other experiments. In addition, the literature review suggested that cladocerans, especially Daphnia, are the most preferred prey of predatory aquatic insects along with trichopteran and dipteran larvae. In our experiment, cladocerans were rarely preferred. This discrepancy most likely reflects different prey composition in the experiments. Previous studies often coupled Daphnia with other zooplankton and thus probably offered a suboptimal prey choice to many predators. Neither have we tested some of the prey reported as non-preferred in previous experiments (Copepoda and Ostracoda; probably invulnerable or non-profitable to most predatory aquatic insects). In addition, some predator-prey combinations were rare in previous experiments. Suggested overall ranking and differences in preferences between the prey (Figure 6) must be therefore taken cautiously as many experiments probably pre-selected prey taxa to test a priori hypotheses.

Prey vulnerability is partly determined by its ability to withstand or avoid predator attacks. Culex larvae are capable of rapid escape movements in the water column, but these were apparently not effective against most predators in our experiment as the predation rates were high. High vulnerability of Culex larvae to various invertebrate predators was repeatedly confirmed in experiments aimed to identify potential anti-malaria control agents [88]. Given their frequently high population densities (e.g., [76]), mosquito larvae possibly serve as keystone prey species that support a large number of different predator species and contribute to the maintenance of high species diversity in small fishless standing water bodies.

Chironomus larvae were exposed and lacked refuge in our experiment, and hence represented highly profitable and easily accessible prey (see Supporting Information File S1 for full discussion of the methodological issues associated with most multiple-choice predation experiments). Cothran & Thorp [89] showed that the presence of a refuge can strongly decrease predation on chironomid larvae as the predator’s selectivity shifts towards alternative prey. Under natural conditions, chironomid larvae burrow in soft sediments to avoid predators and are vulnerable only when migrating [90]. At the same time, predators can successfully specialize on chironomids [90]. Our experimental setup thus corresponds to the presence of a large chironomid population with sufficiently many larvae available to predators. Alternatively, high preference for chironomid larvae in an experiment lacking refuges can indicate preference for benthic prey with high probability of successful attack and high profitability (such as any injured/diseased animals), or tendency towards facultative or obligatory scavenging, which could be the preferred feeding strategy for some predators [78].

On the other hand, prey vulnerability is greatly reduced by reaching a size refuge and/or mechanical defences (such as external hard shell), although some predators may adapt their foraging strategy to overcome the defences (e.g., Dytiscus preying on caddisflies; [58]). In our experiment, Lymnaea snails were the only protected and also the least vulnerable prey. In addition, we used relatively large snails that have apparently reached a size refuge, as small snails are vulnerable to predation [91]. Size refuge, along with species-specific diets or different size and identity of the prey assemblages, could also explain the observed lack of preference of adult Dytiscus marginalis for snails in our experiment, contrary to Dytiscus alascanus [21]. Similarly, Rana tadpoles lack mechanical defence but were the second least consumed prey overall, apparently as they were too large and difficult to handle for most predators except the largest ones (Anax and Dytiscus; compare [31], [57], [59]).

Conclusions

Predatory aquatic insects of standing waters were often seen as generalists in the past [14][16]. This traditional view must be revised: the unfolding story on predator-prey interactions in standing fishless waters is one of complex, challenging patterns. Diets of predatory insects in these habitats vary from highly specialized to broadly general, but hardly any species appear to feed indiscriminately. By combining a simple experiment with a literature survey, we provide a basis for future studies on food webs involving predatory aquatic insects in small standing water bodies.

We found a highly interconnected experimental food web, separable into several modules based on microhabitat use (bottom or water column) and body size of the predators and their prey. We thus suggest that predatory aquatic insects in small standing water bodies are “generalists in a narrow sense”: species with similar size foraging in the same microhabitat have widely overlapping diets. Moreover, ontogenetic diet shifts associated with individual growth in size and changes in foraging microhabitats seem common across all major groups of predatory aquatic insects. That is, food web interactions of predatory aquatic insects might be equally affected by intraspecific and interspecific differences.

All these results have potentially crucial implications for the structuring and stability of food webs. We thus call for further predation experiments using artificially assembled food webs, in combination with other methods such as gut content analyses wherever applicable. Particularly needed are more systematic studies of ontogenetic niche shifts, studies of unlikely and random predator-prey combinations, which could resolve the otherwise overlooked issue of avoided prey or reveal unexpected predator-prey links, and comparative studies of prey preferences in the presence and absence of (semi-realistic) habitat structure. Only such pluralistic approach can map the structure of food webs in small standing waters.

Supporting Information

Table S1.

Summary of the results of our experiments. Initial number of prey, number of prey consumed and Manly’s alpha is shown for all experiments. The results of each experiment, identified by the number in “Predator ID” column, are presented in seven rows, one for each prey.

https://doi.org/10.1371/journal.pone.0037741.s001

(XLS)

Table S2.

Summary of published experiments on selective predation by diving beetles, water bugs and larvae of odonates. Only studies with more than one prey species for a predator are included, and studies which focused primarily on vertebrate prey and studies on cannibalism and intraguild predation are excluded. Each row corresponds to one experiment or a set of experiments with the same set of species, possibly with different size classes/developmental stages or densities. Stage: A = adults, L = larvae (number of instars or size classes given in parentheses). Prey No. = total number of prey types (species and size classes) used (superscripts: a = multiple experiments with prey in different combinations; b = multiple experiments with the same prey combinations but different prey abundances). Multiple choice = more prey species were offered simultaneously; yes/no means that some experiments were conducted with individual prey types separately (typically to study single prey functional responses) and other experiments with a mixture of prey. Test of selectivity: Roger = Roger’s index (Lundkvist et al., 2003) and mortality = comparison of prey mortalities or numbers of prey consumed by the predator. The most preferred prey is highlighted in bold. When the most preferred prey differ between treatments (in experiments which offered prey in different combinations/abundances), all of them are highlighted. Bottom structure presence: yes/no = manipulated presence/absence. Preferred microhabitat of predators and prey was classified into two categories: water column and benthic. References appearing only in Table S2 are listed in File S1.

https://doi.org/10.1371/journal.pone.0037741.s002

(XLS)

File S1.

Multiple-choice experiments: theoretical background. Summary of main conceptual issues related to our and previous experiments (including references) and references appearing only in Table S2.

https://doi.org/10.1371/journal.pone.0037741.s003

(DOC)

Acknowledgments

The structure and presentation of the paper benefited greatly from comments by three anonymous reviewers.

Author Contributions

Conceived and designed the experiments: JK DSB. Performed the experiments: JK. Analyzed the data: JK. Contributed reagents/materials/analysis tools: JK DSB. Wrote the paper: JK DSB. Compiled the literature survey: JK DSB.

References

  1. 1. Woodward G, Hildrew AG (2001) Invasion of a stream food web by a new top predator. J Anim Ecol 70: 273–288.
  2. 2. Woodward G, Speirs DC, Hildrew AG (2005) Quantification and resolution of a complex, size-structured food web. Adv Ecol Res 36: 85–135.
  3. 3. Ings TC, Montoya JM, Bascompte J, Blüthgen N, Brown L, et al. (2009) Ecological networks – beyond food webs. J Anim Ecol 78: 253–269.
  4. 4. Rudolf VHW, Lafferty KD (2011) Stage structure alters how complexity affects stability of ecological networks. Ecol Lett 14: 75–79.
  5. 5. Schoener TW (1971) Theory of feeding strategies. Annu Rev Ecol Syst 2: 369–404.
  6. 6. Pianka ER (1973) The structure of lizard communities. Annu Rev Ecol Syst 4: 53–74.
  7. 7. Brose U, Cushing L, Berlow EL, Jonsson T, Banasek-Richter C, et al. (2005) Body sizes of consumers and their resources. Ecology 86: 2545.
  8. 8. Wellborn GA, Skelly DK, Werner EE (1996) Mechanisms Creating Community Structure Across a Freshwater Habitat Gradient. Annu Rev Ecol Syst 27: 337–363.
  9. 9. Naisbit RE, Kehrli P, Rohr RP, Bersier L-F (2011) Phylogenetic signal in predator – prey body-size relationships. Ecology 92: 2183–2189.
  10. 10. Riede JO, Brose U, Ebenman B, Jacob U, Thompson R, et al. (2011) Stepping in Elton’s footprints: a general scaling model for body masses and trophic levels across ecosystems. Ecol Lett 14: 169–178.
  11. 11. Warren PH (1989) Spatial and Temporal Variation in the Structure of a Freshwater Food Web. Oikos 55: 299.
  12. 12. Hutchinson GE (1993) A Treatise on Limnology. Vol. IV. The Zoobenthos. New York: John Wiley and Sons. 944 p.
  13. 13. Klecka J, Boukal DS (2011) Lazy ecologist’s guide to water beetle diversity: Which sampling methods are the best? Ecol Indic 11: 500–508.
  14. 14. Batzer DP, Wissinger SA (1996) Ecology of insect communities in nontidal wetlands. Annu Rev Entomol 41: 75–100.
  15. 15. Cummins K (1973) Trophic relations of aquatic insects. Annu Rev Entomol 18: 183–206.
  16. 16. Bay E (1974) Predator-prey relationships among aquatic insects. Annu Rev Entomol 19: 441–453.
  17. 17. Kehl S, Dettner K (2003) Predation by pioneer water beetles (Coleoptera, Dytiscidae) from sandpit ponds, based on crop-content analysis and laboratory experiments. Arch Hydrobiol 158: 109–126.
  18. 18. Blois C (1985) The larval diet of three anisopteran (Odonata) species. Freshw Biol 15: 505–514.
  19. 19. Lundkvist E, Landin J, Jackson M, Svensson C (2003) Diving beetles (Dytiscidae) as predators of mosquito larvae (Culicidae) in field experiments and in laboratory tests of prey preference. Bull Entomol Res 93: 219–226.
  20. 20. Cooper SD, Smith DW, Bence JR (1985) Prey selection by freshwater predators with different foraging strategies. Can J Fish Aquat Sci 42: 1720–1732.
  21. 21. Cobbaert D, Bayley SE, Greter JL (2010) Effects of a top invertebrate predator (Dytiscus alaskanus; Coleoptera: Dytiscidae) on fishless pond ecosystems. Hydrobiologia 644: 103–114.
  22. 22. Allesina S, Pascual M (2008) Network structure, predator–prey modules, and stability in large food webs. Theor Ecol 1: 55–64.
  23. 23. Stouffer DB, Bascompte J (2011) Compartmentalization increases food-web persistence. Proc Natl Acad Sci USA 108: 3648–3652.
  24. 24. Woodward G, Hildrew AG (2002) Body-size determinants of niche overlap and intraguild predation within a complex food web. J Anim Ecol 71: 1063–1074.
  25. 25. Cohen JE, Jonsson T, Carpenter SR (2003) Ecological community description using the food web, species abundance, and body size. Proc Natl Acad Sci USA 100: 1781–1786.
  26. 26. Martin ID, Mackay RJ (1982) Interpreting the diet of Rhyacophila larvae (Trichoptera) from gut analyses: an evaluation of techniques. Can J Zool 60: 783–789.
  27. 27. Grey J, Jones RI, Sleep D (2001) Seasonal changes in the importance of the source of organic matter to the diet of zooplankton in Loch Ness, as indicated by stable isotope analysis. Limnol Oceanogr 46: 505–513.
  28. 28. Jennings S, Pinnegar JK, Polunin NVC, Boon TW (2001) Weak cross-species relationships between body size and trophic level belie powerful size-based trophic structuring in fish communities. J Anim Ecol 70: 934–944.
  29. 29. Pace M, Carpenter S, Cole J, Coloso J, Kitchell J, et al. (2007) Does terrestrial organic carbon subsidize the planktonic food web in a clear-water lake? Limnol Oceanogr 52: 2177–2189.
  30. 30. Leclair RJ, Alarie Y, Burassa J-P (1986) Prey choice in larval Dytiscus harrisii Kirby and D. verticalis Say (Coleoptera, Dytiscidae). Entomologica Basiliensa 11: 669–689.
  31. 31. Rubbo MJ, Mirza RS, Belden LK, Falkenbach JJ, Storrs SI, et al. (2006) Evaluating a predator-prey interaction in the field: the interaction between beetle larvae (predator) and tadpoles (prey). J Zool 269: 1–5.
  32. 32. Blunck H (1923) Die Entwicklung des Dytiscus marginalis L. vom Ei bis zur Imago. 2. Teil. Die Metamorphose (B. Das Larven- und das Puppenleben). Zeitschrift für Wissenschaftliche Zoologie 121: 171–391.
  33. 33. Johnson DM, Akre BG, Crowley PH (1975) Modeling arthropod predation: wasteful killing by damselfly naiads. Ecology 56: 1081–1093.
  34. 34. R Core Development Team (2010) R: A language and environment for statistical computing. 1: The R project for statistical computing website Available: http://www.r-project.org. Accessed 2011 November.
  35. 35. Manly B (1974) A model for certain types of selection experiments. Biometrics 30: 281–294.
  36. 36. Chesson J (1983) The estimation and analysis of preference and its relationship to foraging models. Ecology 64: 1297–1304.
  37. 37. Manly BFJ (1995) Measuring selectivity from multiple choice feeding-preference experiments. Biometrics 51: 709–715.
  38. 38. Gotelli NJ, Entsminger G (2001) EcoSim: Null models software for ecology. Version 7. EcoSim website. 1: Available: http://garyentsminger.com/ecosim/index.htm. Accessed 2011 November.
  39. 39. Oksanen J, Blanchet FG, Kindt R, Legendre P, O’Hara R, et al. (2010) vegan: Community Ecology Package. R package version 1.17–2. 1: Available: http://cran.r-project.org/package=vegan. Accessed 2011 November.
  40. 40. Krause AE, Frank KA, Mason DM, Ulanowicz RE, Taylor WW (2003) Compartments revealed in food-web structure. Nature 426: 282–285.
  41. 41. Newman MEJ (2004) Analysis of weighted networks. Phys Rev E 70: 056131.
  42. 42. Dormann CF, Fründ J, Blüthgen N, Gruber B (2009) Indices, graphs and null models: analyzing bipartite ecological networks. Open Ecol J 2: 7–24.
  43. 43. Lim DH, Wolfe DA (2001) An efficient alternative to average ranks for testing with incomplete ranking data. Biom J 43: 187–206.
  44. 44. Hothorn T, Bretz F, Westfall P (2008) Simultaneous inference in general parametric models. Biom J 50: 346–363.
  45. 45. Pianka ER (1974) Niche overlap and diffuse competition. Proc Natl Acad Sci USA 71: 2141–2145.
  46. 46. Kuhara N, Nakano S, Miyasaka H (1999) Interspecific competition between two stream insect grazers mediated by non-feeding predatory fish. Oikos 87: 27–35.
  47. 47. Pimentel CR, Joyeux J-C (2010) Diet and food partitioning between juveniles of mutton Lutjanus analis, dog Lutjanus jocu and lane Lutjanus synagris snappers (Perciformes: Lutjanidae) in a mangrove-fringed estuarine environment. J Fish Biol 76: 2299–2317.
  48. 48. Andheria AP, Karanth KU, Kumar NS (2007) Diet and prey profiles of three sympatric large carnivores in Bandipur Tiger Reserve, India. J Zool 273: 169–175.
  49. 49. Murdoch JD, Munkhzul T, Buyandelger S, Reading RP, Sillero-Zubiri C (2010) Seasonal food habits of corsac and red foxes in Mongolia and the potential for competition. Mammalian Biology - Zeitschrift fur Saugetierkunde 75: 36–44.
  50. 50. Luiselli L (2008) Do lizard communities partition the trophic niche? A worldwide meta-analysis using null models. Oikos 117: 321–330.
  51. 51. Heinrich B (1976) Resource partitioning among some eusocial insects: bumblebees. Ecology 57: 874–889.
  52. 52. Wertheim B, Sevenster JG, Eijs IEM, Van Alphen JJM (2000) Species diversity in a mycophagous insect community: the case of spatial aggregation vs. resource partitioning. J Anim Ecol 69: 335–351.
  53. 53. Takahashi KH, Tuno N, Kagaya T (2005) The relative importance of spatial aggregation and resource partitioning on the coexistence of mycophagous insects. Oikos 109: 125–134.
  54. 54. Ibarra-Nunez G, Garcia J, Lopez J, Lachaud J (2001) Prey analysis in the diet of some ponerine ants (Hymenoptera : Formicidae) and web-building spiders (Araneae) in coffee plantations in Chiapas, Mexico. Sociobiology 37: 723–755.
  55. 55. Siepielski AM, McPeek MA (2010) On the evidence for species coexistence: a critique of the coexistence program. Ecology 91: 3153–3164.
  56. 56. Febria CM, Magnusson AK, Williams DD (2005) Seasonal abundance and prey selection of the nymphs of three sympatric species of Sympetrum (Odonata: Libellulidae) in an intermittent pond. Can Entomol 137: 723–727.
  57. 57. Henrikson B-I (1990) Predation on amphibian eggs and tadpoles by common predators in acidified lakes. Ecography 13: 201–206.
  58. 58. Johansson A, Nilsson AN (1992) Dytiscus latissimus and D. circumcinctus (Coleoptera, Dytiscidae) larvae as predators on three case-making caddis larvae. Hydrobiologia 248: 201–213.
  59. 59. Relyea RA (2001) The relationship between predation risk and antipredator responses in Larval Anurans. Ecology 82: 541.
  60. 60. Johansson F (1993) Diel feeding behavior in larvae of four odonate species. J Insect Behav 6: 253–264.
  61. 61. McPeek MA (1990) Determination of species composition in the Enallagma damselfly assemblages of permanent lakes. Ecology 71: 83–98.
  62. 62. Johansson F (1991) Foraging modes in an assemblage of odonate larvae - effects of prey and interference. Hydrobiologia 209: 79–87.
  63. 63. Yee DA (2010) Behavior and aquatic plants as factors affecting predation by three species of larval predaceous diving beetles (Coleoptera: Dytiscidae). Hydrobiologia 637: 33–43.
  64. 64. Thompson DJ (1978) Prey size selection by larvae of the damselfly, Ischnura elegans (Odonata). J Anim Ecol 47: 769.
  65. 65. Larson DJ (1990) Odonate predation as a factor influencing dytiscid beetle distribution and community structure. Questiones Entomologicae 26: 151–162.
  66. 66. Juliano SA, Lawton JH (1990) Extrinsic vs. intrinsic food shortage and the strength of feeding links: effects of density and food availability on feeding rate of Hyphydrus ovatus. Oecologia 83: 535–540.
  67. 67. Pianka ER, Pianka HD (1976) Comparative ecology of twelve species of nocturnal lizards (Gekkonidae) in the western Australian desert. Copeia. 1976. : 125–142.
  68. 68. Pianka ER, Huey RB (1978) Comparative ecology, resource utilization and niche segregation among gekkonid lizards in the southern Kalahari. Copeia. 1978. : 691–701.
  69. 69. Luz-Agostinho KDG, Agostinho AA, Gomes LC, Júlio HF (2008) Influence of flood pulses on diet composition and trophic relationships among piscivorous fish in the upper Paraná River floodplain. Hydrobiologia 607: 187–198.
  70. 70. Spitz J, Cherel Y, Bertin S, Kiszka J, Dewez A, et al. (2011) Prey preferences among the community of deep-diving odontocetes from the Bay of Biscay, Northeast Atlantic. Deep Sea Res Part I Oceanogr Res Pap 58: 273–282.
  71. 71. Hamer T, Hays D, Senger C, Forsman E (2001) Diets of northern barred owls and northern spotted owls in an area of sympatry. J Raptor Res 35: 221–227.
  72. 72. Wathne J, Haug T, Lydersen C (2000) Prey preference and niche overlap of ringed seals Phoca hispida and harp seals P-groenlandica in the Barents Sea. Mar Ecol Prog Ser 194: 233–239.
  73. 73. de Roos A, Schellekens T, van Kooten T, Persson L (2008) Stage-specific predator species help each other to persist while competing for a single prey. Proc Natl Acad Sci USA 105: 13930–13935.
  74. 74. Rozkošný R (1980) Klíč vodních larev hmyzu [Key to the larvae of aquatic insects]. Prague: Academia. 521 p.
  75. 75. Boukal DS, Boukal M, Fikáček M, Hájek J, Klečka J, et al. (2007) Katalog vodních brouků České republiky/Catalogue of water beetles of the Czech Republic. Klapalekiana 43 Suppl. pp. 1–289.
  76. 76. Chandra G, Mandal SK, Ghosh AK, Das D, Banerjee SS, et al. (2008) Biocontrol of larval mosquitoes by Acilius sulcatus (Coleoptera: Dytiscidae). BMC Infect Dis 8: 138.
  77. 77. Culler LE, Lamp WO (2009) Selective predation by larval Agabus (Coleoptera: Dytiscidae) on mosquitoes: support for conservation-based mosquito suppression in constructed wetlands. Freshw Biol 54: 2003–2014.
  78. 78. Velasco J, Millán V (1998) Feeding habits of two large insects from a desert stream: Abedus herberti (Hemiptera: Belostomatidae) and Thermonectus marmoratus (Coleoptera: Dytiscidae). Aquat Insects 20: 85–96.
  79. 79. Ings TC, Montoya JM, Bascompte J, Blüthgen N, Brown L, et al. (2009) Ecological networks – beyond food webs. J Anim Ecol. pp. 253–269.
  80. 80. Urban MC (2007) The growth-predation risk trade-off under a growing gape-limited predation threat. Ecology 88: 2587–2597.
  81. 81. Werner EE, Gilliam JF (1984) The ontogenetic niche and species interactions in size-structured populations. Annu Rev Ecol Syst 15: 393–425.
  82. 82. Scott MA, Murdoch WW (1983) Selective predation by the backswimmer, Notonecta. Limnol Oceanogr 28: 352–366.
  83. 83. Murdoch WW, Scott MA, Ebsworth P (1984) Effects of the general predator, Notonecta (Hemiptera) upon a freshwater community. J Anim Ecol 53: 791–808.
  84. 84. Orr B, Murdoch W, Bence J (1990) Population regulation, convergence, and cannibalism in Notonecta (Hemiptera). Ecology 71: 68–82.
  85. 85. Hampton SE, Gilbert JJ (2001) Observations of insect predation on rotifers. Hydrobiologia. pp. 115–121.
  86. 86. Walsh EJ, Salazar M, Remirez J, Moldes O, Wallace RL (2006) Predation by invertebrate predators on the colonial rotifer Sinantherina socialis. Invertebr Biol 125: 325–335.
  87. 87. Lawton JH (1970) Feeding and food energy assimilation in larvae of the damselfly Pyrrhosoma nymphula (Sulz.) (Odonata: Zygoptera). J Anim Ecol 39: 669–689.
  88. 88. Shaalan EAS, Canyon DV (2009) Aquatic insect predators and mosquito control. Trop Biomed 26: 223–261.
  89. 89. Cothran ML, Thorp JH (1985) Tests of preference and switching behavior of the dragonfly Celithemis fasciata. Oikos 44: 350–355.
  90. 90. Pajunen V (1983) Prey selection by larvae of Potamonectes griseostriatus (Degeer)(Coleoptera, Dytiscidae). Ann Zool Fenn 20: 31–35.
  91. 91. Turner AM, Chislock MF (2007) Dragonfly predators influence biomass and density of pond snails. Oecologia 153: 407–415.