Gene expression and activity of digestive proteases in Daphnia: effects of cyanobacterial protease inhibitors

Comparative transcriptome analysis of two Daphnia galeata genotypes displaying contrasting phenotypic variation induced by fish kairomones in the same environment of the Han River, Korea

BACKGROUND

Phenotypic plasticity is a crucial adaptive mechanism that enables organisms to modify their traits in response to changes in their environment. Predator-induced defenses are an example of phenotypic plasticity observed across a wide range of organisms, from single-celled organisms to vertebrates. In addition to morphology and behavior, these responses also affect life-history traits. The crustacean Daphnia galeata is a suitable model organism for studying predator-induced defenses, as it exhibits life-history traits changes under predation risk. To get a better overview of their phenotypic plasticity under predation stress, we conducted RNA sequencing on the transcriptomes of two Korean Daphnia galeata genotypes, KE1, and KB11, collected in the same environment.

RESULTS

When exposed to fish kairomones, the two genotypes exhibited phenotypic variations related to reproduction and growth, with opposite patterns in growth-related phenotypic variation. From both genotypes, a total of 135,611 unigenes were analyzed, of which 194 differentially expressed transcripts (DETs) were shared among the two genotypes under predation stress, which showed consistent, or inconsistent expression patterns in both genotypes. Prominent DETs were related to digestion and reproduction and consistently up-regulated in both genotypes, thus associated with changes in life-history traits. Among the inconsistent DETs, transcripts encode vinculin (VINC) and protein obstructor-E (OBST-E), which are associated with growth; these may explain the differences in life-history traits between the two genotypes. In addition, genotype-specific DETs could explain the variation in growth-related life-history traits between genotypes, and could be associated with the increased body length of genotype KE1.

CONCLUSIONS

The current study allows for a better understanding of the adaptation mechanisms related to reproduction and growth of two Korean D. galeata genotypes induced by predation stress. However, further research is necessary to better understand the specific mechanisms by which the uncovered DETs are related with the observed phenotypic variation in each genotype. In the future, we aim to unravel the precise adaptive mechanisms underlying predator-induced responses.

Integrative analysis of microbiome and metabolome reveals the linkage between gut microbiota and carp growth

Gut microbiota and their metabolites are increasingly recognized for their crucial role in regulating the health and growth of the host. The mechanism by which the gut microbiome affects the growth rate of fish (Cyprinus carpio) in the rice-fish coculture system, however, remains unclear. In this study, the gut contents of the fast-growing and slow-growing (FG and SG) carp were collected from the rice-fish coculture system for both the fish gut microbiome and metabolome analyses. High throughput 16 S rRNA gene sequencing showed that the overall gut microbiota of FG group was distinct from that of SG group. For example, the cyanobacteria were highly enriched in the guts of SG carp (18.61%), in contrast, they only represented a minor fraction of gut microbiota for FG group (<0.20%). The liquid chromatography-mass spectrometry (LC-MS)-based metabolomics analysis revealed that 191 identified metabolites mostly located in 18 KEGG pathways were differentially present between the two groups, of which more than 50% of these metabolites were involved in lipid and amino acids metabolism. Compared with the FG group, the gut microbiota of SG group significantly enriched the metabolic pathways involved in the steroid (hormone) biosynthesis, whereas reducing those associated with beta-alanine metabolism, biosynthesis of unsaturated fatty acids and bile secretion. The enrichment and depletion of these metabolic pathways resulted in an increase in steroid metabolites and a decrease in the concentration of spermidine, which may have a major impact on the growth rate of carp. The metabolome results were further supported by the predicated KEGG functions of the gut microbiomes of the two groups, pointing out that the gut microbiota could substantially affect the growth of fish via their unique metabolic functions. Together, our integrated fish gut microbiome and metabolome analysis has substantial implications for the development of engineered microbiome technologies in aquaculture.

Negative Effects of Cyanotoxins and Adaptative Responses of Daphnia

The plethora of cyanobacterial toxins are an enormous threat to whole ecosystems and humans. Due to eutrophication and increases in lake temperatures from global warming, changes in the distribution of cyanobacterial toxins and selection of few highly toxic species/ strains are likely. Globally, one of the most important grazers that controls cyanobacterial blooms is Daphnia, a freshwater model organism in ecology and (eco)toxicology. Daphnia-cyanobacteria interactions have been studied extensively, often focusing on the interference of filamentous cyanobacteria with Daphnia's filtering apparatus, or on different nutritional constraints (the lack of essential amino acids or lipids) and grazer toxicity. For a long time, this toxicity only referred to microcystins. Currently, the focus shifts toward other deleterious cyanotoxins. Still, less than 10% of the total scientific output deals with cyanotoxins that are not microcystins; although these other cyanotoxins can occur just as frequently and at similar concentrations as microcystins in surface water. This review discusses the effects of different cyanobacterial toxins (hepatotoxins, digestive inhibitors, neurotoxins, and cytotoxins) on Daphnia and provides an elaborate and up-to-date overview of specific responses and adaptations of Daphnia. Furthermore, scenarios of what we can expect for the future of Daphnia-cyanobacteria interactions are described by comprising anthropogenic threats that might further increase toxin stress in Daphnia.

Food quantity and quality shapes reproductive strategies of Daphnia

In freshwater environments, one of the challenges aquatic grazers face are periods of suboptimal food quantity and quality. In a life table experiment, the effects of food quantity (a gradient of algae concentration) and quality (a diet of cyanobacteria) on the life histories and resource allocation strategy in Daphnia magna were tested. Growth‐related traits were similarly affected under different food regimes while the reproductive strategies differed in animals exposed to low food quantity and quality. The per‐clutch investment (clutch volume) did not differ between Daphnia fed with cyanobacteria and underfed mothers, but resources were differently allocated; underfed mothers increased their per‐offspring investment by producing fewer, but larger eggs, whereas cyanobacteria‐fed mothers invested in a greater number of eggs of smaller size. I argue that both strategies of resource allocation (number vs. size of eggs) may be adaptive under the given food regime. The results of the study show that the cyanobacteria diet‐driven fitness losses are comparable to losses caused by food quantity, which is only slightly above the growth capability threshold for Daphnia.

Proteome changes in an aquatic invertebrate consumer in response to different nutritional stressors

Nutrient imbalances in zooplankton are caused by the differences in elemental content of producers and the demand for elements in consumers, which alter the life-history traits in consumers. Changes in life-history traits are mediated through metabolic pathways that affect gene expression and the metabolome. However, less is known about proteomic changes to elemental-limitation in zooplankton. Here, we grew Daphnia pulex under high food quantity and quality (HF), low food quantity (LF), and phosphorus (P)-limited (PL) diets for six days and measured growth, elemental composition, and the proteome. Daphnids in both LF and PL diets grew less. Animals in LF diets had less carbon (C), while daphnids in PL diets had less P compared to HF fed animals. In total, we identified 1719 proteins that were used in a partial least squares regression discriminant analysis (PLS-DA). Focusing on a subset of the proteome, the PLS-DA resulted in a clear separation between animals fed HF diets and PL and LF diets. Many proteome changes in nutrient-limited diets are associated with growth, reproduction, lipid metabolism, and nutrient assimilation. Regardless of the limiting nutrient, there were less hemoglobin and small subunit processome component proteins compared to HF fed animals. Daphnids fed LF diets had less vitellogenin fused superoxide dismutase and more lipid-droplet hydrolase, whereas Daphnia fed PL diets had higher abundances of cytochrome P450 and serine protease. Our proteome results compliment other "omic" studies that could be used to study Daphnia physiology in lakes.

Gene expression changes in Daphnia magna following waterborne exposure to cyanobacterial strains from the genus Nostoc

