Variation in behaviour of native prey mediates the impact of an invasive species on plankton communities

Study site and experimental design

From July 7th to August 1st 2014, we conducted a field mesocosm experiment to assess the influence of differences in vertical position of Daphnia on effects of Bythotrephes predation on plankton communities. Bythotrephes impacts on zooplankton community structure and algal production can occur over short periods of time comparable to the length of our experiment (Strecker & Arnott, 2005; Jokela, Arnott & Beisner, 2017; Azan & Arnott, 2017). Mesocosms were set up in Fletcher Lake (45.20.452′N, 78.47.798′W) Haliburton County, Ontario Canada (Table 1), where Bythotrephes was first detected in 2006 (Cairns et al., 2007). Five mil food grade polythene enclosures, 1m in diameter, 13 m deep (Filmtech Plastics, Brampton, Ontario) and filled with 10,990 L of water, were closed at the bottom and suspended from floating wooden frames anchored in the lake. Each mesocosm was filled with water pumped from 1.5 m and filtered through an 80-µm mesh to exclude crustacean zooplankton, but allow most phytoplankton and some small rotifers and nauplii to pass through. Each mesocosm was covered with screen mesh to prevent colonization by aerial insects. Dissolved oxygen (DO) levels in each mesocosm ranged from an average of 8.9–10.6g/L at the start of the experiment. Based on data from other mesocosm experiments conducted in this region, where average DO ranged from 7.9 to 8.1 g/L over the course of 8 weeks, we did not expect DO levels to change significantly during the course of the experiment.

