The Potential Impacts of Invasive Quagga and Zebra Mussels on Macroinvertebrate Communities: An Artificial Stone Substrate Based Field Experiment Using Stable Isotopes

Introduction

The spread of invasive species in freshwater environments is observed as a global phenomenon (Lodge et al., 1998; Francis and Chadwick, 2012; Jones et al., 2021). While investigators have recorded multiple impacts of invasive species on native or established biological communities (Gherardi and Acquistapace, 2007; Stiers et al., 2011; Boltovskoy and Correa, 2015), the prediction of the impact of a newly established non-native species remains a challenge (Williamson, 1999; Roy et al., 2014). This is mainly due to the fact that impacts of a given species may vary over space and time, and might equally depend on the specificity of the invaded ecosystem (Strayer and Malcom, 2006; Karatayev et al., 2021). Moreover, the impact might be potentially measurable only after a considerable delay following invasion (Crooks, 2005; Ricciardi et al., 2013).

Ecologists have, therefore, exerted efforts to understand the potential impacts of invasive species on freshwater systems by utilizing different approaches that include documented knowledge of the literature (Copp et al., 2009; Roy et al., 2014) coupled with mathematical and statistical modeling (Ricciardi, 2003; Kulhanek et al., 2011). While these studies often depend on existing trends, it is more difficult to predict the impacts of the invasion when establishments occur in a novel location with no prior baseline data (Kulhanek et al., 2011). Moreover, some aquatic systems are invaded by more than one bio-physically similar or genetically related taxa, although the important species-specific traits seem to be underestimated (Ricciardi, 2003).

Two highly invasive taxa of particular concern in multiple freshwater lakes and rivers are dreissenids of the Ponto-Caspian origin (Son, 2007; Strayer et al., 2020): quagga mussel (Dreissena rostriformis bugensis) and zebra mussel (D. polymorpha) (Molloy et al., 2007; Evariste et al., 2018). The establishment of these two dreissenids has been extensively reported with synergistic and antagonistic shifts in the community structure of pre-existing freshwater ecosystems (Mörtl and Rothhaupt, 2003; Sousa et al., 2011). Specific changes following invasion include a measurable increase in benthic macroinvertebrate density (Ricciardi, 2003). As such, the physical structural complexity of the mussel bed provides a safe haven or fine-scale refugia against predators (González and Downing, 1999; Ward and Ricciardi, 2007), providing protection against unintended movement of waves that facilitate macroinvertebrates (Ricciardi et al., 1997; Stewart et al., 1998). Evidently, the biofouling of mussels has deleterious effects on certain macroinvertebrates (Sousa et al., 2011; Strayer and Malcom, 2018), and in many lakes, the native bivalves showed a large decrease in abundance after the invasion of dreissenids (Strayer and Malcom, 2007; Burlakova et al., 2014).

Most interestingly, detritivores and grazing herbivorous species may benefit from nutrient release following pseudo-feces excretion and biofilm deposition from dreissenids (Kobak et al., 2013). Following mussel invasion, remarkable increases in macroinvertebrate abundance and biomass of gammarids and snails have been documented (MacIsaac, 1996; Bially and MacIsaac, 2000). Finally, as posited by the invasional meltdown hypothesis, in multiple invaded ecosystems, the positive interactions among invaders initiate positive population-level feedback that intensifies influence and promotes further invasions (Simberloff and von Holle, 1999).

In the past decades, several invasive freshwater mussels have colonized most of the European perialpine lakes. In particular, the zebra and quagga mussels pose continuous challenges in Lake Constance. Adult quagga mussel was first reported in the lake in 2016 (Haltiner et al., 2022), while zebra mussel has colonized the lake since the late 1960s (Walz, 1973). In fact, the quagga mussel recently increased its range and density, and replaced the zebra mussel over most parts of the lake, only 2 years after the first record. These observations, combined with the existing increasing invasion rate and invasive species in Lake Constance, tentatively support the “invasional meltdown” model (Simberloff and von Holle, 1999). It could be argued, at least theoretically, that Lake Constance is experiencing an “easier” invasion as the cumulative number of introduced species increases and as facilitative interactions are expected to exacerbate the effect of these invaders. Indeed, the impacts of zebra mussels on the macroinvertebrate communities of Lake Constance have been well investigated (Mörtl and Rothhaupt, 2003; Gergs and Rothhaupt, 2008, 2015). However, the impacts of the recent invader quagga mussels remain unknown.

Therefore, the first objective of this study was to measure and compare the potential impacts of quagga and zebra mussels on macroinvertebrate communities in Lake Constance. Using an experimental approach, we aimed to observe macroinvertebrate colonization across four series of manipulated stone substrates treated with different specimens: blank stones, dead shell analogs of dreissenids, living adult quagga mussels, or living adult zebra mussels. By simulating these four potential substrate “scenarios” here, we aimed to improve our knowledge regarding these two highly concerning invasive species in the absence of robust baseline data (specifically for quagga mussels) for Lake Constance. Understanding food-web interaction using possible scenarios is critical for accurate prediction of future impacts, which would facilitate undertaking necessary prevention and management approaches (Byers et al., 2002). We expected that substrates with living mussel and shell treatments would cause higher macroinvertebrate abundance compared to the blank treatments.

Unraveling the mechanisms that allow multiple aquatic invasive species to coexist can provide insights into freshwater community ecology. However, dietary information on component species of these assemblages is often difficult to obtain. Therefore, our second objective was to identify potential food sources and estimate the relative trophic position of newly settled macroinvertebrate species (on the stone substrates). Therefore, we measured stable carbon and nitrogen isotopes (δ¹³C and δ¹⁵N) in soft tissue specimens obtained from five macroinvertebrate species (Gergs et al., 2011; Tran et al., 2015). Using the four manipulated substrate setups, we examined whether trophic niche differentiation and microhabitat segregation explain their co-existence in the lake. Both mussels are able to exploit phytoplankton or suspended detritus as significant energy sources. Yet, it is expected that the species might compete with each other in areas of sympatry. Therefore, we hypothesized that the presence of these mussels will influence not only the community composition but also the food resources of newly settled macroinvertebrates. Specifically, we measured and compared the isotopic niche size of each species and that of the whole community, and assessed isotopic niche overlap among the two invasive mussel species using the δ¹³C and δ¹⁵N values.

At present, the vast majority of dreissenids found in the littoral zone of Lake Constance are quagga mussels, while a short time ago, this area was dominated by zebra mussels. Species-specific population size may fluctuate because of the variation in the adult mortality rate, and the survival of adult mussels may be an important component of competitive displacement. Therefore, using this in situ setup, we opportunistically examined the survivorship of adult quagga and zebra mussels in the four treatments.

Materials and Methods

Study Area

The manipulated stone substrate experiment was conducted at the littoral area in Upper Lake Constance (N: 47°41′31.0″, E: 9°12′07.0″), one of the two major basins in Lake Constance, Germany. The littoral area has predominantly stony gravel and pebble substrate. It is surrounded by a semi-natural habitat, with scarce suburban housing and with hardly any extensive artificial land use (Gergs et al., 2011). The stone substrates were deployed at a depth of 0.80 m.