Cyanobacteria can produce highly potent cyanotoxins, however, limited information is provided about their toxicity mechanisms in exposed aquatic invertebrates at the molecular level. In the present study, the effects of cyanobacterial strains from the genus Nostoc (Nostoc Z1 and Nostoc 2S3B) in Daphnia magna after waterborne exposure were investigated. Examined endpoints included immobilization (survival) in acute toxicity tests and selected gene expression changes (cyp314, cyp360A8, gst, p-gp, vtg) analyzed by the quantitative real-time polymerase chain reaction (RT-PCR). In addition, enzyme-linked immunosorbent assay (ELISA) was performed to determine whether the observed changes could be due to the presence of microcystins, the most widespread group of cyanotoxins. The results of acute toxicity tests have shown only minor changes in survival rates, which have not exceeded 20% after 48 h of exposure to either strain. On the other hand, significant changes were recorded in molecular responses of Daphnia to tested strains. Treatment with the aquatic strain Nostoc Z1 altered the expression levels of all analyzed genes. Both strains caused a significant p-glycoprotein (p-gp) induction at 75 µg ml^-1 which suggests the involvement of p-gp mediated multixenobiotic resistance mechanism (MXR) in facilitating excretion of toxic cyanobacterial compounds in daphnids. Additionally, these strains caused an increase in the expression levels of cyp360A8, indicating that genes related to detoxification processes could be sensitive indicators of cyanobacterial toxicity. Statistically significant induction of cyp314, as well as increases in expression of gst and vtg, were observed only after exposure to Nostoc Z1. This study indicates the potential of certain cyanobacterial metabolites to modify the expression of toxicant responsive genes involved in phase I and phase III of the xenobiotic metabolism, as well as possible interference with growth and reproduction in D. magna. Low microcystin concentrations found in both samples suggest that these cyanotoxins were not responsible for the detected toxic effects.

Oscillatoria sp. as a Potent Anti-phytopathogenic Agent and Plant Immune Stimulator Against Root-Knot Nematode of Soybean cv. Giza 111

BACKGROUND

Plant-parasitic nematodes are one of the major constraints to soybean production around the world. Plant-parasitic nematodes cause an estimated $78 billion in annual crop losses worldwide, with a 10-15% crop yield loss on average. Consequently, finding and applying sustainable methods to control diseases associated with soybean is currently in serious need.

METHODS

In this study, we isolated, purified, characterized, and identified a novel cyanobacterial strain Oscillatoria sp. (blue-green alga). Based on its microscopic examination and 16S rRNA gene sequence, the aqueous and methanolic extracts of Oscillatoria were used to test their nematicidal activity against Meloidogyne incognita hatchability of eggs after 72 h of exposure time and juvenile mortality percentage in vitro after 24, 48, and 72 h of exposure time and reduction percentage of galls, eggmass, female number/root, and juveniles/250 soil. Also, the efficacy of the extract on improving the plant growth parameter and chlorophyll content under greenhouse conditions on soybean plant cv. Giza 111 was tested. Finally, the expression of PR-1, PR-2, PR-5, and PR15 (encoding enzymes) genes contributing to plant defense in the case of M. incognita invasion was studied and treated with Oscillatoria extract.

RESULTS

The aqueous and methanolic extracts of Oscillatoria sp. had nematicidal activity against M. incognita. The percentage of mortality and egg hatching of M. incognita were significantly increased with the increase of time exposure to Oscillatoria extract 96.7, 97, and 98 larvae mortality % with S concentration after 24, 48, and 72 h of exposure time. The aqueous extract significantly increased the percentage of Root-Knot nematodes (RKN) of egg hatching, compared with Oxamyl and methanol extract at 96.7 and 97% after 72 h and 1 week, respectively. With the same concentration in the laboratory experiment. Furthermore, water extracts significantly reduced the number of galls in soybean root, egg masses, and female/root by 84.1, 87.5, and 92.2%, respectively, as well as the percentage of J2s/250 g soil by 93.7%. Root, shoot lengths, dry weight, number of pods/plant, and chlorophyll content of soybean treated with Oscillatoria water extract were significantly higher than the control increasing by 70.3, 94.1, 95.5, and 2.02%, respectively. The plant defense system's gene expression was tracked using four important pathogenesis-related genes, PR-1, PR-2, PR-5, and PR15, which encode enzymes involved in plant defense.

CONCLUSIONS

Oscillatoria extract is a potential nematicide against root-knot nematode invasion in soybean.

More Light Please: Daphnia Benefit From Light Pollution by Increased Tolerance Toward Cyanobacterial Chymotrypsin Inhibitors

Cryptochromes are evolutionary ancient blue-light photoreceptors that are part of the circadian clock in the nervous system of many organisms. Cryptochromes transfer information of the predominant light regime to the clock which results in the fast adjustment to photoperiod. Therefore, the clock is sensitive to light changes and can be affected by anthropogenic Artificial Light At Night (ALAN). This in turn has consequences for clock associated behavioral processes, e.g., diel vertical migration (DVM) of zooplankton. In freshwater ecosystems, the zooplankton genus Daphnia performs DVM in order to escape optically hunting predators and to avoid UV light. Concomitantly, Daphnia experience circadian changes in food-supply during DVM. Daphnia play the keystone role in the carbon-transfer to the next trophic level. Therefore, the whole ecosystem is affected during the occurrence of cyanobacteria blooms as cyanobacteria reduce food quality due to their production of digestive inhibitors (e.g., protease inhibitors). In other organisms, digestion is linked to the circadian clock. If this is also the case for Daphnia, the expression of protease genes should show a rhythmic expression following circadian expression of clock genes (e.g., cryptochrome 2). We tested this hypothesis and demonstrated that gene expression of the clock and of proteases was affected by ALAN. Contrary to our expectations, the activity of one type of proteases (chymotrypsins) was increased by ALAN. This indicates that higher protease activity might improve the diet utilization. Therefore, we treated D. magna with a chymotrypsin-inhibitor producing cyanobacterium and found that ALAN actually led to an increase in Daphnia’s growth rate in comparison to growth on the same cyanobacterium in control light conditions. We conclude that this increased tolerance to protease inhibitors putatively enables Daphnia populations to better control cyanobacterial blooms that produce chymotrypsin inhibitors in the Anthropocene, which is defined by light pollution and by an increase of cyanobacterial blooms due to eutrophication.

Daphnia magna Gut-Specific Transcriptomic Responses to Feeding Inhibiting Chemicals and Food Limitation

Transcriptomic responses combined with apical adverse ecologically relevant outcomes have proven to be useful to unravel and anchor molecular mechanisms of action to adverse outcomes. This is the case for feeding inhibition responses in the model ecotoxicological species Daphnia magna. The aim of the present study was to assess the transcriptomic responses in guts dissected from D. magna individuals exposed to concentrations of selected compounds that inhibit feeding and compare them with the responses associated to 2 levels of food restriction (low food and starvation). Chemical treatments included cadmium, copper, fluoranthene, λ-cyhalothrin, and the cyanotoxin anatoxin-a. Although the initial hypothesis was that exposure to chemical feeding inhibitors should elicit similar molecular responses as food limitation, the corresponding gut transcriptomic responses differed significantly. In moderate food limitation conditions, D. magna individuals increased protein and carbohydrate catabolism, likely to be used as energetic sources, whereas under severe starving conditions most metabolism-related pathways appeared down-regulated. Treatment with chemical feeding inhibitors promoted cell turnover-related signaling pathways in the gut, probably to renew tissue damage caused by the reported oxidative stress effects of these compounds, and inhibited the transcription of gut digestive gene enzymes and energetic metabolic pathways. We conclude that chemical feeding inhibitors, rather than mimicking the physiological response to low- or no-food conditions, cause specific toxic effects, preventing Daphnia both from feeding and from adjusting its metabolism to the resulting low energy intake. Environ Toxicol Chem 2021;40:2510-2520. © 2021 SETAC.

Daphnia's Adaptive Molecular Responses to the Cyanobacterial Neurotoxin Anatoxin-α Are Maternally Transferred