To assess if differences in vertical position of Daphnia influence Bythotrephes predation on plankton communities, we manipulated Bythotrephes presence and absence in mesocosms across a gradient of Daphnia vertical position. We used thermal gradients in lieu of light measurements to determine the vertical position of Daphnia in our mesocosms. With 4.6 mg/L of dissolved organic carbon (DOC) present in Fletcher Lake (CAISN 2011), only 1% of ultra-violet and photosynthetically active radiation (PAR) was estimated to be present at 4m (based on average summer solar radiation for south-central Ontario (NASA Langley Atmospheric Science DataCenter, 2016), and relationships between DOC and light attenuation (Morris et al., 1995). Therefore, we expected no light to be present in the hypolimnion (9–13 m) and quantified the vertical position of Daphnia based on their numeric proportion in this region. Bythotrephes predation was expected to be completely absent in the hypolimnion, where light penetration does not occur, as it is a visual predator requiring light to feed (Jokela, Arnott & Beisner, 2013).

After allowing phytoplankton assemblages to increase in biomass without the presence of large zooplankton (>80 µm) for three days, mesocosms were stocked with the entire zooplankton community from two uninvaded lakes; Bonnie Lake (45.17.36′N 79.06.45′W; Bracebridge Municipality) and Echo Lake (45.17.36′N, 79.06.45′W; Lake of Bays Municipality), in south-central Ontario (Table 1, Hasnain & Arnott, 2019). These lakes contained the same Daphnia species, but with contrasting vertical positions (Table 1). Each mesocosm was randomly inoculated with an ambient density of zooplankton sampled from the same volume of water from either Echo or Bonnie Lake, or half volumes for mesocosms stocked from both lakes, resulting in 16 mesocosms stocked from Echo Lake, 16 mesocosms stocked from Bonnie Lake and eight mesocosms stocked from both lakes. This study was part of a larger experiment where an additional treatment was applied to all mesocosms three weeks after the start of the experiment. Only data from the first three weeks of the experiment was included in this study.

Half of the mesocosms were randomly assigned the invasion treatment. Bythotrephes individuals were collected from Fletcher Lake and Lake of Bays, Muskoka, Ontario (45°15.00′N, 79°04.00′W) using an 80 µm mesh net, and stocked at a density of 10 individuals per m³ of epilimnion volume (23 individuals per mesocosm; similar to the average density observed in invaded lakes in this region Jokela, Arnott & Beisner, 2011) at the beginning of the experiment. Zooplankton communities in all mesocosms were acclimated for one week prior to Bythotrephes addition. There were no a priori differences in the proportion of hypo- and non-epilimnetic (hypo- and metalimnetic) Daphnia between invaded and uninvaded mesocosms in week 0 (Fig. S2, Gamma GLMs, hypo: p = 0.911; non-epi: p = 0.792). Daphnia daytime vertical position remained stable across mesocosms throughout the duration of the experiment (see Fig. S4).

There were also significant differences in the abundance of zooplankton between Echo Lake and Bonnie Lake, which influenced initial zooplankton densities in our mesocosms. Overall zooplankton density was greater in enclosures stocked from Echo Lake as compared to Bonnie Lake (Fig. S1); log normally distributed linear model with lake origin (either Echo, Bonnie, or Both) as a predictor variable and zooplankton density as a response variable (p < 0.0001, Echo: 2.490 ±  1.040, Bonnie: 1.497 ± 0.974, Both: 2.910 ± 0.861, see Table S2 for initial densities in each mesocosm). There was no difference in zooplankton starting densities between mesocosms assigned Bythotrephes presence and absence treatments (Fig. S5; Absent:1.954 ± 1.026, Present: 2.287 ± 1.126 linear model, p > 0.05 for all major taxonomic groups). There was very little variation in total algal densities between the mesocosms at the start of the experiment (mean = 6.438 ±  0.110 µg/L, Table S2).

Sampling protocol and identification

All zooplankton samples were collected between 10 AM and 3 PM. Each mesocosm was sampled prior to the addition of Bythotrephes (week 0) to determine abundance and depth distribution of zooplankton and phytoplankton. Samples were also taken at the end of the study (week 3). Epi-, meta- and hypolimnetic boundaries were determined weekly using the thermal profile of Fletcher Lake using a YSI model 600 OMS V2. Zooplankton samples from each mesocosm were collected by towing a closing net with an 80 µm mesh (15 cm diameter) through each thermal layer (starting 20 cm above the enclosure bottom for the hypolimnion). Each layer was sampled in spatial order (i.e., epilimnion followed by metalimnion, and then hypolimnion) to prevent disruption of any layer prior to sampling. Samples were preserved in 90% ethanol. Water samples were collected from the middle of the meta- and hypolimnion using a 2L Van Dorn sampler in weeks 0 and 3. For the epilimnion, water samples were collected from 10 cm below the surface of the mesocosm by submerging the sample container. We determined total algal biomass as well as biomass of green algae, cyanobacteria, diatoms and cryptophytes spectrophotometrically by analyzing a well-mixed 25 ml subsample within 24 h of sample collection (BBE moldaenke Algae Lab Analyser; Schwentinental, Schleswig-Holstein, Germany).

Zooplankton were enumerated using sub-samples of a known volume and identifying all individuals within each subsample until no new species were found three sub-samples in a row. We counted a minimum of seven sub-samples for each thermal layer in each mesocosm and identified all individuals to the species level (Ward & Whipple, 1959; Smith & Fernando, 1978; Melo & Hebert, 1994; Witty, 2004; Haney et al., 2013). Bosmina freyi and Bosmina liederi were grouped as “Bosmina freyi/liederi” and Daphnia pulex and Daphnia pulicaria as “Daphnia pulex/pulicaria” due to morphological similarities between these species. For Daphnia species, only adults were enumerated. Juvenile copepods were identified as either nauplii or copepodites, without distinguishing between cyclopoid and calanoid juveniles as both have similar diets and occupy a similar trophic position (Finlay & Roff, 2004).

Statistical analysis

All analyses were conducted in R v 3.2.4 (R Core Team, 2016) using bbmle v 1.0.17, glmmADMB v 0.7.7, fitdistrplus 1.0.7, piecewiseSEM v 2.1.2, robustbase v 0.93-3 with α = 0.05.

Because of differences in starting densities between mesocosms stocked from Bonnie, Echo or both lakes (Table S2, Fig. S1), we standardized the change in density between weeks 0 and 3 by calculating per capita change in density for each species and functional group per mesocosm; calculated as density in week 3 divided by density in week 0, which may be sensitive to sampling and demographic effects (see Supplemental Information). We separately assessed Bythotrephes impacts on Daphnia, small (<0.85 mm) and large cladocerans (>1.0 mm) excluding Daphnia as well as calanoid, cyclopoid, and juvenile copepods (see Supplemental Information for details on zooplankton categorization). It is possible that the proportion of total hypolimnetic Daphnia may have changed during the experiment due to differences in abiotic and biotic variables between the lakes of origin and Fletcher Lake. We found no significant differences in the proportion of hypolimnetic Daphnia between the beginning (week 0) and the end of the experiment (week 3, see Supplemental Information for statistical assessment details).

We calculated the density of each species in each thermal layer for every mesocosm by estimating the total number of individuals present in sub samples. Species-specific density in each layer was calculated by dividing the total number of individuals by the volume sampled with a single vertical tow. To calculate the total density of a species in a mesocosm, we summed the total number of individuals in the epi- meta- and hypolimnion samples for each mesocosm and divided this by the total volume of the mesocosm that would be sampled in a single vertical tow. We could not conclusively state whether the absence of a species in our sample represented a true absence in our mesocosm as we could only sample 2% of our mesocosm volume (230L) for zooplankton and <1% for phytoplankton. We addressed this by adding minimum detection densities to all taxa assessed (see Supplemental Information).

We used piecewise structural equation modelling (piecewise SEM, Shipley, 2016; Grace, 2006; Lefchek, 2015) to explore causal relationships between per capita change in density of major zooplankton groups (see above), total algae biomass in week 3, Bythotrephes presence, and the proportion of total hypolimnetic Daphnia in week 0 (see Supplemental Information for SEM implementation details). We chose to assess the effect of Daphnia vertical position using the proportion of total hypolimnetic Daphnia in week 0 as it represents the Daphnia vertical position in mesocosms at the start of the experiment, thereby removing the possibility of any direct (i.e., predation) and indirect (i.e., induced response) Bythotrephes effects that may be present in invaded mesocosms after the start of the experiment.

Piecewise SEMs were also fit for the most common species in each zooplankton group with statistically significant causal relationships detected in the full model (Fig. 1), with Bythotrephes presence, and proportion of total hypolimnetic Daphnia in week 0 as exogenous variables, excepting paths for the most common calanoid and cyclopoid species where per capita change in copepodite density was also included. Significant causal relationships in piecewise SEMs for both common species were also re-fit with GLMs to assess if there was a significant interaction between Daphnia vertical position and Bythotrephes presence. We also refit piecewise SEMs with Bythotrephes presence and the proportion of hypolimnetic individuals for Daphnia species whose per capita change in density was identified as significantly impacted by Bythotrephes presence in earlier SEMs to assess if any effect of proportion of total hypolimnetic Daphnia observed was driven by these species. Separate piecewise SEMs were also fit for week 3 biomass of the four major algal groups; green algae, cyanobacteria, diatoms, and cryptophyta using the protocol outlined for the full model, with the nested version also confirmed using AIC.

Paths with significant causal relationships (p < 0.05) between per capita change in density of a zooplankton group or final algal biomass, Bythotrephes presence, and vertical position of Daphnia were re-fit with either normal or gamma distributed generalized linear models to assess the possibility of a significant interaction between the explanatory/exogenous variables. We added stocking lake as additional explanatory variable to our analysis with generalized linear models as it is possible that there are unknown differences between stocking lakes that may be underlying the observed patterns. For all GLMs, we used AICc to assess fit between gamma and log-normal distributions. Model fit was also assessed by visually examining plots of residual versus fitted values and square root of the standard deviance of residuals versus fitted values. Cook’s distance was used to identify influential points (leverage >1.0). Minimum adequate models were chosen using log-likelihood ratio tests based on Crawley’s (2011) procedure. If influential points were detected, gamma or log-normal robust GLMs were fitted using Mallows or Huber type robust estimators (Cantoni & Ronchetti, 2001; Cantoni & Ronchetti, 2006), which down weights the effect of influential points on model fit. For robust GLMs, minimum adequate models were chosen using Wald-type tests. We did not find any density dependent effects of Bythotrephes predation due to differences in starting densities (details about statistical methods used provided in Supplemental Information).

Algal biomass

We observed a significant interaction between Bythotrephes presence and proportion of hypolimnetic Daphnia in week 0 on total algal biomass at the end of study (Gamma Robust Regression, p = 0.0005, df = 314, Table S13). Total algal biomass was greater in invaded mesocosms as compared to uninvaded mesocosms, with the greatest increase observed in mesocosms with greatest proportion of hypolimnetic Daphnia in week 0 (Fig. 3C, Table S13). Total algal biomass at the end of the study was greater in mesocosms stocked from Echo Lake as compared to those stocked from Bonnie Lake or both lakes (p = 0.001, Fig. S8).

Cyanobacteria were the most abundant functional group, followed by green algae, diatoms, and cryptophytes. Green algal biomass was predicted by the proportion of total Daphnia in the hypolimnion, and per capita change in density of large cladocerans, and juvenile copepods (Fig. 2C, vertical position SE: 0.44, Large cladocerans SE: 0.43, Juvenile copepods SE: −0.14, Table S6), and was greater in mesocosms with a greater proportion of hypolimnetic Daphnia in week 0 (Fig. 4E, p = 0.010, df = 314, Table S13). Green algal biomass was greater in mesocosms stocked from Bonnie Lake as compared to those stocked from Echo Lake or both lakes (Fig. S8, p = 0.02 df = 37, Table S13). Diatom biomass was predicted by Bythotrephes presence, the proportion of hypolimnetic Daphnia, and per capita change in density of small cladocerans, large cladocerans, calanoids and juvenile copepods (Fig. 2D, Table S7, Bythotrephes SE: −0.44, vertical position SE:1.48, Small cladoceran SE: −0.62, Large cladoceran SE: 0.65, Calanoid SE: −1.75, Juvenile copepods SE: −0.28), and was greater in mesocosms with a greater proportion of hypolimnetic Daphnia (Fig. 4F). We did not observe any significant interactive effect of Bythotrephes presence and proportion of total hypolimnetic Daphnia (Table S13, p = 0.15) or any effect of stocking lakes on diatom biomass (p = 0.114). There was no effect of Bythotrephes presence, the proportion of total hypolimnetic Daphnia, or per capita change in density in any zooplankton group on cyanobacteria and cryptophyta biomass (Tables S8–S9).

Species-level responses

For Daphnia mendotae, only Bythotrephes presence predicted per capita change in density (piecewise SEM, Fig. 2B, Table S10, Bythotrephes SE: −0.09). Per capita density for D. mendotae was lower in invaded mesocosms as compared to uninvaded mesocosms (Fig. 4A). There was no effect of stocking lake on per capita density for D. mendotae (p = 0.556, Table S13). There was no effect of Bythotrephes presence or proportion of total hypolimnetic Daphnia on per capita change in D. catawba density (Table S10).

Per capita change in density in B. freyi/leideri was only predicted by the proportion of total hypolimnetic Daphnia in week 0 in our SEM analysis (Fig. 2B, Table S10, vertical position SE: 0.01), with no effect of Bythotrephes presence. However, we observed a significant interactive effect of Bythotrephes presence and Daphnia vertical position in our GLMs (p = 0.008, Table S13, Fig. 4B), with greater increase in B. freyi/leideri observed in uninvaded mesocosms with greater proportion of hypolimnetic Daphnia. Greater increase in B. freyi/leideri was observed in mesocosms stocked from Bonnie Lake as compared to those stocked from Echo Lake or both lakes (Fig. S8, p < 0.001). Both Bythotrephes presence and proportion of total hypolimnetic Daphnia in week 0 predicted per capita change in E. longispina density (piecewise SEM, Fig. 2B, Table S10, Bythotrephes SE: −0.03, vertical position SE: −0.05). However, there was no interaction between these variables (Table S13, p = 0.7431). Per capita density was lesser in invaded mesocosms as compared to uninvaded mesocosms and decreased in mesocosms with a greater proportion of hypolimnetic Daphnia (Fig. 4C). There was no effect of stocking lake on per capita change in E. longispina density (p = 0.750). There was no effect of Bythotrephes presence or the proportion of total hypolimnetic Daphnia in week 0 on per capita change in E. tubicen density (piecewise SEM, Table S10).

For copepods, per capita change in C. scutifer (cyclopoid) density was predicted by the proportion of total hypolimnetic Daphnia in week 0 and Bythotrephes presence (piecewise SEM, Fig. 2B, Table S8, Bythotrephes SE: −0.18, vertical position SE: 0.13). Per capita increase in C. scutifer density was smaller in invaded mesocosms as compared to uninvaded mesocosms with a greater proportion of hypolimnetic Daphnia (Fig. 4D, Gamma GLMs, Table S13, p = 0.02, df = 36). There were no differences in per capita change in C. scutifer density between mesocosms stocked from Bonnie, Echo, or both lakes (p = 0.126, Table S13). There was no effect of Bythotrephes presence or proportion of total hypolimnetic Daphnia on per capita change in S. oregonensis (calanoid) density (Table S10).

Indirect impacts on other zooplankton groups

Fewer small cladocerans were present in invaded mesocosms with deeper vertical position of Daphnia (i.e., greater proportion of hypolimnetic Daphnia). This may be linked to Bythotrephes preference for large-bodied cladocerans such as Daphnia (Schulz & Yurista, 1998) in mesocosms where more Daphnia are epilimnetic. In mesocosms with a deeper Daphnia, Bythotrephes predation on small cladocerans led to the smaller increases observed. Bythotrephes predation reduces small cladoceran abundance (Vanderploeg, Liebig & Omair, 1993; Yan & Pawson, 1997; Yan et al., 2001; Barbiero & Tuchman, 2004; Strecker et al., 2006; Strecker & Arnott, 2008; Kerfoot et al., 2016), although this is not consistently observed (Lehman & Cáceres, 1993; Barbiero & Tuchman, 2004; Strecker & Arnott, 2005; Hessen, Bakkestuen & Walseng, 2011). The effect of Bythotrephes predation on small cladocerans is varied, with declines in abundance ranging from 40–126% across invaded lakes (Vanderploeg, Liebig & Omair, 1993; Yan et al., 2001; Kerfoot et al., 2016). Our results suggest that vertical position of Daphnia could be an important factor explaining this variation.

B. freyi/leideri density increased as the proportion of hypolimnetic Daphnia increased. For E. longispina, deeper vertical position of Daphnia was associated with smaller increases in density. These contrasting effects suggest that B. freyi/leideri in our mesocosms could be constrained by competition for grazing with Daphnia (DeMott & Kerfoot, 1982). The greater increase in B. freyi/leideri density observed in mesocosms with deeper Daphnia vertical position may also be a result of the greater increase in density observed in mesocosms stocked from Bonnie Lake as compared to Echo Lake or both lakes. Since Daphnia in Bonnie Lake have a deeper vertical position, Daphnia vertical position in mesocosms stocked from this lake was likely deeper, thereby contributing to this observed pattern. For E. longispina, Bythotrephes predation is likely underlying the smaller increase in density observed in invaded mesocosms. This contradicts observations from literature (Kelly et al., 2013), where E. longispina abundance increased in Bythotrephes-invaded lakes, suggesting that the negative effect of Bythotrephes observed may be mediated by other factors. It is unclear what is driving the lack of increase in E. longispina density with deeper vertical position of Daphnia. A deeper vertical position of Daphnia may increase competition between E. longispina and other bosminid species.

Contrary to our prediction, Bythotrephes presence did not impact change in total calanoid or S. oregonensis density. This is surprising as negative effects of Bythotrephes on calanoid copepods (Strecker & Arnott, 2005; Strecker et al., 2006; Hessen, Bakkestuen & Walseng, 2011; Bourdeau, Pangle & Peacor, 2011; Kelly et al., 2013) have been broadly observed in the literature. Furthermore, the lack of an effect of deeper vertical position of Daphnia on calanoid density also contradicts expected increases due to reduced competition between Daphnia and calanoids (Sommer et al., 2003). Both Daphnia and calanoids primarily feed on algae, therefore less competition for algae between Daphnia and calanoids is expected in mesocosms with greater proportion of hypolimnetic Daphnia, leading to an increase per capita calanoid density. We may not have detected the impacts of these competitive interactions due to the short experiment time frame.

We observed interactive negative effects of Bythotrephes and the vertical position of Daphnia on change in total cyclopoid density. A larger increase in C. scutifer density was observed in uninvaded mesocosms with more hypolimnetic Daphnia as compared to invaded mesocosms, suggesting that Bythotrephes presence negatively impacts C. scutifer density. Kelly et al. (2013) observed similar Bythotrephes impacts on C. scutifer abundance in Canadian and Norwegian lakes. In contrast, no Bythotrephes effect was observed by Hessen, Bakkestuen & Walseng (2011), while a positive effect was observed by Walseng, Andersen & Hessen (2015). Although Bythotrephes effects on C. scutifer remain unresolved in observational studies, our results provide the first experimental evidence of negative Bythotrephes impacts on C. scutifer abundance. The negative effect of vertical position of Daphnia on total cyclopoid and C. scutifer density may be a consequence of indirect interactions. Deeper vertical position of Daphnia was linked to increased juvenile copepod density, which had a positive effect on total cyclopoid density. The negative effect of vertical position of Daphnia on cyclopoids may be a result of competitive interactions which were not included in our analysis.

It is unclear why Daphnia vertical position altered in our mesocosms as compared to their source lakes after inoculation. We expected three clearly differentiated groups of mesocosms where vertical position aligned with the lakes used for inoculation, i.e., high (stocked from Echo Lake, very low proportion of hypolimnetic Daphnia), low (stocked from Bonnie Lake, very high proportion of hypolimnetic Daphnia), and mix (stocked from both lakes, with the median proportion of hypolimnetic Daphnia between Echo Lake and Bonnie Lake), instead of the gradient that we observed. Daytime vertical position in Daphnia species is influenced by many ecological forces, including predator presence and location in the water column, predator type, food availability, exposure to ultra-violet radiation, and temperature-related metabolic costs (Leibold, 1990; Boeing, Leech & Williamson, 2004; Williamson et al., 2011; Larsson & Lampert, 2012). In addition, there is also a strong genetic component, with large intra-population variation in vertical position observed in response to predator cues (De Meester, 1993; De Meester, 1996). Similar zooplanktivorous fish species were present in all three lakes (Ontario Ministry of Natural Resources and Forestry, 2020), making it unlikely that these changes were driven by different predator cues. Mesocosms were filled with water from Fletcher Lake, which has similar physicochemical properties as Echo Lake and Bonnie Lake, with the exception of DOC, which was highest in Echo, followed by Fletcher, and Bonnie Lake (Fletcher: 4.6 mg/L, Bonnie: 2.3 mg/L, Echo: 6.6 mg/L). It is possible that differences in light availability and UV radiation as a consequence of differences in DOC influenced vertical position of Daphnia.

We also did not observe any change in the proportion of hypolimnetic Daphnia between week 0 and week 3 in invaded mesocosms. Bythotrephes can apply strong predation pressure on many Daphnia species (Boudreau & Yan, 2003; Strecker et al., 2011, reviewed in Azan, Arnott & Yan, 2015) and has been shown to induce a deeper vertical position in some D. mendotae populations (Bourdeau et al., 2013). Our observations match those of Kiehnau & Weider (2019), where no differences in phototactic behaviour in response to Bythotrephes kairomone were observed for D. ambigua, D. mendoate, and D. pulicaria clones from Lake Mendota, Wisconsin prior to and after Bythotrephes invasion.

Primary production and trophic cascades

Daphnia vertical position mediated the impacts of Bythotrephes on primary production by altering grazing pressure. Total algal biomass was greater in invaded mesocosms with deeper vertical position of Daphnia as compared to uninvaded mesocosms, likely a result of reduced small cladoceran abundance due to Bythotrephes predation and reduced epilimnetic grazing by Daphnia. Interestingly, we did not observe any effect of vertical position of Daphnia on total algal biomass in invaded mesocosms. For invaded mesocosms with shallower Daphnia, reduced grazing pressure on algae due to Bythotrephes preferential predation on Daphnia (Schulz & Yurista, 1998) may have resulted in the high algal biomass observed. In invaded mesocosms with deeper vertical position of Daphnia, Bythotrephes predation on smaller cladocerans in the absence of preferred Daphnia prey likely reduced grazing pressure on algal biomass resulting in a similar outcome. The invasion of Bythotrephes in some North American temperate lakes is associated with trophic cascades due to reduced grazer biomass (Walsh, Carpenter & Vander Zanden, 2016; Martin, Walsh & Vander Zanden, 2022) which has been observed in some lakes but not others (Strecker & Arnott, 2008). Our results suggest that in addition to the nutrient status of these lakes (Walsh, Carpenter & Vander Zanden, 2016), vertical position of Daphnia could be an important factor explaining the varied Bythotrephes impacts on primary production in these studies.

Total algal biomass increased across most mesocosms (Fig. S7) regardless of grazing by Daphnia and other cladocerans. Daphnia grazing was reduced in mesocosms with more hypolimnetic Daphnia, likely increasing green algal and diatom biomass. Bythotrephes presence negatively impacted diatom biomass. We also observed a negative effect of increasing small cladoceran density on diatom biomass. Apparent competition due to Bythotrephes preferential predation on larger-bodied Daphnia could lead to increased small cladoceran grazing. Furthermore, bosminids (i.e., B. freyi/leideri, E. longispina) were the most abundant small cladoceran taxa in our mesocosms. Their demonstrated selectivity on diatoms and green algae (Fulton III, 1988; Tõnno et al., 2016) likely resulted in the decline observed.