Expression of target and reference genes in Daphnia magna exposed to ibuprofen
© Heckmann et al; licensee BioMed Central Ltd. 2006
Received: 06 March 2006
Accepted: 07 July 2006
Published: 07 July 2006
Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e.g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application.
We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l-1. Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was up-regulated following 24 h exposure to 63–81 mg IB l-1.
As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.
Toxicogenomics is an evolving discipline investigating stressor impact at the genome level . Microarrays are state-of-the-art tools for global gene expression profiling of the response of an organism in a particular biological context. However, expression levels of key genes responding on the microarray need to be validated . The preferred technique of validation is by real-time quantitative polymerase chain reaction (QPCR) [3, 4]. Although QPCR is a fairly reliable technique, amplification can vary depending on factors such as RNA integrity, reverse transcriptase (RT) efficiencies, sample-to-sample variations in amplification efficiency, and variation in cDNA sample loading. Using the same sample sizes, assessing RNA integrity and equalising RNA concentrations prior to RT are some of the basic normalisation steps in QPCR . However, normalisation to some internal control is essential for accurate QPCR to balance sample-to-sample variations within the RT and PCR reactions. Internal control is generally achieved using reference genes, also known as housekeeping genes.
Some of the most commonly used reference genes, such as β-actin and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), are historical carryovers from northern blotting, a "predecessor" of QPCR . It is good practice to rigorously validate the suitability of these reference genes under the specific experimental condition to which they are applied [5–7]. Total RNA and rRNA have also been suggested as internal controls, although again differences between tissue, individuals and experimental conditions apply . Moreover, using RNA as an internal control has the disadvantage of not controlling for variation inherent in the RT and PCR reactions . None of these genes or RNA products are completely invariable, and thus each needs validation before use [5–7]. A further problem of using only one reference gene for normalisation, is that it theoretically should have been normalised itself prior to target gene normalisation . Various methods have been proposed to overcome these problems. One involves the use of a normalisation factor (NF) based on the geometric mean of multiple carefully selected reference genes. This NF can be calculated by geNorm (see Vandesompele et al., 2002  and the geNorm website ) or other freely available Excel (Microsoft) based software applications, e.g. BestKeeper  (see Pfaffl et al., 2004 ). Vandesompele et al., 2002  recommend that NF s should be estimated from a minimum of three reference genes. The advantage of this approach is that it allows for the cumulative error of the entire process from RNA extraction to the QPCR.
Candidate reference genes and target gene investigated
GenBank accession no.
wFleaBase EST no.
Forward primer (5'-3')
Reverse primer (5'-3')
Amplicon size (bp)
Structural constituent of cytoskeleton (F) GO:0005200
CCA CAC TGT CCC CAT TTA TGA A
CGC GAC CAG CCA AAT CC
Microtubule (C) GO:0005874
TGG AGG TGG TGA CGA CT
CCA AGT CGA CAA AGA CAG CA
Protein folding (P) GO:0006457
GAC TTT CCA CCA GTG CCA TT
AAC TTT CCA TCG CAT CAT CC
Glyceraldehyde-3- phosphate dehydrogenase
Glycolysis (P) GO:0006096
GGC AAG CTA GTT GTC AAT GG
TAT TCA GCT CCA GCA GTT CC
Tricarboxylic acid cycle (P) GO:0006099
TGC CAT TTA GTC GCA CTC AG
GTG AGC TTG TCC TCC TTT GC
TATA-box binding protein
Transcription initiation (P) GO:0006367
GCA GGG AAG TTT AGT TTC TGG A
TGG TAT GCA CAG GAG CAA AG
Tryptophanyl- tRNA synthetase
Tryptophanyl-tRNA ligase activity (F) GO:0004830
GCA CTG TGA TCA GCT CTT GC
TTT CGA TTG GAG GAA ATT CG
Ubiquitin Conjugating Enzyme
Ubiquitin-dependent protein catabolism (P) GO:0006511
TCA CCT GCA CTC ACC ATT TC
AAT CTC CGG AAC CAA AGG AT
18S ribosomal RNA
Protein biosynthesis (P) GO:0006412
CGC TCT GAA TCA AGG GTG TT
TGT CCG ACC GTG AAG AGA GT
28S ribosomal RNA
Protein biosynthesis (P) GO:0006412
GAG GCG CAA TGA AAG TGA AG
TGT TCG AGA CGG GAT CA
Leukotriene B4 12-hydroxy- dehydrogenase b
Response to toxin (P) GO:0009636
AAC CTA CAC CGA GGG TTT CG
TCC AAC ATT AAC GCC ATT AAG C
We use the geNorm algorithm  to estimate the variability of the reference genes, and to discover an optimal normalisation factor (NF), based on the geometric mean of three of them. geNorm estimates reference gene variability by calculating an expression level ratio for any two candidate reference genes. The variation between the expression ratios is then used as an inverse measure to estimate the variability of the analysed reference genes (please see Methods or Vandesompele et al., (2002)  for further information). The approach is illustrated using as target gene the D. magna ortholog of leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was identified as being up-regulated in a suppressive subtractive hybridisation following 24 h exposure to 63–81 mg IB l-1 . Ltb4dh is of special interest because NSAIDs are known to inhibit the mammalian biosynthesis of various eicosanoids (e.g. prostaglandins, epoxyeicosatrienoic acids and leukotriene) that play important regulatory and signalling functions, for instance regulation of ion flux [14, 15]. As anticipated, use of the NF clarified the response of Ltb4dh to ibuprofen.
Water chemistry and quantification of ibuprofen
The chemical parameters pH and conductivity remained stable throughout the exposure, ranging between 7.8–7.9 and 422–452 μS cm-1, respectively. Water temperatures ranged from 20.8 to 21.0°C. Ibuprofen-sodium increased measured pH by approximately 9% and conductivity by approximately 7% in the highest treatment (80 mg IB l-1) compared with the control (data not shown).
Quantification of IB revealed that nominal concentrations were within ± 10% of the measured concentration; except for 20 mg IB l-1 replicate three that was 10.6% lower than the nominal (data not shown). No mortality was observed in D. magna in any of the treatments following exposure.
Relative expression and variability of candidate reference genes
Relative gene expression and variability of candidate reference genes
20 mg IB l-1
40 mg IB l-1
80 mg IB l-1
Amplification efficiency (E)a
1.000 ± 0.078
0.963 ± 0.035
1.030 ± 0.044
1.001 ± 0.040
0.795 ± 0.066
1.000 ± 0.129
0.896 ± 0.126
0.885 ± 0.095
0.613 ± 0.045 c
1.010 ± 0.016
1.000 ± 0.040
0.949 ± 0.103
0.900 ± 0.015
0.817 ± 0.067 d
1.038 ± 0.016
1.000 ± 0.093
0.914 ± 0.059
0.919 ± 0.047
0.879 ± 0.067
0.946 ± 0.029
1.000 ± 0.145
0.889 ± 0.170
0.950 ± 0.180
1.027 ± 0.026
0.992 ± 0.015
1.000 ± 0.096
0.786 ± 0.176
0.788 ± 0.153
1.038 ± 0.084
0.925 ± 0.016
1.000 ± 0.117
0.954 ± 0.098
0.906 ± 0.060
0.955 ± 0.081
0.922 ± 0.064
1.000 ± 0.040
0.995 ± 0.095
1.024 ± 0.016
0.888 ± 0.097
1.034 ± 0.013
1.000 ± 0.279
2.096 ± 1.074
2.506 ± 1.640
2.204 ± 0.910
0.853 ± 0.089
1.000 ± 0.261
1.247 ± 0.445
1.410 ± 0.560
1.457 ± 0.616
0.846 ± 0.020
Variability of normalisation factors
Relative expression of reference genes following normalisation
There is a circular problem with normalising a target gene with reference genes that should be normalised themselves prior to target gene normalisation . Thus, in addition to normalising the target gene, we also normalised all 10 reference genes by NF3. Atb had the same significant (p < 0.05) differences as prior to normalisation, and should therefore be regarded as a target gene. None of the other reference genes, including Cyp, were significantly (p > 0.05) affected by IB following normalisation (data not shown).
Relative expression of the target gene Ltb4dh following normalisation
We investigated the variability of 10 candidate reference genes in D. magna following a 24 h exposure to IB to discover the least variable internal control(s) for QPCR normalisation. A comparison of the reference genes, using the geNorm software, ranked the least variable genes as; UBC = GAPDH <Act <WARS <SDH <TBP <18S <28S. Atb and Cyp were not included in the geNorm analysis as they were down-regulated by IB. Furthermore, geNorm identified the optimal normalisation factor as NF3 based on the geometric mean expression of UBC, GAPDH and Act. This NF was based on the lowest recommended number of reference genes with the lowest level of variation . The response to IB of an example target gene, Ltb4dh, was little changed by application of NF3 (Fig. 2), although at the highest concentration of IB, there was a slight increase in the response, and a diminution of the variation between replicates. Such reduction in the variation between replicates is exactly what one hopes to achieve through normalisation.
QPCR confirmed that Ltb4dh was up-regulated following exposure to IB . NSAIDs are known to inhibit the biosynthesis of various eicosanoids that play important regulatory and signalling functions, e.g. regulation of ion flux. In mammals, Ltb4dh is involved in the metabolism of leukotriene B4 an eicosanoid that is formed in the lipoxygenase pathway . There is evidence of eicosanoid pathways in invertebrates being similar to the mammalian pathways . Further experimentation, involving a global expression profile based on cDNA microarrays, is underway to support this hypothesis and reveal the overall molecular stress response of D. magna to IB. To our knowledge, this is one of the first genomic studies to look at the impact of a NSAID on an aquatic invertebrate. The only non-mammalian studies are on yeast Saccharomyces cerevisiae , a Drosophila cell line  and Bacillus megaterium .
Following normalisation to NF3 it transpired that Atb, assayed as a candidate reference gene, was affected by IB, and so should be regarded as an IB target gene. In a review on the molecular mechanism of NSAIDs in mammals, Tarnawski and Jones, 2003  reveal that it is likely that NSAIDs exert many inhibitory effects, one of which is the inhibition of cell division. The observed suppression of Atb following exposure to IB could be linked to this, since alpha-tubulin helps form microtubules that are involved in cell division, cell structure, and transportation of vesicles and organelles.
Of the remaining candidate reference genes, Act, GAPDH, UBC, WARS, SDH and TBP had a constant expression following exposure to IB. But although, GAPDH and UBC were the most stable reference genes they both had a weakly significant (p < 0.10) linear trend of down-regulation. Thus, relying solely on either of these two reference genes for normalisation could obscure the output.
Using reference genes as an internal control in QPCR is currently the recommended approach . However, Bustin and Nolan, 2004  advocate the introduction of more standard analysis and reporting procedures in QPCR similar to the MIAME guidelines  established for microarray technology. In union with other authors [5, 10, 21, 23] we would like to recommend some basic applications for QPCR: (i) apply a dissociation curve in every QPCR, and omit any samples diverging from the curve; (ii) use a minimum of two non-template controls (NTCs) in every QPCR, and omit any samples having a cycle threshold (Ct) value close to or higher than the mean NTC Ct-value; (iii) avoid inter-assay comparison if possible, as it is more appropriate to duplicate the samples in one QPCR than using triplicates of each sample, but having to run two QPCRs (S.N. Peirson, personal communication); (iv) use a NF based on a minimum of three validated reference genes for target gene normalisation ; (v) integrating the DART-PCR approach , which verifies sample-to-sample variations in amplification efficiency with a normalisation factor software package, e.g. geNorm , creates a powerful tool for analysing relative gene expression of target genes.
Investigators engaging in toxicogenomic research using invertebrates, or vertebrate cell lines, exposed to NSAIDs may be able to apply our findings directly under similar experimental conditions, e.g. acute exposure. However, we strongly recommend validating reference genes prior to commencing target gene expression profiling under experimental conditions different than those described herein. Our findings stress the importance of performing fundamental research, by validating the most invariant internal QPCR control(s) of a particular experimental condition, prior to investigating the expression of target genes examined by QPCR only, or as QPCR validation of microarray data.
Daphnia magna Straus (Clone Type 5 – IRCHA), originally obtained from the Water Research Centre (WRc), Medmenham, UK, were investigated. For further information on culturing conditions of D. magna please see Hooper et al., 2006 .
Experimental design and analytical chemistry
The experimental design consisted of three replicates of a control and three treatments with ibuprofen-sodium (Sigma-Aldrich, CAS no. 31121-93-4) containing 20, 40 and 80 mg IB l-1, respectively. Each replicate consisted of 100 third brood neonates (< 24 h old) placed in a 1000 ml glass beaker containing one litre of reconstituted freshwater with or without the addition of IB. Neonates were exposed to IB for 24 h under the same conditions as stock cultures (see Hooper et al., 2006 ), but without feeding. To quantify IB, 1.5 ml was sampled from each replicate at 0 h and 24 h and stored at -20°C. Subsequently, the free base concentration of IB was measured at 217 nm by UV-spectrophotometry using an Ultrospec3000 (Biochrom) after the method of Pascoe et al., 2003 . Conductivity, pH and water temperature were checked at the beginning and at the end of the exposure to verify stable water chemistry.
Tissue preparation, RNA extraction and reverse transcription
Following exposure, the neonates were immediately placed in 0.2 ml of RNAlater® (Ambion) and stored at -80°C. Subsequently, total RNA was extracted using the RNeasy Mini kit with on-column DNase treatment (Qiagen) to remove any traces of genomic DNA following the manufacturer's instructions. RNA concentrations were determined by spectrophotometry using GeneQuant Pro (Biochrom). The integrity of the RNA was verified by 1% agarose gel electrophoresis. Four micrograms of DNase-treated total RNA was reverse transcribed with Oligo(dT)12–18 primers (Invitrogen) using the Omniscript Reverse Transcriptase kit (Qiagen) following the manufacturer's instructions. Once synthesised, cDNA was diluted 10-fold resulting in total RNA concentrations of 10 ng/μl and stored at -20°C.
Real-time quantitative PCR
Ten expressed sequence tags (ESTs) were obtained from the NCBI website  and the Daphnia Genomics Consortium website  (see Colbourne et al., 2005 ). ESTs were verified by BLASTN and BLASTX analysis . Primers were designed using Primer3 , and synthesised by MWG. The investigated ESTs, accession numbers, primer sequences and gene ontology are shown in Table 1.
QPCR was conducted on the GeneAmp 5700 Sequence Detection System (SDS) (Applied Biosystems) using ABsolute™ QPCR SYBR® Green ROX (500 nM) mix (ABgene). Each reaction was run in triplicate and contained 2.5 μl of cDNA template (equivalent to 25 ng total RNA) along with 900 nM primers in a final volume of 25 μl. Cycling parameters were 95°C for 15 min to activate the DNA polymerase, then 40 cycles of 95°C for 15 s and 60°C for 1 min. Melting curves were performed by using dissociation curve SDS software version 1.3 (Applied Biosystems) to verify that only a single product with no primer-dimers were amplified.
Data analysis and statistics
Following QPCR, the raw data was exported to an Excel workbook (Microsoft), entitled Data Analysis for Real-Time PCR (DART-PCR) . This programmed workbook enables calculation of threshold cycles and amplification efficiencies for every sample. Differences in amplification efficiency were assessed using one-way analysis of variance (ANOVA), based upon the null hypotheses: (i) that amplification efficiency is comparable within sample groups (outlier detection) and (ii) that amplification efficiency is comparable between sample groups (amplification equivalence) . Outliers identified by DART-PCR and samples diverging from the dissociation curve were omitted prior to further analysis.
Differences in relative expression were assessed on log-transformed data using ANOVA and Tukey's honestly significant difference (HSD) for post hoc comparisons. Equality of variance was tested using Levene's test. The relative expression of the reference genes calculated by DART-PCR was exported as an Excel workbook (Microsoft) to geNorm version 3.4 . The "geNorm" algorithm is based on the proposition that, in any set of assays, the expression levels of suitable reference genes should be perfectly correlated. However, because of technical difficulties in assessing correlation in circumstances when there is little variation, a modified procedure is adopted as follows . First, the ratio of expression levels is calculated for any two candidate reference genes. Assuming that the expression ratios should not vary between assays, the variation between the expression ratios is used as an inverse measure, here referred to as gene variability, of the adequacy of the reference genes to be used for normalisation. Thus, the lower the variation, the better the reference genes. Variation is here measured as standard deviation. Given a set of assays that covers the treatments of interest, this procedure can be used to identify the best two reference genes from all possible pairs of reference genes available to the experimenter. A normalisation factor (NF2) is then produced for these two genes, calculated as the geometric mean of their expression levels. The procedure is then extended to assess whether stepwise inclusion of additional reference genes – producing NF3 to NF n – will reduce the average pairwise variation between them. The optimal number of reference genes to be included in a NF is estimated by comparing the pairwise variation between sequential normalisation factors , e.g. NF3 vs. NF4. All statistical tests were performed in SPSS 12.0.1 for Windows (SPSS, Chicago, IL, USA) applying a significant level of 5%, except reference gene variability and NF analysis, which was computed in geNorm version 3.4 .
We gratefully acknowledge the financial support of AstraZeneca, Syngenta, NERC (project NER/D/S/2002/00413 "The population and molecular stress responses of an ecotoxicology indicator species") and The Research Endowment Trust Fund of The University of Reading. Furthermore, we are thankful to Dr Stuart N. Peirson and Dr Godfrina Mckoy for valuable discussions on QPCR, and three anonymous reviewers for their valuable comments on the manuscript.
- Orphanides G: Toxicogenomics: challenges and opportunities. Toxicol Lett. 2003, 140–141: 145-148. 10.1016/S0378-4274(02)00500-3.PubMedView ArticleGoogle Scholar
- Dallas PB, Gottardo NG, Firth MJ, Beesley AH, Hoffmann K, Terry PA, Treitas JR, Boag JM, Cummings AJ, Kees UR: Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR – how well do they correlate?. BMC Genomics. 2005, 6: 59-10.1186/1471-2164-6-59.PubMedPubMed CentralView ArticleGoogle Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000, 25: 169-193. 10.1677/jme.0.0250169.PubMedView ArticleGoogle Scholar
- Williams TD, Gensberg K, Minchin SD, Chipman JK: A DNA expression array to detect toxic stress response in European flounder (Platichthys flesus). Aquat Toxicol. 2003, 65: 141-157. 10.1016/S0166-445X(03)00119-X.PubMedView ArticleGoogle Scholar
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.6364190.PubMedView ArticleGoogle Scholar
- Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotech. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7.View ArticleGoogle Scholar
- Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004, 313: 856-862. 10.1016/j.bbrc.2003.11.177.PubMedView ArticleGoogle Scholar
- Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002, 29: 23-39. 10.1677/jme.0.0290023.PubMedView ArticleGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper – Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47.PubMedView ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: research0034.1-0034.11. 10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- geNorm software. [http://medgen.ugent.be/~jvdesomp/genorm/]
- BestKeeper software. [http://www.gene-quantification.info/]
- Heckmann LH, Connon R, Hooper HL, Maund SJ, Hutchinson TH, Sibly RM, Callaghan A: Molecular and population stress responses of Daphnia magna exposed to ibuprofen [abstract]. Society of Environmental Toxicology and Chemistry. Europe 15th Annual Meeting, Lille, France, 308-22–26 May 2005Google Scholar
- Coulson CJ: Molecular Mechanisms of Drug Action. 1994, London: Taylor and FrancisGoogle Scholar
- Stanley-Samuelson DW: Assessing the significance of prostaglandins and other eicosanoids in insect physiology. J Insect Physiol. 1994, 40: 3-11. 10.1016/0022-1910(94)90106-6.View ArticleGoogle Scholar
- Rowley AE, Vogan CL, Taylor GW, Clare AS: Prostaglandins in non-insectan invertebrates: recent insights and unsolved problems. J Exp Biol. 2005, 208: 3-14. 10.1242/jeb.01275.PubMedView ArticleGoogle Scholar
- Tucker CL, Fields S: Quantitative genome-wide analysis of yeast deletion strain sensitivities to oxidative and chemical stress. Comp Func Genomics. 2004, 5: 216-224. 10.1002/cfg.391.View ArticleGoogle Scholar
- Kim KS, Yoon JH, Kim JK, Baek SJ, Eling TE, Lee WJ, Ryu JH, Lee JG, Lee JH, Yoo JB: Cyclooxygenase inhibitors induce apoptosis in oral cavity cancer cells by increased expression of nonsteroidal anti-inflammatory drug-activated gene. Biochem Biophys Res Commun. 2004, 325: 1298-1303. 10.1016/j.bbrc.2004.10.176.PubMedView ArticleGoogle Scholar
- English N, Hughes V, Wolf GR: Induction of cytochrome P-450(BM-3) (CYP 102) by non-steroidal anti-inflammatory drugs in Bacillus megaterium. Biochem J. 1996, 316: 279-283.PubMedPubMed CentralView ArticleGoogle Scholar
- Tarnawski AS, Jones MK: Inhibition of angiogenesis by NSAIDs: molecular mechanisms and clinical implications. J Mol Med. 2003, 81: 627-636. 10.1007/s00109-003-0479-y.PubMedView ArticleGoogle Scholar
- Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004, 15: 155-166.PubMedPubMed CentralGoogle Scholar
- MIAME. [http://www.mged.org/miame]
- Peirson SN, Butler JN, Foster RG: Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis. Nucl Acids Res. 2003, 31: e73-10.1093/nar/gng073.PubMedPubMed CentralView ArticleGoogle Scholar
- Hooper HL, Connon R, Callaghan A, Maund SJ, Liess M, Duquesne S, Hutchinson TH, Moggs J, Sibly RM: The use of image analysis methods to estimate population growth rate in Daphnia magna. J Appl Ecol. 2006,Google Scholar
- Pascoe D, Karntanut W, Muller CT: Do pharmaceuticals affect freshwater invertebrates? A study with the cnidarian Hydra vulgaris. Chemosphere. 2003, 51: 521-528. 10.1016/S0045-6535(02)00860-3.PubMedView ArticleGoogle Scholar
- NCBI. [http://www.ncbi.nlm.nih.gov/]
- wFleaBase. [http://wfleabase.org/]
- Colbourne JK, Singan VR, Gilbert DG: wFleaBase: The Daphnia genome database. BMC Bioinformatics. 2005, 6: 45-10.1186/1471-2105-6-45.PubMedPubMed CentralView ArticleGoogle Scholar
- Primer3 software. [http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi]
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.