Cyanobacterial blooms are an omnipresent and well-known result of eutrophication and climate change in aquatic systems. Cyanobacteria produce a plethora of toxic secondary metabolites that affect humans, animals and ecosystems. Many cyanotoxins primarily affect the grazers of phytoplankton, e.g., Daphnia. The neurotoxin anatoxin-α has been reported world-wide; despite its potency, anatoxin-α and its effects on Daphnia have not been thoroughly investigated. Here, we investigated the effects of the anatoxin-α-producing Tychonema on life-history parameters and gene expression of nicotine-acetylcholine receptors (NAR), the direct targets of anatoxin-α, using several D. magna clones. We used juvenile somatic growth rates as a measure of fitness and analyzed gene expression by qPCR. Exposure to 100% Tychonema reduced the clones' growth rates and caused an up-regulation of NAR gene expression. When 50% of the food consisted of Tychonema, none of the clones were reduced in growth and only one of them showed an increase in NAR gene expression. We demonstrate that this increased NAR gene expression can be maternally transferred and that offspring from experienced mothers show a higher growth rate when treated with 50% Tychonema compared with control offspring. However, the addition of further (anthropogenic) stressors might impair Daphnia's adaptive responses to anatoxin-α. Especially the presence of certain pollutants (i.e., neonicotinoids), which also target NARs, might reduce Daphnia's capability to cope with anatoxin-α.

Cyanobacterial anabaenopeptin-B, microcystins and their mixture cause toxic effects on the behavior of the freshwater crustacean Daphnia magna (Cladocera)

Comparison of the toxic effects caused by the pure cyanobacterial cyclic hexapeptide anabaenopeptin-B (AN-B), the heptapeptides: microcystin-LR (MC-LR) and MC-LF as well as a binary mixture of AN-B with MC-LR on the swimming speed and hopping frequency - essential activities of Daphnia, was experimentally determined. Till now, no information on behavioral effects of AN-B and its mixture with microcystins, commonly produced by cyanobacteria, was available. Also MC-LF effect on aquatic crustaceans was determined for the first time. The results showed that AN-B exerted considerable inhibition of D. magna swimming speed and hopping frequency similar to MC-LR and MC-LF. The mixture of AN-B and MC-LR caused stronger toxic effects, than the individual oligopeptides used at the same concentration. The much lower 48 h- EC₅₀ value of the AN-B and MC-LR mixture (0.95 ± 0.12 μg/mL) than those of individual oligopeptides AN-B (6.3 ± 0.63 μg/mL), MC-LR (4.0 ± 0.27 μg/mL), MC-LF (3.9 ± 0.20 μg/mL) that caused swimming speed inhibition explains the commonly observed stronger toxicity of complex crude cyanobacterial extracts to daphnids than individual microcystins. The obtained results indicated that AN-B, microcystins and their mixture exerted time- and concentration-dependent motility disturbances of crustaceans and they can be good candidates for evaluation of toxicity in early warning systems. Other cyanobacterial oligopeptides beyond microcystins should be considered as a real threat for aquatic organisms.

Cell free Microcystis aeruginosa spent medium affects Daphnia magna survival and stress response

Primary consumers in freshwater ecosystems, such as the zooplankton organism Daphnia magna, are highly affected by cyanobacteria, both as they may use it as a food source but also by cyanobacterial metabolites present in the water. Here, we investigate the impacts of cyanobacterial metabolites focussing on the environmental realistic scenario of the naturally released mixture without crushing cyanobacterial cells or their uptake as food. Therefore, D. magna were exposed to two concentrations of cell free cyanobacterial spent medium from Microcystis aeruginosa PCC 7806 to represent higher and lower ecologically-relevant concentrations of cyanobacterial metabolites. Including microcystin-LR, 11 metabolites have been detected of which 5 were quantified. Hypothesising concentration and time dependent negative impact, survival, gene expression marking digestion and metabolism, oxidative stress response, cell cycle and molting as well as activities of detoxification and antioxidant enzymes were followed for 7 days. D. magna suffered from oxidative stress as both catalase and glutathione S-transferase enzyme activities significantly decreased, suggesting enzyme exhaustibility after 3 and 7 days. Moreover, gene-expressions of the 4 stress markers (glutathione S-transferase, glutathione peroxidase, catalase and thioredoxin) were merely downregulated after 7 days of exposure. Energy allocation (expression of glyceraldehyde-3-phosphate dehydrogenase) was increased after 3 days but decreased as well after 7 days exposure. Cell cycle was impacted time dependently but differently by the two concentrations, along with an increasing downregulation of myosin heavy chain responsible for cell arrangement and muscular movements. Deregulation of nuclear hormone receptor genes indicate that D. magna hormonal steering including molting seemed impaired despite no detection of microviridin J in the extracts. As a consequence of all those responses and presumably of more than investigated molecular and physiological changes, D. magna survival was impaired over time, in a concentration dependent manner. Our results confirm that besides microcystin-LR, other secondary metabolites contribute to negative impact on D. magna survival and stress response.

Locally adapted gut microbiomes mediate host stress tolerance

While evidence for the role of the microbiome in shaping host stress tolerance is becoming well-established, to what extent this depends on the interaction between the host and its local microbiome is less clear. Therefore, we investigated whether locally adapted gut microbiomes affect host stress tolerance. In the water flea Daphnia magna, we studied if the host performs better when receiving a microbiome from their source region than from another region when facing a stressful condition, more in particular exposure to the toxic cyanobacteria Microcystis aeruginosa. Therefore, a reciprocal transplant experiment was performed in which recipient, germ-free D. magna, isolated from different ponds, received a donor microbiome from sympatric or allopatric D. magna that were pre-exposed to toxic cyanobacteria or not. We tested for effects on host life history traits and gut microbiome composition. Our data indicate that Daphnia interact with particular microbial strains mediating local adaptation in host stress tolerance. Most recipient D. magna individuals performed better when inoculated with sympatric than with allopatric microbiomes. This effect was most pronounced when the donors were pre-exposed to the toxic cyanobacteria, but this effect was also pond and genotype dependent. We discuss how this host fitness benefit is associated with microbiome diversity patterns.

Next-generation sequencing of DNA from resting eggs: signatures of eutrophication in a lake's sediment

Hatching resting stages of ecologically important organisms such as Daphnia from lake sediments, referred to as resurrection ecology, is a powerful approach to assess changes in alleles and traits over time. However, the utility of the approach is constrained by a few obstacles, including low and/or biased hatching among genotypes. Here, we eliminated such bottlenecks by investigating DNA sequences isolated directly (i.e. without hatching) from resting eggs found in the sediments of Lake Constance spanning pre-, peri-, and post-eutrophication. While we expected genome-wide changes, we specifically expected changes in alleles related to pathways involved in mitigating effects of cyanobacterial toxins. We used pairwise F_ST-analyses to identify transcripts that showed strongest divergence among the four different populations and a clustering analysis to identify correlations between allele frequency shifts and changes in abiotic and biotic lake parameters. In a cluster that correlated with the increased abundance of cyanobacteria in Lake Constance we find genes that have been reported earlier to be differentially expressed in response to the cyanobacterial toxin microcystin and to microcystin-free cyanobacteria. We further reveal the enrichment of gene ontology terms that have been shown to be involved in microcystin-related responses in other organisms but not yet in Daphnia and as such are candidate loci for adaptation of natural Daphnia populations to increased cyanobacterial abundances. In conclusion this approach of investigating DNA extracted from Daphnia resting stages allowed to determine frequency changes of loci in a natural population over time.

How Copepods Can Eat Toxins Without Getting Sick: Gut Bacteria Help Zooplankton to Feed in Cyanobacteria Blooms