Field Experimental Design

The experimental design followed the former study by Mörtl and Rothhaupt (2003). The standardized stone substrates (commercially available concrete cobblestones) had an upper surface area of 121 cm². The upper surface area of each stone substrate was manipulated to obtain four treatments with eight replicates each. The first control treatment (B treatment) consisted of blank and clean stones with no shells and no mussels on them. Treatment 2 (S treatment) consisted of stones with eviscerated shells of mussels that were glued together to mimic intact dead mussels (Botts et al., 1996). Shells were glued to the substrate at their posterior side and positioned in a similar direction, simulating field observations of adult dead mussels at the lakeshore. Treatment 3 (Q treatment) and treatment 4 (Z treatment) consisted of stones with upper surfaces covered by living adult quagga and zebra mussels, respectively. Living mussels of individual species were glued to the substrate at their posterior side. Care was taken to ensure that individual mussels could open and close their shells as in natural cases to ensure survival. The abundance of shells or living mussels on the upper surfaces of stones was equivalent to 3,719 ± 413 ind./m². The substrate appearance resembled that of a natural mussel bed or druse, simulating field observations of living mussels. Caution was taken to ensure that glued individual mussels or shells were approximately equidistant.

To prepare the stone treatments, adult quagga and zebra mussels of ca. 20 mm shell length had been collected by hand at a littoral site of the lake where the two dreissenids coexisted by then (i.e., June 2018). Based on the morphological traits, mussels had been identified and separated as quagga or zebra mussels (Pathy and Mackie, 1993; Kerambrun et al., 2018). To prepare the shell treatments, shells of eviscerated mussels were washed and dried in the laboratory at room temperature for 72 h. Stones were deployed at the study site for the experiment between June and July 2018 (28 days). Elsewhere, artificial substrates have been reported to reach stabilized newly settled macroinvertebrate communities within 20–60 days (Roby et al., 1978; Wise and Molles, 1979). Therefore, through the experimental period of 28 days, we envisioned to capture the invasion processes and evaluate their effects on the density and taxonomic richness of newly settled macroinvertebrates within the first 4 weeks (short time duration).

On the last day of the experiment, all the stones were collected from the lake, thoroughly washed, and carefully inspected for attached (newly settled) macroinvertebrates. The newly settled macroinvertebrates were gently removed with a brush, washed, and collected from each stone after filtering through a 200 μm mesh. We took all the newly settled macroinvertebrates in the entire area of the respective stone into account. Afterward, living adult quagga and zebra mussels that had been glued on the stones were gently removed from their respective stones, and their survival was assessed by observing shell gape and soft tissue. If no soft tissue was present inside the shells, then mortality was recorded. One stone in the Z treatment was upside down at the end of the experiment, and this sample was omitted from further analysis.

Macroinvertebrate Counting

In the laboratory, all newly settled macroinvertebrate samples were sorted and counted under a dissecting microscope (× 10 magnification, Zeiss Stemi 2000-C). Identification was made to species or genus level, except for Diptera, which mainly included Chironomidae. Chironomidae were identified at the subfamily level due to morphological ambiguity at their smallest size ranges. The size of macroinvertebrates was measured using a microscope equipped with an ocular micrometer. Newly settled dreissenids larger than 1 mm in length were separated into quagga or zebra mussels based on morphological traits. We pooled the small dreissenids (≤ 1 mm in length) into one group.

Stable Isotope Analysis

The five common newly settled macroinvertebrate species recorded from the experiment, which included the two mussel species (quagga mussel D. rostriformis bugensis and zebra mussel D. polymorpha), two snail species (New Zealand mud snail Potamopyrgus antipodarum and faucet snail Bithynia tentaculata), and one gammarid species (killer shrimp Dikerogammarus villosus), were used for stable isotope measurement. The soft tissues of the two mussel species and two snail species were eviscerated for stable isotope analysis. Individuals of the gammarid species were used for stable isotope analysis.

Lipid content was removed from dried and pulverized samples (Post et al., 2007; Rothhaupt et al., 2014) and weighed (ca. 0.85 mg) in tin cups to the nearest 0.001 mg, using a micro-analytical balance. The ¹³C/¹²C or ¹⁵N/¹⁴N ratios of samples were determined using two instruments. Half of the samples were measured at the stable isotope lab of the University of Konstanz in Germany. Here, samples were combusted in a vario PYRO cube elemental analyzer (Elementar Analysensysteme, Germany) at 1,120°C. Resulting CO₂ and N₂ were separated by gas chromatography and passed into an IsoPrime isotope ratio mass spectrometer (IRMS; Isoprime Ltd., Manchester, United Kingdom). The rest of the samples were measured at the stable isotope laboratory of the University of Erlangen-Nürnberg, Germany. Here samples were combusted in a Costech Elemental Analyzer (ECS 4010; Costech International, Pioltello, Italy; now NC Technologies, Bussero, Italy). The resulting gases were separated and measured using Thermo Scientific Delta V plus IRMS (Thermo Fisher Scientific, Bremen, Germany). The repeatability and comparability of results obtained from the two laboratories were first ensured by measuring multiple macroinvertebrate samples (“dummy” samples) in both laboratories before proceeding with the samples used for this experiment.

All stable isotope values (δ¹³C or δ¹⁵N) were reported in the δ notation (per mill) where δ = [1,000 × (Rsample/Rstandard) –1] ‰, relative to the Pee Dee Belemnite (PDB) for carbon and atmospheric N₂ for nitrogen in parts per thousand deviations (‰).

Data Analysis

Non-parametric Kruskal–Wallis test was used to analyze the differences in the abundances of newly settled macroinvertebrates among four different treatments, as the normality of the whole dataset was not reached (Shapiro–Wilk test). Dunn’s test was used for post hoc pairwise comparisons. Non-parametric Mann–Whitney Wilcoxon test was used to compare the abundances of quagga and zebra mussels in each treatment.

To identify differences in community composition across all treatments, we examined dissimilarities among macroinvertebrate assemblages by applying non-metric multidimensional scaling (NMDS) using Bray–Curtis dissimilarities. Mean abundance data (ind./m² of newly settled macroinvertebrate per taxon) were used in the analysis, and data were standardized by the decostand function using the method “total.” We then assessed whether the composition of the macroinvertebrate community differed among treatments using analysis of similarity (ANOSIM). In order to reduce distortion of macroinvertebrate assemblage in the community, taxa accounting for less than 0.5% of the total mean abundance per treatment were excluded from the analysis.

Graphical summaries of the mean abundance of abundant newly settled macroinvertebrates (ind./m² + SE) per treatment were made for taxa that can be identified up to the genus or species level, which included D. rostriformis bugensis, D. polymorpha, Dreissena (≤ 1 mm), P. antipodarum, B. tentaculata, D. villosus, Katamysis warpachowskyi, and Caenis (Figure 1). A summation of abundances for all the newly settled macroinvertebrates present in each sample is expressed as the total abundance.

