Effects of mussels on nutrient cycling and bioseston in two contrasting tropical freshwater habitats

Freshwater mussels (Unionida) can strongly affect nutrient cycling in temperate ecosystems but data from the tropics is lacking. We quantified the effects of mussel filtration, excretion and biodeposition on nutrient and photosynthetic pigment concentrations in a tropical eutrophic lake and mesotrophic river, featuring one non-native and two native species, respectively. Changes in nutrient and pigment concentrations were measured over a 3 h period to assess effects on (1) the water column in field enclosures, and (2) water column and benthos combined in controlled laboratory experiments. In field enclosures in both systems, mussel density and biomass were significantly correlated with the magnitude of reduction in sestonic pigment concentrations. In laboratory experiments, presence of mussels led to reduced PO4 and increased TAN concentrations in both systems, lower combined sestonic and deposited pigment concentrations in the river but increases in the same in the lake. We conclude that excretion by mussels probably accelerated bioseston growth in both systems due to N-fertilisation, an effect that may be particularly common in tropical freshwaters, which are frequently N-limited. However, whilst river mussels reduced bioseston concentrations through rapid filtration, higher rates of N-excretion and/or deposition of undigested bioseston by lake mussels apparently resulted in a net increase of pigment concentrations.


Introduction
Freshwater mussels (Bivalvia: Unionida, hereafter mussels) are sedentary, benthic molluscs, that feed primarily through filter-feeding and can occur in high densities in various freshwater habitats globally (Graf & Cummings, 2015). Their filtering activity, which has been assessed to attain between about 0.1-3 l h -1 mussel -1 (Kryger & Riisgård, 1988;Chowdhury et al., 2016;Cyr et al., 2017 and references therein), affects ecosystems at various trophic levels and compartments . Firstly, mussels affect seston abundance and composition in the water column by removing organic and inorganic particles, including phytoplankton, zooplankton and bacteria (Christian et al., 2004;Howard & Cuffey, 2006;Spooner & Vaughn, 2006 and references therein). High mussel densities can result in increased water clarity, affecting growth and community composition of plankton and macrophytes (Welker & Walz, 1998;Chowdhury et al., 2016;Vaughn, 2018). At the same time, mussels excrete nutrients as solutes (predominantly NH 3 and PO 4 ) to the water column, thereby providing nutrients for primary producers (Nalepa et al., 1991). Mussel excretion generally increases the N:P ratio, which can alleviate strict N-limitation in streams and can lead to a subsequent change in algal communities (Atkinson et al., 2013).
Mussels also affect nutrient concentrations and composition of the sediment by biodeposition. Whilst ingested edible material and part of the nutrients are assimilated and converted into soft tissue and shell, much of the material taken up from the water column is deposited to the sediment as non-assimilated pseudofaeces and egested faeces (Howard & Cuffey, 2006;Spooner & Vaughn, 2006;Saraiva et al., 2011). Biodeposition of faeces and pseudofaeces results in an increased availability of nutrients (particularly N and P) and other resources for benthic organisms, which in turn enhances benthic production and biodiversity (Howard & Cuffey, 2006). In addition, at least some species and populations of freshwater mussels have been shown to not exclusively feed on seston but exhibit a much more diverse diet including benthic algae and allochthonous organic material of terrestrial origin (James, 1987;Nichols & Garling, 2000;Raikow & Hamilton, 2001;Collier et al., 2017;Zhang et al., 2017).
Much of the research on the effect of mussel filtration, excretion and deposition on suspended algae (measured as pigments) and nutrient concentrations has been conducted in controlled ex situ experiments. Clearance rates have been estimated for a number of species by measuring the decline in chlorophyll a (chl a) concentrations after mussel filtration as a proxy for rates of phytoplankton removal from the water column (e.g., Kryger & Riisgård, 1988;McIvor, 2004;Chowdhury et al., 2016). Clearance rates, excretion rates, deposition rates and quality of excretion and deposition products (e.g., N:P ratios) vary with species, size and reproductive stage of mussels, food, temperature and other environmental conditions (see Vaughn et al., 2008 for a comprehensive list of references). For example, mussel populations from fast flowing rivers have been found to be more efficient in clearing phytoplankton from the water column than those from slow flowing rivers or lentic systems (Byllaardt & Ackerman, 2014).
Whilst controlled laboratory experiments are useful in determining mussel clearance, excretion and deposition rates, field-or semi-natural experiments are needed to understand the in situ effects of freshwater mussels on their environment. Vaughn et al. (2004) used controlled mesocosm experiments to assess how mussel density and species composition affects TAN (= total ammonia nitrogen = NH 3 ? NH 4 ) and chl a concentrations in the water column. TAN concentrations increased with mussel density as expected, but so did chl a due to a TAN-fertilisation effect. In contrast, Soto & Mena (1999) showed that Diplodon chilensis (Gray, 1852) significantly reduced water columnconcentrations of chl a (* tenfold in 18 days), TAN (* fourfold), SRP (soluble reactive phosphorus; * fourfold) and TP (total phosphorus; * tenfold) in outdoor fish tanks compared to tanks without mussels. Depending on the ecosystem and characteristics of dominant species present, the effects of freshwater mussels on nutrient cycling and phytoplankton production can thus be highly variable, from reducing nutrient loadings to stimulating primary production where nutrients are limiting by converting suspended material to dissolved nutrients (Vaughn et al., 2007). For example, the effect of mussels on nutrient cycling of lotic systems has been shown to be much stronger under low discharge compared to high discharge conditions (Vaughn, 2018).
Data on freshwater mussel filtration, excretion and biodeposition published to date are almost exclusively restricted to temperate systems. Nutrient cycling in tropical freshwaters differs considerably from that in temperate ones, owing to consistently high temperatures and strongly seasonal patterns in river discharge, water depth, conductivity, and concentrations of oxygen and nutrients (Winemiller & Jepsen, 2005). In comparison with temperate rivers, tropical freshwaters are characterised by much higher sediment load and are more frequently N-limited, resulting in greatly increased primary production when pristine tropical lands are disturbed (Downing et al., 1999). In small tropical streams, benthic organisms and food resources are frequently diminished by flash floods that scour the substrate (Winemiller & Jepsen, 2005). Filtration rates and effects of mussel feeding on tropical ecosystems may therefore be expected to differ from those in temperate systems. At the same time, understanding the functional role and importance of biota in tropical ecosystems is particularly important when considering the high rates and extent of species loss in these systems (Dudgeon et al., 2006;Cumberlidge et al., 2009). Despite these facts, we are aware of only a single study on the effects of tropical freshwater mussels, which determined ex situ chl a clearance rates of two Unionidae species collected from an urban lake in Bangladesh (Chowdhury et al., 2016).
The present study aims to assess the effects of mussel filtration, excretion and deposition on nutrient and bioseston concentrations in two tropical ecosystems; a natural river and an artificial lake in Peninsular Malaysia.