Toxin-producing cyanobacteria can be harmful to aquatic biota, although some grazers utilize them with often beneficial effects on their growth and reproduction. It is commonly assumed that gut microbiota facilitates host adaptation to the diet; however, the evidence for adaptation mechanisms is scarce. Here, we investigated the abundance of mlrA genes in the gut of the Baltic copepods Acartia bifilosa and Eurytemora affinis during cyanobacteria bloom season (August) and outside it (February). The mlrA genes are unique to microcystin and nodularin degraders, thus indicating the capacity to break down these toxins by the microbiota. The mlrA genes were expressed in the copepod gut year-round, being >10-fold higher in the summer than in the winter populations. Moreover, they were significantly more abundant in Eurytemora than Acartia. To understand the ecological implications of this variability, we conducted feeding experiments using summer- and winter-collected copepods to examine if/how the mlrA abundance in the microbiota affect: (1) uptake of toxic Nodularia spumigena, (2) uptake of a non-toxic algal food offered in mixtures with N. spumigena, and (3) concomitant growth potential in the copepods. The findings provide empirical evidence that the occurrence of mlrA genes in the copepod microbiome facilitates nutrient uptake and growth when feeding on phytoplankton mixtures containing nodularin-producing cyanobacteria; thus, providing an adaptation mechanism to the cyanobacteria blooms.

The complexity of co-limitation: nutrigenomics reveal non-additive interactions of calcium and phosphorus on gene expression in Daphnia pulex

Many lakes across Canada and northern Europe have experienced declines in ambient phosphorus (P) and calcium (Ca) supply for over 20 years. While these declines might create or exacerbate nutrient limitation in aquatic food webs, our ability to detect and quantify different types of nutrient stress on zooplankton remains rudimentary. Here, we used growth bioassay experiments and whole transcriptome RNAseq, collectively nutrigenomics, to examine the nutritional phenotypes produced by low supplies of P and Ca separately and together in the freshwater zooplankter Daphnia pulex. We found that daphniids in all three nutrient-deficient categories grew slower and differed in their elemental composition. Our RNAseq results show distinct responses in singly limited treatments (Ca or P) and largely a mix of these responses in animals under low Ca and P conditions. Deeper investigation of effect magnitude and gene functional annotations reveals this patchwork of responses to cumulatively represent a co-limited nutritional phenotype. Linear discriminant analysis identified a significant separation between nutritional treatments based upon gene expression patterns with the expression patterns of just five genes needed to predict animal nutritional status with 99% accuracy. These data reveal how nutritional phenotypes are altered by individual and co-limitation of two highly important nutritional elements (Ca and P) and provide evidence that aquatic consumers can respond to limitation by more than one nutrient at a time by differentially altering their metabolism. This use of nutrigenomics demonstrates its potential to address many of the inherent complexities in studying interactions between multiple nutritional stressors in ecology and beyond.

Insights into the genetic basis of predator-induced response in Daphnia galeata

Phenotypic plastic responses allow organisms to rapidly adjust when facing environmental challenges-these responses comprise morphological, behavioral but also life-history changes. Alteration of life-history traits when exposed to predation risk have been reported often in the ecological and genomic model organism Daphnia. However, the molecular basis of this response is not well understood, especially in the context of fish predation. Here, we characterized the transcriptional profiles of two Daphnia galeata clonal lines with opposed life histories when exposed to fish kairomones. First, we conducted a differential gene expression, identifying a total of 125 candidate transcripts involved in the predator-induced response, uncovering substantial intraspecific variation. Second, we applied a gene coexpression network analysis to find clusters of tightly linked transcripts revealing the functional relations of transcripts underlying the predator-induced response. Our results showed that transcripts involved in remodeling of the cuticle, growth, and digestion correlated with the response to environmental change in D. galeata. Furthermore, we used an orthology-based approach to gain functional information for transcripts lacking gene ontology (GO) information, as well as insights into the evolutionary conservation of transcripts. We could show that our candidate transcripts have orthologs in other Daphnia species but almost none in other arthropods. The unique combination of methods allowed us to identify candidate transcripts, their putative functions, and evolutionary history associated with predator-induced responses in Daphnia. Our study opens up to the question as to whether the same molecular signature is associated with fish kairomones-mediated life-history changes in other Daphnia species.

Toward Disentangling the Multiple Nutritional Constraints Imposed by Planktothrix: The Significance of Harmful Secondary Metabolites and Sterol Limitation

The harmful bloom-forming cyanobacterium Planktothrix is commonly considered to be nutritionally inadequate for zooplankton grazers, resulting in limited top-down control. However, interactions between Planktothrix and zooplankton grazers are poorly understood. The food quality of Planktothrix is potentially constrained by morphological properties (i.e., filament formation), the production of harmful secondary metabolites, and a deficiency in essential lipids (i.e., primarily sterols). Here, we investigated the relative significance of toxin production (microcystins, carboxypeptidase A inhibitors, protease inhibitors) and sterol limitation for the performance of Daphnia feeding on one Planktothrix rubescens and one P. agardhii wild-type/microcystin knock-out mutant pair. Our data suggest that the poor food quality of both Planktothrix spp. is due to deleterious effects mediated by various harmful secondary metabolites and that the impact of sterol limitation is partially or completely superimposed by toxicity. The significance of the different factors seems to depend on the metabolite profile of the considered Planktothrix strain and the Daphnia clone that is used for the experiments. The toxin-responsive gene expression (transporter genes, gpx, and trypsin) and enzyme activity patterns revealed strain-specific food quality constraints and that Daphnia is capable of modulating its physiological responses according to the ingested Planktothrix strain. Future studies need to consider that Planktothrix-grazer interactions are simultaneously modulated by multiple factors to improve our understanding of top-down influences on Planktothrix bloom formation.

Disentangling of the ecotoxicological signal using "omics" analyses, a lesson from the survey of the impact of cyanobacterial proliferations on fishes

Omics technologies offer unprecedented perspectives for the rational investigation of complex biological systems. Indeed, omics present the ability of offering an extensive perception of the biochemistry and physiology of the cell and of any perturbing consequences of contaminants through the joint investigation of thousands of molecular responses simultaneously; then it has recently conducted to a fervent attention by research ecotoxicologists. Beyond the presentation of latest advances, exemplified here by omics investigation of cyanobacterial deleterious effects on various fishes (at various experimental and biological scales and with various analytical tools and pipeline), the present review paper re-explores the promising perspectives and also the pitfalls of such holistic investigations of the ecotoxicological response of organisms for environmental assessment.

Microviridin 1777: A Toxic Chymotrypsin Inhibitor Discovered by a Metabologenomic Approach

The toxicity of the cyanobacterium Microcystis aeruginosa EAWAG 127a was evaluated against the sensitive grazer Thamnocephalus platyurus, and the extract possessed strong activity. To investigate the compounds responsible for cytotoxicity, a series of peptides from this cyanobacterium were studied using a combined genomic and molecular networking approach. The results led to the isolation, structure elucidation, and biological evaluation of microviridin 1777, which represents the most potent chymotrypsin inhibitor characterized from this family of peptides to date. Furthermore, the biosynthetic gene clusters of microviridin, anabaenopeptin, aeruginosin, and piricyclamide were located in the producing organism, and six additional natural products were identified by tandem mass spectrometry analyses. These results highlight the potential of modern techniques for the identification of natural products, demonstrate the ecological role of protease inhibitors produced by cyanobacteria, and raise ramifications concerning the presence of novel, yet uncharacterized, toxin families in cyanobacteria beyond microcystin.

Phosphate Limitation Increases Content of Protease Inhibitors in the Cyanobacterium Microcystis aeruginosa

Increased anthropogenic nutrient input has led to eutrophication of lakes and ponds, resulting worldwide in more frequent and severe cyanobacterial blooms. In particular, enhanced availability of phosphorus (P) can promote cyanobacterial mass developments and may affect the content of secondary metabolites in cyanobacteria, such as protease inhibitors (PIs). PIs are common among cyanobacteria and have been shown to negatively affect herbivorous zooplankton. Here, we test the hypothesis that P-limitation reduces the growth of Microcystis, but increases the content of PIs. In batch culture experiments with eight different initial phosphate concentrations (5-75 µM) we determined growth, stoichiometry, and PI content of Microcystis aeruginosa NIVA Cya 43. This strain produces the protease inhibitor BN920 that is converted by chlorination to CP954, which constitutes the major PI in this strain. C:N:P-ratios of the biomass indicated variation of P-limitation with treatment and time. When normalized to biomass, the PI content varied up to nearly nineteen-fold with treatment and time and was highest in the low-P treatments, especially during the mid-exponential growth phase. However, these effects were alleviated under nitrogen co-limitation. The content of CP954 showed an inverse u-shaped response to growth rate and C:N-ratio of the cyanobacterial biomass, whereas it increased with cyanobacterial C:P. The results indicate that P-limitation supports a higher content of defensive PIs and may indirectly foster cyanobacterial blooms by increasing the negative interference of cyanobacteria with their consumers.

