Pressure-Induced Shifts in Trophic Linkages in a Simplified Aquatic Food Web

Introduction

Aquatic ecosystems harbor a diverse array of organisms that interact across trophic levels (Wallace et al., 1997). As species in natural communities are strongly interdependent (Pimm et al., 1991), effects of a pressure on a given species can be expected to resonate through the entire food web (Gessner et al., 2010; Vörösmarty et al., 2010). Hence, to generate realistic predictions on how existing and emerging anthropogenic stressors can affect major ecosystem properties and impede ecosystem functioning, an understanding of their effects on trophic linkages is required (Moore et al., 2005; Jabiol et al., 2013).

A major pressure on aquatic communities is the use of agricultural chemicals (nutrients and pesticides) which are used to increase and protect crop yield (Helweg et al., 2003). Aquatic animals can react to such pollutants by altering their food source, leading to changes in predation efficiency (Smith and Weis, 1997; Schulz and Dabrowski, 2001), bioturbation (Hinkle-Conn et al., 1998; Hunting et al., 2013) and nutrient cycling (Fleeger et al., 2003). Such alterations in trophic linkages can be regarded as rippling effects in food webs which may result in ecosystem-wide effects and subsequently affect ecosystem functioning and the services it provides (Fleeger et al., 2003; Corvalan et al., 2005). However, how anthropogenic pressures propagate through trophic linkages in aquatic food webs and inherent effects on ecosystem processes remains poorly understood.

To date, the vast majority of studies focusing on the effect of pesticides have evaluated the effect of a single chemical substance on survival of individuals, as this is the common approach in ecotoxicological risk assessment (van den Brink, 2006; Gessner and Tlili, 2016). Most ecosystems however, particularly aquatic ecosystems bordering agricultural land, are exposed to a myriad of interacting stressors (Ormerod et al., 2010; Gagic et al., 2017), making it important to predict their direct and indirect effects on networks of trophic linkages. Little doubt exists that the effects of these substances (e.g., of bioaccumulation and biomagnification, effects on trophic linkages) are not necessarily additive (Côté et al., 2016), which is traditionally assumed by models that assess effects of mixtures of contaminants (Altenburger et al., 2000; Backhaus et al., 2004). Another glaring issue is that most traditional studies to date have typically ignored sublethal effects such as changes in resource preference that further mediate the effects of abiotic factor on species interactions and identities (Rohr et al., 2006; van den Brink, 2006). The traditional approach thus presents a considerable oversimplification of real environmental conditions (Cairns, 1983; Kimball and Levin, 1985; Calow and Forbes, 2003; Rohr et al., 2006; van den Brink, 2006; Vijver et al., 2017), thereby warranting experimental studies that assess the effect of multiple stressors on trophic linkages in aquatic food webs.

Pesticides and herbicides affect different components of aquatic ecosystems, in which herbicides are more likely to indirectly affect herbivorous organisms, and insecticides to directly affect (sensitive) invertebrates. It is thus conceivable that herbicides and pesticides differentially affect the structure and trophic linkages of aquatic food webs. This study therefore aims to evaluate whether (1) environmentally relevant concentrations of key agricultural pesticides and nutrients, and interactions between these, induce shifts in trophic linkages, and (2) observed changes in trophic linkages are qualitatively different between the green (algal-based) and brown (detritus-based) part of the food web. To this end, we carried out a full factorial outdoor mesocosm experiment to investigate the interactive effects of three chemical stressors (herbicide, insecticide, and excess nutrients) at environmentally realistic concentrations on the trophic linkages within a simplified, model aquatic food web.

Methods

To adequately separate the effect of agricultural pressures we used a full factorial design containing an aquatic community composed of commonly co-occurring species, in which we mimicked all combinations of three different realistic pressures: of a single herbicide and a single insecticide, in presence and absence of agricultural nutrients and the interactions between all three factors. We applied the following eight different treatments in a randomized fashion, each with a sample size of six: (I) imidacloprid, (II) terbuthylazine, (III) nutrients, (IV) imidacloprid and terbuthylazine, (V) imidacloprid and nutrients, (VI) terbuthylazine and nutrients, (VII) imidacloprid, terbuthylazine and nutrients, and (VIII) a control with no added substances. The insecticide and herbicide of choice were imidacloprid (CAS No. 138261-41-3, purity 99.7%) and terbuthylazine (CAS No. 5915-41-3, purity 99.4%) respectively, purchased from Sigma Aldrich (Zwijndrecht, The Netherlands). These chemical substances were selected for the following reasons: First, they are representative for a group of insecticides (nitroguanidines, which affect the nervous system) and herbicides (triazines, which are photosystem II -inhibitors), which are commonly used groups of pesticides. Second, both types of pesticides are used across a range of intensively managed agricultural systems and present in high concentrations in surface water in the Netherlands (Pesticide Atlas, 2015).

Experimental Setup