FIGURE 1

Figure 1. The abundance of abundant newly settled macroinvertebrates in four different treatments. Taxa that can be identified to species or genus level are presented. Error bars indicate standard errors.

The stable isotope values of newly settled macroinvertebrates were not normally distributed (Shapiro–Wilk test). Kruskal–Wallis test was used to analyze the differences in δ¹³C and δ¹⁵N values for the five main newly settled macroinvertebrate species among four treatments, and Dunn’s test was used for post hoc pairwise comparisons. Mann–Whitney Wilcoxon test was used to compare the δ¹³C and δ¹⁵N values of quagga and zebra mussels in each treatment. We used Bayesian standard ellipse area (SEA) and size-corrected standard ellipse area for small samples (SEAc) to describe species and community isotopic niche width of newly settled macroinvertebrates among different treatments. We also calculated and reported the 95% confidence intervals of SEA using 10,000 posterior Markov chain Monte Carlo draws (Jackson et al., 2011; Catry et al., 2016).

All calculations and figures were conducted using R (R Core Team, 2021) with the package readxl, tidyverse, DescTools, ggplot2, ggpubr, vegan, and SIBER, in order to import data, reshape data, conduct Dunn’s test, make figures, arranges figures, conduct NMDS analysis, and analyze Bayesian standard ellipse areas, respectively.

Results

The Abundance of Newly Settled Macroinvertebrates

The most abundant taxa in all the four treatments were mussels [D. rostriformis bugensis, D. polymorpha, and Dreissena spp. (≤ 1 mm)], snails (P. antipodarum and B. tentaculata), gammarids (D. villosus), mysids (K. warpachowskyi), mayfly nymphs (Caenis spp.), and chironomid larvae (Chironominae and Orthocladiinae) (Figures 1, 2). Other taxa, such as Trichoptera, that occur infrequently were used only for total abundance analysis.

FIGURE 2

Figure 2. Non-metric multidimensional scaling (NMDS) ordination plot of newly settled macroinvertebrate abundances, representing the community structure in four different treatments. Ellipses draw standard errors of the treatments and show a 95% confidence interval. Db means D. rostriformis bugensis, Dp means D. polymorpha, Dr means Dreissena (≤ 1 mm), Pa means P. antipodarum, Bt means B. tentaculata, Dv means D. villosus, Kw means K. warpachowskyi, Ca means Caenis, Ch means Chironominae, and Or means Orthocladiinae.

Differences in the abundances of newly settled macroinvertebrates among the four treatments are presented in Table 1. The abundances of newly settled quagga mussels were not significantly different. The abundances of newly settled zebra mussels were also not significantly different. The abundances of newly settled snail P. antipodarum and mayfly nymphs Caenis were significantly different (p < 0.01), with the B treatment having the lowest numbers. The abundances of newly settled gammarid D. villosus were less in B treatments than in Q and Z treatments (p < 0.05). The abundances of newly settled chironomid larvae of subfamily Chironominae were higher in Z and B treatments than in S and Q treatments (p < 0.05). Total abundances of all newly settled macroinvertebrates were significantly higher in Z treatments than in B treatments (p = 0.002).

TABLE 1

Table 1. The differences in abundances of newly settled macroinvertebrates among four different treatments.

The total abundance of quagga mussels in the entire experiment area was higher than the total abundance of zebra mussels. In Z treatments, the abundance of newly settled quagga mussels was significantly higher than that of zebra mussels (p = 0.021), and additional information is provided in Supplementary Table 1.

In the NMDS plot (stress 0.09), the average composition of colonizing macroinvertebrate community across treatments appeared to be most segregated by mussel treatments (ANOSIM global test R = 0.473, p = 0.001) (Figure 2). For this 28-day experiment, Q and S treatments were closely grouped and separated from the Z treatment, while the B treatment was more isolated on the plot.

The mortality rate of glued adult zebra mussels (24.00 ± 2.28%) was significantly higher than that of glued adult quagga mussels (7.33 ± 1.33%) (p < 0.01; Figure 3). The abundance of newly settled chironomid larvae (subfamily Chironominae) showed a positive correlation with the mortality rate of adult dreissenids (Figure 3).

FIGURE 3

Figure 3. (A) Mortality rates of adult living quagga and zebra mussels glued on the stone substrates. D.rostriformis bugensis is the adult mussel used in Q treatments and D. ploymorpha is the adult mussel used in Z treatments. (B) Linear regression (solid line) and the 95% confidence interval of the regression (dotted line) for abundances of newly settled chironomid larvae (subfamily Chironominae) and mortality rates of adult dreissenids.

The Stable Isotope Values of Newly Settled Macroinvertebrates

The stable isotope values of five main macroinvertebrate species are shown in the stable isotope (δ¹³C and δ¹⁵N) bi-plots (Figure 4). Confidence intervals of the standard ellipse areas for different treatments and the whole experimental area are shown in the density plot (Figure 5 and Supplementary Figure 1).

FIGURE 4

Figure 4. The stable isotope bi-plots (δ¹³C and δ¹⁵N) of five main newly settled macroinvertebrates in four different treatments. Error bars indicate standard errors.

FIGURE 5

Figure 5. Density plot showing the confidence intervals of the standard ellipse area for five main newly settled macroinvertebrates and the whole community in four different treatments (B, S, Q, and Z). Red crosses correspond to the sample size-corrected standard ellipse areas (SEAc). Boxes shaded from dark to light gray represent the 50, 75, and 95% confidence intervals.

Among four treatments, the δ¹³C values of newly settled quagga mussels were significantly higher in the B treatment than in the Z treatment (p = 0.016; Table 2). The δ¹³C values of newly settled snail P. antipodarum were significantly higher in Z treatment than in S treatment (p = 0.014; Table 2). The other three species showed no significant differences in δ¹³C values among the four treatments. For all five main species, there were no significant differences in δ¹⁵N values among the four treatments.

TABLE 2

Table 2. The differences in stable isotope values of newly settled macroinvertebrates among four different treatments.

The SEAc results indicated that the newly settled quagga mussels showed a smaller isotopic niche in Z treatment than in B, S, and Q treatments (Figure 5 and Table 3). The newly settled zebra mussels also exhibited a smaller isotopic niche in Z treatment. The remaining three newly settled (frequent) macroinvertebrate species showed equivalent isotopic niches among the different treatments. Taken together, SEAc for the whole community was higher in B treatment.

TABLE 3

Table 3. Sample size-corrected standard ellipse areas (SEAc) of five newly settled macro-invertebrates and whole community from the four different treatments.

For the two dreissenids, the δ¹³C values of newly settled quagga mussels were higher than the δ¹³C values of newly settled zebra mussels. This pattern was evident in all four treatments (p ≤ 0.01; Supplementary Table 2). In each treatment, the δ¹⁵N values of all newly settled quagga mussels showed no significant differences from zebra mussels. Based on SEAc results, there were no overlaps in the isotopic niches between quagga and zebra mussels in each treatment.