Timing of maternal exposure to toxic cyanobacteria and offspring fitness in Daphnia magna: Implications for the evolution of anticipatory maternal effects

Organisms that regularly encounter stressful environments are expected to use cues to develop an appropriate phenotype. Water fleas (Daphnia spp.) are exposed to toxic cyanobacteria during seasonal algal blooms, which reduce growth and reproductive investment. Because generation time is typically shorter than the exposure to cyanobacteria, maternal effects provide information about the local conditions subsequent generations will experience. Here, we evaluate if maternal effects in response to microcystin, a toxin produced by cyanobacteria, represent an inheritance system evolved to transmit information in Daphnia magna. We exposed mothers as juveniles and/or as adults, and tested the offspring's fitness in toxic and non-toxic environments. Maternal exposure until reproduction reduced offspring fitness, both in the presence and in the absence of toxic cyanobacteria. However, this effect was accompanied by a small positive fitness effect, relative to offspring from unexposed mothers, in the presence of toxic cyanobacteria. This effect was mainly elicited in response to maternal exposure to toxic cyanobacteria early in life and less so during reproduction. None of these effects were explained by changes in egg size. A meta-analysis using our and others' experiments suggests that the adaptive value of maternal effects to cyanobacteria exposure is weak at best. We suggest that the beneficial maternal effect in our study is an example of phenotypic accommodation spanning generations, rather than a mechanism evolved to transmit information about cyanobacteria presence between generations.

Heterologous expression and characterization of a novel serine protease from Daphnia magna: A possible role in susceptibility to toxic cyanobacteria

Mass developments of toxin-producing cyanobacteria are frequently observed in freshwater ecosystems due to eutrophication and global warming. These mass developments can partly be attributed to cyanobacterial toxins, such as protease inhibitors (PIs), which inhibit digestive serine proteases of Daphnia, the major herbivore of phytoplankton and cyanobacteria. To date, mechanisms of this inhibition in the gut of the crustacean Daphnia magna are not known. Here, we characterize a single serine protease, chymotrypsin 448 (CT448), which is present in the gut of the crustacean D. magna. Sequence alignments with human serine proteases revealed that CT448 has a putative N-terminal pro-peptide which is extended compared to the mammalian homologs and within this pro-peptide two N-linked glycosylation motifs were found. CT448 was heterologously expressed in Sf21 insect cells using a baculovirus expression system for optimized protein production and secretion into the medium. The protein was purified via a one-step affinity chromatography, which resulted in a protein yield of 3.45 mg/l medium. The inactive precursor (zymogen) could be activated by tryptic digestion. This is the first example of a recombinant expression of an active crustacean serine protease, which functions in the gut of Daphnia. Proteomic identification of protease cleavage sites (PICS) and hydrolysation of various synthetic substrates showed that CT448 is a chymotrypsin-like elastase. In this study, we confirm that CT448 is a target of cyanobacterial protease inhibitors. Local evolutionary modifications of CT448 might render this proteolytic enzyme less susceptible against cyanobacterial secondary metabolites and might improve the fitness of Daphnia during cyanobacterial blooms.

Spatial analysis of toxic or otherwise bioactive cyanobacterial peptides in Green Bay, Lake Michigan

Cyanobacterial harmful algal blooms (cyanoHABs) are a growing problem in freshwater systems worldwide. CyanoHABs are well documented in Green Bay, Lake Michigan but little is known about cyanoHAB toxicity. This study characterized the diversity and spatial distribution of toxic or otherwise bioactive cyanobacterial peptides (TBPs) in Green Bay. Samples were collected in 2014 and 2015 during three cruises at sites spanning the mouth of the Fox River north to Chambers Island. Nineteen TBPs were analyzed including 11 microcystin (MC) variants, nodularin, three anabaenopeptins, three cyanopeptolins and microginin-690. Of the 19 TBPs, 12 were detected in at least one sample, and 94% of samples had detectable TBPs. The most prevalent TBPs were MCRR and MCLR, present in 94% and 65% of samples. The mean concentration of all TBPs was highest in the Fox River and lower bay, however, the maximum concentration of all TBPs occurred in the same sample north of the lower bay. MCs were positively correlated with chlorophyll and negatively correlated with distance to the Fox River in all cruises along a well-established south-to-north trophic gradient in Green Bay. The mean concentration of MC in the lower bay across all cruises was 3.0 +/- 2.3 μg/L. Cyanopeptolins and anabaenopeptins did not trend with the south-north trophic gradient or varied by cruise suggesting their occurrence is driven by different environmental factors. Results from this study provides evidence that trends in TBP concentration differ by congener type over a trophic gradient.

In silico analysis and effects of environmental salinity in the expression and activity of digestive α-amylase and trypsins from the euryhaline crab Neohelice granulata

Studies on molecular characteristics and modulation of expression of α-amylase and trypsin in the hepatopancreas of intertidal euryhaline crabs are lacking. In this work, we cloned and studied by in silico approaches the characteristics of cDNA sequences for α-amylase and two trypsins isoforms, as well as the effect of environmental salinity, on gene expression and protein activities in the hepatopancreas of Neohelice granulata (Dana, 1851), which is a good invertebrate model species. The cDNA sequence of α-amylase is 1637 bp long, encoding 459 amino acid residues. Trypsin 1 and 2 are 689 and 1174 bp long, encoding 204 and 151 amino acid residues, respectively. Multiple sequence alignment of deduced protein sequences revealed the presence of conserved motifs found in other invertebrates. In crabs acclimated at 37 psu (hyporegulation), α-amylase mRNA level and total pancreatic amylase activity were higher than at 10 psu (hyperregulation) and 35 psu (osmoconformation). Trypsin 1 mRNA levels increased at 37 psu, while trypsin 2 levels decreased at 10 and 37 psu. Total trypsin activity was similar in all salinities. Our results showed a differential modulation of α-amylase and trypsin expression and total amylase activity by salinity acclimation, suggesting the occurrence of distinct mechanisms of regulation at different levels that could lead to digestive adjustments in relation to hyperregulation and (or) hyporegulation.

Cyanobacteria facilitate parasite epidemics in Daphnia

The seasonal dominance of cyanobacteria in the phytoplankton community of lake ecosystems can have severe implications for higher trophic levels. For herbivorous zooplankton such as Daphnia, cyanobacteria have poor nutritional value and some species can produce toxins affecting zooplankton survival and reproduction. Here we present another, hitherto largely unexplored aspect of cyanobacteria, namely that they can increase Daphnia susceptibility to parasites. In a 12-yr monthly time-series analysis of the Daphnia community in Greifensee (Switzerland), we observed that cyanobacteria density correlated significantly with the epidemics of a common gut parasite of Daphnia, Caullerya mesnili, regardless of what cyanobacteria species was present or whether it was colonial or filamentous. The temperature from the previous month also affected the occurrence of Caullerya epidemics, either directly or indirectly by the promotion of cyanobacterial growth. A laboratory experiment confirmed that cyanobacteria increase the susceptibility of Daphnia to Caullerya, and suggested a possible involvement of cyanotoxins or other chemical traits of cyanobacteria in this process. These findings expand our understanding of the consequences of toxic cyanobacterial blooms for lake ecosystems and might be relevant for epidemics experienced by other aquatic species.

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

High dietary quality of non-toxic cyanobacteria for a benthic grazer and its implications for the control of cyanobacterial biofilms

BACKGROUND