The experiment was set up in the experimental garden of the Hortus Botanicus in Leiden, The Netherlands and was conducted in 48 65 L polyethylene tubs, hereafter called mesocosms. To mimic natural pond conditions and to prevent the mesocosms from unnatural heating, tubs were placed 30 cm into the soil. The 8 cm of the mesocosms that was raised above the ground was covered in reflective tape in order to reflect sunlight to prevent heating. Additionally, mesocosms were tightly covered with 50% shadow cloth nets to prevent colonization, escape and excessive heating. In order to provide the animals with oxygen and to homogenize the water column, each mesocosm was aerated using a 5 watt aquarium pump.

In order to set up an adequate sediment layer, a combination of finely ground hay and fine-grained, ignited quartz sand as mineral substrate (grain size: 0.1–0.5 mm) was used. First, 500 g of finely cut dry hay (5 cm) was pre-soaked and mixed with 12.5 kg moist sand, which was placed at the bottom of the mesocosms. After placement of the sediment layer, ~36 L of copper-free tap water was added. Copper removal was conducted by filtering tap water with active carbon for 14 days. Subsequently, each mesocosm received 4l of filtered ditch water (mesh width 0.5 mm) that was mixed first in a large volume and was added to inoculate the mesocosms with a natural species composition of unicellular organisms. The water column and sediment layer were left for 1 week in order to allow the microbial community to establish.

The food web composition of the species in the experiment consisted of 10 commonly co-occurring species in natural and agricultural freshwater ecosystems in the Netherlands (Ieromina et al., 2016) of a variety of feeding types; (I) grazer, (II) shredder, (III) filter-collector, (IV) predator, and (V) gatherer/collector (Ieromina et al., 2016). Species were assigned to trophic groups based on literature. This resulted in the following three trophic groups: (1) herbivorous diet: Cyclops sp., Cloeon dipterum and Daphnia magna (Ahlgren et al., 1990; Santer and Van Den Bosch, 1994) (2) Detritivore diet: Tubifex sp., Chironomus riparius, Culex pipiens, and Asellus aquaticus (Graça et al., 1993; David et al., 2000), (3) Mixed diet: Sigara striata, Lymnaea stagnalis (Klecka and Boukal, 2012; Groendahl and Fink, 2016), (4) Predator group consisting of a single species (Notonecta glauca). Notonecta species are generalist predators that can potentially have large impacts on pond community structure (Cockrell, 1984; Wilbur and Dec, 1997; Klecka and Boukal, 2012) and are known to prey on all the species in the introduced community except for Tubifex sp. and only very rarely on L. stagnalis.

Species were either grown in lab culture or were collected from non-polluted field sites. To be able to trace the change in stable isotope levels during the experiment we added individuals of most species as juveniles. A full description of species and their numbers is presented in Table 1. Seven days after setting up the sediments, non-predacious animals were placed in the mesocosms. One day later, we added the predators. Herbicide, insecticide and fertilizer treatments were installed 1 day after predator additions.

TABLE 1

Table 1. Characteristics of species and numbers used in the mesocosm experiment.