Discussion

The Community Composition of Newly Settled Macroinvertebrates

The reaction of newly settled macroinvertebrates toward the stone substrate treatments complimented the hypothesis of facilitative associations between invasive dreissenids (Ward and Ricciardi, 2007; Ozersky et al., 2011), wherein one species has a positive effect on the establishment of another species. For example, the highest mean density of macroinvertebrates was consistently found in substrate treatments with shells or living mussels. By implication, the presence of mussels or shells was driving significant variation across experiments. Given the uniform in situ physicochemical conditions for stone substrate, this variation could be attributed to differences in the presence or absence of zebra and quagga mussels and their respective shells. Indeed, such effects were anticipated, since macroinvertebrates have shown a positive reaction to the physical structure of mussel shells (Ricciardi et al., 1997; Burlakova et al., 2012).

The most frequent newly settled macroinvertebrate taxa in our study, P. antipodarum, D. villosus, and Caenis spp., were mainly present at comparatively higher densities on stone substrate shells (S) and living mussel treatments (Q and Z). The rarity of settling macroinvertebrates on blank stone treatments (B) clearly indicates the feasible benefit of macroinvertebrates from the effects provided. In all three cases, each taxon could indeed benefit from the increased complexity of the substratum and surface area associated with natural dreissenid beds (Stewart et al., 1998; Ward and Ricciardi, 2007). Briefly, for P. antipodarum, shells may provide a bed substratum on which tubular refuges with protrusion could be constructed for grazing and scraping on detritus and sediments (Broekhuizen et al., 2001). In particular, the zebra mussel has shown mutualistic interactions with killer shrimp D. villosus and can benefit them in various ways. As shown during ex situ laboratory experiments on dreissenids, killer shrimp evidently utilizes mussel beds as protection from fish and allopatric predators more efficiently than other gammarids (Kobak et al., 2014). Similarly, juvenile zebra mussels can attach to the hard cover of the killer shrimp which could facilitate their dispersal (Yohannes et al., 2017), while the production of feces and pseudo-feces by zebra mussels may provide food (Klerks et al., 1996; Gergs and Rothhaupt, 2008). For Caenis, the taxon’s size might have allowed the use of the physical structure of dreissenids as nearby refuges through the exploitation of small interstices that are typical characteristics of dreissenid beds (Werner and Rothhaupt, 2007).

Like most of the three numerically dominant species, the abundances of newly settled chironomid larvae of subfamily Chironominae were higher in response to the presence of living zebra mussels (Z treatments). However, this species was one of the few species that showed higher density in the absence of mussels (B treatment). Almost all the newly settled chironomid larvae in B treatment were rather small with the average width of head measuring less than 180 μm. The observed pattern is presumably linked to the natural temporal or seasonal life history patterns (in growth, size, and abundance of the species), and perhaps this response was evident due to the 4 weeks of methodological effect. The overall pattern of chironomid larvae abundance and richness from our experimental study reinforces previous findings that zebra mussels facilitate the growth of benthic macroinvertebrate communities. This is in complement with previous studies that showed that the abundance of chironomid larvae increased after the colonization of zebra mussels (Botts et al., 1996; Gergs and Rothhaupt, 2008). Nevertheless, the abundance of chironomid larvae in Z treatments was higher than in Q treatments (Table 1). The high mortality rate of glued adult zebra mussels than quagga mussels might contribute to the differences, as the animal remains are part of the diet for chironomid larvae (Johnson, 1987).

The results of our analysis using analog shell and live quagga mussels should well be noted again, wherein no significant difference in the abundance of all taxa between S treatments and Q treatments was observed. This finding suggests that macroinvertebrates might primarily respond to the physical structure of quagga mussel shells, as discussed in previous investigations (Botts et al., 1996; Burlakova et al., 2012).

Although zebra mussels invaded the lake first, the rapid expansion in the range of the invasive quagga mussel in Lake Constance has caused a dominance shift from zebra to quagga mussels. However, specific effects of live mussels could offer additive benefits for multiple taxa, including newly settled conspecific zebra and quagga mussels. In fact, taking the whole experimental setup into consideration, the entire abundance of newly settled quagga mussels was higher than the entire abundance of newly settled zebra mussels. These data coupled with the current field observation in Lake Constance indicate that this and/or other similar facilitative invasions do cause major changes in biodiversity and ecosystem. The higher mortality rate of adult zebra mussels than adult quagga mussels may also contribute to the dominance shift. We also observed, based on a laboratory study (own unpublished data), that the quagga mussel has a faster growth rate than the zebra mussel.

Our experiment was conducted in 2018, only 2 years after the initial record of quagga mussels in the lake. Nevertheless, as shown by prior studies in other invaded systems, such initial reactions of macroinvertebrate communities toward dreissenids might reduce after several years (Strayer and Malcom, 2006; Karatayev et al., 2015). At the current particular population levels (5,708 ind./m²), we may conclude that quagga mussels would significantly impact the structure of the benthic community in this system. Indeed, this initial establishment of the species at this littoral site could be considered surprising; however, the highest densities of invasive populations have been recorded in higher depths of North American lakes, such as 342,000, 75,000, and 16,000 ind./m² in Lakes Erie, Huron, and Michigan, respectively (Nalepa, 1995; Howell et al., 1996; Nalepa et al., 2009).

Food Resources of Newly Settled Macroinvertebrates

We were able to note significant differences in the food resources of newly settled macroinvertebrates using mainly δ¹³C values (Layman et al., 2012). Among the four treatments applied, the newly settled quagga mussels in B treatments showed higher δ¹³C values than in Z treatments, which was probably related to the smaller size of newly settled quagga mussels in B treatments (p < 0.001). The newly settled snail P. antipodarum is known as a grazer and feeds on periphyton (Larson and Ross Black, 2016), and as expected, showed higher δ¹³C values in Z treatments than in S treatments, thus the higher mortality of the adult zebra mussels in Z treatments may have generated a different nutrient resource. Comparing Q treatments and the other treatments, there were no significant differences in δ¹³C and δ¹⁵N values of newly settled macroinvertebrates. Previous studies by Gergs et al. (2011) found that the stable isotope values of gammarids were significantly and positively related to the biodeposition of zebra mussels. Here, the presence of the adult quagga mussels did not show a strong influence on the isotope values of the newly settled macroinvertebrates, which might be related to the smaller study area (experiment) and shorter study (invasion) period than the previous study.