Mass occurrences of cyanobacteria frequently cause detrimental effects to the functioning of aquatic ecosystems. Consequently, attempts haven been made to control cyanobacterial blooms through naturally co-occurring herbivores. Control of cyanobacteria through herbivores often appears to be constrained by their low dietary quality, rather than by the possession of toxins, as also non-toxic cyanobacteria are hardly consumed by many herbivores. It was thus hypothesized that the consumption of non-toxic cyanobacteria may be improved when complemented with other high quality prey. We conducted a laboratory experiment in which we fed the herbivorous freshwater gastropod Lymnaea stagnalis single non-toxic cyanobacterial and unialgal diets or a mixed diet to test if diet-mixing may enable these herbivores to control non-toxic cyanobacterial mass abundances.

RESULTS

The treatments where L. stagnalis were fed non-toxic cyanobacteria and a mixed diet provided a significantly higher shell and soft-body growth rate than the average of all single algal, but not the non-toxic cyanobacterial diets. However, the increase in growth provided by the non-toxic cyanobacteria diets could not be related to typical determinants of dietary quality such as toxicity, nutrient stoichiometry or essential fatty acid content.

CONCLUSIONS

These results strongly contradict previous research which describes non-toxic cyanobacteria as a low quality food resource for freshwater herbivores in general. Our findings thus have strong implications to gastropod-cyanobacteria relationships and suggest that freshwater gastropods may be able to control mass occurrences of benthic non-toxic cyanobacteria, frequently observed in eutrophied water bodies worldwide.

Thermal tolerance in the keystone species Daphnia magna-a candidate gene and an outlier analysis approach

Changes in temperature have occurred throughout Earth's history. However, current warming trends exacerbated by human activities impose severe and rapid loss of biodiversity. Although understanding the mechanisms orchestrating organismal response to climate change is important, remarkably few studies document their role in nature. This is because only few systems enable the combined analysis of genetic and plastic responses to environmental change over long time spans. Here, we characterize genetic and plastic responses to temperature increase in the aquatic keystone grazer Daphnia magna combining a candidate gene and an outlier analysis approach. We capitalize on the short generation time of our species, facilitating experimental evolution, and the production of dormant eggs enabling the analysis of long-term response to environmental change through a resurrection ecology approach. We quantify plasticity in the expression of 35 candidate genes in D. magna populations resurrected from a lake that experienced changes in average temperature over the past century and from experimental populations differing in thermal tolerance isolated from a selection experiment. By measuring expression in multiple genotypes from each of these populations in control and heat treatments, we assess plastic responses to extreme temperature events. By measuring evolutionary changes in gene expression between warm- and cold-adapted populations, we assess evolutionary response to temperature changes. Evolutionary response to temperature increase is also assessed via an outlier analysis using EST-linked microsatellite loci. This study provides the first insights into the role of plasticity and genetic adaptation in orchestrating adaptive responses to environmental change in D. magna.

Bisulfite Sequencing with Daphnia Highlights a Role for Epigenetics in Regulating Stress Response to Microcystis through Preferential Differential Methylation of Serine and Threonine Amino Acids

Little is known about the influence that environmental stressors may have on genome-wide methylation patterns, and to what extent epigenetics may be involved in environmental stress response. Yet, studies of methylation patterns under stress could provide crucial insights on stress response and toxicity pathways. Here, we focus on genome-wide methylation patterns in the microcrustacean Daphnia magna, a model organism in ecotoxicology and risk assessment, exposed to the toxic cyanobacterium Microcystis aeruginosa. Bisulfite sequencing of exposed and control animals highlighted differential methylation patterns in Daphnia upon exposure to Microcystis primarily in exonic regions. These patterns are enriched for serine/threonine amino acid codons and genes related to protein synthesis, transport and degradation. Furthermore, we observed that genes with differential methylation corresponded well with genes susceptible to alternative splicing in response to Microcystis stress. Overall, our results suggest a complex mechanistic response in Daphnia characterized by interactions between DNA methylation and gene regulation mechanisms. These results underscore that DNA methylation is modulated by environmental stress and can also be an integral part of the toxicity response in our study species.

Anti-Inflammatory Activity of Cyanobacterial Serine Protease Inhibitors Aeruginosin 828A and Cyanopeptolin 1020 in Human Hepatoma Cell Line Huh7 and Effects in Zebrafish (Danio rerio)

Intensive growth of cyanobacteria in freshwater promoted by eutrophication can lead to release of toxic secondary metabolites that may harm aquatic organisms and humans. The serine protease inhibitor aeruginosin 828A was isolated from a microcystin-deficient Planktothrix strain. We assessed potential molecular effects of aeruginosin 828A in comparison to another cyanobacterial serine protease inhibitor, cyanopeptolin 1020, in human hepatoma cell line Huh7, in zebrafish embryos and liver organ cultures. Aeruginosin 828A and cyanopeptolin 1020 promoted anti-inflammatory activity, as indicated by transcriptional down-regulation of interleukin 8 and tumor necrosis factor α in stimulated cells at concentrations of 50 and 100 µmol·L⁻¹ aeruginosin 828A, and 100 µmol·L⁻¹ cyanopeptolin 1020. Aeruginosin 828A induced the expression of CYP1A in Huh7 cells but did not affect enzyme activity. Furthermore, hatched zebrafish embryos and zebrafish liver organ cultures were exposed to aeruginosin 828A. The transcriptional responses were compared to those of cyanopeptolin 1020 and microcystin-LR. Aeruginosin 828A had only minimal effects on endoplasmic reticulum stress. In comparison to cyanopeptolin 1020 our data indicate that transcriptional effects of aeruginosin 828A in zebrafish are very minor. The data further demonstrate that pathways that are influenced by microcystin-LR are not affected by aeruginosin 828A.

Changes in iTRAQ-Based Proteomic Profiling of the Cladoceran Daphnia magna Exposed to Microcystin-Producing and Microcystin-Free Microcystis aeruginosa

Global warming and increased nutrient fluxes cause cyanobacterial blooms in freshwater ecosystems. These phenomena have increased the concern for human health and ecosystem services. The mass occurrences of toxic cyanobacteria strongly affect freshwater zooplankton communities, especially the unselective filter feeder Daphnia. However, the molecular mechanisms of cyanobacterial toxicity remain poorly understood. This study is the first to combine the established body growth rate (BGR), which is an indicator of life-history fitness, with differential peptide labeling (iTRAQ)-based proteomics in Daphnia magna influenced by microcystin-producing (MP) and microcystin-free (MF) Microcystis aeruginosa. A significant decrease in BGR was detected when D. magna was exposed to MP or MF M. aeruginosa. Conducting iTRAQ proteomic analyses, we successfully identified and quantified 211 proteins with significant changes in expression. A cluster of orthologous groups revealed that M. aeruginosa-affected differential proteins were strongly associated with lipid, carbohydrate, amino acid, and energy metabolism. These parameters could potentially explain the reduced fitness based on the cost of the substance metabolism.

Cyanobacteria Affect Fitness and Genetic Structure of Experimental Daphnia Populations

Zooplankton communities can be strongly affected by cyanobacterial blooms, especially species of genus Daphnia, which are key-species in lake ecosystems. Here, we explored the effect of microcystin/nonmicrocystin (MC/non-MC) producing cyanobacteria in the diet of experimental Daphnia galeata populations composed of eight genotypes. We used D. galeata clones hatched from ephippia 10 to 60 years old, which were first tested in monocultures, and then exposed for 10 weeks as mixed populations to three food treatments consisting of green algae combined with cyanobacteria able/unable of producing MC. We measured the expression of nine genes potentially involved in Daphnia acclimation to cyanobacteria: six protease genes, one ubiquitin-conjugating enzyme gene, and two rRNA genes, and then we tracked the dynamics of the genotypes in mixed populations. The expression pattern of one protease and the ubiquitin-conjugating enzyme genes was positively correlated with the increased fitness of competing clones in the presence of cyanobacteria, suggesting physiological plasticity. The genotype dynamics in mixed populations was only partially related to the growth rates of clones in monocultures and varied strongly with the food. Our results revealed strong intraspecific differences in the tolerance of D. galeata clones to MC/non-MC-producing cyanobacteria in their diet, suggesting microevolutionary effects.

