1. Introduction
Factors such as intensified industrialization, urban growth, and agricultural activities lead to increased inputs of nitrogen and phosphate nutrients in aquatic ecosystems, causing eutrophication (Chorus & Welker, 2021). This process has favored the occurrence of cyanobacterial harmful algal blooms (CyanoHABs) which have been boosted by climate changes (Huisman et al., 2018). Bloom-forming cyanobacteria can cause negative impacts in aquatic ecosystems by cyanotoxins production, oxygen depletion, shading of benthic habitats and trophic decoupling (Havens & Hudnel, 2008; Pearl & Huisman, 2008; Moustaka-Gouni & Sommer, 2020).
Cyanobacteria are photosynthetic prokaryotes that usually compose marine and freshwater phytoplankton, and contribute to the primary production in aquatic food webs as a significant food source for zooplankton consumers, especially in eutrophic water bodies (Lürling & Verschoor, 2003; Nogueira et al., 2004; Ger et al., 2016). However, dietary exposure to cyanobacteria can affect zooplankton due to low nutritional quality, such as the lack of polyunsaturated fatty acids (PUFAs) and sterols (Martin-Creuzeburg et al., 2008); difficulted digestibility (DeMott & Van Donk, 2013) or by the ability of some cyanobacterial genera on producing bioactive secondary metabolites (cyanotoxins) which can affect zooplankton motility, feeding behavior, survival and life cycle (Ferrão-Filho & Kozlowsky-Suzuki, 2011; Vilar, Ferrão-Filho & Azevedo, 2022).
Among cyanobacteria, Raphidiopsis raciborskii (formerly Cylindrospermopsis raciborskii) (Aguilera et al., 2018) stands out as a N-fixing and bloom-forming species which has been recently reported for higher latitudes as an invasive aquatic microbe. The ability of R. raciborskii on spreading over several regions is attributed to its high phenotypic plasticity regarding tolerance to light, salinity, ability to storage phosphorus, toxin production, and resistance to herbivory (Antunes et al. 2015; Burford et al., 2016; Vilar et al., 2021). Toxin-producing strains of this species can produce cytotoxic (cylindrospermopsins) (Hawkins et al., 1997; Wood & Stirling, 2003; Lei et al., 2014) and neurotoxic (saxitoxins, goniautoxins, and C toxins, hereafter STXs) alkaloids; being the latter toxins produced up to date only by south American strains (Lagos et al.,1999; Molica et al., 2005; Pomati et al., 2003; Ferrão-Filho et al., 2007; Piccini et al., 2011). R. raciborskii blooms are increasingly frequent in tropical lakes and reservoirs which makes this species a representative component of phytoplankton communities in these water bodies (Soares et al., 2013; Rzymski et al., 2014; Yang et al., 2017; Pagni et al., 2020).
Overall, the establishment of cyanoHABs usually represents a change from diverse and edible phytoplankton to a less diverse and inedible one (large colonies and filaments, and production of harmful metabolites) (Hansson et al., 2007; Moustaka-Gouni & Sommer, 2020). In turn, it can affect zooplankton composition, which shifts from a community dominated by large generalist filter feeders such as cladocerans (e.g., Daphnia), toward small-bodied selective organisms such as copepods, rotifers, and planktonic protozoans (Przytulska et al., 2015; Moustaka-Gouni & Sommer, 2020). However, although there is a consensus that cyanobacteria play a role in the trophic decoupling of freshwater food webs (Heathcote et al., 2016), mainly due to defensive chemical traits against phytoplankton consumers (Rangel et al., 2020; Vilar et al., 2021), little is known whether the increased cyanobacterial biomass affects the zooplankton due to a relatively low abundance of nutritious algae or if cyanobacterial toxicity overrides nutritional constraint regardless of high quality food availability at a non-limiting biomass concentration. Therefore, we evaluated whether the negative impacts of dietary exposure of Daphnia to toxic cyanobacteria are mostly ruled out by nutritional constraints due to nutritious algae availability. Therefore, we tested the hypothesis that the toxicity of STXs-producing R. raciborskii impairs Daphnia feeding behavior and fitness regardless of the high availability of nutritious algae.
2. Methods
2.1. Phytoplankton and zooplankton culture conditions
The saxitoxin-producing strain Raphidiopsis raciborskii CYLCAM-01 (Mesquita et al., 2019) was kindly provided by Dr. Marcelo M. Marinho (Laboratory of Ecology and Physiology of Phytoplankton – Rio de Janeiro State University). The cyanobacterium was isolated from the Camorim reservoir; a eutrophic water body located at the Pedra Branca State Park (PEPB) – Rio de Janeiro, Brazil. The strain was cultivated in a WC medium (Guillard, 1975) with initial pH adjusted to 8.0 and maintained at a temperature of 23±1 °C, the light intensity of 50 μmol m-2s-1 and a 12:12 h light/dark photoperiod.
The tested zooplankton consisted of neotropical Daphnia species isolated from different water bodies: D. laevis (Lagoa do Ibirité-Minas Gerais), D. laevis (Lagoa Jacaré-Minas Gerais) and D. gessneri (Reservatório do Apertado-Bahia). Animals were kept in 1000 mL beakers filled with artificial RT medium (Tollrian, 1993) enriched with 30% filtered (filter 1-2 µm pore size, Sartorius) natural water from decantation tanks of the CEDAE treatment plant (State Company for Water and Sewage, Camorim Reservoir) and supplied with commercial 0.05% humic acid extract (~1.125 mg C L-1, Microbe Lift® Amazon Black & Soft Water Conditioner, USA) at an initial pH 7.6. The animals were incubated at 23.5 ± 1 °C under low light intensity (20 μmol m-2s-1) and a 12:12h light/dark photoperiod. To feed the animals, the green microalgae Monoraphidium capricornutum and Ankistrodesmus stiptatus were cultivated in oligo medium (Afnor, 1980) with initial pH adjusted to 7.0. For cultures maintenance, the animals are usually fed at each two days with a total carbon biomass of 0.4 mg L-1 of both algal species (1:1).
2.2. Phytoplankton biomass estimative and saxitoxin analysis
Cyanobacterial filaments and microalgal cells were counted in a Fuchs-Rosenthal chamber under an optical microscope to determine the cell/filament density (unit mL-1). Subsequently, the population density was converted to biovolume (mm3 L-1) from the average cell/filament volume (µm3) (Hillebrand et al., 2002; Sun & Liu, 2003) obtained by cells and filaments dimensions (length and width) measurements under an optical microscope with the aid of CellB software (Soft Imaging System GmbH, Münster, Germany). Finally, the carbon biomass (mgC L-1) was estimated from cell volume measurements as described by Rocha & Duncan (1985).
STX was analyzed by ELISA (Enzyme-Linked Immunoabsorbent Assay) using a Beacon saxitoxin kit (Beacon Analytical Systems, Inc.) according to the manufacturer’s recommendation.
2.3. Life table experiments
Life-table experiments were run to evaluate whether the amount of saxitoxin-producing cyanobacteria in the diet could affect zooplankton fitness, biomass acquisition, and survivorship, despite the availability of nutritious green microalgae (hereafter nutritious food) as follows:
2.4. Dietary exposure to STX-producing cyanobacteria under constant nutritious food
In the first experiment, neonates (<24h) of the species D. laevis (RD clone), D. laevis (IB clone) e D. gessneri were exposed to increased concentrations of the toxic R. raciborskii CYLCAM-01 (0, 0.5, 1.0 and 1.5 mg C L-1) and a constant biomass concentration of nutritious food (0.4 mg C L-1 of M. capricornutum) during 15 days. The experiment was run with one individual per test tube (n=15) filled with 30 mL of the phytoplankton cell suspensions at the different biomass concentrations (experimental diets) and incubated over 15 days. Furthermore, we assumed to keep a constant nutritious food concentration in an attempt to separate the effects of nutritional inadequacy and toxicity (Jónasdóttir et al. 1998).
2.5. Dietary exposure to STX-producing cyanobacteria under variable nutritious food
In the second experiment, neonates (<24h) from each Daphnia species/clone were exposed to increased cyanobacterial proportions (25-90%) relative to nutritious food (75-10%). For the nutritious food (green algae), we kept a higher proportion of A. stiptatus (90%) relative to M. capricornutum since this microalgae seems to be more suitable for zooplankton nutrition regarding its lipid composition (as well as phosphate content) which consists mainly of unsaturated fatty acids (C18 chains such as alpha-linolenic acid, oleic acid, and linoleic acid) (Miranda et al., 2016). Thus, the Daphnia dietary exposure to toxic R. raciborskii CYLCAM-01 was set into three experimental diets: i) 25% CYLCAM-01 + 75% nutritious food; ii) 50% CYLCAM-01 + 50% nutritious food; and iii) 90% CYLCAM-01 + 10% nutritious food. The control group consisted of 100% nutritious food. All proportions were relative to a total carbon biomass concentration of 1.0 mg C L-1. The experiment was run in test tubes (n=15) filled with 30 mL of phytoplankton cell suspensions and incubated over 15 days. For all experimental diets, the medium was daily renewed over the experiment duration.
For the life table parameters, neonates individually exposed to the different diets were assessed from newborn to maturity during 15 days or up to at least the third brood, to assess the (
i) age at first reproduction (primiparous) and the number of neonates produced per female in each brood. The survival (
lx) and mean fecundity (
mx) were used to estimate the mean intrinsic rate of population increase (
r) by bootstrap analysis (Software Rm 2.0, Taberner
et al., 1993) using the following equation:
where
lx is age-specific survival,
mx is age-specific fecundity and
x is the age in days.
At the beginning and the end of the experiments, the organisms of each treatment were measured for the body length from the top of the head to the base of the caudal spine under an Olympus magnifying glass (SZ61) with 20× magnification to estimate the somatic growth rate by the following formula:
where
Li is the initial length,
Lf is the final length and
t is the incubation time in days.
2.6. Grazing assays
Grazing assays was performed to evaluate whether an increased proportion of saxitoxin-producing cyanobacteria in the diet, at a constant non-limiting biomass concentration, affect or not zooplankton feeding behavior.
Adult individuals (~6 days) of all Daphnia clones were dietary exposed to phytoplankton cell suspensions at increased proportions of R. raciborskii CYLCAM-01 (25, 50, and 90%) relative to the nutritious food (=green microalgae M. capricornutum), at a total carbon biomass concentration of 1.0 mg C L-1. The animals were starved (~12h) before the beginning of the experiment and then conditioned in 10 mL test tubes (n=4; 2 individuals/test tube) filled with the different experimental diets. A cellular suspension composed of 100% nutritious food was assumed as the negative control for cyanobacterial impact by dietary exposure. Alongside, all experimental diets were also established without animals as the negative control for cell loss by Daphnia grazing.
The
Daphnia clones were incubated in the dark for 3 hours and subsequently, the clearance rate (CR, in mL ind
-1 h
-1) was estimated from chlorophyll fluorescence (ChlF0) emitted from phytoplankton cells in the experimental diets and detected by a PHYTO-PAM fluorimeter (Walz, Germany) (Lürling & Verschoor, 2003). The clearance rate was estimated using the following equation:
with
Chlacontrol and
Chlatreatment as the difference in the algal concentration (as chlorophyll-
a specific for Cyanobacteria and green algae; µg L
-1) at the initial and final time in the control (without animals) and treatments (with animals), respectively. Δ
t is the incubation time (h),
V is the volume of culture medium (mL) and
N is the number of animals.
The calculation of clearance rate was based on the fluorescence signals given: green (for green algae), blue (for cyanobacteria) and total chlorophyll (for the sum of both signals).
2.7. Data analysis
CL50 (Lethal Concentration to 50% of individuals) and CE50 (Effective Concentration where 50% of the population exhibits a response) were estimated in the online platform MOSAIC (Charles et al., 2018).
Population parameters (reproduction parameters and intrinsic rate of population increase (r)), as well as feeding behavior were examined for normality and variance homoscedasticity. Once the parametric premises were assumed, the data were compared using one-way ANOVA, and differences between treatments were analyzed by Tukey HSD test (p<0.05) using GraphPad Prism 5.0 software. To verify differences regarding the population parameters between Daphnia species/clones and treatments we used the two-way ANOVA, with species/clone and treatment as categorical factors. The intrinsic rate of population increase (r) was analyzed by Student's t-test (p<0.05), using Primer 1.0 program.
4. Discussion
The results corroborated our main hypothesis that toxicity of STXs-producing R. raciborskii impairs Daphnia feeding behavior and fitness regardless of the high availability of nutritious algae. This effect was, however, species/clone dependent. In both life table assays, there was a reduction in the somatic growth and population increase rates of all species exposed to cyanobacterial concentrations relative to the control. Probably, the energetic cost invested in survival may have decreased the fitness of these animals, causing reduced fecundity and growth. The quality of food is determined by a variety of substances in natural environments. Cyanobacteria are generally considered a poor food due to lack of specific micronutrients, such as polyunsaturated fatty acids and sterols (Martin-Creuzburg et al., 2008; Von Elert et al., 2003). If diets are deficient in essential substances necessary for growth, such as sterols and fatty acids, growth rates may be affected (Sarrazin & Sperfeld, 2022). However, species/clones may respond differently to the variable diets, according to their sensitivity and nutritional requirements (Fileto et al., 2007; Ferrão-Filho et al., 2019).
The two experimental designs had some important differences that should be pointed out. In the first one we had the purpose to exclude que nutritional effect by providing a non-limiting amount of nutritious food (green algae) at a constant ratio while varying the concentration of cyanobacteria (0 to 1.5 mg C L
-1), thus the total food concentration varied (from 0.4 to 1.9 mg C L
-1) and so the proportion of nutritious food in the diet varied (100 to 21%). The second one was performed with variable proportions between cyanobacteria and green algae so the total food concentration remained constant at 1.0 mg C L
-1. Therefore, the equivalent concentration of cyanobacteria in both assays do not represent the same proportion of the total food in each assay. For example, at 1.0 mg C L
-1 of cyanobacteria in assay 1 there was ~30% of nutritious food, while at 0.9 mg C L
-1 in assay 2 there was only 10% of nutritious food (see
Table 1). Thus, at equivalent cyanobacterial concentrations, nutritional restriction was much higher in assay 2.
The negative effects on reproductive parameters and somatic growth of D. laevis (RD) were more pronounced in the experiment without nutritional restriction. Therefore, only the low nutritional value of cyanobacteria does not explain the reduction in fecundity and somatic growth, since these organisms received a sufficient amount of food with high nutritional quality, probably evidencing a toxic effect due to the presence of STX. Also, the fact that at equivalent concentrations of cyanobacteria there was a higher concentration of nutritious food in assay 1 than in assay 2 corroborates this hypothesis. In a previous study, a D. laevis clone isolated from Lagoa da Pampulha-MG showed reduced fecundity when chronically exposed to CYRF-1 strain (STXs producer) and NPCS-1 strain (non-STXs producer), both strains of the species R. raciborskii. However, the effects of the CYRF-1 strain were more pronounced than those of the NPCS-1 strain, showing that besides toxicity, nutritional constraints must play a role on the reduced fecundity (Restani & Fonseca, 2014).
On the other hand, in the life table assay with variable nutritional food availability, D. laevis clones did not show much difference in reproductive parameters, both being equally affected only in high share of cyanobacteria (90%). At lower share of R. raciborskii (25 – 50%) D. laevis seems to cope with the presence of this poor food resource, showing good growth rates. This species can coexist for long periods with cyanobacterial blooms, but its total lipid reserves decrease (Pinto-Coelho et al., 2003), which could explain the reduced fitness of these animals when they were exposed to treatments with higher proportions of cyanobacteria. In studies with cylindrospermopsin-producing (R. raciborskii and Aphanizomenon ovalisporum) and non-producing (R. raciborskii) strains, effects on survival and somatic growth of D. magna, as well as damage to the digestive epithelium were also observed (Nogueira et al., 2006).
Among the species, D. gessneri showed the most contrasting responses between the two life table assays performed. While in the assay without nutritional restriction D. gessneri was the most tolerant among the species tested, in the trial with the nutritional restriction it suffered the most deleterious effects on reproduction, having extremely low fecundity, total offspring production and negative rate of population increase in the highest proportion (90%) of cyanobacteria in the diet.
Regarding survival, in both life table assays, the results did not show resistance of D. gessneri as in previous studies where this species showed high survival. Using the R. raciborskii T3 strain (STXs producer), authors found that it did not interfere with the survival or reproduction of D. gessneri, while the cladocerans D. pulex and M. micrura, had reductions in survivorship and population increase rate (Costa et al., 2013). Therefore, it is possible that D. gessneri is more tolerant to toxins and sensitive to the lack of nutritious food.
Using mixtures of cyanobacteria and green algae (Scenedesmus obliquus), effects of the CYRF-1 strain on D. magna, such as reductions in survival, body growth, population growth rate, and clearance rate, have been reported (Soares et al., 2009). As for D. laevis in our experiment, a proportion of 1:1 (50:50%) between green algae and cyanobacteria was enough to reverse the effects of CYRF strain on D. magna. In an assay with D. pulicaria being exposed to treatments containing S. obliquus or filaments of a non-toxic strain of R. raciborskii, the presence of a non-toxic strain affected the reproductive success of Daphnia, not only through reduced fecundity but also through an increase in the number of aborted eggs (Bednarska & Slusarczyk, 2013). This study, as well as that of Restani and Fonseca (2014) with the non-toxic strain NPCS-1, reinforce the fact that the observed effects on reproduction are not explained only by the toxins and nutritional issues are involved.
In a study with another R. raciborskii strain (CYLCAM-2) isolated from the same reservoir and with the same experimental design of our study, but with more limited food condition (total food = 200 g C L-1), the four cladocerans tested showed decreased fitness even in the lower proportion of cyanobacteria (25%), with differences in sensitivity between species/clones (Ferrão-Filho et al., 2019). Apparently, D. laevis (IB clone) was one of the most resistant species being more affected only at the highest proportion of cyanobacteria in the diet (100%). Thus, this clone seems to be more affected by nutritional deficiency than toxicity. In the present study, D. gessneri showed to be very sensitive to lack of nutritious food, having its fitness severely depressed in the highest proportion of cyanobacteria (90% of the diet).
The grazing assay showed decreased total clearance rate only for D. laevis (IB) clone when exposed to a diet with cyanobacteria. When the clearance rate of green algae was analyzed, there was a significant reduction for both D. laevis clones. Probably, those clones reduced its clearance rate due to the presence of toxins, poor nutritional quality, and/or morphology of R. raciborskii. In a previous study, D. magna, even presenting feeding inhibition, continued with a positive population increase rate (r), indicating that these organisms can thrive in eutrophic environments dominated by cyanobacteria (Soares et al., 2009). Similarly, it is possible that despite the reduction in clearance rate D. laevis maintained its fitness, since in the life table assay with variable food availability both clones were less affected in reproductive parameters. Curiously, the clearance rate of cyanobacteria for the clone RD was higher in the 90% treatment, showing that this clone is capable of ingesting high amounts of cyanobacterial filaments. This may have compensated the decreased ingestion of green algae.
It has been suggested that filaments (with an average length of 426.7 ± 208.4 µm in this study) could impair the ingestion of food by Daphnia (Wilson et al., 2006; Panosso & Lürling, 2010; Bednarska et al., 2014). Although the D. laevis clones are medium-sized (1.8 mm), ingestion of large filaments has been observed in other studies in cladocerans of similar size. D. pulex (0.6 mm) when exposed to STXs-producing strains, was significantly more affected in its clearance rate by ingesting long filaments of MVCC19 (863.5 ± 483.3 µm) and CYRF-01 (338.7 ± 190.3 µm) than the copepod Notodiaptomus iheringi (1.1 mm) (Fabre et al., 2017). Copepods seem to cope better with the presence of filamentous cyanobacteria (Ger et al., 2011). They are able to cut long filaments (Bouvy et al., 2000) and can avoid toxic strains (Panosso et al., 2003). Cladocerans, however, are generalists filter-feeders and filamentous cyanobacteria can clog their filtering apparatus (DeMott et al., 2001). Nevertheless, D. laevis (RD) seems to ingest the filaments of R. raciborskii, becoming more exposed to the toxins of the strain,. This is the reason it may have been more affected in the life table experiment with no food restriction. In addition, toxin production of R. raciborskii can also reduce the fitness of copepods, especially when exposed to nutritionally restricted diets, indicating that toxins interfere more with herbivory than filament size (Rangel et al., 2016; 2020).
In the grazing assay, D. gessneri showed a tendency to increase the total clearance rate in the treatments with cyanobacteria, which may be attributed to an increase in the ingestion of green algae. By the other hand, this species showed a reduced clearance rate of cyanobacteria in the 50% treatment, but a comparable rate to green algae in the 90% treatment, which may explain its reduced fitness, as it may be able to ingest the poor food in an environment with low nutritious food availability. DeMott (1999) has shown that Daphnia spp. exposed to toxic Microcystis exhibited inhibition of feeding after 1 h of exposure to toxic food, but had substantially recovered after 24 h in the same mixtures (“hunger effect”). Thus, it is possible that D. gessneri elevated its feeding rate in 90% cyanobacterial diet to decrease the risk of starvation, as the availability of nutritious food was too low.
Some studies have already demonstrated selectivity in cladocerans. Although in general they are considered non-selective filter-feeders, experiments with D. ambigua revealed that they can avoid toxic strains of M. aeruginosa (Tillmanns et al., 2011). Using video recording, D. pulicaria was able to discriminate between toxic Microcystis strains and the green algae Scenedesmus, slowing its mandible movement in the presence of this cyanobacterium (Ghadouani et al., 2004). Higher clearance rates of the green algae S. capricornutum were observed in treatments with a mixture of Microcystis aeruginosa (MIRF-1) and R. raciborskii (CYRF-1) than in all other treatments for D. laevis and D. similis, suggesting no feeding inhibition of these cladocerans and perhaps selective clearance (Ferrão-Filho et al., 2017). Therefore, the possibility that D. gessneri has the ability to discriminate toxic cyanobacteria at the same time that it increases it is clearance rate of green algae cannot be discarded.
Differences between Daphnia clones may be dependent on the genotype and life history of these organisms (Tillmanns et al., 2008; Bednarska et al., 2014). The D. gessneri clone was isolated from the Apertado Reservoir, located in Mucugê (BA), where there are no reported occurrences of cyanobacterial blooms. The D. laevis clone (IB) was isolated from Ibirité Lagoon (MG), where phytoplankton is dominated by Microcystis spp. (Garcia et al., 2009). The D. laevis clone (RD) was isolated from Lake Jacaré (Rio Doce State Park, MG) where R. raciborskii has already been reported to occur (Miranda et al., 2013). Cyanobacteria dominating environments for long periods can act as a natural selection force on Daphnia populations, with genotypes with higher growth and reproduction rates being favored (Mello & Maia-Barbosa, 2015). Thus, differences in sensitivity among Daphnia clones to cyanobacteria may be a selective factor, determining the composition of zooplankton, favoring the competitive ability of some tolerant species over other sensitive species (Ferrão-Filho et al., 2007).
In conclusion, we found adverse effects of toxic R. raciborskii CYLCAM-1 on Daphnia fitness, even when exposed to a diet without food restriction. Particularly, in the life table assay with constant nutritious food, a reduction in the reproduction parameters for all species was observed, evidencing that the toxic effect of CYLCAM-1 overrides the nutritional effect. In the assay with variable nutritious food, D. gessneri presented as the most sensitive species in the diet with the highest proportion of cyanobacteria, which suggests a strongest nutritional effect in that clone. Daphnia clones that showed a high clearance rate (i.e. continued feeding even in diets with high cyanobacterial ratios) appeared to be less tolerant to exposure to CYLCAM-01. While clones with low clearance rate that reduced their clearance rate in the presence of cyanobacteria, were more tolerant and seem to avoid ingesting filaments, being able to coexist with R. raciborskii in conditions where dominance of this cyanobacteria occurs.
In a scenario with climate change and increased dominance of cyanobacteria, the zooplankton community may suffer alterations of its composition due to different levels of sensitivity among organisms. In a natural environment, reduced predation of zooplankton on phytoplankton may favor the dominance of toxic cyanobacteria, causing the extinction of sensitive species and reduction in the diversity of a natural ecosystem. Studies investigating behavioral resistance (such as rejecting filaments) or physiological tolerance to STX-producing strains are needed for more elucidation on the feeding selectivity and ecophysiology of different Daphnia clones.