Nonetheless, the newly settled quagga mussels had significantly higher δ¹³C values than the newly settled zebra mussels in each treatment. Quagga and zebra mussels attained from shallow depth of Lake Constance in 2018 also showed higher δ¹³C values of quagga mussels than that of zebra mussels (own unpublished data). The diet of dreissenids can be very flexible, and contains phytoplankton and suspended detritus (Garton et al., 2005). The benthic primary producers normally have a higher δ¹³C value compared to pelagic primary producers (more negative) (France, 1995; Mitchell et al., 1996). In our study period, the newly settled quagga mussels might consume a higher percentage of benthic food or a lower percentage of pelagic food, in comparison to newly settled zebra mussels. A similar result showed that after the invasion by smallmouth bass and rock bass, lake trout switched to consuming a lower percentage of littoral prey and had a lower δ¹³C value than the reference lakes (vander Zanden et al., 1999). One study found that the δ¹³C values of quagga mussels are more negative than those of zebra mussels (Verhofstad et al., 2013), and the explanation was that quagga mussels consumed more chemoautotrophs which have highly reduced δ¹³C values. The difference between our study and the study of Verhofstad et al. (2013) may be due to the differences in the substrates, depth ranges, or measurements. Some other studies did not find any difference in the stable isotope values between quagga and zebra mussels (Garton et al., 2005).

No overlaps in the isotopic niche of quagga and zebra mussels were observed in our field experiment. The lengths of newly settled quagga and zebra mussels in S, Q, and Z treatments did not show significant differences, thus the differences in stable isotope values of quagga and zebra mussels were not caused by the differences in the size of dreissenids. At our study site where quagga and zebra mussels coexisted, the dietary segregation or partition of colonized space may be the mechanism for the co-existence of quagga and zebra mussels (Diggins et al., 2004; Rothhaupt et al., 2014). Our experiment was conducted in 2018 in Lake Constance, only 2 years after the first record of quagga mussels. During the early colonization period, quagga mussels might occupy a different niche from zebra mussels to reduce resource competition with zebra mussels (Jackson and Britton, 2014).

Conclusion

In our study, the abundance of gastropods, gammarids, and some insect larvae increased with the presence of dreissenid shells in the littoral region. The macroinvertebrate communities might change after quagga mussels replace zebra mussels and occupy a high percentage in the total biomass of benthic macroinvertebrates, particularly in deep regions where adult quagga mussels can also live (Mills et al., 1993; Nalepa et al., 2020). After the invasion of dreissenids in the Great Lakes, the abundance of dressenids increased in the profundal region; however, the abundance of non-dreissenid profundal macroinvertebrates decreased, which is different from the changes in the littoral region (Burlakova et al., 2018). In our study, there was apparent niche partitioning between newly settled quagga and zebra mussels. The interesting patterns of dietary segregation between quagga and zebra mussels might be a common phenomenon in the early stage of quagga mussel invasions when the two dreissenids coexisted in lakes. These results provide a benchmark by which future changes in trophic ecology and invasion dynamics can be measured across the ecosystem. However, for future experiments, we recommend a longer stone substrate deployment duration encompassing standardized seasonal periods for individual species to achieve stabilized assemblages of wider coverage of macroinvertebrates. Yet, this study offers the first overview of the quagga mussel’s progressive invasion and the reaction of the zebra mussels and other newly settled macroinvertebrates, and contributes to the hypothesis of facilitative associations between invasive mussel species.

In summary, following the occurrence of quagga mussels in Lake Constance, significant changes in the benthic macroinvertebrate community are observed. Only about 5 years after the first detection of quagga mussel, it has become the most dominant macroinvertebrate in the littoral area of the lake. A recent study that analyzed the stomach content showed that benthic whitefish (Coregonus macrophthalmus), roach (Rutilus rutilus), and tench (Tinca tinca) showed high levels of consumption of quagga mussels (Baer et al., 2022). Although equivalent and empirical data on the actual impact of the invasive mussels on the accompanying macroinvertebrates are not available, these dietary shifts observed in multiple benthic dwelling fish species imply alternation in the dietary sources following quagga mussel colonization in Lake Constance.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

K-OR, EY, and HZ designed the field experiment. K-OR secured the funding for the project. HZ conducted the field experiment, counted the macroinvertebrate samples, and prepared the samples for stable isotope analysis. HZ analyzed the data with guidance from EY and K-OR. HZ and EY wrote the draft. All authors revised the draft, read, and approved the submitted version.

Funding

This project was supported by funding from the Limnological Institute, University of Konstanz (AFF grant to K-OR). HZ received funding from the Chinese Scholarship Council (grant no. 201704910903) for a living stipend.

Conflict of Interest

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank Christian Fiek for his great help with conducting the experiment and counting macroinvertebrates. Ioanna Salvarina also provided guidance for counting macro-invertebrates. We would also like to thank Wolfgang Kornberger, Yao Zhiyijun, Christian Hanke, and Robert van Geldern for their assistance in the stable isotope analysis of the macroinvertebrates and the members of SeeWandel: Life in Lake Constance – the past, present and future for valuable discussions about the experiment. Martin Wolf also helped with the preparation of stone substrates for the experiment.

Supplementary Material

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

References

Baer, J., Spiessl, C., Auerswald, K., Geist, J., and Brinker, A. (2022). Signs of the times: isotopic signature changes in several fish species following invasion of Lake Constance by quagga mussels. J. Great Lakes Res. 48, 746–755. doi: 10.1016/j.jglr.2022.03.010

CrossRef Full Text | Google Scholar

Bially, A., and MacIsaac, H. J. (2000). Fouling mussels (Dreissena spp.) colonize soft sediments in Lake Erie and facilitate benthic invertebrates. Freshwater Biol. 43, 85–97. doi: 10.1046/j.1365-2427.2000.00526.x

CrossRef Full Text | Google Scholar

Boltovskoy, D., and Correa, N. (2015). Ecosystem impacts of the invasive bivalve Limnoperna fortunei (golden mussel) in South America. Hydrobiologia 746, 81–95. doi: 10.1007/s10750-014-1882-1889

CrossRef Full Text | Google Scholar

Botts, P. S., Patterson, B. A., and Schloesser, D. W. (1996). Zebra mussel effects on benthic invertebrates: physical or biotic? J. North Am. Benthol. Soc. 15, 179–184. doi: 10.2307/1467947

CrossRef Full Text | Google Scholar

Broekhuizen, N., Parkyn, S., and Miller, D. (2001). Fine sediment effects on feeding and growth in the invertebrate grazers Potamopyrgus antipodarum (Gastropoda, Hydrobiidae) and Deleatidium sp. (Ephemeroptera, Leptophlebiidae). Hydrobiologia 457, 125–132. doi: 10.1023/A:1012223332472

CrossRef Full Text | Google Scholar

Burlakova, L. E., Barbiero, R. P., Karatayev, A. Y., Daniel, S. E., Hinchey, E. K., and Warren, G. J. (2018). The benthic community of the Laurentian Great Lakes: analysis of spatial gradients and temporal trends from 1998 to 2014. J. Great Lakes Res. 44, 600–617. doi: 10.1016/j.jglr.2018.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Burlakova, L. E., Karatayev, A. Y., and Karatayev, V. A. (2012). Invasive mussels induce community changes by increasing habitat complexity. Hydrobiologia 685, 121–134. doi: 10.1007/s10750-011-0791-794

CrossRef Full Text | Google Scholar