The interaction between cyanobacteria and zooplankton in a more eutrophic world

As blooms of cyanobacteria expand and intensify in freshwater systems globally, there is increasing interest in their ecological effects. In addition to being public health hazards, cyanobacteria have long been considered a poor quality food for key zooplankton grazers that link phytoplankton to higher trophic levels. While past laboratory studies have found negative effects of nutritional constraints and defensive traits (i.e., toxicity and colonial or filamentous morphology) on the fitness of large generalist grazers (i.e., Daphnia), cyanobacterial blooms often co-exist with high biomass of small-bodied zooplankton in nature. Indeed, recent studies highlight the remarkable diversity and flexibility in zooplankton responses to cyanobacterial prey. Reviewed here are results from a wide range of laboratory and field experiments examining the interaction of cyanobacteria and a diverse zooplankton taxa including cladocerans, copepods, and heterotrophic protists from temperate to tropical freshwater systems. This synthesis shows that longer exposure to cyanobacteria can shift zooplankton communities toward better-adapted species, select for more tolerant genotypes within a species, and induce traits within the lifetime of individual zooplankton. In turn, the function of bloom-dominated plankton ecosystems, the coupling between primary producers and grazers, the stability of blooms, and the potential to use top down biomanipulation for controlling cyanobacteria depend largely on the species, abundance, and traits of interacting cyanobacteria and zooplankton. Understanding the drivers and consequences of zooplankton traits, such as physiological detoxification and selective vs. generalist grazing behavior, are therefore of major importance for future studies. Ultimately, co-evolutionary dynamics between cyanobacteria and their grazers may emerge as a critical regulator of blooms.

Conserved transcriptional responses to cyanobacterial stressors are mediated by alternate regulation of paralogous genes in Daphnia

Despite a significant increase in genomic data, our knowledge of gene functions and their transcriptional responses to environmental stimuli remains limited. Here, we use the model keystone species Daphnia pulex to study environmental responses of genes in the context of their gene family history to better understand the relationship between genome structure and gene function in response to environmental stimuli. Daphnia were exposed to five different treatments, each consisting of a diet supplemented with one of five cyanobacterial species, and a control treatment consisting of a diet of only green algae. Differential gene expression profiles of Daphnia exposed to each of these five cyanobacterial species showed that genes with known functions are more likely to be shared by different expression profiles, whereas genes specific to the lineage of Daphnia are more likely to be unique to a given expression profile. Furthermore, while only a small number of nonlineage-specific genes were conserved across treatment type, there was a high degree of overlap in expression profiles at the functional level. The conservation of functional responses across the different cyanobacterial treatments can be attributed to the treatment-specific expression of different paralogous genes within the same gene family. Comparison with available gene expression data in the literature suggests differences in nutritional composition in diets with cyanobacterial species compared to diets of green algae as a primary driver for cyanobacterial effects on Daphnia. We conclude that conserved functional responses in Daphnia across different cyanobacterial treatments are mediated through alternate regulation of paralogous gene families.

Physiological interaction of Daphnia and Microcystis with regard to cyanobacterial secondary metabolites

Cyanobacterial blooms in freshwater ecosystems are a matter of high concern with respect to human health and ecosystem services. Investigations on the role of cyanobacterial secondary metabolites have largely been confined to microcystins, although cyanobacteria produce a huge variety of toxic or inhibitory secondary metabolites. Mass occurrences of toxic cyanobacteria strongly impact freshwater zooplankton communities; especially the unselective filter feeder Daphnia. Daphnids have been shown to successfully suppress bloom formation. However, the opposite situation, i.e. the suppression of Daphnia populations by cyanobacteria can be observed as well. To understand these contradictory findings the elucidation of the underlying physiological mechanisms that help daphnids to cope with cyanotoxins is crucial. We fed Daphnia magna with the cyanobacterium Microcystis aeruginosa PCC7806 for 24h and used high-resolution LCMS analytics to analyze the Microcystis cells, the Daphnia tissue and the surrounding medium in order to investigate the fate of seven investigated cyanobacterial compounds (cyanopeptolins A-C, microcyclamide 7806A and aerucyclamides B-D). For none of these bioactive compounds evidence for biotransformation or biodegradation by Daphnia were found. Instead feeding and subsequent release experiments point at the importance of transport mechanisms in Daphnia with regard to the cyanopeptolins A and C and microcyclamide 7806A. In addition we found hints for new inducible defense mechanism in Microcystis against predation by Daphnia. These putative defense mechanisms include the elevated production of toxic compounds other than microcystins, as could be demonstrated here for aerucyclamide B and D, cyanopoeptolin B and microcyclamide 7806A. Moreover, our data demonstrate the elevated active export of at least one cyanobacterial compound (microcyclamide 7806A) into the surrounding medium as a response to grazer presence, which might constitute an entirely new not yet described cyanobacterial defense mechanism.

Trypsin isozymes in the lobster Panulirus argus (Latreille, 1804): from molecules to physiology

Trypsin enzymes have been studied in a wide variety of animal taxa due to their central role in protein digestion as well as in other important physiological and biotechnological processes. Crustacean trypsins exhibit a high number of isoforms. However, while differences in properties of isoenzymes are known to play important roles in regulating different physiological processes, there is little information on this aspect for decapod trypsins. The aim of this review is to integrate recent findings at the molecular level on trypsin enzymes of the spiny lobster Panulirus argus, into higher levels of organization (biochemical, organism) and to interpret those findings in relation to the feeding ecology of these crustaceans. Trypsin in lobster is a polymorphic enzyme, showing isoforms that differ in their biochemical features and catalytic efficiencies. Molecular studies suggest that polymorphism in lobster trypsins may be non-neutral. Trypsin isoenzymes are differentially regulated by dietary proteins, and it seems that some isoenzymes have undergone adaptive evolution coupled with a divergence in expression rate to increase fitness. This review highlights important but poorly studied issues in crustaceans in general, such as the relation among trypsin polymorphism, phenotypic (digestive) flexibility, digestion efficiency, and feeding ecology.

Proteomic analysis of Daphnia magna hints at molecular pathways involved in defensive plastic responses

Background

Phenotypic plasticity in defensive traits occurs in many species when facing heterogeneous predator regimes. The waterflea Daphnia is well-known for showing a variety of these so called inducible defences. However, molecular mechanisms underlying this plasticity are poorly understood so far. We performed proteomic analysis on Daphnia magna exposed to chemical cues of the predator Triops cancriformis. D. magna develops an array of morphological changes in the presence of Triops including changes of carapace morphology and cuticle hardening.

Results

Using the 2D-DIGE technique, 1500 protein spots could be matched and quantified. We discovered 179 protein spots with altered intensity when comparing Triops exposed animals to a control group, and 69 spots were identified using nano-LC MS/MS. Kairomone exposure increased the intensity of spots containing muscle proteins, cuticle proteins and chitin-modifying enzymes as well as enzymes of carbohydrate and energy metabolism. The yolk precursor protein vitellogenin decreased in abundance in 41 of 43 spots.

Conclusion

Identified proteins may be either directly involved in carapace stability or reflect changes in energy demand and allocation costs in animals exposed to predator kairomones. Our results present promising candidate proteins involved in the expression of inducible defences in Daphnia and enable further in depth analysis of this phenomenon.

Genome-wide transcription profiles reveal genotype-dependent responses of biological pathways and gene-families in Daphnia exposed to single and mixed stressors

