Target gene approaches: Gene expression in Daphnia magna exposed to predator-borne kairomones or to microcystin-producing and microcystin-free Microcystis aeruginosa
© Schwarzenberger et al; licensee BioMed Central Ltd. 2009
Received: 26 March 2008
Accepted: 16 November 2009
Published: 16 November 2009
Two major biological stressors of freshwater zooplankton of the genus Daphnia are predation and fluctuations in food quality. Here we use kairomones released from a planktivorous fish (Leucaspius delineatus) and from an invertebrate predator (larvae of Chaoborus flavicans) to simulate predation pressure; a microcystin-producing culture of the cyanobacterium Microcystis aeruginosa and a microcystin-deficient mutant are used to investigate effects of low food quality. Real-time quantitative polymerase chain reaction (QPCR) allows quantification of the impact of biotic stressors on differential gene activity. The draft genome sequence for Daphnia pulex facilitates the use of candidate genes by precisely identifying orthologs to functionally characterized genes in other model species. This information is obtained by constructing phylogenetic trees of candidate genes with the knowledge that the Daphnia genome is composed of many expanded gene families.
We evaluated seven candidate reference genes for QPCR in Daphnia magna after exposure to kairomones. As a robust approach, a combination normalisation factor (NF) was calculated based on the geometric mean of three of these seven reference genes: glyceraldehyde-3-phosphate dehydrogenase, TATA-box binding protein and succinate dehydrogenase. Using this NF, expression of the target genes actin and alpha-tubulin were revealed to be unchanged in the presence of the tested kairomones. The presence of fish kairomone up-regulated one gene (cyclophilin) involved in the folding of proteins, whereas Chaoborus kairomone down-regulated the same gene.
We evaluated the same set of candidate reference genes for QPCR in Daphnia magna after exposure to a microcystin-producing and a microcystin-free strain of the cyanobacterium Microcystis aeruginosa. The NF was calculated based on the reference genes 18S ribosomal RNA, alpha-tubulin and TATA-box binding protein. We found glyceraldehyde-3-phosphate dehydrogenase and ubiquitin conjugating enzyme to be up-regulated in the presence of microcystins in the food of D. magna. These findings demonstrate that certain enzymes of glycolysis and protein catabolism are significantly upgregulated when daphnids ingest microcystins. Each differentially regulated gene is a member of an expanded gene family in the D. pulex genome. The cyclophilin, GapDH and UBC genes show moderately large sequence divergence from their closest paralogs. Yet actin and alpha-tubulin genes targeteted by our study have nearly identical paralogs at the amino acid level.
Gene expression analysis using a normalisation factor based on three reference genes showed that transcription levels of actin and alpha-tubulin were not substantially changed by predator-borne chemical cues from fishes or invertebrates, although changes in expression on the protein level were shown elsewhere. These changes in protein level could be caused by others than the investigated paralogs, showing the importance of the construction of phylogenetic trees for candidate gene approaches. However, fish kairomones caused an up-regulation, and Chaoborus kairomone caused a down-regulation of cyclophylin, which proved to be a potential target gene for further analysis of kairomone effects on the life history of daphnids. Changes in food quality required a different set of reference genes compared to the kairomone experiment. The presence of dietary microcystins led to an up-regulation of two genes involved in the basic metabolism of D. magna, i.e. glyceraldehyde-3-phosphate dehydrogenase and ubiquitin conjugating enzyme, which suggests that microcystins in cyanobacteria have more general effects on the metabolism of D. magna than previously thought. Phylogenetic trees resolving relationships among paralogs that share the same gene name are shown to be important for determining the identity of the candidate genes under investigation.
Notwithstanding other so called 'model organisms', of which whole genome sequences have been obtained, the ecology of the model organism Daphnia sp. is outstandingly well known. Therefore it is a challenge to investigate gene/environment interactions for major ecological interactions of this cladoceran. In most freshwater lakes and ponds, Daphnia sp. is the major consumer of algae and cyanobacteria and is also the most important prey for predatory invertebrates and planktivorous fishes. In line with its intermediate position in the food chain, Daphnia sp. populations can be controlled by predation (top-down) or resources (bottom-up). However, the impact of top-down and bottom-up factors on Daphnia sp. population dynamics shows a pronounced seasonality . Predation pressure is low in spring, but peaks with the appearance of young-of-the-year fishes and fourth-instar larvae of Chaoborus water midges in early summer and remains moderate until fall [2, 3]. Bottom-up factors become a major constraint on Daphnia sp. population growth, particularly in eutrophic lakes in the summer when mass developments of toxic cyanobacteria lead to a suppression of Daphnia sp. biomass [4, 5]. The low predictability of intensity and seasonality of both predation pressure and dominance of toxic cyanobacteria should lead to the evolution of plastic instead of fixed adaptations . Indeed, adaptive phenotypic plasticity in Daphnia magna to both toxic cyanobacteria and predator-borne chemical cues has been reported [7–14]. Daphnia magna has been shown to reduce size at first reproduction (SFR) in response to kairomones from fish whereas chemical cues from larvae of Chaoborus flavicans led to increased SFR; both responses have been proven to be adaptive as fishes and Chaoborus differ in size-selectivity of their prey [10–14]. A different kind of phenotypic plasticity constitutes the enhanced tolerance of D. magna against cyanobacterial toxins upon exposure to a toxin producing cyanobacterium. The recent release of the Daphnia pulex genome sequence (wFleaBase: http://wFleaBase.org, JGI Genome Portal: http://www.Jgi.doe.gov/Daphnia/) creates the opportunity to precisely identify candidate genes that differ in their expression in response to predator-borne chemical cues (i.e. kairomones) and to a toxic cyanobacterium as a first step to decipher the underlying molecular mechanisms of adaptive phenotypic plasticity in D. magna.
Toxin production is a characteristic feature of several strains of the bloom-forming freshwater cyanobacterium Microcystis aeruginosa. A large variety of cyclic heptapeptides, termed microcystins (which can become a health hazard to humans and livestock ), have been identified in M. aeruginosa. Wild type M. aeruginosa PCC7806 produces relatively large amounts of two microcystin variants (LR and RR); experiments with a microcystin-free mutant of this strain  have led to the conclusion that microcystins contribute to the daphnid poisoning by M. aeruginosa [17, 18].
The microcystins of M. aeruginosa PCC7806 are known to inhibit protein phosphatases 1 and 2A in warm-blooded animals and in Daphnia sp ., which suggests that the poisoning effect of microcystin-LR in Daphnia sp. is due to the inhibition of these two protein phosphatases. However, it remains entirely unclear which major physiological pathways in Daphnia sp. are affected by the binding of microcystin to protein phosphatases 1 and 2A.
For single clones of D. magna it has been shown that they develop tolerance against a microcystin producing cyanobacterium [20, 21]. Although the physiological mechanisms have not been elucidated, this increased tolerance can be transferred to the offspring via maternal effects . For a better understanding of the process of physiological adaptation of daphnids to toxic cyanobacteria, it is important to examine the genes that are differentially regulated in the presence of microcystins. To achieve this goal we quantified the expression of a set of genes involved in the basic metabolism in D. magna when cultured on a M. aeruginosa PCC7806 mutant in which the production of microcystins had been knocked out, or on the microcystin-synthesizing wild type strain. As a reference cyanobacterium, we used a strain of Synchecoccus elongatus which is non-toxic to daphnids  and easily ingested; a strain of the green alga Scenedesmus obliquus which is widely used as standard food for daphnids was used as reference for high quality food .
Much of the recent evidence for inducible defences in freshwater ecology which has contributed to the general understanding of predator-prey interactions has been derived from experimental studies with daphnids. Despite considerable progress in the understanding of inducible defences, the underlying plasticity has rarely been studied at a molecular level. Only recently has the response of D. magna to kairomones from fish and invertebrates been investigated on the protein level, and it has been shown that a clone of D. magna that was isolated from a habitat where it coexists with fishes and invertebrate predators responds with a decrease in the proteins actin and alpha-tubulin . We hypothesized that this decrease might be due to a change in transcription of the actin and alpha-tubulin genes and that these genes might function as targets for predator-borne chemical cues. We therefore exposed the same clone of D. magna as Pijanowska & Kloc (2004)  to similar kairomones and quantified the transcription levels of putative target genes.
Kairomone experiment: Effects on the life-history of Daphnia magna
Kairomone experiment: Identification of reference genes and normalisation factors
Relative gene expression of the reference genes in D. magna after normalisation in the fish-kairomone experiment.
1.0000 ± 0.0051
1.123 ± 0.0179
1.0000 ± 0.0078
1.145 ± 0.0073
1.0000 ± 0.0063
0.7778 ± 0.0164
Relative gene expression of the reference genes in D. magna after normalisation in the Chaoborus-kairomone experiment.
1.0000 ± 0.0024
1.4451 ± 0.0059
1.0000 ± 0.0035
0.8881 ± 0.0113
1.0000 ± 0.0018
0.7792 ± 0.0058
Kairomone experiment: Relative normalized expression of the target genes actin and alpha-tubulin
A different response was found in the mean relative expression of alpha-tubulin. There was a significant 1.7-fold increase between fish-kairomone treatment and control (Tukey's HSD post-hoc, p < 0.001 after one-way ANOVA: F9;0.00026 = 64420.31, p < 0.001; Fig. 3C), and also a significant 1.07-fold increase of the mean relative expression of alpha-tubulin between control and Chaoborus-kairomone treatment (Tukey's HSD post-hoc, p < 0.001 after one-way ANOVA: F9;0.00026 = 64420.31, p < 0.001; Fig. 3D). The standard deviation of the relative expression of alpha-tubulin across fish-kairomone treatments was 0.50 and was thus around four orders of magnitude higher than SD in the reference genes (Tab. 1). The standard deviation of the relative expression of alpha-tubulin across Chaoborus-kairomone treatments was 0.05 and had thus almost the same SD across treatments as the reference genes (Tab. 2).
Kairomone experiment: Relative normalized expression of non-reference genes
Relative gene expression of the non-reference genes in D. magna after normalisation in the fish-kairomone experiment.
1.0000 ± 0.0123
2.2449 ± 0.0278
1.0000 ± 0.0519
1.8555 ± 0.0873
1.0000 ± 0.0012
2.9216 ± 0.0020
1.0000 ± 0.0046
2.0671 ± 0.0251
Relative gene expression of the non-reference genes in D. magna after normalisation in the Chaoborus kairomone experiment.
1.0000 ± 0.0033
1.8183 ± 0.0069
1.0000 ± 0.0065
1.2141 ± 0.0013
1.0000 ± 0.0012
0.3877 ± 0.0022
1.0000 ± 0.0042
1.7394 ± 0.0042
Microcystin experiment: Effects on the growth of Daphnia magna
Four different treatments were analysed in the microcystin experiment. The green alga S. obliquus, a good food-alga, served as calibrator. To account for potential general cyanobacterial effects, S. elongatus was included in the experimental design. The aim was to be able to differentiate between the wild type and the mutant of M. aeruginosa PCC 7806.
Microcystin experiment: Normalisation factors and identification of reference genes
Relative gene expression of the reference genes in D. magna after normalisation in the microcystin experiment. D. magna was fed either the green alga S. obliquus or microcystin-free cyanobacteria (S. elongatus or the mutant of M. aeruginosa PCC 7806 [Mut]) or the microcystin-producing wild type of M. aeruginosa PCC 7806 (WT).
M. aeruginosa WT
M. aeruginosa Mut
1.0000 ± 0.0030
0.3937 ± 0.0021
0.4077 ± 0.0006
0.3028 ± 0.0013
1.0000 ± 0.0094
1.0186 ± 0.0056
1.3176 ± 0.0075
3.0758 ± 0.2210
1.0000 ± 0.0015
2.4935 ± 0.0020
1.8613 ± 0.0073
1.0736 ± 0.0002
Microcystin experiment: Relative expression of the target genes GapDH, SucDH and UBC following normalisation
Resolving gene identities by homology to the Daphnia pulex genome
The occurrence of cyclic heptapeptides, termed microcystins, is widespread in cyanobacteria and not restricted to the genus Microcystis; mass development of microcystin-producing cyanobacteria can constitute a high risk for intoxication of humans and livestock . Numerous studies have been carried out in order to determine the ecological significance of microcystin production in cyanobacteria. The availability of the mcy- mutant of PCC7806, which is genetically identical to the wild type except for its inability to synthesize microcystins , made it possible to more precisely analyse the role of microcystins in daphnid poisoning. Life-table experiments with the wild type and the mcy- mutant of PCC7806 with Daphnia galeata have shown that the wild type was poisonous to D. galeata, whereas the mutant strain had no lethal effects [17, 18]. These findings and similar results for D. magna  suggest that microcystins play a role in the defence of M. aeruginosa against zooplankton grazing.
Feeding on the cyanobacterium M. aeruginosa led to significantly reduced growth of D. magna compared to animals grown on the high quality food alga S. obliquus  or the non-toxic cyanobacterium S. elongatus . Although there was reduced growth in D. magna feeding on M. aeruginosa compared to the reference cyanobacterium, there was no difference between the wild type and the mcy- mutant treatment. However, a specific microcystin effect became evident on day five, when D. magna raised on the wild type strain died, whereas no mortality was observed in D. magna raised on the mutant strain.
In an in-vitro system, microcystin-LR has been shown to inhibit protein phosphatases 1 and 2A in crude extracts of Daphnia sp . However protein phosphatase 1 and 2A each comprise a family of protein serine/threonine phosphatases with a wide range of different specificities that are mediated by different interactors  and regulatory subunits [29, 30]. Hence it remains entirely unclear which specific physiological pathways in daphnids are affected by the binding of microcystin to protein phosphatases 1 and 2A.
Here for the first time in-situ effects of dietary microcystins on gene expression of daphnids were investigated. The experiments presented in this paper were designed to identify genes involved in the general metabolism in D. magna in which the expression level responds to the presence of microcystins. We therefore compared the effects of the microcystin-producing wild type M. aeruginosa PCC7806 and the mcy- mutant of this strain on the relative expression of genes involved in basic metabolism. We found substantial up-regulation of GapDH (Dappu-302823) and UBC (Dappu-120690) in response to the presence of microcystins in the food of D. magna, which demonstrates that certain enzymes of glycolysis and protein catabolism are significantly up-regulated when daphnids ingest microcystins. For the first time a specific gene regulation in response to dietary microcystins has been demonstrated in daphnids. This up-regulation might have enabled D. magna to avoid a microcystin-specific depression of growth until day four but could not prevent mortality on day five of the growth experiment.
Upon exposure to the microcystin-producing wild type of M. aeruginosa PCC7806, D. magna has been shown to develop a tolerance against this toxic strain within an individual's lifespan and to transfer this tolerance to the next generation through maternal effects, a fact that has been interpreted as an inducible defence against microcystin . It remains to be tested which role the observed up-regulation of GapDH and UBC plays in the inducible tolerance of D. magna to microcystins. Furthermore, clones of D. magna have been shown to differ in their tolerance to M. aeruginosa PCC7806 , which suggests a genetic basis for increased toxin tolerance. It remains to be investigated whether the up-regulation of GapDH and UBC contributes to the tolerance to M. aeruginosa PCC7806.
In addition to the microcystins in PCC7806 wild type, both the wild type and mcy- mutant PCC7806 produce other classes of secondary metabolites of unknown biological activity [32, 33]. D. magna feeding on either of these two strains revealed a substantial up-regulation of SucDH, and it remains to be seen which cyanobacterial compounds induce this up-regulation of a key enzyme of the tricarboxylic acid cycle. In order to account for possible general effects of cyanobacteria on expression of the investigated genes, we fed Synechococcus elongatus to D. magna. This cyanobacterium is easily ingested by daphnids and does not contain toxins or inhibitors . The effects of S. elongatus on GapDH, UBC and SucDH were negligible compared to the afore mentioned effects of M. aeruginosa, which indicates that the up-regulation of the tested loci of GapDH, UBC and SucDH in D. magna is a specific and not a general response to cyanobacterial secondary metabolites. It would be interesting to see, if this holds true for all different paralogs of the affected genes or if the up-regulation is restricted to specific clusters or single paralogs of these highly variable genes (Fig. 7E-F).
Predation is an important stressor in aquatic communities, and many studies using Daphnia sp. have contributed to an understanding of the adaptive value of inducible anti-predator defences in the genus Daphnia. Achieving a better understanding of the mechanisms and constraints of the evolution of inducible anti-predator defences requires more research on the mechanisms of inducible defences at the molecular level. Only recently has this field been started to be explored. Our work was stimulated by the paper of Pijanowska & Kloc, (2004) , who used a clone of D. magna which has been shown to be plastic with regard to life-history traits and behaviour [10, 11, 13, 14, 34] in response to kairomones from fish and Chaoborus. Pijanowska & Kloc (2004)  have shown a dramatic decrease of the proteins actin and alpha-tubulin in this clone of D. magna when it was exposed to kairomones from planktivorous fish or the larvae of Chaoborus water midges. These identical effects of vertebrate and invertebrate kairomones suggested that actin might play a major role in anti-predator responses in D. magna in general. Using the same clone of D. magna, we here demonstrate that an exposure to chemical cues from both invertebrate and vertebrate predators results in a change in actin expression. However, although significant, the 1.75-fold (fish) increase and 0.94-fold (invertebrate) decrease in actin expression was rather moderate and did not reflect the dramatic decrease of the protein actin reported by Pijanowska & Kloc . The same holds true for the weak although significant increase in the gene alpha-tubulin in the fish (1.71) and the Chaoborus treatments (1.07). Since we found two possible alpha-tubulin orthologous protein sequences in D. pulex, which were very similar to each other (Fig. 7C), we concluded that the effect on the expression holds true for all paralogs in their cluster. Therefore, the substantial decrease of actin and alpha-tubulin on the protein level reported by Pijanowska & Kloc  could be a posttranslational process, e.g. miRNA-mediated regulation or increased degradation, as has been suggested by the authors . We conclude that these loci of actin and alpha-tubulin are no strong target genes for anti-predator defences. However, construction of phylogenetic trees reveals very high variability between the different paralogs of actin and alpha-tubulin (Fig. 7A-C). It remains to be tested if the decrease of actin and alpha-tubulin on the protein level reported by Pijanowska & Kloc  is caused by another paralogous sequence sharing the same gene name.
Following normalisation to NF, it turned out that the expression of 28S, UBC, 18S and cyclophilin was affected by the type of kairomone. Genes involved in protein biosynthesis (18S, 28S) and protein catabolism (UBC) were up-regulated by kairomone. These effects were considerably stronger for fish kairomone. The expression of cyclophylin (Dappu-92663), a gene involved in protein folding, was up-regulated in the presence of kairomones from vertebrate and down-regulated by kairomones from invertebrate predators. The finding that the two kairomones differ in their effect on cyclophylin in D. magna is in accord with the observation that the life-history response of this clone of D. magna differs between kairomones released from fish or Chaoborus . Cyclophilin, could serve as a potential target gene for further analysis of kairomone effects on daphnids. It remains to be seen how cyclophilin is involved in mediating kairomone effects on life history of daphnids and if this is specific to the orthologous sequence and to related paralogous sequences of cyclophilin.
Our study is the first detailed study that investigates effects of kairomones from vertebrate and invertebrate predators and of microcystin on gene expression of genes involved in different basic metabolic processes in D. magna. Kairomones from both vertebrate and invertebrate predators led to the well-established adaptive shifts in SFR in D. magna giving evidence for biologically active incubation water from either predator. Similarly, evidence for specific effects of microcystin comes from the higher mortality of D. magna on the wild type strain than on the mutant of M. aeruginosa PCC 7806. Calculating a combination normalisation factor based on the geometric mean of three genes for the kairomone experiment and for the growth experiment, stressor-specific regulation of some of the genes involved in basic metabolism is demonstrated.
All target genes in Daphnia show a surprisingly high variability between paralogs. If such a high variability holds true for other genes in D. magna, this could hint at a highly plastic genome, which might be adaptive for an animal that living in a very complex aquatic environment and therefore has to maintain a high potential for adaptations.
Three (GapDH, TBP, SucDH) of the seven genes investigated (GapDH, TBP SucDH, 28S, UBC, 18S, cyclophilin) were found to be stable across the kairomone treatments and were used as reference genes for normalization. Although significant, no substantial kairomone-dependent regulation of actin and alpha-tubulin was found, indicating that the dramatic decrease of actin and alpha-tubulin at the protein level in response to kairomones (reported earlier) was not due to a regulation of the transcription of the actin and alpha-tubulin loci investigated here. Therefore these gene loci cannot serve as target genes in the analysis of kairomone effects on D. magna. If this holds true for the other paralogs sharing the same gene name remains to be tested. However, the expression of other genes involved in protein biosynthesis, protein catabolism and protein folding, especially the regulation of cyclophilin by kairomones, indicated major effects on protein folding. These genes have the potential to serve as target genes in further analysis of kairomone effects on the life history of daphnids.
Three genes (i.e. 18S, TBP, alpha-tubulin) proved to be stable across microcystin-free and microcystin-containing cyanobacterial food treatments and were used for normalization. Two of the candidate genes (UBC and GapDH) were shown to have toxin-specific regulation and were clearly up-regulated when microcystins were present in the food. This indicates that microcystins strongly affect protein catabolism and glycolysis in D. magna when the animals ingest microcystins via the natural route of exposure, i.e. uptake of microcystin-containing food items; it remains to be seen which role the observed up-regulation of GapDH and UBC plays in the inducible tolerance of D. magna to microcystins.
The construction of phylogenetic trees is an essential step in target gene analysis in Daphnia in order to account for the high variability between different paralogs. Phylogenetic trees of the different paralogs are indispensable for clustering the utilized loci with similar ones and to delineate them from others. Related paralogs might have a similar relevance within the genome. This approach is especially important in an organism like Daphnia sp. with p-distances showing a very high variability between different paralogs.
Test species and cultures
A clone of Daphnia magna originating from Lake Binnensee, Germany, inhabited by fishes and various invertebrates, was cultured at 20°C in membrane-filtered tap water (conductivity: 740 μS/l; pH 7.2; major ions: Ca2+ (110 mg/l) and HCO3- (270 mg/l)).
Fifteen animals per litre were kept under non-limiting food concentrations (2 mg Cpart/l) with Scenedesmus obliquus (SAG-276-3a) (Stammsammlung für Algen, Göttingen, Germany) as food alga. Only third clutch neonates which had been born within 12 h were used for the experiments.
The green alga Scenedesmus obliquus, the cyanobacteria Synechococcus elongatus (SAG 89.79) and Microcystis aeruginosa (UTEX LB 2063 and PCC 7806), and a genetically engineered microcystin synthetase knock out mutant of Microcystis aeruginosa (PCC 7806 mcy- ) were cultivated semi-continuously in cyanophyceen medium  at 20°C, with half of the medium exchanged weekly. The medium consisted of 0.6 mM CaCl2 × 2 H20, 0.8 mM NaNO3, 0.4 mM K2HPO4 × 3 H2O, 0.4 mM MgSO4 × 7 H2O, 0.01 mM NaFeEDTA, 0.8 mM KCl, 100 μM H3BO3 and 20 μM Na2MoO4 × 2 H2O and had a pH of 8.5. Cyanobacteria were cultivated with constant light at 95 μE, S. obliquus at 130 μE. Carbon concentrations of the autotrophic food suspensions were estimated from photometric light extinction (800 nm) and from carbon-extinction equations previously determined.
Fish-conditioned water was prepared by exposing four sunbleaks (Leucaspius delineatus) in 4 l of tap water at 20°C for 24 h without feeding during the whole experiment. The fishes were stopped being fed for 24 h prior to the experiment; the water thus did not contain any faeces.
The Chaoborus incubation water was prepared by exposing 60 fourth-instar larvae of Chaoborus flavicans (which had previously been allowed to feed on zooplankton for 4 h each day of the experiment) in 1 l of tap water at 15°C for 20 h.
Predator-conditioned water was filtered (GFF), and in case of fish-conditioned water diluted 1:4, prior to introducing D. magna. From a cohort of D. magna neonates that had been born within 12 h from synchronized mothers, five animals each were exposed in 250 ml glass containers to predator-conditioned or control water from birth until maturity. All treatments were run in triplicate. D. magna were fed the green alga S. obliquus (4 mg Cpart/l); the media were changed daily. At the day when the 1st clutch was visible size at first reproduction (SFR) was determined as the size of the egg-bearing D. magna. For each replicate a mean SFR was calculated, and these mean values were used to calculate the respective mean value and the variance for the treatment.
From a cohort of new born D. magna, 8-10 animals each were transferred to 1 l of tap water with a food concentration of 2 mg Cpart/l. The animals were either fed the green alga S. obliquus as a control for high quality food or one of the three cyanobacteria. Each day the medium and the food were exchanged. The experiment took place under low light conditions at 20°C and lasted for four days for the real-time PCR analysis. All food treatments were run in triplicate, and somatic growth rates of D. magna were determined from dry weight of animals collected at the start and at day four of the experiment. according to .
RNA extraction and reverse transcription
At the end of the experiments the animals' RNA was extracted using the RNeasy Mini Kit (Qiagen). In order to remove any traces of genomic DNA, the RNA was treated with Desoxyribonuclease I (Fermentas) following the manufacturer's instructions. The integrity of the RNA was verified with 1.5% agarose gel electrophoreses. RNA concentrations were determined with a Qubit Fluorometer (Invitrogen). 1 μg of RNA was reverse transcribed using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). The cDNA was diluted 50-fold resulting in total RNA concentrations of 1 ng/μl. The cDNA was stored at -20°C.
Quantitative real-time PCR (QPCR)
Nine different housekeeping genes recently introduced for QPCR in D. magna by Heckmann et al. (2001)  were used in QPCR analysis: actin, alpha-tubulin, cyclophilin, glyceraldehyde-3-phosphate dehydrogenase (GapDH), succinate dehydrogenase (SucDH), TATA-box binding protein (TBP), ubiquitin conjugating enzyme (UBC), 18S ribosomal RNA (18S), 28S ribosomal RNA (28S).
QPCR was conducted on the 7300 real time PCR system (Applied Biosystems). Each reaction contained 5 μl of cDNA template, 10 μl Power SYBR® Green PCR Master Mix (Applied Biosystems) and 2.5 μM of each primer in a final volume of 20 μl. Each reaction was conducted in triplicate. Cycling parameters were 95°C for 10 min to activate the DNA polymerase followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. After the actual analysis, dissociation curves were performed to verify that no primer-dimers had been amplified. Outliers and samples diverging from the dissociation curve were omitted.
Data analysis and statistics
The raw data were analysed after QPCR. Because of the differing amplification efficiencies of the primer pairs , the relative expressions were calculated as quantities using the formula WR = (E+1)(CTmin-CTx) (a modification of the 2-ΔΔCT Method ), in which WR is a quantity for the relative expression of one sample, E is the amplification efficiency of its assay, CTmin is the lowest threshold cycle of all samples of this assay, and CTx is the threshold cycle of the analysed sample. The quantities could then be imported into geNorm version 3.4 , an Excel (Microsoft) based tool which calculates the minimum required number and best-suited combination from a given set of reference genes and from that generates a "normalisation factor" to be used for stable normalisation of QPCR measurements.
After determination of a normalisation factor from the set of reference genes, the raw data of the QPCR runs were imported into qBase version 1.3.5 , an Excel (Microsoft) based tool which calculates relative gene expression normalised using the normalisation factor generated with geNorm. D. magna fed entirely with S. obliquus served as calibrator which was always set as 1.
The growth rates and the body length were (×2)-1 transformed when needed to ensure homogeneity of variances and analysed with ANOVA and Tukey's honestly significant difference (HSD) for post hoc comparisons to assess differences in relative expression.
The data generated with qBase were log-transformed when needed to ensure homogeneity of variances and analysed with ANOVA and Tukey's honestly significant difference (HSD) for post hoc comparisons to assess differences in relative expression.
The statistics were performed with Statistica 6.0.
Database search and construction of phylogenetic trees
To compare the Daphnia magna sequences with the Daphnia pulex database http://wfleabase.org and to discover homologs in the genome for the genes of interest, the D. magna sequences of the target genes actin, alpha-tubulin, cyclophilin, GapDH and UBC  were used as queries for sequence similarity searches using BLASTn against the Dappu v1.1 draft genome assembly in wFleabase. The best alignments with the highest score was taken as the ortholog for the sequence of interest in D. pulex. Its protein sequence was subsequently used to query for additional homologs using BLASTp (blastp; e-value cut off limit = 0.01) against the the v1.1 gene builds (July, 2007) archived in wFleaBase. All protein sequences with a significant hit were taken and aligned (BioEdit v.126.96.36.199 ) with the orthologous protein sequence. Using the program MEGA 4  Neighbor-Joining trees with Bootstrap Test of Phylogeny were constructed and p-distances were calculated. Hypothetically, the D. pulex trees and the D. magna trees are approximately identical.
particulate organic carbon
TATA-box binding protein
ubiquitin conjugating enzyme
18S ribosomal RNA
28S ribosomal RNA.
The authors would like to thank Prof. Dr. M. Deckert from the University Hospital Cologne, Department of Neuropathology, for kindly providing access to real time PCR.
The sequencing and portions of the analyses were performed at the DOE Joint Genome Institute under the auspices of the U.S. Department of Energy's Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48, Lawrence Berkeley National Laboratory under Contract No. DE-AC02-05CH11231, Los Alamos National Laboratory under Contract No. W-7405-ENG-36 and in collaboration with the Daphnia Genomics Consortium (DGC), http://daphnia.cgb.indiana.edu. Additional analyses were performed by wFleaBase, developed at the Genome Informatics Lab of Indiana University with support to Don Gilbert from the National Science Foundation and the National Institutes of Health. Coordination infrastructure for the DGC is provided by The Center for Genomics and Bioinformatics at Indiana University, which is supported in part by the METACyt Initiative of Indiana University, funded in part through a major grant from the Lilly Endowment, Inc. Our work benefits from and contributes to the Daphnia Genomics Consortium.
The authors thank Frederic Bartlett for the correction of the manuscript. Thanks to John Colbourne and Frank Nitsche for the assistance in constructing the phylogenetic trees.
- Sommer U, Gliwicz ZM, Lampert W, Duncan A: The PEG-model of seasonal succession of planktonic events in fresh waters. Arch Hydrobiol. 1986, 106: 433-471.Google Scholar
- Ringelberg J, Flik BJG, Lindenaar D, Royackers K: Diel vertical migration of Daphnia hyalina (sensu latiori) in Lake Maarsseveen: Part 1. Aspects of seasonal and daily timing. ArchHydrobiol. 1991, 121: 129-145.Google Scholar
- Voss S, Mumm H: Where to stay by night and day: Size-specific and seasonal differences in horizontal and vertical distribution of Chaoborus flavicans larvae. Freshw Biol. 1999, 42: 201-213. 10.1046/j.1365-2427.1999.444444.x.View ArticleGoogle Scholar
- de Bernardi R, Giussani G: Are blue-green algae a suitable food for zooplankton? An overview. Hydrobiologia. 1990, 200/201: 29-41. 10.1007/BF02530326.View ArticleGoogle Scholar
- Gilbert JJ: Differential effects of Anabaena affinis on cladocerans and rotifers: Mechanisms and implications. Ecology. 1990, 71: 1727-1740. 10.2307/1937581.View ArticleGoogle Scholar
- Miner BG, Sultan SE, Morgan SG, Padilla DK, Reylea RA: Ecological consequences of phenotypic plasticity. TREE. 2005, 20: 685-692.PubMedGoogle Scholar
- Gustafsson S, Hansson LA: Development of tolerance against toxic cyanobacteria in Daphnia. Aqu Ecol. 2004, 38: 37-44. 10.1023/B:AECO.0000020985.47348.5e.View ArticleGoogle Scholar
- Gustafsson S, Rengefors K, Hansson LA: Increased consumer fitness following transfer of toxin tolerance to offspring via maternal effects. Ecology. 2005, 86: 2561-2567. 10.1890/04-1710.View ArticleGoogle Scholar
- Hansson LA, Gustafsson S, Rengefors K, Bomark L: Cyanobacterial chemical warfare affects zooplankton community composition. Freshw Biol. 2007, 52: 1290-1301. 10.1111/j.1365-2427.2007.01765.x.View ArticleGoogle Scholar
- Weider LJ, Pijanowska J: Plasticity of Daphnia life histories in response to chemical cues from predators. Oikos. 1993, 67: 385-392. 10.2307/3545351.View ArticleGoogle Scholar
- Mikulski A, Czernik M, Pijanowska J: Induction time and reversibility of changes in Daphnia life history caused by the presence of fish. J Plankton Res. 2005, 27: 757-762. 10.1093/plankt/fbi048.View ArticleGoogle Scholar
- Loose CJ, Von Elert E, Dawidowicz P: Chemically induced diel vertical migration in Daphnia: a new bioassay for kairomones exuded by fish. Arch Hydrobiol. 1993, 126: 329-337.Google Scholar
- Von Elert E, Loose CJ: Predator-induced diel vertical migration in Daphnia - enrichment and preliminary chemical characterization of a kairomone exuded by fish. J Chem Ecol. 1996, 22: 885-895. 10.1007/BF02029942.View ArticlePubMedGoogle Scholar
- Von Elert E, Pohnert G: Predator specificity of kairomones in diel vertical migration of Daphnia: a chemical approach. Oikos. 2000, 88: 119-128. 10.1034/j.1600-0706.2000.880114.x.View ArticleGoogle Scholar
- Carmichael WW: The toxins of cyanobacteria. Scientific American. 1994, 270: 64-72.View ArticleGoogle Scholar
- Dittmann E, Neilan BA, Erhard M, Von Doehren H, Börner T: Insertional mutagenesis of a peptide synthetase gene that is responsible for hepatotoxin production in the cyanobacterium Microcystis aeruginosa PCC 7806. Mol Microbiol. 1997, 26: 779-787. 10.1046/j.1365-2958.1997.6131982.x.View ArticlePubMedGoogle Scholar
- Kaebernick M, Rohrlack T, Christoffersen K, Neilan BA: A spontaneous mutant of microcystin biosynthesis: genetic characterization and effect on Daphnia. Environ Microbiol. 2001, 3: 669-679. 10.1046/j.1462-2920.2001.00241.x.View ArticlePubMedGoogle Scholar
- Rohrlack T, Dittmann E, Henning M, Boerner T, Kohl J: Role of microcystins in poisoning and food ingestion inhibition of Daphnia galeata caused by the cyanobacterium Microcystis aeruginosa. Applied & Environmental Microbiology. 1999, 65: 737-739.Google Scholar
- DeMott WR, Dhawale S: Inhibition of in-vitro protein phosphatase-activity in three zooplankton species by microcystin-lr, a toxin from cyanobacteria. Archiv für Hydrobiologie. 1995, 134: 417-424.Google Scholar
- Gustafsson S, Rengefors K, Hansson LA: Increases consumer fitness following transfer of toxin tolerance to offspring via maternal effects. Ecology. 2005, 86: 2561-2567. 10.1890/04-1710.View ArticleGoogle Scholar
- Hansson LA, Gustafsson S, Rengefors K, Bomark L: Cyanoabacterial chemical warfare affects zooplankton community composition. Freshwat Biol. 2007, 52: 1290-1301. 10.1111/j.1365-2427.2007.01765.x.View ArticleGoogle Scholar
- Von Elert E, Wolffrom T: Supplementation of cyanobacterial food with polyunsaturated fatty acids does not improve growth of Daphnia. Limnology and Oceanography. 2001, 46: 1552-1558.View ArticleGoogle Scholar
- Von Elert E: Determination of limiting polyunsaturated fatty acids in Daphnia galeata using a new method to enrich food algae with single fatty acids. Limnology and Oceanography. 2002, 47: 1764-1773.View ArticleGoogle Scholar
- Pijanowska J, Kloc M: Daphnia response to predation threat involves heat-shock proteins and the actin and tubulin cytoskeleton. Genesis: the Journal of Genetics & Development. 2004, 38: 81-86.View ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology. 2002, 3: research0034.I-0034.II. 10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- Heckmann LH, Connon R, Hutchinson TH, Maund SJ, Sibly RM, Callaghan A: Expression of target and reference genes in Daphnia magna exposed to ibuprofen. BMC Genomics. 2006, 7: 175-182. 10.1186/1471-2164-7-175.PubMed CentralView ArticlePubMedGoogle Scholar
- Chorus I: Cyanotoxins: Occurrence, causes, consequences. 2001, Berlin, Heidelberg: SpringerView ArticleGoogle Scholar
- Martin-Creuzburg D, Von Elert E, Hoffmann KH: Nutritional constraints at the cyanobacteria--Daphnia interface: the role of essential lipids. Limnology and Oceanography. 2008, 53: 456-468.View ArticleGoogle Scholar
- Ceulemans H, Bollen M: Functional diversity of protein phosphatase-1, a cellular economizer and reset button. Physiological Reviews. 2004, 84: 1-39. 10.1152/physrev.00013.2003.View ArticlePubMedGoogle Scholar
- Janssens V, Goris J: Protein phosphatase 2A: a highly regulated family of serine/threonine phosphatases implicated in cell growth and signalling. Biochemical Journal. 2001, 353: 417-439. 10.1042/0264-6021:3530417.PubMed CentralView ArticlePubMedGoogle Scholar
- Gustafsson S: Inter- and intra-clonal variability in Daphnia in response to toxic cyanobacteria. subm. 2008, 49-57.Google Scholar
- Moore RE: Cyclic-peptides and depsipeptides from cyanobacteria - a review. Journal of Industrial Microbiology. 1996, 16: 134-143. 10.1007/BF01570074.View ArticlePubMedGoogle Scholar
- Rao PVL, Gupta N, Bhaskar ASB, Jayaraj R: Toxins and bioactive compounds from cyanobacteria and their implications on human health. Journal of Environmental Biology. 2002, 23: 215-224.PubMedGoogle Scholar
- Loose CJ: Daphnia diel vertical migration behavior: Response to vertebrate predator abundance. Archiv für Hydrobiologie Beihefte. 1993, 39: 29-36.Google Scholar
- Jüttner F, Leonhardt J, Möhren S: Environmental factors affecting the formation of mesityloxid, dimethylallylic alcohol and other volatile compounds excreted by Anabaena cylindrica. J o General Microbiology. 1983, 129: 407-412.Google Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 408-10.1006/meth.2001.1262.View ArticleGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biology. 2007, 8: R19-10.1186/gb-2007-8-2-r19.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.