Burlakova, L. E., Tulumello, B. L., Karatayev, A. Y., Krebs, R. A., Schloesser, D. W., Paterson, W. L., et al. (2014). Competitive replacement of invasive congeners may relax impact on native species: interactions among zebra, quagga, and native unionid mussels. PLoS One 9:e114926. doi: 10.1371/journal.pone.0114926

PubMed Abstract | CrossRef Full Text | Google Scholar

Byers, J. E., Reichard, S., Randall, J. M., Parker, I. M., Smith, C. S., Lonsdale, W. M., et al. (2002). Directing research to reduce the impacts of nonindigenous species. Conservation Biol. 16, 630–640. doi: 10.1046/j.1523-1739.2002.01057.x

CrossRef Full Text | Google Scholar

Catry, T., Lourenço, P. M., Lopes, R. J., Carneiro, C., Alves, J. A., Costa, J., et al. (2016). Structure and functioning of intertidal food webs along an avian flyway: a comparative approach using stable isotopes. Funct. Ecol. 30, 468–478. doi: 10.1111/1365-2435.12506

CrossRef Full Text | Google Scholar

Copp, G. H., Vilizzi, L., Mumford, J., Fenwick, G. V., Godard, M. J., and Gozlan, R. E. (2009). Calibration of FISK, an invasiveness screening tool for nonnative freshwater fishes. Risk Anal. 29, 457–467. doi: 10.1111/j.1539-6924.2008.01159.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Crooks, J. A. (2005). Lag times and exotic species: the ecology and management of biological invasions in slow-motion. Ecoscience 12, 316–329. doi: 10.2980/i1195-6860-12-3-316.1

CrossRef Full Text | Google Scholar

Diggins, T. P., Weimer, M., Stewart, K. M., Baier, R. E., Meyer, A. E., Forsberg, R. F., et al. (2004). Epiphytic refugium: are two species of invading freshwater bivalves partitioning spatial resources? Biol. Invasions 6, 83–88.

Google Scholar

Evariste, L., David, E., Cloutier, P. L., Brousseau, P., Auffret, M., Desrosiers, M., et al. (2018). Field biomonitoring using the zebra mussel Dreissena polymorpha and the quagga mussel Dreissena bugensis following immunotoxic reponses. is there a need to separate the two species? Environ. Pollut. 238, 706–716. doi: 10.1016/j.envpol.2018.03.098

PubMed Abstract | CrossRef Full Text | Google Scholar

France, R. L. (1995). Differentiation between littoral and pelagic food webs in lakes using stable carbon isotopes. Limnol. Oceanography 40, 1310–1313. doi: 10.4319/lo.1995.40.7.1310

CrossRef Full Text | Google Scholar

Francis, R. A., and Chadwick, M. A. (2012). “Invasive alien species in freshwater ecosystems: a brief overview,” in A Handbook of Global Freshwater Invasive Species, ed. R. A. Francis (Abingdon: Earthscan), 3–21. doi: 10.1093/icb/icaa023

PubMed Abstract | CrossRef Full Text | Google Scholar

Garton, D. W., Payne, C. D., and Montoya, J. P. (2005). Flexible diet and trophic position of dreissenid mussels as inferred from stable isotopes of carbon and nitrogen. Canadian J. Fisheries Aquatic Sci. 62, 1119–1129. doi: 10.1139/f05-025

CrossRef Full Text | Google Scholar

Gergs, R., and Rothhaupt, K. O. (2008). Effects of zebra mussels on a native amphipod and the invasive Dikerogammarus villosus: the influence of biodeposition and structural complexity. J. North Am. Benthol. Soc. 27, 541–548. doi: 10.1899/07-151.1

CrossRef Full Text | Google Scholar

Gergs, R., and Rothhaupt, K. O. (2015). Invasive species as driving factors for the structure of benthic communities in Lake Constance, Germany. Hydrobiologia 746, 245–254. doi: 10.1007/s10750-014-1931-1934

CrossRef Full Text | Google Scholar

Gergs, R., Grey, J., and Rothhaupt, K. O. (2011). Temporal variation in zebra mussel (Dreissena polymorpha) density structure the benthic food web and community composition on hard substrates in Lake Constance, Germany. Biol. Invasions 13, 2727–2738. doi: 10.1007/s10530-011-9943-9948

CrossRef Full Text | Google Scholar

Gherardi, F., and Acquistapace, P. (2007). Invasive crayfish in Europe: the impact of Procambarus clarkii on the littoral community of a Mediterranean lake. Freshwater Biol. 52, 1249–1259. doi: 10.1111/j.1365-2427.2007.01760.x

CrossRef Full Text | Google Scholar

González, M. J., and Downing, A. (1999). Mechanisms underlying amphipod responses to zebra mussel (Dreissena polymorpha) invasion and implications for fish-amphipod interactions. Canadian J. Fisheries Aquatic Sci. 56, 679–685. doi: 10.1139/f98-211

CrossRef Full Text | Google Scholar

Haltiner, L., Zhang, H., Anneville, O., de Ventura, L., DeWeber, T. J., Hesselschwerdt, J., et al. (2022). The distribution and spread of quagga mussels in perialpine lakes north of the Alps. Aquatic Invasions 17, 153–173. doi: 10.3391/ai.2022.17.2.02

CrossRef Full Text | Google Scholar

Howell, E. T., Marvin, C. H., Bilyea, R. W., Kauss, P. B., and Somers, K. (1996). Changes in environmental conditions during Dreissena colonization of a monitoring station in eastern Lake Erie. J. Great Lakes Res. 22, 744–756. doi: 10.1016/S0380-1330(96)70993-70990

CrossRef Full Text | Google Scholar

Jackson, A. L., Inger, R., Parnell, A. C., and Bearhop, S. (2011). Comparing isotopic niche widths among and within communities: SIBER - Stable Isotope Bayesian Ellipses in R. J. Animal Ecol. 80, 595–602. doi: 10.1111/j.1365-2656.2011.01806.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, M. C., and Britton, J. R. (2014). Divergence in the trophic niche of sympatric freshwater invaders. Biol. Invasions 16, 1095–1103. doi: 10.1007/s10530-013-0563-563

CrossRef Full Text | Google Scholar

Johnson, R. K. (1987). Seasonal variation in diet of Chironomus plumosus (L.) and C.anthracinus Zett. (Diptera-Chironomidae) in mesotrophic Lake Erken. Freshwater Biol. 17, 525–532. doi: 10.1111/j.1365-2427.1987.tb01073.x

CrossRef Full Text | Google Scholar

Jones, P. E., Tummers, J. S., Galib, S. M., Woodford, D. J., Hume, J. B., Silva, L. G. M., et al. (2021). The use of barriers to limit the spread of aquatic invasive animal species: a global review. Front. Ecol. Evol. 9:611631. doi: 10.3389/fevo.2021.611631

CrossRef Full Text | Google Scholar

Karatayev, A. Y., Burlakova, L. E., and Padilla, D. K. (2015). Zebra versus quagga mussels: a review of their spread, population dynamics, and ecosystem impacts. Hydrobiologia 746, 97–112. doi: 10.1007/s10750-014-1901-x