The present study investigated the possibilities and limitations of implementing a genome-wide transcription-based approach that takes into account genetic and environmental variation to better understand the response of natural populations to stressors. When exposing two different Daphnia pulex genotypes (a cadmium-sensitive and a cadmium-tolerant one) to cadmium, the toxic cyanobacteria Microcystis aeruginosa, and their mixture, we found that observations at the transcriptomic level do not always explain observations at a higher level (growth, reproduction). For example, although cadmium elicited an adverse effect at the organismal level, almost no genes were differentially expressed after cadmium exposure. In addition, we identified oxidative stress and polyunsaturated fatty acid metabolism-related pathways, as well as trypsin and neurexin IV gene-families as candidates for the underlying causes of genotypic differences in tolerance to Microcystis. Furthermore, the whole-genome transcriptomic data of a stressor mixture allowed a better understanding of mixture responses by evaluating interactions between two stressors at the gene-expression level against the independent action baseline model. This approach has indicated that ubiquinone pathway and the MAPK serine-threonine protein kinase and collagens gene-families were enriched with genes showing an interactive effect in expression response to exposure to the mixture of the stressors, while transcription and translation-related pathways and gene-families were mostly related with genotypic differences in interactive responses to this mixture. Collectively, our results indicate that the methods we employed may improve further characterization of the possibilities and limitations of transcriptomics approaches in the adverse outcome pathway framework and in predictions of multistressor effects on natural populations.

An investigation of the inter-clonal variation of the interactive effects of cadmium and Microcystis aeruginosa on the reproductive performance of Daphnia magna

Interactive effects between chemical and natural stressors as well as genetically determined variation in stress tolerance among individuals may complicate risk assessment and management of chemical pollutants in natural ecosystems. Although genetic variation in tolerance to single stressors has been described extensively, genetic variation in interactive effects between two stressors has only rarely been investigated. Here, we examined the interactive effects between a chemical stressor (Cd) and a natural stressor (the cyanobacteria Microcystis aeruginosa) on the reproduction of Daphnia magna in 20 genetically different clones using a full-factorial experimental design and with the independent action model of joint stressor action as the reference theoretical framework. Across all clones, the reduction of 21-day reproduction compared to the control treatment (no Cd, no M. aeruginosa) ranged from -10% to 98% following Cd exposure alone, from 44% to 89% for Microcystis exposure alone, and from 61% to 98% after exposure to Cd+Microcystis combined. Three-way ANOVA on log-transformed reproduction data of all clones together did not detect a statistically significant Cd×Microcystis interaction term (F-test, p=0.11), meaning that on average both stressors do not interact in inhibiting reproductive performance of D. magna. This finding contrasted expectations based on some known shared mechanisms of toxicity of Cd and Microcystis and therefore cautions against making predictions of interactive chemical+natural stressor effects from incomplete knowledge on affected biological processes and pathways. Further, still based on three-way ANOVA, we did not find statistically significant clone×Cd×Microcystis interaction when data for all clones were analyzed together (F-test, p=0.07), suggesting no inter-clonal variation of the interactive effect between Cd and Microcystis. However, when the same data were quantitatively analyzed on a clone-by-clone scale, we found a relatively wide range of deviations between observed and IA-model-predicted reproduction in combined Cd+Microcystis treatments (both in direction and magnitude), suggesting some biological significance of inter-clonal variation of interactive effects. In one of the twenty clones this deviation was statistically significant (two-way ANOVA, F-test, p=0.005), indicating an interactive Cd×Microcystis effect in this clone. Together, these two observations caution against the extrapolation of conclusions about mixed stressor data obtained with single clones to the level of the entire species and to the level of natural, genetically diverse populations.

Clonal variation in growth plasticity within a Bosmina longirostris population: the potential for resistance to toxic cyanobacteria

Many aquatic organisms respond phenotypically, through morphological, behavioral, and physiological plasticity, to environmental changes. The small-size cladoceran Bosminalongirostris, a dominant zooplankter in eutrophic waters, displayed reduced growth rates in response to the presence of a toxic cyanobacterium, Microcystisaeruginosa, in their diets. The magnitude of growth reduction differed among 15 clones recently isolated from a single population. A significant interaction between clone and food type indicated a genetic basis for the difference in growth plasticity. The variation in phenotypic plasticity was visualized by plotting reaction norms with two diets. The resistance of each clone to dietary cyanobacteria was measured as the relative change in growth rates on the "poor" diet compared with the "good" diet. The enhanced resistance to M. aeruginosa in B. longirostris was derived from both the reduced slope of reaction norms and the increased mean growth rates with two diets. The large clonal variation within a B. longirostris population may contribute to local adaptation to toxic cyanobacteria and influence ecosystem function via clonal succession.

Interspecific differences between D. pulex and D. magna in tolerance to cyanobacteria with protease inhibitors

It is known that cyanobacteria negatively affect herbivores due to their production of toxins such as protease inhibitors. In the present study we investigated potential interspecific differences between two major herbivores, Daphnia magna and Daphnia pulex, in terms of their tolerance to cyanobacteria with protease inhibitors. Seven clones each of D. magna and of D. pulex were isolated from different habitats in Europe and North America. To test for interspecific differences in the daphnids’ tolerance to cyanobacteria, their somatic and population growth rates were determined for each D. magna and D. pulex clone after exposure to varying concentrations of two Microcystis aeruginosa strains. The M. aeruginosa strains NIVA and PCC⁻ contained either chymotrypsin or trypsin inhibitors, but no microcystins. Mean somatic and population growth rates on a diet with 20% NIVA were significantly more reduced in D. pulex than in D. magna. On a diet with 10% PCC⁻, the population growth of D. pulex was significantly more reduced than that of D. magna. This indicates that D. magna is more tolerant to cyanobacteria with protease inhibitors than D. pulex. The reduction of growth rates was possibly caused by an interference of cyanobacterial inhibitors with proteases in the gut of Daphnia, as many other conceivable factors, which might have been able to explain the reduced growth, could be excluded as causal factors. Protease assays revealed that the sensitivities of chymotrypsins and trypsins to cyanobacterial protease inhibitors did not differ between D. magna and D. pulex. However, D. magna exhibited a 2.3-fold higher specific chymotrypsin activity than D. pulex, which explains the observed higher tolerance to cyanobacterial protease inhibitors of D. magna. The present study suggests that D. magna may control the development of cyanobacterial blooms more efficiently than D. pulex due to differences in their tolerance to cyanobacteria with protease inhibitors.

Effect of nutrient limitation of cyanobacteria on protease inhibitor production and fitness of Daphnia magna

Herbivore-plant interactions have been well studied in both terrestrial and aquatic ecosystems as they are crucial for the trophic transfer of energy and matter. In nutrient-rich freshwater ecosystems, the interaction between primary producers and herbivores is to a large extent represented by Daphnia and cyanobacteria. The occurrence of cyanobacterial blooms in lakes and ponds has, at least partly, been attributed to cyanotoxins, which negatively affect the major grazer of planktonic cyanobacteria, i.e. Daphnia. Among these cyanotoxins are the wide-spread protease inhibitors. These inhibitors have been shown (both in vitro and in situ) to inhibit the most important group of digestive proteases in the gut of Daphnia, i.e. trypsins and chymotrypsins, and to reduce Daphnia growth. In this study we grew cultures of the cyanobacterium Microcystis aeruginosa strain BM25 on nutrient replete, N-depleted or P-depleted medium. We identified three different micropeptins to be the cause for the inhibitory activity of BM25 against chymotrypsins. The micropeptin content depended on nutrient availability: Whereas N-limitation led to a lower concentration of micropeptins per biomass, P-limitation resulted in a higher production of these chymotrypsin inhibitors. The altered micropeptin content of BM25 was accompanied by changed effects on the fitness of Daphnia: A higher content of micropeptins led to lower IC50 values for Daphnia gut proteases and vice-versa. Following expectations, the lower content of micropeptins in the N-depleted BM25 caused higher somatic growth of Daphnia. Therefore, protease inhibitors can be regarded as a nutrient-dependent defence against grazers. Interestingly, although the P-limitation of the cyanobacterium led to a higher micropeptin content, high growth of D. magna was observed when they were fed with P-depleted BM25. This might be due to reduced digestibility of P-depleted cells with putatively thick mucilaginous sheaths. These findings indicate that both the grazer and the cyanobacterium benefit from P-reduction in light of digestibility and growth inhibition, which is an interesting starting point for further studies.