The following nominal concentrations were used: 3.5 μg/L terbuthylazine and 4.0 μg/L imidacloprid, which were based on realistic concentrations in Dutch surface waters (Pesticide Atlas, 2015). The substances were dissolved in demineralized water in glass bottles and mixed by magnetic stirring. In addition, we added 6.16 mL of soluble plant feed in nutrient treatments (Nitrogen: Phosphorus: Potassium 7:4:7 (Pokon Naturado), in order to reach nominal concentrations of soluble phosphorus of 3.8 mg/L, which are realistic for Dutch surface waters and were previously shown to lead to high algal production (Ieromina et al., 2016). A second addition of insecticides, herbicides and nutrients, in equal concentrations to the first application, was conducted halfway through the experiment (t = 18, t = time in days), thus mimicking a biweekly application of substances. In order to further promote homogenization of the water column (next to continuous aerating of the water), all substances were added to the mesocosms in large volumes. Presence of either terbuthylazine (Omegam Laboratoria BV (Amsterdam, The Netherlands) and imidacloprid (Wageningen University, Wageningen, The Netherlands) was measured using liquid chromatography-tandem mass spectrometry (LC-MS/MS) on a weekly basis. The measured concentrations showed a gradual decrease in pesticide concentrations, confirming expectations (Electronic Appendix Figure S1). The time weighted average for terbuthylazine was 2.04 μg l⁻¹ and for imidacloprid 2.9μg l⁻¹. Also for nutrients we found that actual concentrations were lower than nominal concentrations, indicating uptake of nutrients.

The experiment was harvested after 35 days. To retrieve all macrofauna, the entire mesocosm was filtered through two sieves, one for the larger specimen (mesh size 2 mm) and one for the smaller specimen (mesh size 0.106 mm). All animals except for D. magna and Cyclops sp. were collected by hand with forceps and pipettes. Individuals of a single species were pooled per mesocosm in a single sample tube. In addition, a 20 × 10 cm plot of periphyton was scraped from the side of the mesocosm interior that received the most sun hours per day. Algal quantity in the mesocosm was insufficient for stable isotope analyses.

Stable Isotope Analyses

To investigate effects of the three different pressures on food web properties, we used carbon and nitrogen stable isotopes. Stable isotopes of carbon are excellent tools to study such changes in food web organization because there is a difference between the δ¹³C/δ¹²C of C3 detritus and freshwater algae (Hamilton and Lewis, 1992). Changes in δ¹³C/δ¹²C thus reflect changes in the utilization of different carbon sources by the respective trophic groups and taxa. Shifts from algal feeding to detritus feeding are expected to result in an increase in δC¹³ levels. Heavy nitrogen isotopes (δ¹⁵N) accumulate in the trophic chain (Hobson and Welch, 1992; Post, 2002), and can thus inform about effects of treatments on trophic structure and position. A change in δ¹⁵N can expected to happen when there is a shift in N source to the base of the food web, for example when algae start using different N sources or as a result of shift in carnivore feeding preference, which could for instance result from weakened prey/intermediate trophic levels.

After the harvesting the experiment, all sampling tubes were stored at −20°C prior to further processing. For the largest species (N. glauca and L. stagnalis), we used only part of the body for stable isotope analysis. For N. glauca, we used the pronotum, which was separated from the remainder of the body manually, and for L. stagnalis, we used only the foot of the snail. We expected that, as a result of growth, these parts would most readily reflect treatment-induced changes in carbon and nitrogen isotopes. Samples were freeze-dried, ground and weighed in tin cups prior to stable isotope analysis. Large samples of N. glauca, L. stagnalis, and S. striata were ground for 5 min. Smaller samples of C. dipterum, A. aquaticus, and periphyton were ground for 1.5 min, using a ball mill (Retsch Technology GmbH, Haan, Germany). Samples of the remaining species were small enough to be used entirely. We used 1.0 (±0.2) mg of each sample of a single species from a single mesocosm, which were weighed into 5 × 8 mm tin capsules. The stable isotopes were measured using an elemental analyser (NC2500; ThermoQuest Italia, Rodana, Italy) coupled with an isotope ratio mass spectrometer (Delta Plus; Thermo-Quest Finnigan, Bremen, Germany) at the VU University Amsterdam. For calibration, USGS 40 and USGS 41 were used. The reproducibility of the δ¹³C and δ¹⁵N analysis determined by repeated analysis of an internal standard (Bovine liver, NIST 1577c) was within 0·15 ‰ (n = 3). Stable isotope abundance was expressed using the δ notation:

$\begin{array}{l}\text{δ}X(‰)=R_{\text{sample}}-R_{\text{standard}}/R_{\text{standard}}\times 1000,\hfill \end{array}$

where X represents the ¹³C content and R_sample and R_standard denote the δ¹³C/δ¹²C ratios of the sample and the standard Vienna PeeDee Belemnite (VPDB), respectively. For N, X represents the δ¹⁵N content and R_sample and R_standard denote the δ¹⁵N/δ¹⁴N ratios of the sample and the standard air, respectively.

Statistics

Statistical analyses were performed with Rstudio, version 0.99.902 (R version 3.3.0) and Statistica 10 (Stat Soft Inc., Tulsa, USA). To illustrate food web relationships between single species, linear least square regressions were carried out for both C and N isotopes. A multiple regression analysis of the nitrogen isotopes was carried out to detect the likely food sources of A. aquaticus in the various treatments (species within 3‰ N from A. aquaticus: S. striata, C. dipterum, periphyton). For the analysis of changes in stable isotope composition within green and brown parts of the food web, estimations for species were averaged within trophic groups, by taking the average of the species in each of the trophic groups. We tested the single and combined effects of the agricultural chemicals on nitrogen and carbon isotopes of each species and of the entire trophic group by means of factorial ANOVA (type-II sum of squares). The factors tested were imidacloprid (2 levels; yes/no), terbuthylazine (2 levels; yes/no), nutrients (2 levels; yes/no) and all possible combinations of these factors. Stable isotopes of the different species/trophic groups were always (and only) used as response variables in these analyses. Assumptions for homogeneity of variances as well as normality of residuals were met.

Results

Effects on Trophic Structure

All species survived until the end of the experiment, with the exception of dipteran larvae C. pipiens, and C. riparius, most likely as a result of predation by N. glauca (Klecka and Boukal, 2012). These are therefore excluded from the remainder of the results and discussion.

As expected, we found that, across all treatments, δ¹³C/δ¹²C ratios differ markedly between species, thus indicating different food preference between species (Figure 1). Similarly, we find different δ¹⁵N/ δ¹⁴N ratios between species, indicating different trophic levels (Figure 1). For both N and C stable isotopes, positive correlations were observed between a number of species (C. dipterum, D. magna, Cyclops sp., L. stagnalis, and periphyton) (Tables 2, 3), suggesting a similar reaction to changing environmental conditions or trophic interactions. In contrast, a negative correlation was found between S. striata and N. glauca, for both nitrogen and carbon isotopes (Tables 2, 3), thus indicating possible competition between these species.

FIGURE 1

Figure 1. Eight panels showing the most important changes in food webs. (A–D) show effects of the nutrient and herbicide treatment, and interaction between these the on herbivore part of the food web. (E–H) Show food webs in nutrient and insecticide treatment and interaction between these, and the effect on the detritivore part of the food web. Green circles in (A–D) indicate the herbivore energy channel; brown circles in (E–H) indicate the detritivore energy channel. Arrows indicate significant effects of treatments on the herbivore and detritivore channel.

TABLE 2

Table 2. Pairwise correlation coefficients between nitrogen isotopes of all species in the experiment.

TABLE 3

Table 3. Pairwise correlation coefficients between carbon isotopes of all species in the experiment.

Effect of Treatments on Base of Food Web

The nutrients that were added to the treatments had a δ¹⁵N/δ¹⁴N ratio that was >3‰ lower than the periphyton, indicating that the periphyton used this food source. It was not possible to determine the stable carbon isotope ratio of added inorganic nutrients, because they contained only trace amounts of carbon. Periphyton δ¹⁵N/δ¹⁴N was positively affected by insecticide additions [F_{(1, 26)}4.3; P = 0.05], and herbicide additions [F_{(1, 26)}9.2; P = 0.005] and was negatively affected by nutrient additions [F_{(1, 26)}8.5; P = 0.007]. Carbon isotope ratios for periphyton were positively affected by the nutrient treatment [F_{(1, 26)} 9.7, P = 0.004].

Effects of Treatments on the Green Web

When averaged over all species, herbivores exhibited a decrease in δ¹⁵N/δ¹⁴N ratio when exposed to nutrients [F_{(1, 40)}25.6, P < 0.0001; Figure 2A]. The other treatments (herbicides, insecticides) had no significant main effect on the δ¹⁵N/δ¹⁴N ratio of herbivores. Carbon isotopes for the herbivores showed a significant main effect of nutrients [F_{(1, 40)} 16.4 P < 0.001] and a significant interaction effect of herbicides and nutrients [F_{(1, 40)} 7.6, P = 0.008; Figure 2B]. Insecticides had a marginally significant positive effect on the carbon isotope ratio of the herbivore channel [F_{(1, 40)} 3.2, P = 0.08]. Herbivore species exhibited changes in nitrogen isotopes as a result of herbicide and nutrient additions. For Cyclops sp. and D. magna, nutrient additions lead to significantly decreased δ¹⁵N/δ¹⁴N ratios [Cyclops: F_{(1, 38)} 10.0; P = 0.003; D. magna: F_{(1, 28)} 5.2, P = 0.03]. Carbon isotope ratios of herbivore species also changed in response to the various treatments. For Cyclops sp, we found a significantly negative effect of nutrients on δ¹³C/δ¹²C ratio [F_{(1, 38)}9.8, P = 0.003] and an interaction effect between herbicides and nutrients [F_{(1, 38)} 9.8, P = 0.004], where the effect of nutrients was larger when herbicides were added [F_{(1, 38)} 4.0, P = 0.05]. For D. magna, no effect of any single stressor was found, only an interaction effect between herbicides and nutrients [F_{(1, 27)} 5.2, P = 0.03], as well as an interaction effect between nutrients, herbicides and insecticides [F_{(1, 27)} 6.1, P = 0.02].

FIGURE 2

Figure 2. Figure showing the significant treatment effect on the green and brown parts of the food web: upper two (A,B) showing effects of treatments on trophic level (δ¹⁵N/¹⁴N) and lower two (C,D) illustrate the effect on food source (δ¹³C/¹²C) for green and brown.

Effects on the Brown Web, N Isotopes

The δ¹⁵N/¹⁴N ratio of the entire group of detritivores group was significantly affected by the insecticide * nutrients treatment [F_{(1, 40)}5.6290, P = 0.02; Figure 2C]; where insecticides increase the δ¹⁵N/δ¹⁴N ratio, but only in absence of nutrients. Carbon isotopes in the detritivore group were only affected by insecticide presence, where δ¹³C/δ¹²C levels decreased when insecticides were present [F_{(1, 35)} 4.7, P = 0.03; Figure 2D].

Within the detritivore guild, nitrogen isotope ratios were only significantly affected by the treatments for L. stagnalis, whereas A. aquaticus and Tubifex sp. showed no significant response. A significant, negative main effect of nutrients on the δ¹⁵N/δ¹⁴N ratio of L. stagnalis was found [F_{(1, 40)}20.3, P < 0.0001] whereas insecticide additions had a positive effect [F_{(1, 40)}3.5820; P = 0.05] on δ¹⁵N/δ¹⁴N-levels of L. stagnalis. Trophic position was also affected by the interaction between herbicides, insecticides and nutrients [F_{(1, 40)}8.1, P = 0.007].

Effects on the Brown Web, C Isotopes

Overall, detritivore species had 3.3‰ higher δ¹³C/δ¹²C ratios than herbivore species (Figure 1), thus confirming that food preferences were different between these functional groups. However, this was not the case for A. aquaticus; carbon isotopes for this species were significantly lower and the δ¹⁵N/ δ¹⁴ N ratio of A. aquaticus was almost an entire trophic level [2.6–3.1‰; one trophic level generally amounts to 3.4‰ (Post, 2002)] higher than the other primary consumers (Figure 1). This indicates that this species may have consumed food sources other than the hay debris. A multiple regression analysis of the nitrogen isotopes of the likely food sources of A. aquaticus (species within 3‰ from A. aquaticus: S. striata, C. dipterum, periphyton) indicated that C. dipterum may have composed a likely food source for this species [Electronic Appendix, Table 1; F_(1.21) 9.4 P 0.02]. A pairwise correlation between A. aquaticus and C. dipterum showed a significant positive relationship between C. dipterum and A. aquaticus for both nitrogen isotopes and carbon isotopes (Figures 3A,B), thus suggesting that A. aquaticus might have had a diet partly composed of dead individuals of C. dipterum.

FIGURE 3

Figure 3. Linear least square regression of nitrogen (A) and carbon (B) isotopes between A. aquaticus and C. dipterum. Line fit, R² and significance level are depicted in both panels. Dotted line indicates hypothetical 1:1 relationship between A. aquaticus and C. dipterum.

Mixed Feeders and Predators

We found no significant effect of nutrient and pesticide treatment on the average δ¹⁵N/δ¹⁴ N level of mixed feeders. Within this guild, nutrients had a negative effect on carbon isotope ratio of C. dipterum [F_{(1, 18)}8.9; P = 0.008], whereas herbicides had a positive effect [F_{(1, 18)}5.0; P = 0.04]. Moreover, we found an interaction effect between herbicides and nutrients [F_{(1, 18)} 6.7, P = 0.02], where the effect of nutrients was larger when herbicides were added. Nutrient additions led to a significantly decreased δ¹⁵N/δ¹⁴N ratio in C. dipterum [F_{(1, 18)}10.084; P = 0.005], but we found no effect of insecticide or herbicide additions on nitrogen isotope ratio.

As indicated by the high δ¹⁵N/ δ¹⁴ N levels, N. glauca was the top predator (Figure 1). We found no significant effect of any of the treatments on both Hemiptera species (S. striata, N. glauca) (Tables 1, 2).

Discussion

By studying the single and joint effects of stressors at environmentally realistic concentrations to an aquatic food web, we show that a number of key species in the community can alter their trophic position and food source as a result of the different stressors. Herbivorous species reacted mainly to the nutrients and herbicides and the synergistic interaction between these, while species in the detritivore part of the food web were observed to be affected by insecticide applications and interactions with nutrients. This suggests that there are qualitatively different effects of stressors on the different functional parts (green and brown) of the food web. This may have strong implications for our understanding of the effects of pesticides and nutrients on key ecosystem processes (Gessner and Tlili, 2016). However, the observed strong interactions between different pesticide classes and nutrients for the different parts of the aquatic food web also present us with a puzzle: we have no clear explanations for the second and third-order effects that go beyond mere speculation. These effects observed in a quasi-natural setting differ markedly and qualitatively from laboratory incubations which show isolated effects on individual species (Salminen et al., 1996; Song et al., 1997; Stoughton et al., 2008), and indicate that effects on species can feed through to higher levels of organization in aquatic ecosystems, and likely inherent ecosystem processes, as described in Mills et al. (1993) and Relyea and Hoverman (2006).

The observed shifts in trophic position and food source suggest that effects of pesticides can cause a rippling effect in the trophic linkages of the food web, depending whether or not excess nutrients are present. In this study, this effect was most pronounced in A. aquaticus and L. stagnalis. Our results show that the trophic level of A. aquaticus was considerably higher than expected based on literature (Graça et al., 1993). While the other detritivores (C. riparius, Tubifex sp.) had δ¹⁵N/δ ¹⁴N-levels of around 4.5-5, δ¹⁵N/δ ¹⁴N-levels of A. aquaticus were almost a full trophic level higher. We suggest that this unexpected finding may be explained by the shift in diet from detritus-feeding to feeding on animal remains. Indeed, C. dipterum (and S. striata) exhibited a great population decline in all mesocosms. Similarly, we find that the trophic position of L. stagnalis depended strongly on input of insecticides. Presence of insecticides leads to an increase in trophic level for this species, which also indicates a shift in diet toward feeding on animal remains. The observed changes in food source and trophic level of A. aquaticus and L. stagnalis suggest a more general ecosystem phenomenon: that pesticides, through their effects on trophic linkages (Moore et al., 2005; Jabiol et al., 2013), can induce food web wide effects by rippling through the food web. Moreover, as these animals represent a significant part of temperate freshwater food webs, this indicates a potential for trophic transfer of toxicants (e.g., bioaccumulation, biomagnification), which in turn might affect species at higher trophic levels.

This study used a simplified community and setting and therefore it remains questionable whether the observed patterns also hold true for natural environments. The mesocosm approach that we chose for this study lacks a number of crucial elements of real ecosystems, most notably, it lacks actual sediments, has relatively low species densities and lacks dispersal: influx and outflux of species (Rohr et al., 2006; van den Brink, 2006). Species replacement is a powerful way to deal with changes in pressures, and could potentially dampen some of the observed effects. Not accounting for this complexity may have resulted in an overestimated toxicity in our setup, yet the emerging pattern that combined toxicity can differentially affect trophic linkages is likely equally relevant in natural environments, suggesting that a consideration of ecosystem complexity is warranted to obtain a more reliable understanding of how stressors can affect species and propagate through the environment.

Conclusion

This study shows that agricultural pressures and interactions between these pressures can induce shifts in trophic linkages in a simplified, model aquatic food web, and highlights that interactions between stressors can variably affect trophic linkages and can cascade toward ecosystem attributes that are essentially not targeted by the chemical pressures used in this study. These types of antagonistic and synergistic interactions are increasingly found to be potentially important explanations for the apparently unpredictable responses of ecosystem properties and processes to anthropogenic pressures, and have important implications for our understanding of how chemical stressors can propagate through food webs (e.g., bio-accumulation and bio-magnification) with inherent risks for their functioning and the services they provide.

Author Contribution

The ideas for this study were conceived by MS, PvB, MV, SB and EH; MS, SB and RvL carried out data collection; data analysis was performed by RvL, SB, MS, and EH, and manuscript writing was carried out by MS, SB, and EH with structural comments from MV and PvB.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer, VG, and handling Editor declared their shared affiliation.

Acknowledgments

We are grateful to the Hortus Botanicus in Leiden for providing the space and help to run this experiment, most notably Prof Paul Keβler and Theo Houthoff. We would like to thank Rody Blom, Arjen van Putten, and Erik Gertenaar for practical assistance and colleagues at the CML for their help with collecting the data. SB, EH, and MV were funded by NWO-ASPASIA and NWO-VIDI 864.13.010 granted to MV.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvs.2017.00075/full#supplementary-material

References

Ahlgren, G., Lundstedt, L., Brett, M., and Forsberg, C. (1990). Lipid composition and food quality of some freshwater phytoplankton for cladoceran zooplankters. J. Plankton Res. 12, 809–818. doi: 10.1093/plankt/12.4.809

CrossRef Full Text | Google Scholar

Altenburger, R., Backhaus, T., Boedeker, W., Faust, M., Scholze, M., and Grimme, L. H. (2000). Predictability of the toxicity of multiple chemical mixtures to Vibrio fischeri: mixtures composed of similarly acting chemicals. Environ. Toxicol. Chem.19, 2341–2347. doi: 10.1002/etc.5620190926

CrossRef Full Text | Google Scholar

Backhaus, T., Faust, M., Scholze, M., Gramatica, P., Vighi, M., and Grimme, L. H. (2004). Joint algal toxicity of phenylurea herbicides is equally predictable by concentration addition and independent action. Environ. Toxicol. Chem. 23, 258–264. doi: 10.1897/02-497

PubMed Abstract | CrossRef Full Text | Google Scholar

Cairns, J. Jr. (1983). Are single species toxicity tests alone adequate for estimating environmental hazard? Hydrobiologia 100, 47–57. doi: 10.1007/BF00027421

CrossRef Full Text | Google Scholar

Calow, P., and Forbes, V. E. (2003). Peer reviewed: does ecotoxicology inform ecological risk assessment? Environ. Sci. Technol. 37, 146A−151A. doi: 10.1021/es0324003

CrossRef Full Text | Google Scholar

Cockrell, B. J. (1984). Effects of water depth on choice of spatially separated prey by Notonecta glauca L. Oecologia 62, 256–261. doi: 10.1007/BF00379023

PubMed Abstract | CrossRef Full Text | Google Scholar

Corvalan, C., Hales, S., and McMichael, A. (2005). Millenium Ecosystem Assessment. Ecosystems and Human well-being: Health Synthesis. Geneve: World Health Organization.

Côté, I. M., Darling, E. S., and Brown, C. J. (2016). Interactions among ecosystem stressors and their importance in conservation. Proc. R. Soc. B Biol. Sci. 283, 20152592. doi: 10.1098/rspb.2015.2592

PubMed Abstract | CrossRef Full Text | Google Scholar

David, J.-P., Rey, D., Pautou, M.-P., and Meyran, J.-C. (2000). Differential toxicity of leaf litter to dipteran larvae of mosquito developmental sites. J. Invertebr. Pathol. 75, 9–18. doi: 10.1006/jipa.1999.4886

PubMed Abstract | CrossRef Full Text | Google Scholar

Fleeger, J. W., Carman, K. R., and Nisbet, R. M. (2003). Indirect effects of contaminants in aquatic ecosystems. Sci. Total Environ. 317, 207–233. doi: 10.1016/S0048-9697(03)00141-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Gagic, V., Kleijn, D., Báldi, A., Boros, G., Jørgensen, H. B., Elek, Z., et al. (2017). Combined effects of agrochemicals and ecosystem services on crop yield across Europe. Ecol. Lett. 20, 1427–1436. doi: 10.1111/ele.12850

PubMed Abstract | CrossRef Full Text | Google Scholar

Gessner, M. O., Swan, C. M., Dang, C. K., McKie, B. G., Bardgett, R. D., Wall, D. H., et al. (2010). Diversity meets decomposition. Trends Ecol. Evol. 25, 372–380. doi: 10.1016/j.tree.2010.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Gessner, M. O., and Tlili, A. (2016). Fostering integration of freshwater ecology with ecotoxicology. Freshw. Biol. 61, 1991–2001. doi: 10.1111/fwb.12852

CrossRef Full Text | Google Scholar

Graça, M. A., Maltby, L., and Calow, P. (1993). Importance of fungi in the diet of Gammarus pulex and Asellus aquaticus I: feeding strategies. Oecologia 93, 139–144. doi: 10.1007/BF00321203

PubMed Abstract | CrossRef Full Text | Google Scholar

Groendahl, S., and Fink, P. (2016). The effect of diet mixing on a nonselective herbivore. PLoS ONE 11:e0158924. doi: 10.1371/journal.pone.0158924

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamilton, S. K., and Lewis, W. M. (1992). Stable carbon and nitrogen isotopes in algae and detritus from the Orinoco River floodplain, Venezuela. Geochim. Cosmochim. Acta 56, 4237–4246. doi: 10.1016/0016-7037(92)90264-J

CrossRef Full Text | Google Scholar

Helweg, C., Mogensen, B. B., and Sørensen, P. B. (2003). Fate of Pesticides in Surface Waters, Laboraty and Field Experiments. Kopenhagen: Ministry of the Environment, Danish Environmental Protection Agency.

Hinkle-Conn, C., Fleeger, J. W., Gregg, J. C., and Carman, K. R. (1998). Effects of sediment-bound polycyclic aromatic hydrocarbons on feeding behavior in juvenile spot (Leiostomus xanthurus Lacepede: Pisces). J. Exp. Mar. Biol. Ecol. 227, 113–132. doi: 10.1016/S0022-0981(97)00265-7

CrossRef Full Text | Google Scholar

Hobson, K. A., and Welch, H. E. (1992). Determination of trophic relationships within a high Arctic marine food web using δ¹³C and δ¹⁵N analysis. Mar. Ecol. Prog. Ser. 84, 9–18.

Google Scholar

Hunting, E. R., Mulder, C., Kraak, M. H. S., Breure, A. M., and Admiraal, W. (2013). Effects of copper on invertebrate-sediment interactions. Environ. Pollut. 180, 131–135. doi: 10.1016/j.envpol.2013.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Ieromina, O., Peijnenburg, W. J. G. M., Musters, C. J. M., and Vijver, M. G. (2016). The effect of pesticides on the composition of aquatic macrofauna communities in field ditches. Basic Appl. Ecol. 17, 125–133. doi: 10.1016/j.baae.2015.08.002

CrossRef Full Text | Google Scholar

Jabiol, J., Mckie, B. G., Bruder, A., Bernadet, C., Gessner, M. O., and Chauvet, E. (2013). Trophic complexity enhances ecosystem functioning in an aquatic detritus-based model system. J. Anim. Ecol. 82, 1042–1051. doi: 10.1111/1365-2656.12079

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimball, K. D., and Levin, S. A. (1985). Limitations of laboratory bioassays: the need for ecosystem-level testing. Bioscience 35, 165–171. doi: 10.2307/1309866

CrossRef Full Text | Google Scholar

Klecka, J., and Boukal, D. S. (2012). Who eats whom in a pool? a comparative study of prey selectivity by predatory aquatic insects. PLoS ONE 7:e37741. doi: 10.1371/journal.pone.0037741

PubMed Abstract | CrossRef Full Text | Google Scholar

Mills, L. S., Soule, M. E., and Doak, D. F. (1993). The keystone-species concept in ecology and conservation. BioScience 43, 219–224. doi: 10.2307/1312122

CrossRef Full Text | Google Scholar

Moore, J. C., McCann, K., and De Ruiter, P. C. (2005). Modeling trophic pathways, nutrient cycling, and dynamic stability in soils. Pedobiologia 49, 499–510. doi: 10.1016/j.pedobi.2005.05.008

CrossRef Full Text | Google Scholar

Ormerod, S. J., Dobson, M., Hildrew, A. G., and Townsend, C. R. (2010). Multiple stressors in freshwater ecosystems. Freshw. Biol. 55, 1–4. doi: 10.1111/j.1365-2427.2009.02395.x

CrossRef Full Text | Google Scholar

Pesticide Atlas, version 2.0. (2015). Leiden University (CML) and Rijkswaterstaat-WVL. Leiden. Available online at: www.bestrijdingsmiddelenatlas.nl (Accessed Feb 2, 2016).

Pimm, S. L., Lawton, J. H., and Cohen, J. E. (1991). Food web patterns and their consequences. Nature 350, 669–674. doi: 10.1038/350669a0

CrossRef Full Text | Google Scholar

Post, D. M. (2002). Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology 83, 703–718. doi: 10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2

CrossRef Full Text | Google Scholar

Relyea, R., and Hoverman, J. (2006). Assessing the ecology in ecotoxicology: a review and synthesis in freshwater systems. Ecol. Lett. 9, 1157–1171. doi: 10.1111/j.1461-0248.2006.00966.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rohr, J. R., Kerby, J. L., and Sih, A. (2006). Community ecology as a framework for predicting contaminant effects. Trends Ecol. Evol. 21, 606–613. doi: 10.1016/j.tree.2006.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Salminen, J., Eriksson, I., and Haimi, J. (1996). Effects of terbuthylazine on soil fauna and decomposition processes. Ecotoxicol. Environ. Saf. 34, 184–189. doi: 10.1006/eesa.1996.0062

PubMed Abstract | CrossRef Full Text | Google Scholar

Santer, B., and Van Den Bosch, F. (1994). Herbivorous nutrition of Cyclops vicinus: the effect of a pure algal diet on feeding, development, reproduction and life cycle. J. Plankton Res. 16, 171–195. doi: 10.1093/plankt/16.2.171

CrossRef Full Text | Google Scholar

Schulz, R., and Dabrowski, J. M. (2001). Combined effects of predatory fish and sublethal pesticide contamination on the behavior and mortality of mayfly nymphs. Environ. Toxicol. Chem. 20, 2537–2543. doi: 10.1002/etc.5620201120

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, G. M., and Weis, J. S. (1997). Predator-prey relationships in mummichogs (Fundulus heteroclitus (L.)): effects of living in a polluted environment. J. Exp. Mar. Biol. Ecol. 209, 75–87. doi: 10.1016/S0022-0981(96)02590-7

CrossRef Full Text | Google Scholar

Song, M. Y., Stark, J. D., and Brown, J. J. (1997). Comparative toxicity of four insecticides, including imidacloprid and tebufenozide, to four aquatic arthropods. Environ. Toxicol. Chem. 16, 2494–2500. doi: 10.1002/etc.5620161209

CrossRef Full Text | Google Scholar

Stoughton, S. J., Liber, K., Culp, J., and Cessna, A. (2008). Acute and chronic toxicity of imidacloprid to the aquatic invertebrates Chironomus tentans and Hyalella azteca under constant-and pulse-exposure conditions. Arch. Environ. Contam. Toxicol. 54, 662–673. doi: 10.1007/s00244-007-9073-6

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Brink, P. J. (2006). Response to recent criticism on aquatic semifield experiments: opportunities for new developments in ecological risk assessment of pesticides. Integr. Environ. Assess. Manag. 2, 202. doi: 10.1002/ieam.5630020302

PubMed Abstract | CrossRef Full Text | Google Scholar

Vijver, M. G., Hunting, E. R., Nederstigt, T. A., Tamis, W. L., van den Brink, P. J., and van Bodegom, P. M. (2017). Postregistration monitoring of pesticides is urgently required to protect ecosystems. Environ. Toxicol. Chem. 36, 860–865. doi: 10.1002/etc.3721

PubMed Abstract | CrossRef Full Text | Google Scholar

Vörösmarty, C. J., McIntyre, P. B., Gessner, M. O., Dudgeon, D., Prusevich, A., Green, P., et al. (2010). Global threats to human water security and river biodiversity. Nature 467, 555–561. doi: 10.1038/nature09440

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallace, J. B., Eggert, S. L., Meyer, J. L., and Webster, J. R. (1997). Multiple trophic levels of a forest stream linked to terrestrial litter inputs. Science 277, 102–104. doi: 10.1126/science.277.5322.102

CrossRef Full Text | Google Scholar

Wilbur, H. M., and Dec, N. (1997). Experimental ecology of food webs: complex systems in temporary ponds. Ecology 78, 2279–2302. doi: 10.1890/0012-9658(1997)078[2279:EEOFWC]2.0.CO2

CrossRef Full Text | Google Scholar