CrossRef Full Text | Google Scholar

Karatayev, A. Y., Karatayev, V. A., Burlakova, L. E., Mehler, K., Rowe, M. D., Elgin, A. K., et al. (2021). Lake morphometry determines Dreissena invasion dynamics. Biol. Invasions 23, 2489–2514. doi: 10.1007/s10530-021-02518-2513

CrossRef Full Text | Google Scholar

Kerambrun, E., Delahaut, L., Geffard, A., and David, E. (2018). Differentiation of sympatric zebra and quagga mussels in ecotoxicological studies: a comparison of morphometric data, gene expression, and body metal concentrations. Ecotoxicol. Environ. Saf. 154, 321–328. doi: 10.1016/j.ecoenv.2018.02.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Klerks, P., Fraleigh, P. C., and Lawniczak, E. (1996). Effects of zebra mussels (Dreissena polymorpha) on seston levels and sediment deposition in western Lake Erie. Canadian J. Fisheries Aquatic Sci. 53, 2284–2291. doi: 10.1139/f96-190

CrossRef Full Text | Google Scholar

Kobak, J., Jermacz, Ł, and Płąchocki, D. (2014). Effectiveness of zebra mussels to act as shelters from fish predators differs between native and invasive amphipod prey. Aquatic Ecol. 48, 397–408. doi: 10.1007/s10452-014-9492-9491

CrossRef Full Text | Google Scholar

Kobak, J., Kakareko, T., Jermacz, Ł, and Poznańska, M. (2013). The impact of zebra mussel (Dreissena polymorpha) periostracum and biofilm cues on habitat selection by a Ponto-Caspian amphipod Dikerogammarus haemobaphes. Hydrobiologia 702, 215–226. doi: 10.1007/s10750-012-1322-1327

CrossRef Full Text | Google Scholar

Kulhanek, S. A., Leung, B., and Ricciardi, A. (2011). Using ecological niche models to predict the abundance and impact of invasive species: application to the common carp. Ecol. Appl. 21, 203–213. doi: 10.1890/09-1639.1

CrossRef Full Text | Google Scholar

Larson, M. D., and Ross Black, A. (2016). Assessing interactions among native snails and the invasive New Zealand mud snail, Potamopyrgus antipodarum, using grazing experiments and stable isotope analysis. Hydrobiologia 763, 147–159. doi: 10.1007/s10750-015-2369-z

CrossRef Full Text | Google Scholar

Layman, C. A., Araujo, M. S., Boucek, R., Hammerschlag-Peyer, C. M., Harrison, E., Jud, Z. R., et al. (2012). Applying stable isotopes to examine food-web structure: an overview of analytical tools. Biol. Rev. 87, 545–562. doi: 10.1111/j.1469-185X.2011.00208.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lodge, D. M., Stein, R. A., Brown, K. M., Covich, A. R., Christer, B., Garvey, J. E., et al. (1998). Predicting impact of freshwater exotic species on native biodiversity: challenges in spatial scaling. Australian J. Ecol. 23, 53–67. doi: 10.1111/j.1442-9993.1998.tb00705.x

CrossRef Full Text | Google Scholar

MacIsaac, H. J. (1996). Potential abiotic and biotic impacts of zebra mussels on the inland waters of North America. Am. Zool. 36, 287–299. doi: 10.1093/icb/36.3.287

PubMed Abstract | CrossRef Full Text | Google Scholar

Mills, E. L., Dermott, R. M., Roseman, E. F., Dustin, D., Mellina, E., Conn, D. B., et al. (1993). Colonization, ecology, and population structure of the “quagga” mussel (Bivalvia: Dreissenidae) in the Lower Great Lakes. Canadian J. Fisheries Aquatic Sci. 50, 2305–2314. doi: 10.1139/f93-255

CrossRef Full Text | Google Scholar

Mitchell, M. J., Mills, E. L., Idrisi, N., and Michener, R. (1996). Stable isotopes of nitrogen and carbon in an aquatic food web recently invaded by Dreissena polymorpha (Pallas). Canadian J. Fisheries Aquatic Sci. 53, 1445–1450. doi: 10.1139/f96-053

CrossRef Full Text | Google Scholar

Molloy, D. P., Bij, De Vaate, A., Wilke, T., and Giamberini, L. (2007). Discovery of Dreissena rostriformis bugensis (Andrusov 1897) in Western Europe. Biol. Invasions 9, 871–874. doi: 10.1007/s10530-006-9078-9075

CrossRef Full Text | Google Scholar

Mörtl, M., and Rothhaupt, K. O. (2003). Effects of adult Dreissena polymorpha on settling juveniles and associated macroinvertebrates. Int. Rev. Hydrobiol. 88, 561–569. doi: 10.1002/iroh.200310640

CrossRef Full Text | Google Scholar

Nalepa, T. F. (1995). Initial colonization of the zebra mussel (Dreissena polymorpha) in Sagniaw Bay, Lake Huron: population recruitment, density, and size structure. J. Great Lakes Res. 21, 417–434.

Google Scholar

Nalepa, T. F., Burlakova, L. E., Elgin, A. K., and Karatayev, A. Y. (2020). Abundance and biomass of benthic macroinvertebrates in Lake Michigan in 2015, with a summary of temporal trends. NOAA Techn. Memorandum GLERL 175, 1–41. doi: 10.25923/g0d3-3v41

CrossRef Full Text | Google Scholar

Nalepa, T. F., Fanslow, D. L., and Lang, G. A. (2009). Transformation of the offshore benthic community in Lake Michigan: recent shift from the native amphipod Diporeia spp. to the invasive mussel Dreissena rostriformis bugensis. Freshwater Biol. 54, 466–479. doi: 10.1111/j.1365-2427.2008.02123.x

CrossRef Full Text | Google Scholar

Ozersky, T., Barton, D. R., and Evans, D. O. (2011). Fourteen years of dreissenid presence in the rocky littoral zone of a large lake: effects on macroinvertebrate abundance and diversity. J. North Am. Benthol. Soc. 30, 913–922. doi: 10.1899/10-122.1

CrossRef Full Text | Google Scholar

Pathy, D. A., and Mackie, G. L. (1993). Comparative shell morphology of Dreissena polymorpha, Mytilopsis leucophaeata, and the “quagga” mussel (Bivalvia: Dreissenidae) in North America. Canadian J. Zool. 71, 1012–1023. doi: 10.1139/z93-135

CrossRef Full Text | Google Scholar

Post, D. M., Layman, C. A., Arrington, D. A., Takimoto, G., Quattrochi, J., and Montaña, C. G. (2007). Getting to the fat of the matter: models, methods and assumptions for dealing with lipids in stable isotope analyses. Oecologia 152, 179–189. doi: 10.1007/s00442-006-0630-x

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2021). R: A Language and Environmental for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Ricciardi, A. (2003). Predicting the impacts of an introduced species from its invasion history: an empirical approach applied to zebra mussel invasions. Freshwater Biol. 48, 972–981. doi: 10.1046/j.1365-2427.2003.01071.x