Study sites
Semenyih Lake is a small, shallow, eutrophic artificial lake (max depth * 2 m, total area * 0.05 km 2 ) located close to a residential area and used for fishing and other recreational activities (Table 1). Muar River is a mesotrophic river (* 7-11 m wide) that flows through a mosaic of oil palm plantation and secondary forest, several km from the nearest human settlement ( Table 1). The only mussel present in the lake is the non-native Sinanodonta woodiana Lea, 1834, which is native to China and Russia, and has been widely introduced in ponds and rivers across Malaysia (Watters, 1997;Zieritz et al., 2016Zieritz et al., , 2018b. The river features two native mussel species, Contradens contradens (Lea, 1838) and Monodontina vondembuschiana (Lea, 1840) .

Experimental setup
We conducted two experiments that were designed to assess the effects of freshwater mussels on their environment in semi-natural conditions: (1) The relationship between mussel density/biomass and changes in nutrient and photosynthetic pigment concentrations in the water column was assessed in field experiments at the two study sites at Semenyih Lake and Muar River (Table 1). (2) Effects of mussel filtration, excretion and deposition on the water column and benthos combined were assessed in laboratory experiments using mussels and water collected from the two study sites.

Field experiments
Field experiments were carried out over four days in the dry season in July 2015 (Muar River) and March 2016 (Semenyih Lake), respectively. Each day, at each of the two study sites, three black, bottomless, open plastic cylinders (diameter = 43 cm, height = 51 cm) were placed at randomly selected spots across a 100-m stretch of river/lake where mussels were present, at a maximum of 3-m distance from the shore where the water level was \ 40 cm and the tops of the cylinders were therefore above water column. Mussels were not disturbed, and mussel density was not manipulated to provide as natural conditions as possible. Water depth was measured in each cylinder. Water samples (1 l) were taken from each of the cylinders at 50% of water depth at the beginning of the experiment and again after 3 h. This corresponds to water before mussel filter-feeding activity and after filter-feeding activity takes place, respectively. After the filtration period, all the mussels from the cylinders were collected, counted and their shell length (SL) measured to ± 0.1 cm accuracy using a sliding calliper. In addition, two replicate water samples were taken per day from outside the cylinders for description of general characteristics of the study sites.

Laboratory experiments
The laboratory filtration experiments were carried out at 25°C air temperature at the University of Nottingham Malaysia Campus during the dry season in August (Semenyih Lake) and September 2016 (Muar River). In the morning of the day of the experiment, a total of 20 (10 C. contradens and 10 M. vondembuschiana) and 15 (for S. woodiana) randomly selected mussels, and approximately 80 l of river/lake water were taken from each of the two sites and transported to the laboratory in darkened containers. In the laboratory, 23 and 18 (for Muar and Semenyih experiment, respectively) transparent, cylindrical, 3.6 l plastic tanks were set up with 3.5 l of lake water and an air stone connected to an air pump in each; no substrate was added to the tanks. A 0.5 l water sample was taken from each tank for analysis. Immediately afterwards, one mussel that was previously scrubbed to remove epibionts was put inside each tank with the exception of three randomly selected tanks, which served as controls. After 3 h, mussels were removed from the tanks and water was collected after homogenisation. SL of each mussel was measured as above. Mussel soft tissue was dissected from each specimen, dried at 80°C for 48 h and weighed to ± 0.01 g accuracy.

Water chemistry analysis
Water samples were refrigerated and processed as soon as possible on the day of the experiment. A known volume of each water sample (0.2-0.4 l, depending on turbidity) was filtered through Whatman Ò GF/C-filters with a vacuum pump. Filters were frozen, and filtered and unfiltered water samples were kept cool for subsequent determination of following parameters in the laboratory: Concentrations of chlorophyll a (chl a), carotenoids, total phosphorus (TP), soluble reactive phosphorus (SRP), total ammoniacal nitrogen (TAN), and nitrate (NO 3 -N) were assessed using standard spectrophotometric/colorimetric methods (Lorenzen, 1967;Mackereth et al., 1989). Particulate phosphorus (PP) was subsequently estimated as PP = TP -SRP.

Data analysis
Comparison of study sites Differences in baseline conditions at the two study sites were assessed using Mann-Whitney U tests on the six water quality parameters measured from eight replicate water samples taken across the four study days, as well as mussel density, biomass and SL in the 12 replicate field cylinders. Mussel density and biomass were calculated both per area and water volume using the surface area covered by the cylinder (i.e., 1452 cm 2 ) and water depth in the cylinder. Biomass per cylinder was estimated by (1) deriving the power equation of SL versus tissue dry weight for each habitat from the specimens used in the laboratory experiment; (2) applying these equations to estimate tissue dry weight for each specimen found in field cylinders based on SL measurements; and (3) chl a (lg l -1 ) 1 6 ± 10 (5-32) (N = 8) 1.6 ± 0.8 (0.6-2.7) (N = 8)** carotenoids (lSPU l -1 ) 20 ± 11 (5-34) (N = 8) 6 ± 6 (1-13) (N = 8)* Data are presented as averages ± standard deviation (range). Asterisks indicate significant differences between the study sites determined by Mann-Whitney U test (*P \ 0.05, **P \ 0.01, ***P \ 0.001; n.s. not significant) chl a chlorophyll a, DW mussel soft tissue dry weight, NO 3 -N nitrate, SL shell length, SRP soluble reactive phosphorus (PO 4 -P), TAN total ammonia nitrogen (NH 4 -N ? NH 3 -N), TP total phosphorus summing together the estimated tissue dry weights of all specimens in a cylinder. The change in concentration of the six parameters (Dchla, Dcarotenoids etc.) after 3 h incubation was determined for each replicate cylinder (field experiments) or tank (laboratory experiments) by subtracting concentrations at the beginning of the experiment from respective concentrations at the end of the experiment.

Field experiments
We assessed the effect of mussel density and biomass in field experiments using General Linear Models (GLMs) on each of the six measured parameters, fitting mussel density or biomass (tissue dry weight) as a covariate, day as a factor with three levels, and the interaction factor. Non-significant factors/covariates were sequentially removed from the model (backwards selection).

Laboratory experiments
We tested for differences in each of the six parameters between mussel (n = 20 and 30 for lake and river experiments, respectively) and control treatments (n = 3) in laboratory experiments using Welch's unequal variances t tests, which are more reliable than Students t tests when sample sizes are unequal (Ruxton, 2006). For the river-laboratory experiment, these were preceded by Students t tests on the mussel dataset only to test for significant differences between the two river mussel species (n = 15 for both species). We then assessed the effect of mussel size and species (for river-laboratory experiments only) on each of the six parameters using GLMs, fitting SL or biomass (tissue dry weight) as a covariate, species as a factor with two levels (for river-laboratory experiment only), and the interaction factor. Non-significant factors/covariates were sequentially removed from the model (backwards selection).
We calculated clearance rates of chl a (CR chla ) for each mussel using the formula ln (C 0 /C t ) 9 V/T, and excretion rates for TAN (R TAN ) using the formula (C t -C 0 ) 9 V/T, where C t and C 0 are final and initial nutrient or pigment concentrations, respectively, V is the volume of water in the tank (3 l), and T is incubation time (3 h) (following Riisgård (2001a, b) and Cyr et al. (2017)), and deducting respective rates obtained in control tanks. All statistical analyses were conducted in R.

Differences between study sites
The two study sites differed significantly in water quality parameters and other characteristics (Table 1). Lake water concentrations of chl a and TAN were significantly higher than in river water, whilst SRP and NO 3 -N concentrations were significantly higher in the river than the lake ( Table 1). Ratios of available N/P were * 5 in the lake and \ 1 in the river.
Mussel densities in the river were significantly higher than the lake (Table 1). However, due to significantly different mussel sizes at the two habitats, mussel biomass (soft tissue dry weight) per area and volume of water did not differ significantly between the two sites (Table 1). SL within each habitat was similar in laboratory and field experiments (t tests: River: t = 0.359, df = 117, P = 0.720; Lake: t = 1.828, df = 63, P = 0.072). The relationship between SL and dry weight of mussel tissue (DW) was DW [g] = 0.043 9 SL [cm] 2.158 (R 2 = 0.67) for lake mussels and DW [g] = 0.050 9 SL [cm] 1.977 for river mussels (R 2 = 0.64); this relationship was not different between the two river species (ANCOVA: dry weight: F = 82.195, P \ 0.0001; Species: F = 2.080, P = 0.161; dry weight 9 Species: F = 0.068, P = 0.797).

Field experiments (effects on the water column)
In field experiments, mussel density had a clear effect on sestonic photosynthetic pigment concentrations, with Dchl a and Dcarotenoid being significantly negatively correlated with mussel density in the river and lake, respectively (Fig. 1a, b). Although changes in sestonic nutrient concentrations in the water column were not significantly correlated with mussel density, R 2 values (C 0.15) indicated associations between mussel density and DTAN at both sites, with the relationship being positive at the lake, but negative at the river (Fig. 1e). Similar but slightly weaker correlations were found with mussel biomass (Lake: Dcarotenoid: P = 0.028, R 2 = 0.40; River: Dchl a: P = 0.040, R 2 = 0.36; all other correlations not

Laboratory experiments (combined effects on the water column and benthos)
In laboratory experiments conducted on the lake system components, combined sestonic and deposited Dchl a, Dcarotenoids, DPP, DTAN and DNO 3 -N was significantly larger in mussel treatments compared to controls (Table 2, Fig. 2). Mussel presence also appeared to negatively affect DSRP though this effect was not statistically significant (Table 2, Fig. 2d). In addition, within mussel treatments, Dchl a, Dcarotenoids, DPP and DTAN significantly increased with SL, whilst no significant effect was found in this respect for DSRP and DNO 3 -N (Fig. 3).
Excretion rates of TAN per mussel differed considerably between systems. One lake mussel on average excreted 17 9 as much TAN as a river mussel (lake: 51 ± 21 lg N mussel -1 h -1 ; river: 3 ± 7 lg N mussel -1 h -1 ), and about 5 9 as much when accounting for differences in mussel size between the two ecosystems (i.e., weight-corrected excretion rates calculated based on average DW of 5.4 g and 1.5 g DW for lake and rivers mussels, respectively: lake: 9.4 lg N g DW -1 h -1 ; river: 2.0 lg N g DW -1 h -1 ). Considering that average densities were 0.47 and 0.28 g DW l -1 in the lake and b Fig. 1 Relationship between mussel density and change in pigment and nutrient concentrations after 3-h enclosure experiments in Semenyih Lake and Muar River, Malaysia. Bold values are significant at 95% level   (Table 1), mussels excrete about 4.4 lg N l -1 h -1 in the lake and 0.56 lg N l -1 h -1 in the river system. CR chla averaged 0.8 l mussel -1 h -1 for C. contradens and 1.0 l mussel -1 h -1 for M. vondembuschiana, reaching a maximum of 3.2 l mussel -1 h -1 in the river. CR chla for S. woodiana was negative, averaging -0.9 l mussel -1 h -1 .

Increasing TAN concentrations by excretion
Laboratory experiments showed that excretion by both lake and river mussels increased TAN concentrations. Compared to previously reported values, which range from excretion of about 4-80 lg N mussel -1 h -1 or 5-80 lg N g DW -1 h -1 (Nalepa et al., 1991;Baker & Hornbach, 2001;Vaughn et al., 2004;Cyr et al., 2017), TAN excretion rates reported in the present study were moderate to low for Semenyih Lake and very low for Muar River. Lake mussels excreted 17 9 more TAN on average compared to river mussels, and almost 5 9 more when accounting for size differences. TAN excretion rates of lake mussels are probably underestimations, as suggested by the significant increase in NO 3 -N concentrations in tanks with mussels compared to controls. It is likely that this increase arises from nitrification of the NH 4 excreted by S. woodiana, which may have been amplified by aeration during the experiment and explains the lack of correlation between NO 3 -N and mussel size (Lewis, 2001). S. woodiana in the lake thus exhibits much higher TAN excretion rates than the two native study species in the river, even after accounting for size differences between the species (Fig. 3e). We cannot distinguish whether the observed difference in TAN excretion rates was caused by differences in seston quantity and quality between sites, by a species-effect, or a combination of the two, but mass-specific TANexcretion rates have been observed to differ by up to [ 3 times even between sympatric freshwater mussel species .

N-fertilisation accelerates bioseston growth
Both study systems were N-limited, with soluble N/P ratios considerably lower than the Redfield ratio. Data from our laboratory experiments indicate that TAN excretion by mussels alleviated this N-limitation and stimulated bioseston growth in both systems. That is, a pronounced decrease in SRP concentrations was observed in mussel tanks compared to control tanks for both the lake and the river system (Fig. 2d). This decrease was unexpected as freshwater mussels are known to excrete SRP (Nalepa et al., 1991;Vaughn et al., 2004) but may be explained by increased growth of bioseston and consequently, increased SRP uptake, in mussel tanks due to increased TAN availability. A similar TAN-fertilisation effect by mussels leading to increased bioseston growth (detected as increased chl a concentrations) has been observed previously in a mesocosm study by Vaughn et al. (2004).
Despite N-fertilisation, for the river-derived laboratory experiments, chl a concentrations decreased significantly in the presence of mussels, suggesting that phytoplankton was rapidly taken up by mussels and passed into their digestive system. The lack of a significant reduction in carotenoid concentrations in mussel treatments is consistent with the idea that carotenoids are generally more likely to pass through invertebrates undigested than chl a and are therefore more likely to be deposited (McLeroy-Etheridge & McManus, 1999).
In contrast, in mussel tanks from the lake system, pigment concentrations increased about 10-20-fold, an effect that was positively correlated with mussel size. We suggest that this vast increase in phytoplankton concentration could be due to a particularly strong N-fertilisation effect by mussel excretion, which averaged about 50 lg TAN l -1 h -1 in the lake laboratory experiment, and which was also positively correlated with mussel size. In the experimental density of 1 mussel 3 l -1 , this might lead to phytoplankton growth rates that are faster than clearance rates of mussels. An alternative or additional b Fig. 2 Boxplots of change in pigment and nutrient concentrations after 3 h in 3 l of lake/river water with single mussels (M; n = 20 for lake, n = 30 for river) and without mussels (C; n = 3). Box limits represent third and first quartile; centre horizontal line represents median; vertical lines represent maximum and minimum; dots represent outliers. Asterisks indicate significant differences between control and mussel tanks determined by Welch Two-Sample t tests (*P \ 0.05, **P \ 0.01, ***P \ 0.001; n.s. not significant). Note that boxplots and y-axis on the left refer to the lake habitat, whereas those on the right refer to the river 123 explanation for the observed increase in photosynthetic pigments in lake-mussel tanks could be bioseston material that was present in the digestive system of lake mussels at the time of collection, and was released as (partly) undigested bioseston during the experimental period. Selective feeding and ejection of different phytoplankton taxa has been observed in Sinanodonta calipygos (Kobelt, 1879) from Lake Biwa, which was shown to preferentially feed on small algae \ 10 lm in size (e.g., Scenedesmus, Selenastrum and Synechococcus), whilst large species were rejected (e.g., Aphanocapsa, Pediastrum spp.) (Miura & Yamashiro, 1990). Similar processes might be important in the closely related S. woodiana in the eutrophic Lake Semenyih. As our experimental design did not allow for discrimination between bioseston in the water column and deposited bioseston, we are unfortunately unable to accurately quantify the relative contribution of each of these components.

Clearance rates
In the field experiments, seston Dchl a and Dcarotenoids were significantly negatively correlated with mussel density and biomass at both sites. This indicates that freshwater mussels are clearing phytoplankton from the water column, with more mussels clearing a larger proportion of the water volume, as has been observed in a number of previous experiments (e.g., Soto & Mena, 1999;Chowdhury et al., 2016;Douda & Č adková, 2018). Clearance rates (CRs) calculated from laboratory experiments for both river species were relatively high, averaging 0.8-1 l mussel -1 h -1 and reaching over 3 l mussel -1 h -1 in M. vondembuschiana. These values exceed most previously published freshwater mussel laboratory filtration rates that were measured in tanks under static conditions, without sediment and using cultured algae or water taken from a different source than the study population. For example, previously reported CRs from other freshwater mussel populations of similar average sizes (* 5 cm SL) averaged 0.1-0.7 l mussel -1 h -1 in different lentic and lotic species of North America (Byllaardt & Ackerman 2014); 0.3 l mussel -1 h -1 in Parreysia caerulea Lea, 1831/Lamellidens marginalis (Lamarck, 1819) from a eutrophic lake in Bangladesh (Chowdhury et al., 2016), and 0.5-0.7 l mussel -1 h -1 in Echyridella menziesi (Gray, 1843) from six oligo-and mesotrophic lakes in New Zealand (Cyr et al., 2017). Our values were, however, comparable to those measured by Kryger & Riisgård (1988) in a more natural circumstances using sediment-containing aquaria, with a single-species algal culture. The reason for the comparatively high filtration rates observed in the present study might be that we used water from the mussels' natural habitat and conducted experiments at a slightly higher temperature (25°C) compared to most previous studies (19-21°C). Nevertheless, we argue that observed filtration rates were still likely underestimations for a number of reasons: (1) mussels are likely to have been disturbed by the sampling process; (2) water temperatures in the lake and river site usually exceed 25°C (A. Zieritz, unpublished data); (3) our experimental design does not distinguish between chl a not ingested from the water column and chl a ingested and deposited to the water column. Finally, (4) Byllaardt & Ackerman (2014) showed that compared to static conditions, freshwater mussel CR (of algae measured with a fluorometer) under continuous supply of algae (algal flux) in a recirculating flowchamber increases by an average of 20 and up to over 40 times, especially for riverine species. High CRs in flow-through chambers were confirmed by Douda & Č adková (2018) on S. woodiana from Kyjovka River, Czech Republic, measuring an average filtration rate of about 2.5 l per mussel. Future studies on tropical river mussels should thus be conducted under more natural flow conditions using river mesocosms or flow chambers.
Due to the observed increase in chl a in lake laboratory experiments (see above for a discussion of potential reasons for this observation), clearance rates could not be calculated for lake mussels.

Conclusions
Our study indicates that the effects of mussels on nutrient cycling and bioseston growth can differ considerably between tropical freshwater ecosystems, depending on nutrient availability in the system, and excretion rates and filtration ability of the mussel b Fig. 3 Relationship between mussel shell length and change in pigment and nutrient concentrations per species after 3 h in 3 l of lake or river water. Bold values are significant at 95% level community. Excretion by mussels probably accelerated bioseston growth in both systems due to N-fertilisation, an effect that may be particularly common in tropical freshwaters, which are frequently N-limited. However, whilst river mussels reduced bioseston concentrations through rapid filtration, higher rates of N-excretion and/or deposition of undigested bioseston by lake mussels apparently resulted in a netincrease of pigment concentrations. Further experiments are required to assess whether these results can be generalized across tropical river and lake systems, and whether the observed differences are triggered by differences in the environmental conditions at the sites, a species-effect or a combination of the two. Answering the question whether non-native and native species play different functional roles in the same ecosystem will be particularly important considering the fast and ongoing spread of S. woodiana across Southeast Asia and beyond (Zieritz et al., , 2018a.