Research Article |
|
Corresponding author: Astrid N. Schwalb ( schwalb@txstate.edu ) Academic editor: David Wong
© 2025 Emily Lorkovic, Jason P. Martina, Monica E. McGarrity, Astrid N. Schwalb.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Lorkovic E, Martina JP, McGarrity ME, Schwalb AN (2025) Two invaders, one reservoir: Hydrilla shapes the distribution of zebra mussels and may facilitate their growth. Aquatic Invasions 20(3): 273-290. https://doi.org/10.3391/ai.2025.20.3.156675
|
Aquatic invasive species can alter ecosystem processes, detrimentally affect native species, and facilitate the invasion of other species. One infamous aquatic invader, the zebra mussel (Dreissena polymorpha), is known to cause declines in phytoplankton through their filtering activity and facilitate the subsequent growth of macrophytes by increasing water clarity. In turn, submerged macrophytes may provide substrate for settlement of zebra mussels. The goal of this study was to examine variation in the distribution of zebra mussels and hydrilla (Hydrilla verticillata subsp. verticillata) in relation to sediment composition, each other (including potential facilitation), and with season (summer vs. fall) in a subtropical reservoir. Surveys of zebra mussels and hydrilla showed that zebra mussel densities tended to be higher in rocky habitats where they were found on hydrilla and rocks (gravel and cobble), compared to muddy habitats where they were found only on hydrilla. Within the rocky habitat, zebra mussels attached to hydrilla had significantly higher densities and a smaller size than those attached to rocks. However, spring populations may be largely transient because only a small fraction of zebra mussels remained on hydrilla in early fall, almost exclusively representing a new settlement cohort based on their size distribution. Nevertheless, hydrilla may directly facilitate zebra mussel dispersal, especially in spring, as mussels attached to plant fragments can be transported downstream by currents or by human activities, such as entanglement in boat propellers and trailers. Laboratory experiments did not detect any significant impact of zebra mussels on the growth, biomass, or nutrient content of hydrilla. However, zebra mussel biomass was higher when hydrilla was present, suggesting that hydrilla may facilitate zebra mussel growth, although the difference was only statistically significant at low hydrilla densities. This study illustrates the complexities of interactions between multiple introduced species which can lead to facilitation of invasion of aquatic ecosystems.
Invader facilitation, invasional meltdown hypothesis, invasive species distribution, plant-mussel interactions, reservoir ecology
Modern human globalization has increased the rate of species introductions, leading to aquatic ecosystems being colonized by multiple non-indigenous and potentially invasive species (
Zebra mussels (Dreissena polymorpha) are considered one of the most problematic invasive species on multiple continents because of their large-scale negative ecosystem effects coupled with their ability to encrust submerged infrastructure (e.g., water supply and control), causing large economic costs (
Macrophytes can provide biotic surfaces for the attachment of zebra mussels, which may be especially important in areas with soft substrates, where zebra mussels otherwise would not be able to attach (
Research on interactions between zebra mussels and hydrilla is lacking. Hydrilla is an obligate submerged aquatic plant native to Southeast Asia that has spread to multiple continents, including North America, where established populations persist in eastern and southern river basins (
The goal of this study was to examine variation in the distribution of both zebra mussels and hydrilla in relation to substrate, each other, and season (i.e., summer vs. fall) and to examine potential interactive effects in laboratory experiments. Our study had four main objectives. First, we examined how substrate composition affects the distribution of zebra mussels and hydrilla within a large reservoir. Second, we examined differences in density and size distribution of mussels between those attached to rocks versus hydrilla in field surveys. Third, we examined these same differences over time between summer and fall. Zebra mussels experience high mortality when water temperatures are high in summer (
This study was conducted in Canyon Lake, a eutrophic reservoir with over 24 kilometers of shoreline located in the Guadalupe River basin in south-central Texas (Fig.
During each field survey, temperature (°C), dissolved oxygen (mg/L-1 and %), and specific conductance (μS cm-1) were measured with a multisonde (YSI model ProDDS). Water temperature averaged 29.3 ± 0.7 °C (mean ± SD) in June and 28.7 ± 0.9 °C in September (Table
Nested analysis of variance for benthic composition with site as a covariate for surveyed hydrilla and zebra mussel densities in June 2022.
| Source | df | F-value | p-value | Percent explained |
|---|---|---|---|---|
| Hydrilla | ||||
| Substrate | 1 | 24.36 | 2.39 e-05 | 0.13 |
| Substrate: Site | 6 | 21.85 | 4.90 e-10 | 0.70 |
| Residuals | 32 | 0.17 | ||
| Zebra Mussels | ||||
| Substrate | 1 | 19.62 | 1.03 e-04 | 0.22 |
| Substrate: Site | 6 | 6.549 | 1.39 e-04 | 0.43 |
| Residuals | 32 | 0.35 | ||
Hydrilla and zebra mussel surveys were carried out by scuba diving using five sampling quadrats (25 × 25 cm), each placed at two-meter intervals along 10-meter transects located parallel to the shore at a water depth between 2.7 and 4.0 m in June and 1.7 and 2.7 m in September. Inside each sampling quadrat, all hydrilla tissues present in the water column (i.e., from surface to lake bottom) were removed by hand using garden shears. Immediately after collection, hydrilla tissue samples were placed in 0.5-millimeter mesh bags. Rocks colonized by zebra mussels were sampled by collecting from 1/4 of the sampling quadrat (156 cm2) for rocky sites in June only due to time constraints. Plant and mussel samples were transported in aerated lake water, to reduce stress, to the laboratory and immediately frozen at -12 °C.
Hydrilla samples were thawed at room temperature for processing, after which zebra mussels were carefully removed from plant tissue and washed with deionized water. Subsampling was performed by removing zebra mussels from a random 1/3 of each plant sample (derived from one sampling quadrat). If fewer than 50 mussels were found, an extra 1/3 of the plant sample was inspected for more individuals. If a total of 50 mussels was still not reached in the subsamples, then the entire plant sample was inspected for individuals. All zebra mussels removed from hydrilla and collected from rocks and hydrilla samples were oven-dried at 60 °C for 24 and 48 hours, respectively, until a constant weight was achieved. Zebra mussel individuals were counted and measured with calipers to the nearest 0.1 mm, and hydrilla samples were cooled to room temperature and further weighed to determine their dry biomass.
A laboratory experiment was conducted to test the effects of zebra mussels on hydrilla growth and vice versa under controlled conditions. The experiment was carried out using a 3 × 2 fully crossed design which included three levels of hydrilla density (without hydrilla, low, or high) and two levels of zebra mussels (without mussels or with mussels) with five replicates of the six treatment combinations, except for the water only control, which had three replicates. Treatment combinations were: hydrilla high density + zebra mussels, hydrilla low density + zebra mussels, zebra mussels only, hydrilla only in low or high density, plus water only to control for water parameters. Each independent experimental unit consisted of a 1.5 L tank filled with filtered (50 μm mesh) lake water from Canyon Lake. The experiment was carried out for 16 days, from November 7th to November 23rd, 2022. Zebra mussels and hydrilla were collected from Canyon Lake at the start of the experiment near boat Ramp 23 (Fig.
Temperature was continuously measured using temperature loggers (HOBO pendant Temperature/Light 64K data logger) in four randomly selected tanks across the treatments. The average water temperature during lab experiments, measured continuously at 1-hour increments, was 20.6±1.4 °C over both light and dark periods. The experiment was conducted at ambient air temperature, selected as a moderate compromise to support the metabolic requirements and health of both zebra mussels and hydrilla. Photosynthetic active radiation (PAR) was calculated by regression analysis (y = 0.0138x, R2 = 0.677)from lux (lumens per m2) measured in the tanks using a PAR meter (LI-190R Quantum Sensor, LI-COR, Lincoln, NE, USA). Three Hygger Programmable Aquarium LED Lights (HG-957, 36 watts) were placed between each set of two rows and set to a 16:8 day/night schedule. Light intensity (PAR) at the center of the light fixture averaged 61.4±3.86 μmol/m²/s during daytime, which is well above sufficient light (30 μmol/m²/s) for photosynthesis in similar submerged macrophytes (
Normality of data was assessed using Shapiro-Wilk’s test, and homogeneity of variances was evaluated using Levene’s test. Several response variables exhibited normal distributions, including hydrilla biomass from surveys, hydrilla C/N ratios, hydrilla growth, and zebra mussel AFDM from the lab experiment. Square-root, logarithmic, or cube-root transformations were applied to non-normal data to meet the normality assumptions. Zebra mussel densities were cube-root transformed for both total datasets and subsets (i.e., June vs September, substrate rocks vs hydrilla). Pooled zebra mussel densities and lengths as well as experimental hydrilla biomass did not respond to transformation and were analyzed using non-parametric tests. To examine the effect of benthic composition type on hydrilla biomass and the densities of zebra mussels attached to hydrilla, a nested Analysis of Variance (ANOVA) was conducted, with site serving as a covariate. One-way ANOVA was employed to evaluate differences in hydrilla C/N ratios for zebra mussel AFDM for lab experiments. Tukey’s post-hoc tests were performed on all ANOVAs to test for significant differences among groups. A Pearson correlation coefficient was computed to examine whether there was a significant correlation between total zebra mussel densities and hydrilla biomass at rocky and muddy sites. Pairwise t-tests were used to compare summer and fall surveyed hydrilla biomass, assess differences between hydrilla-bound and rock-bound mussels at rocky sites, and evaluate treatment effects of zebra mussels on hydrilla change in biomass and growth in the field and laboratory experiment. Non-parametric Mann-Whitney U tests were employed to test for zebra mussel density changes between summer and fall, size differences between mussels found on hydrilla tissues and rocks, and pairwise comparisons between controls and associated treatments plus low-density and high-density hydrilla biomass in the laboratory experiment. All statistical analyses were performed in R version 4.1.0 (2021).
In the field surveys, substrate and site (nested) influenced the density of zebra mussels and biomass of hydrilla (Fig.
Hydrilla biomass (a, c) and zebra mussel density on hydrilla (b, d) at different sites (a, b) and mean by treatment at muddy vs. rocky sites (c, d). The lower and upper lines of the boxplots represent the 25th and 75th percentiles and the bold line represents the median. Data falling outside the percentile range are plotted as outliers. Different letters denote significant differences (Tukey HSD pairwise comparison, p-values < 0.05). Differences between muddy vs. rocky sites (c, d) were not significant.
Zebra mussel densities and size differed between those attached to rocks and those attached to hydrilla. At rocky sites in June, zebra mussel densities were significantly higher attached to hydrilla compared to those attached to rocks at all sites (Fig.
Zebra mussel density on rocks (white boxplots) and on hydrilla (grey boxplots) and mean by treatment for Hydrilla vs Rock as substrates in June 2022 for sites where rocky substrates were dominant. The lower and upper lines of the boxplots represent the 25th and 75th percentiles and the bold line represents the median. Data falling outside the percentile range are plotted as outliers. Asterisks above the diagram denote statistically significant differences between attachment surfaces (rocks vs. hydrilla) at each site (paired student’s t-test, * p-value < 0.05, ** p-value < 0.01, *** p-value < 0.001).
At all rocky sites, zebra mussels on rocks were significantly larger (6.1 ± 3.1 mm, mean ± SD) compared to mussels on hydrilla (3.4 ± 0.8 mm) (Fig.
Length of zebra mussels in June attached to rocks (white boxplots) and hydrilla (light grey) and attached to hydrilla in September (dark grey). The lower and upper lines of the boxplots represent the 25th and 75th percentiles and the bold line represents the median. Data falling outside the percentile range are plotted as outliers. Asterisks above the boxplots denote statistically significant differences between mussels on rocks vs. hydrilla in June (Mann Whitney U tests with Bonferroni corrections, * p-value < 0.05, ** p-value < 0.01. There were no significant differences between June and September except for Crane A and BR1 Cove G (denoted by a plus, pairwise t-test, p-value < 0.05).
Zebra mussel presence did not have a significant effect on hydrilla biomass (Low Density, W = 20, p = 0.14; High Density, W = 9, p = 0.53; Fig.
Mussel biomass was higher when hydrilla was present (F2,11 = 4.3, p = 0.04; Fig.
Ash Free Dry Mass (AFDM) of zebra mussels in low density (LD) and high density (HD) of hydrilla treatments (grey boxplots) and control (white boxplot). The lower (25th) and upper (75th) lines of the boxplots represent respective percentiles and the thick bold line indicates median values. Different letters denote significant differences (Tukey HSD pairwise comparison, p-values < 0.05).
This study examines variation in the distribution of both zebra mussels and hydrilla in relation to substrate, each other, and season (i.e., summer vs. fall). It also quantifies densities of zebra mussels attached to hydrilla. Our results support findings by previous studies that macrophytes support transient populations of seasonally high numbers of small zebra mussels (
The lower numbers of zebra mussels found on hydrilla in September were most likely the result of loss of the spring cohort and lower settlement occurring during summer months (
Other factors may have influenced population declines such as predation by fish (e.g., catfish) or dissolved oxygen depletion. Respiration by bacteria and macrophytes at night can cause oxygen levels in the water column to plummet due to the absence of photosynthetic activities (
Although zebra mussels were found attached to hydrilla, densities of mussels in this study were moderate compared to attachment on other macrophytes found in other studies. For example, we found up to 46,729 ind. m-2 located on hydrilla dry mass, whereas densities up to 750,000 ind. m-2 were found on clasping-leaf pondweed (Potamogeton perifoliatus) and Eurasian watermilfoil (Myriophyllum spicatum) in Lake Balaton, Hungary (
Higher densities of hydrilla in muddy compared to rocky substrates were observed. However, the presence of hydrilla observed in rocky areas further corroborates past evidence that hydrilla is often found in rocky areas despite limited sediment nutrients (
In Lake Lewisville, Texas, hydrilla above-ground biomass peaked in early August (
Zebra mussels have been shown to improve the physiological condition of Eurasian Watermilfoil and Eelgrass (
Zebra mussels may benefit physiologically from the presence of hydrilla as shown in our laboratory experiment where mussels increased biomass when hydrilla was present, particularly at low density. One possible explanation is that macrophytes can provide attachment surface for biofilms and other bacteria which can be utilized as a food source for zebra mussels (
Reservoirs, albeit man-made ecosystems, provide habitat for plants and animals as well as sources of drinking water and hydroelectric power, but abundant recreational access points may result in species introductions, making it crucial to fully understand and effectively manage invasive species to preserve biodiversity and prevent economic impacts. This study shows that managing hydrilla is important as it may directly facilitate zebra mussel dispersal through transport of mussel-fouled hydrilla fragments downstream by currents or to new water bodies by human vectors, as plants with attached zebra mussels are transported on boat propellers or trailers (
We would like to thank students and technicians in the Texas State Stream Ecology Lab for their assistance with field and lab work. We also extend our gratitude to the anonymous reviewers for their valuable feedback, which helped improve the clarity and quality of this manuscript. This work was funded in part by the US Army Corps of Engineers’ Aquatic Nuisance Species Program’s focus on Next Generation Ecological Modeling. The findings and conclusions in this article are those of the authors and do not necessarily represent the views of the United States Army Corps of Engineers.
Supplementary data
Data type: docx
Explanation note: table S1: Physical characteristics of the eight survey locations arranged in order (Top to bottom) from far to near distance from the dam, stars denote sites with sediments dominated by rock; fig. S1: Relationship between total densities of zebra mussels (benthic and hydrilla attached) and biomass of hydrilla for each quadrat in June 2022.