CrossRef Full Text | Google Scholar

Ricciardi, A., Hoopes, M. F., Marchetti, M. P., and Lockwood, J. L. (2013). Progress toward understanding the ecological impacts of nonnative species. Ecol. Monographs 83, 263–282. doi: 10.1890/13-0183.1

CrossRef Full Text | Google Scholar

Ricciardi, A., Whoriskey, F. G., and Rasmussen, J. B. (1997). The role of the zebra mussel (Dreissena polymorpha) in structuring macroinvertebrate communities on hard substrata. Canadian J. Fisheries Aquatic Sci. 54, 2596–2608. doi: 10.1139/f97-174

CrossRef Full Text | Google Scholar

Roby, K. B., Newbold, J. D., and Erman, D. C. (1978). Effectiveness of an artificial substrate for sampling macroinvertebrates in small streams. Freshwater Biol. 8, 1–8. doi: 10.1111/j.1365-2427.1978.tb01420.x

CrossRef Full Text | Google Scholar

Rothhaupt, K. O., Hanselmann, A. J., and Yohannes, E. (2014). Niche differentiation between sympatric alien aquatic crustaceans: an isotopic evidence. Basic Appl. Ecol. 15, 453–463. doi: 10.1016/j.baae.2014.07.002

CrossRef Full Text | Google Scholar

Roy, H. E., Peyton, J., Aldridge, D. C., Bantock, T., Blackburn, T. M., Britton, R., et al. (2014). Horizon scanning for invasive alien species with the potential to threaten biodiversity in Great Britain. Global Change Biol. 20, 3859–3871. doi: 10.1111/gcb.12603

PubMed Abstract | CrossRef Full Text | Google Scholar

Simberloff, D., and von Holle, B. (1999). Positive interactions of nonindigenous species: invasional meltdown? Biol. Invasions 1, 21–32. doi: 10.1023/A:1010086329619

CrossRef Full Text | Google Scholar

Son, M. O. (2007). Native range of the zebra mussel and quagga mussel and new data on their invasions within the Ponto-Caspian Region. Aquatic Invasions 2, 174–184. doi: 10.3391/ai.2007.2.3.4

CrossRef Full Text | Google Scholar

Sousa, R., Pilotto, F., and Aldridge, D. C. (2011). Fouling of European freshwater bivalves (Unionidae) by the invasive zebra mussel (Dreissena polymorpha). Freshwater Biol. 56, 867–876. doi: 10.1111/j.1365-2427.2010.02532.x

CrossRef Full Text | Google Scholar

Stewart, T. W., Miner, J. G., and Lowe, R. L. (1998). Quantifying mechanisms for zebra mussel effects on benthic macroinvertebrates: organic matter production and shell-generated habitat. J. North Am. Benthol Soc. 17, 81–94. doi: 10.2307/1468053

CrossRef Full Text | Google Scholar

Stiers, I., Crohain, N., Josens, G., and Triest, L. (2011). Impact of three aquatic invasive species on native plants and macroinvertebrates in temperate ponds. Biol. Invasions 13, 2715–2726. doi: 10.1007/s10530-011-9942-9949

CrossRef Full Text | Google Scholar

Strayer, D. L., and Malcom, H. M. (2006). Long-term demography of a zebra mussel (Dreissena polymorpha) population. Freshwater Biol. 51, 117–130. doi: 10.1111/j.1365-2427.2005.01482.x

CrossRef Full Text | Google Scholar

Strayer, D. L., and Malcom, H. M. (2007). Effects of zebra mussels (Dreissena polymorpha) on native bivalves: the beginning of the end or the end of the beginning? J. North Am. Benthol Soc. 26, 111–122.

Google Scholar

Strayer, D. L., and Malcom, H. M. (2018). Long-term responses of native bivalves (Unionidae and Sphaeriidae) to a Dreissena invasion. Freshwater Sci. 37, 697–711. doi: 10.1086/700571

CrossRef Full Text | Google Scholar

Strayer, D. L., Fischer, D. T., Hamilton, S. H., Malcom, H. M., Pace, M. L., and Solomon, C. T. (2020). Long-term variability and density dependence in Hudson River Dreissena populations. Freshwater Biol. 65, 474–489. doi: 10.1111/fwb.13444

CrossRef Full Text | Google Scholar

Tran, T. N. Q., Jackson, M. C., Sheath, D., Verreycken, H., and Britton, J. R. (2015). Patterns of trophic niche divergence between invasive and native fishes in wild communities are predictable from mesocosm studies. J. Animal Ecol. 84, 1071–1080. doi: 10.1111/1365-2656.12360

PubMed Abstract | CrossRef Full Text | Google Scholar

vander Zanden, M. J., Casselman, J. M., and Rasmussen, J. B. (1999). Stable isotope evidence for the food web consequences of species invasions in lakes. Nature 401, 464–467. doi: 10.1038/46762

CrossRef Full Text | Google Scholar

Verhofstad, M. J., Grutters, B. M. C., van der Velde, G., and Leuven, R. S. E. W. (2013). Effects of water depth on survival, condition and stable isotope values of three invasive dreissenid species in a deep freshwater lake. Aquatic Invasions 8, 157–169. doi: 10.3391/ai.2013.8.2.04

CrossRef Full Text | Google Scholar

Walz, N. (1973). Studies on the biology of Dreissena polymorpha in Lake Constance. Arch. Hydrobiol. Suppl. 42, 452–482.

Google Scholar

Ward, J. M., and Ricciardi, A. (2007). Impacts of Dreissena invasions on benthic macroinvertebrate communities: a meta-analysis. Diversity Distrib. 13, 155–165. doi: 10.1111/j.1472-4642.2007.00336.x

CrossRef Full Text | Google Scholar

Werner, S., and Rothhaupt, K.-O. (2007). Effects of the invasive bivalve Corbicula fluminea on settling juveniles and other benthic taxa. J. North Am. Benthol Soc. 26, 673–680. doi: 10.1899/07-017R.1

CrossRef Full Text | Google Scholar

Williamson, M. (1999). Invasions. Ecography 22, 5–12. doi: 10.1111/j.1600-0587.1999.tb00449.x

CrossRef Full Text | Google Scholar

Wise, D. H., and Molles, M. C. (1979). Colonization of artificial substrates by stream insects: influence of substrate size and diversity. Hydrobiologia 65, 69–74. doi: 10.1007/BF00032721

CrossRef Full Text | Google Scholar

Yohannes, E., Ragg, R. B., Armbruster, J. P., and Rothhaupt, K. O. (2017). Physical attachment of the invasive zebra mussel Dreissena polymorpha to the invasive gammarid Dikerogammarus villosus: supplementary path for invasion and expansion? Fundamental Appl. Limnol. 191, 79–85. doi: 10.1127/fal/2017/1063

CrossRef Full Text | Google Scholar