- Research article
- Open Access
Reference genes for QRT-PCR tested under various stress conditions in Folsomia candida and Orchesella cincta (Insecta, Collembola)
BMC Molecular Biologyvolume 10, Article number: 54 (2009)
Genomic studies measuring transcriptional responses to changing environments and stress currently make their way into the field of evolutionary ecology and ecotoxicology. To investigate a small to medium number of genes or to confirm large scale microarray studies, Quantitative Reverse Transcriptase PCR (QRT-PCR) can achieve high accuracy of quantification when key standards, such as normalization, are carefully set. In this study, we validated potential reference genes for their use as endogenous controls under different chemical and physical stresses in two species of soil-living Collembola, Folsomia candida and Orchesella cincta. Treatments for F. candida were cadmium exposure, phenanthrene exposure, desiccation, heat shock and pH stress, and for O. cincta cadmium, desiccation, heat shock and starvation.
Eight potential reference genes for F. candida and seven for O. cincta were ranked by their stability per stress factor using the programs geNorm and Normfinder. For F. candida the succinate dehydrogenase (SDHA) and eukaryotic transcription initiation factor 1A (ETIF) genes were found the most stable over the different treatments, while for O. cincta, the beta actin (ACTb) and tyrosine 3-monooxygenase (YWHAZ) genes were the most stable.
We present a panel of reference genes for two emerging ecological genomic model species tested under a variety of treatments. Within each species, different treatments resulted in differences in the top stable reference genes. Moreover, the two species differed in suitable reference genes even when exposed to similar stresses. This might be attributed to dissimilarity of physiology. It is vital to rigorously test a panel of reference genes for each species and treatment, in advance of relative quantification of QRT-PCR gene expression measurements.
Genomic techniques have undergone major developments in the last two decades. As a result, they have become conducive for evolutionary and ecotoxicological studies, which generally use non genomic-model organisms. Quantitative Reverse-Transcriptase Polymerase Chain Reaction (QRT-PCR) is a technique to estimate gene expression levels. This technique is often used to confirm high throughput systems like microarrays. Its application has mainly been limited to small numbers of genes per experiment due to constraints of low throughput coinciding with relatively high costs per assay. This is about to change, as high throughput QRT-PCR systems using small volume (capillary) PCR are becoming available . QRT-PCR is a valuable tool for ecological studies as it provides a relatively straightforward way to measure the direct transcriptional response of an organism exposed to different treatments . QRT-PCR has been applied to study adaptive evolution at the transcriptional level [3–6]. For instance, Roelofs et al.  conducted a QRT-PCR study to assess the relevance of transcriptional regulation in the adaptive evolution of stress tolerance.
The conditions that have to be met for a successful QRT-PCR experiment are reviewed by Bustin . Different strategies have been developed for quantifying gene expression with QRT-PCR data. The most widely used method is to quantify the relative amount of target mRNA between samples, using for instance the comparative CT (2-ΔΔCT) method , or the more comprehensive Pfaffl method where relative quantities are adjusted for amplification efficiencies. Relative quantification methods depend on reference genes for normalization . QRT-PCR reference genes, sometimes called 'housekeeping' genes, can either be internal, when measured in the same reaction tube as the target gene, or external, when measured in a different reaction tube. The use of reference genes is necessary to correct for factors such as RNA input differences and reverse transcriptase efficiency variation .
An essential requirement for QRT-PCR reference genes is stability of their transcriptional level across the various conditions to which an organism is exposed during an experiment. Classically, the most used reference genes are carry-overs from the Northern-blot days  such as beta actin (ACTb), glyceraldehyde-3P-dehydrogenase (GAPDH) and ribosomal RNA genes (18S, 28S). Unfortunately, these classical reference genes are not always found to be suitable for this use. For example, mRNA levels of ACTb and GAPDH can fluctuate widely in human T-cells exposed to different treatments . Also, commonly used housekeeping genes like ACTb, GAPDH, cyclophilin A (CYP) and 28S are up or down regulated in cell lines exposed to hypoxia stress . Therefore, it is vital to validate the stability of a panel of reference genes in order to select the most suitable ones for each new treatment and species of choice.
Currently two Visual Basic Applications for Microsoft Excel are widely used to determine reference gene suitability: geNorm  and Normfinder . GeNorm is based on the principle that the expression ratio of two ideal control genes should be identical in all samples and experimental conditions. It calculates gene expression stability (M), which is the mean pair-wise variation between an individual gene and all other tested reference genes. Subsequently, the least stable reference gene (highest M value) is excluded from the set and the calculation is reiterated until the two most stable reference genes remain. Normfinder  estimates intergroup and intragroup variation to calculate reference gene stability and to rank them. As with geNorm it calculates a stability value for each potential reference gene but it uses a variance model based approach (mixed linear effect modeling), instead of the reiterative approach used by geNorm. Normfinder calculates intragroup variability for the genes in each of the groups, and the intergroup variability or bias among the groups . Since the program can differentiate between groups, Normfinder is best suited when the stability of reference genes needs to be assessed over multiple treatments. When reference gene stability only needs to be calculated for samples exposed to a single treatment the two methods are similar and should give the same results .
In this paper we develop a panel of reference genes for two species of springtails (Collembola), and study their stability across various treatments. Collembola are important model organisms in evolutionary ecology and ecotoxicological studies. The soil-dwelling collembolan Folsomia candida is a standard test animal used in a soil pollution test that is certified by the International Organization for Standardization . Its parthenogenetic mode of reproduction and the availability of a recently sequenced Expressed Sequence Tag (EST) database, make F. candida also a proper test animal for genomic studies on the effects of soil pollution [19–21]. An F. candida microarray is currently in use for testing soil toxicants  and physical conditions (Timmermans et al., unpubl. data; De Boer et al., unpubl. data). Therefore, confirmation of microarray results by QRT-PCR is abound. Orchesella cincta is a sexually reproducing collembolan that lives in the litter layer rather than in the soil and is generally used to study adaptation and phenotypic plasticity [23, 24]. Furthermore, it is an emerging genomic model to study adaptive evolution in polluted environments [20, 21, 25].
Here we study the stability of potential QRT-PCR reference genes across various treatments. We exposed the two collembolan species to several stressors (F. candida: cadmium, phenanthrene, desiccation, temperature and pH stress. O. cincta: cadmium, desiccation, temperature stress and starvation) that are currently under investigation in gene expression studies. Based on previous studies (e.g. [12, 26]), a panel of potential QRT-PCR reference genes was developed (Table 1) using the Collembola EST database Collembase http://www.collembase.org. In addition, a target gene (expected to be differential) was measured for each treatment, to validate the impact of the treatment at the transcriptional level of the organism. The stability of the potential reference genes was determined using both geNorm and Normfinder.
Expression levels under different conditions
Average cycle threshold (Ct) values varied widely between conditions and treatments, and ranged between 9.7 (28S, O. cincta) and 28.6 (CYP, F. candida, see Additional file 1 for an overview). An example of the different expression levels in treatments and conditions is given in Figure 1 for ACTb and the differential metallothionein gene (MT) of O. cincta.
In O. cincta, the target gene MT showed a 7 fold up regulation (FR) in the cadmium treatment (P < 0.05) and HSP70 was 5 fold up regulated after heat shock (P < 0.001). Starvation and desiccation treatments showed no significant changes in either of the two target genes, implying that those treatments had fairly small effects, at least on the transcriptional response of these two genes. In F. candida, all but one treatment showed a significant effect on the target genes (cadmium: mitochondrial chaperone BCS-1 (BCS1), FR = 20, P < 0.0001; phenanthrene: BCS1, FR = 6.4, P = 0.001; temperature: HSP70, FR > 2, P < 0.001 (for all conditions); desiccation: Cuticle protein (CP), FR = 4, P < 0.05). The pH treatment showed no significant differences in the target V-type ATPase (ATPase); in fact the ATPase gene used in the pH treatment was more stable than most reference genes, again implying that this treatment had only a modest effect on the transcriptional level of this gene.
Ranking of the reference genes
Stability rankings were established using geNorm and Normfinder. As described in the methods section, the Normfinder analysis should be executed with genes that show no substantial treatment-specific response . Therefore, differential genes were excluded in most cases as well as some of the other genes that showed a systematic response to the treatment (see Table 2 for an overview).
The stability rankings generated by geNorm or Normfinder were largely similar, even though the ranking order of the genes differed to some extent. An overview of which genes to use for each treatment is found in Table 3. From the three treatments tested for both species, only the cadmium treatment gave similar results between F. candida (tyrosine 3-monooxygenase (YWHAZ), succinate dehydrogenase (SDHA) and GAPDH) and O. cincta (alpha-tubulin (TBa) – not tested in F. candida – SDHA and YWHAZ). To compare interspecies parallels, an additional geNorm analysis was done including only the reference genes available for both species (see Additional file 2). This did not change the fact that most rankings did not correspond between species. In F. candida, overall analyses show the same outcomes for the cadmium treatment and the phenanthrene treatment: YWHAZ, SDHA, GAPDH. For the temperature and desiccation treatments (F. candida) results also overlapped with SDHA ETIF and elongation factor 1α (EF1a) being the best suited reference genes. The results for O. cincta seem to be quite variable. The commonly used reference gene ACTb is placed in all top rankings, but not in the cadmium treatment. 28S, a well known but also controversial reference gene to use, showed a high stability in the temperature treatment. Also remarkable is the stability of HSP70, expected to be differential, in the desiccation and starvation treatments. This result has been reported previously for Drosophila melanogaster, which did not show an increase in HSP70Aa mRNA levels in response to starvation and desiccation . Despite the fact that the ranking of reference genes shows no typical overall uniformity over all different treatments, the most generally applicable reference genes are suggested to be: ACTb, YWHAZ and TBa for O. cincta and SDHA, ETIF and YWHAZ for F. candida.
Optimum number of reference genes
The optimum numbers of reference genes are shown in Table 3. In only two out of nine treatments the two methods agreed on the number of genes to use for normalization (two genes in these cases). However, in nearly all cases either geNorm or Normfinder recommended the use of two reference genes. Comparing the two sets proposed by either program, the differences in significance levels and relative expression levels of a differential gene were found to be very small.
Effect of method of normalization on relative expression of HSP70
We selected the temperature treatment with the differential HSP70 to illustrate the effects of reference gene selection on the calculated relative expression level of a gene of interest. This treatment was chosen because there were four different conditions and HSP70 expression clearly responded differently to each of these conditions.
We normalized HSP70 expression in the F. candida temperature treatment with four different sets of reference genes: i) only with a commonly used single reference gene (ACTb); ii) with the appropriate number of selected genes according to the geNorm and iii) the Normfinder analysis and iv) with all available reference genes in the F. candida temperature dataset. Each set of reference genes showed significant upregulation of HSP70 expression in the animals exposed to 0°C and 30°C as compared to those exposed to 10°C and 20°C (P < 0.05), as well as a significant difference between 10°C and 20°C (Figure 2). However, the difference in HSP70 expression between 0°C and 30°C was only significant when using the reference genes proposed by the geNorm and Normfinder analysis. This indicates that selection of reference genes can influence the resolution with which differences in gene expression between two samples can be detected.
Ranking O. cincta reference genes over all treatments
Making comparisons between experiments is difficult when no standardization of reaction conditions has taken place . Reverse transcription is the most crucial step in the introduction of technical sampling variation in the QRT-PCR procedure . In the O. cincta dataset we attempted to meet the standards needed for a multi-treatment analysis, which allowed us to analyze the whole O. cincta dataset similarly to the single treatment analyses using Normfinder. We also included a non-normalized situation ('total RNA input') which was equal for all samples. The differential MT was omitted from ranking due to bias beyond the bias threshold as well as 'total RNA input' (Table 2, see materials section and Additional file 3. Although geNorm is less suited to handle heterogeneity than Normfinder, both algorithms produced the same top three of generally applicable reference genes: ACTb TBa, and YWHAZ.
Discussion and Conclusion
In this study we identified appropriate reference genes of two collembolan species that show invariant expression across experimental treatments. The most important result of our study is the lack of similarity in stability of genes i) within the species upon exposure to different treatments, and ii) between species that have undergone similar treatments. The two collembolan species have a distinct physiology and ecology, but they belong to phylogenetically related families . An earlier study on two flatfish species with a comparable phylogenetic relatedness as the two collembolan species used here, found congruent results in reference gene stabilities between the two species . In our study, differences were most pronounced in the temperature and desiccation treatments, while the cadmium treatments show a greater congruency between species. It may be argued that the physiological response to xenobiotic stress is more similar in the two species than responses to the environmental factors. Still, experimental methods also differed between species and therefore it cannot be excluded that part of the observed interspecific differences was due to differences in experimental set-up.
Previously, it has been suggested that the results of reference gene selection studies might serve as a resource for future gene expression studies in the same or related species [31, 32]. Indeed, some genes were stable across different treatments. ACTb and YWHAZ were the dominant genes found in the top of the rankings of three out of four O. cincta treatments, but this general pattern did not hold up for ACTb in the cadmium treatment and YWHAZ in the desiccation treatment. In F. candida, the two chemical treatments (cadmium and phenanthrene) and the environmental treatments (temperature, desiccation and pH-stress) differed in selected reference genes, indicating that different classes of physiological responses, for instance detoxification and acclimation, can be expected to cause fluctuations in different 'housekeeping' genes to maintain cellular homeostasis. Therefore one cannot assume that two species responding similarly to a certain treatment will also respond in a similar way to another treatment. In addition, differences in priming strategy and QRT-PCR assay characteristics can also introduce variation in stability comparisons. In the light of our present results we must caution against the use of literature data from related species to select a normalization standard, unless careful notion is taken of the parallels and differences in the technical context of the experiment as well as the internal processes of the organisms that are being studied.
Optimizing a set of reference genes not only requires making the right choice of genes, it also implies choosing the right number of genes. The difference in calculated levels of HSP70 expression, using different sets of reference genes, exemplifies the importance of the number of reference genes used. Resolution was too low to detect small differences in expression, when either a single reference gene was used for normalization or when all available reference genes were included. In fact the small difference between 0°C and 30°C was only statistically significant when normalization was done with the sets proposed by geNorm and Normfinder.
In the geNorm analysis, inclusion of a differential gene gives insight into the variability of all other genes . Hypothetically, this differential gene should be the most unstable gene in the dataset. In both our species and all treatments (except for the pH treatment) at least one of the selected differential genes was indeed the most variable of the set tested. This provides evidence that the applied treatments caused an effect at the transcriptional level, which strengthens the validation of the reference genes that remained stable under the changed regimes. The only exception was observed for the pH treatment, where the pre-selected differential gene was positioned among the housekeeping genes. Most likely, the treatment may not have been severe enough or it may indicate that some exposure types do not initiate an effect on the transcriptional level.
Normfinder performs relative comparisons between the potential reference genes; hence this method is more sensitive to the presence of differentially expressed genes in the dataset than geNorm. In this study, we set a threshold for the maximal allowable bias at 0.13 times the standard deviation of the treatments' intergroup variation as calculated with Normfinder. After the pre-selection, the reference genes selected by geNorm and Normfinder were remarkably similar. In all treatments, except for the F. candida pH treatment and the O. cincta starvation treatment, the top three consisted of the same three genes even though the order differed.
QRT-PCR has proven its value in many areas of genetic and genomic research. Knowledge on genetic pathways and molecular responses to external environmental cues and chemical factors now also make a significant contribution to evolutionary ecology and ecotoxicology. Experiments in these scientific disciplines often focus on non-model organisms, such as Collembola, and routinely use large sample sizes. Now that high-throughput systems QRT-PCR are becoming available , the technique will become even more valuable for molecular ecological research. The reference genes presented in this study can therefore act as a starting point for scientists who use identical collembolan species for ecological, evolutionary and ecotoxicological research. Nonetheless, as previously stated by Stürzenbaum and Kille , it should always be kept in mind that techniqually successful QRT-PCR depends not only on the genes of interest alone, and therefore reference genes should be carefully validated prior to experimentation.
Folsomia candida was kept in plastic containers with a water-saturated plaster of Paris bottom containing 10% charcoal at 20°C in a 12:12 light dark regime. The animals were fed dried baker's yeast (Dr. Oetker) ad libitum. For all experiments animals of at least 20 days old were used.
Orchesella cincta was held comparably, but fed algae (Desmococcus sp) growing on twigs of pine trees. For all the experiments animals of 4–5 weeks old were used, with a maximum age difference of seven days.
Folsomia candida treatments
pH & temperature treatment
Animals of 20 days old were exposed in 100 ml jars with 30 grams of OECD artificial soil. OECD soil pH was adjusted with CaCO3 (J.T. Baker) to four different values (3.5, 4.5, 5.5 & 6.5) according to OECD guideline 207 . Pools of thirty animals per jar were exposed to the different pH values for three days at 75% humidity and 20°C. For the temperature treatment OECD soil at a pH of 5.5 was used. Pools of thirty animals per jar were exposed to 0, 10, 20 and 30°C for 3 days. For each condition two biological replicates were used (separately) for RNA extraction.
Drought exposure was performed as described by Bayley and Holmstrup . Per replicate, a pool of twenty-five to thirty animals was exposed in plastic containers to a relative air humidity controlled at 98.2% by placing a NaCl solution of 31.6 g L-1 inside the closed container. Animals were sacrificed after 8, 27, 53 and 174 hours of exposure for RNA isolation. For each condition two biological replicates were used (separately) for RNA extraction.
Cadmium and phenanthrene treatment
For the cadmium and phenanthrene treatments animals were exposed in 100 ml jars on a compressed layer of 10 g wet weight of LUFA 2.2 soils (for details see ). Procedures of the standard ISO protocol 11267  were followed to spike the soils to nominal concentrations equivalent to the LC50 28 days for cadmium (6.86 mmol kg-1 dry soil ) and phenanthrene (422 μmol kg-1 dry soil ). For the cadmium spiking a solution of hydrated CdCl2 (purity 99%; J.T. Baker, The Netherlands) was used, while phenanthrene (purity 96%; Sigma-Aldrich Chemie, Germany) was dissolved in acetone (Riedel-de Haën, Seelze, Germany), followed by an overnight evaporation time for the acetone. Both soils did not undergo a period of aging. Pools of 15 animals per jar were exposed for a period of 48 hr (cadmium) or 96 hr (phenanthrene) to spiked and clean control soils. For each condition four biological replicates were used (separately) for RNA extraction.
O. cincta treatments
Cadmium exposure was performed as described in Roelofs et al. . Animals were exposed individually to a nominal concentration of 1 μmole cadmium per gram dry weight food (algal paste). The exposure was started immediately after moulting to exclude hormonal effects on gene expression and lasted for three days. Five individuals were pooled per replicate; a control was used with the same set-up but fed clean algae. For each condition three biological replicates were used (separately) for RNA extraction.
Springtails were exposed in glass vials containing slightly moistened foam at the bottom and moistened foam stoppers. Five individuals per vial were exposed to three temperature treatments: cold (10°C), control (20°C) and heat (35°C). Temperature treatments consisted of four hours at 10°C respectively 20°C in a climate room, or one hour placement in a water bath of 35°C with a one hour recovery period at 20°C. For each condition three biological replicates of five pooled animals were used (separately) for RNA extraction.
The desiccation treatment followed the protocol described by Bahrndorff et al. . Springtails were exposed to 97.2% relative humidity in a tightly sealed container containing a NaCl solution of 50.66 g L-1 at 20°C for five days. The control treatment followed the same protocol, but instead of the NaCl solution demineralized water was used. We used five individuals per vial and three vials per treatment. Animals from each vial were pooled for RNA isolation.
We used five individuals per vial and three vials for the starvation treatment and its control. Springtails were transferred to glass vials containing slightly moistened foam at the bottom and moistened foam stoppers. In the control treatment food was made available by adding a piece of bark overgrown with green algae to the animals, while the animals in the starvation treatment did not have access to food. The experimental vials were kept at 20°C in a 12:12 light dark regime for 8 days. Animals from each vial were pooled for RNA isolation.
RNA isolation and reverse transcription
After exposure animals were snap-frozen in liquid nitrogen and total RNA was isolated with the SV Total RNA isolation system (Promega) according to manufacturer's instructions. Genomic DNA was removed via a DNAse treatment supplied with the kit. RNA integrity was confirmed on a 1% agarose gel and RNA quantities were assessed with a nanodrop ND-1000 spectrophotometer (Nanodrop Technologies) and ranged between 30 and 100 ng μL-1 of total RNA. As indicated by 260/280 and 260/230 nm ratios, all samples used in this study were assumed free from protein contamination and (organic) salts. Absence of amplicons after PCR with Taq-polymerase (MRC Holland, The Netherlands) and 1 μL RNA solution confirmed that no trace DNA contamination was present in the samples used in the further analyses.
Synthesis of cDNA was performed using the reverse-transcriptase system of Promega with the M-MLV reverse transcriptase enzyme and an oligo-dT primer (F. candida samples) or random hexamer primers (O. cincta). Random hexamer primers were used in the case of O. cincta because 28S ribosomal RNA was used as one of the potential reference genes. In the case of the O. cincta samples, reverse transcription input amounts were equalized by diluting the total RNA concentrations to 0.5 μg μL-1, and samples were reverse transcribed together in a single run. cDNA samples were diluted 4 times before QPCR was carried out. All treatment conditions were reverse transcribed in triplicate, except for the pH treated samples which were performed in duplicate.
Besides reference genes, one or two differentially expressed genes were included in order to observe transcriptional effects of the treatment. These differentials were previously assessed for their response to the majority of the treatments ([25, 39, 40], TE de Boer, MJTN Timmermans, unpubl. data). For the desiccation and starvation treatments in O. cincta, lack of prior knowledge of characteristic responsive genes made us look to a fixed set of two differentials for all four treatments. QPCR assays for seven candidate reference genes (Table 1) and two genes differential for a range of stressful conditions, HSP70 and MT, [41, 42] were developed. For F. candida, eight candidate reference genes were analyzed (Table 1) together with one differential gene for each treatment. Primers were based on sequences present in the Collembola EST database Collembase http://www.collembase.org and generated with Primer Express 1.5 (Applied Biosystems) with the following settings: Melting temperatures were kept between 59°C and 60°C and the amplified fragment length was kept between 90 and 120 base pairs with an optimum amplicon melting temperature of 80°C. For primer sequences, GC content and melting temperatures (see Additional file 4). Details on the positions of the QPCR amplicons in the full coding sequences are given in (Additional file 5).
Reaction efficiency of each of the QPCR assays was determined by means of a standard curve consisting of 5 samples each four fold diluted, from an initial cDNA pool. Each reaction was carried out in a total volume of 20 μL, using 2 μL cDNA template, 10 μL SYBR Green I master mix (Applied Biosystems) and 20 pmole of each gene specific primer (supplied by Isogen). QPCR Cycling was performed on an DNA engine Opticon1 (Biorad), with three technical replicates per sample. Cycling conditions were kept constant for all assays. For details on PCR mix and program see Roelofs et al. ). Specificity of the PCR products was confirmed after each amplification by analysis of the melting curve; 60–90°C with a heating rate of 0.1°C per second and one fluorescence measurement per second. Each run included a non-template control for each assay.
For two of the O. cincta assays (GAPDH and SDHA) no sequence information was available. Therefore degenerate primers were developed based on sequences from multiple organisms that were taken from Genbank. Generated PCR templates were subsequently cloned, sequenced and used as a template for QPCR primers as described above. Proper QPCR data were not retrieved from one F. candida pH treatment replicate (pH 3.5). Therefore this sample was discarded from further analysis.
Ct values were calculated with the Opticon Monitor 3 software (Biorad), using a manually set cycle threshold. At a level of 0.01 raw fluorescent units which in all assayed plates fell within the exponential phases of the QPCR reactions. Averages of the three technical replicates were used in case of a standard deviation ≤ 0.5. When standard deviation exceeded this number, fluorescence curves were evaluated. The analysis was always based on at least two replicates.
Ranking of reference genes was determined using geNorm and Normfinder applications, as implemented in the Genex Light software package . Optimum numbers of reference genes are based on a Vn/Vn+1 value of > 0.15 for geNorm  and the minimum of accumulated standard deviations for Normfinder . The original geNorm VBA applet for Excel was used for automated calculation of Vn/Vn+1.
Relative gene expressions of all differential genes were calculated with the Pfaffl method  and normalized with both of the optimum sets of reference genes proposed by geNorm and Normfinder. Significance levels were tested by Student's t-test, for both algorithms. To be conservative we report only the larger P-value in the results. The temperature exposures of F. candida (Figure 2) were normalized as described in the results section, and tested for their significance by one-way ANOVA analysis with a Bonferroni corrected post-hoc test between the different temperatures for each normalization method.
Bias threshold definition
In Normfinder the total stability ('variability') of a gene is defined by the magnitude of the intragroup variation relative to the intergroup variation . Differential genes, which are responding to one or more treatments, will greatly increase the intergroup variation and hence have a disproportionally large effect on the calculated variability and the ranking order of the other candidate genes. To avoid such biases, differential genes should be excluded from the NormFinder analysis using a pre-selection procedure. The pre-selection procedure consisted of an initial NormFinder analysis with all genes, from which we calculated a bias threshold for the amount of total absolute bias a gene was allowed to have, compared to the total mean absolute bias of the group. Assuming a normal distribution of biases and genes, we considered the 10% most stable genes to be suitable as potential reference genes, which sets the bias threshold between -0.13*SD and 0.13*SD of the normal distribution . In the final NormFinder analysis only those genes that met the criterion of being among the 10% most stable genes were included. Ranking of the remaining genes was subsequently done by determining the standard deviation in a Normfinder analysis with the software settings set to not taking groups into account (M. Kubista, pers. comm.).
Morrison T, Hurley J, Garcia J, Yoder K, Katz A, Roberts D, Cho J, Kanigan T, Ilyin SE, Horowitz D, Dixon JM, Brenan C: Nanoliter high throughput quantitative PCR. Nucleic Acids Res. 2006, 34 (18):
Snell TW, Brogdon SE, Morgan MB: Gene Expression Profiling in Ecotoxicology. Ecotoxicology. 2003, 12 (6): 475-483. 10.1023/B:ECTX.0000003033.09923.a8
Ellers J, Mariën J, Driessen G, Van Straalen NM: Temperature-induced gene expression associated with different thermal reaction norms for growth rate. J Exp Zool. 2008, 310B (2): 137-147. 10.1002/jez.b.21194. 10.1002/jez.b.21194
Muller WEG, Grebenjuk VA, Thakur NL, Thakur AN, Batel R, Krasko A, Muller IM, Breter HJ: Oxygen-controlled bacterial growth in the sponge Suberites domuncula: toward a molecular understanding of the symbiotic relationships between sponge and bacteria. Appl Environ Microbiol. 2004, 70 (4): 2332-2341. 10.1128/AEM.70.4.2332-2341.2004
Scharf ME, Wu-Scharf D, Pittendrigh BR, Bennett GW: Caste- and development-associated gene expression in a lower termite. Genome Biol. 2003, 4 (10): 11-10.1186/gb-2003-4-10-r62. 10.1186/gb-2003-4-10-r62
Zientz E, Beyaert N, Gross R, Feldhaar H: Relevance of the endosymbiosis of Blochmannia floridanus and carpenter ants at different stages of the life cycle of the host. Appl Environ Microbiol. 2006, 72 (9): 6027-6033. 10.1128/AEM.00933-06
Roelofs D, Mariën J, Van Straalen NM: Differential gene expression profiles associated with heavy metal tolerance in the soil insect Orchesella cincta. Insect Biochem Mol Biol. 2007, 37 (4): 287-295. 10.1016/j.ibmb.2006.11.013
Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002, 29 (1): 23-39. 10.1677/jme.0.0290023
Applied_Biosystems: Comparative CT Method (Separate Tubes). User Bulletin 2: Rev. B ABI Prism 7700 Sequence Detection System. 2001, 2: 11-15.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29 (9): 2002-2007. 10.1093/nar/29.9.e45. 10.1093/nar/29.9.e45
Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6 (4): 279-284. 10.1038/sj.gene.6364190
Heckmann L-H, Connon R, Hutchinson T, Maund S, Sibly R, Callaghan A: Expression of target and reference genes in Daphnia magna exposed to ibuprofen. BMC Genomics. 2006, 7 (1): 175- 10.1186/1471-2164-7-175
Bas A, Forsberg G, Hammarstrom S, Hammarstrom ML: Utility of the Housekeeping Genes 18S rRNA, beta-Actin and Glyceraldehyde-3-Phosphate-Dehydrogenase for Normalization in Real-Time Quantitative Reverse Transcriptase-Polymerase Chain Reaction Analysis of Gene Expression in Human T Lymphocytes. Scand J Immunol. 2004, 59 (6): 566-573. 10.1111/j.0300-9475.2004.01440.x
Zhong H, Simons JW: Direct Comparison of GAPDH, beta-Actin, Cyclophilin, and 28S rRNA as Internal Standards for Quantifying RNA Levels under Hypoxia. Biochem Biophys Res Commun. 1999, 259 (3): 523-526. 10.1006/bbrc.1999.0815
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 (7): research0034.0031-research0034.0011. 10.1186/gb-2002-3-7-research0034. 10.1186/gb-2002-3-7-research0034
Andersen CL, Jensen JL, Orntoft TF: Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: A Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and Colon Cancer Data Sets. Cancer Res. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496
Bergkvist A, Forootan A, Zoric N, Strombom L, Sjoback R, Kubista M: Choosing a Normalization Strategy for RT-PCR. Genet Eng Biotechn N. 2008, 28 (13): 1-
, : ISO, Soil Quality. Inhibition of Reproduction of Collembola (Folsomia candida). ISO Guideline 11267. 1999, 16-International Standardization Organization. Switzerland,
Fountain MT, Hopkin SP: Folsomia candida (Collembola): A "standard" soil arthropod. Annu Rev Entomol. 2005, 50: 201-222. 10.1146/annurev.ento.50.071803.130331
Timmermans MJTN, De Boer ME, Nota B, De Boer TE, Mariën J, Klein-Lankhorst RM, Van Straalen NM, Roelofs D: Collembase: a repository for springtail genomics and soil quality assessment. BMC Genomics. 2007, 8 (341): 10-
Timmermans MJTN, Ellers J, Van Straalen NM: Allelic diversity of metallothionein in Orchesella cincta (L.): traces of natural selection by environmental pollution. Heredity. 2007, 98 (5): 311-319. 10.1038/sj.hdy.6800942
Nota B, Timmermans MJTN, Franken O, Montagne-Wajer K, Mariën J, De Boer ME, De Boer TE, Ylstra B, Van Straalen NM, Roelofs D: Gene Expression Analysis of Collembola in Cadmium Containing Soil. Environ Sci Technol. 2008, 42 (21): 8152-8157. 10.1021/es801472r
Bahrndorff S, Holmstrup M, Petersen H, Loeschcke V: Geographic variation for climatic stress resistance traits in the springtail Orchesella cincta. J Insect Physiol. 2006, 52 (9): 951-959. 10.1016/j.jinsphys.2006.06.005
Liefting M, Ellers J: Habitat-specific differences in thermal plasticity in natural populations of a soil arthropod. Biol J Linn Soc. 2008, 94 (2): 265-271. 10.1111/j.1095-8312.2008.00969.x. 10.1111/j.1095-8312.2008.00969.x
Roelofs D, Overhein L, De Boer ME, Janssens TKS, Van Straalen NM: Additive genetic variation of transcriptional regulation: metallothionein expression in the soil insect Orchesella cincta. Heredity. 2006, 96 (1): 85-92.
Spinsanti G, Panti C, Lazzeri E, Marsili L, Casini S, Frati F, Fossi C: Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies. BMC Mol Biol. 2006, 7 (1): 32- 10.1186/1471-2199-7-32
Sinclair BJ, Gibbs AG, Roberts SP: Gene transcription during exposure to, and recovery from, cold and desiccation stress in Drosophila melanogaster. Insect Mol Biol. 2007, 16 (4): 435-443. 10.1111/j.1365-2583.2007.00739.x
Bustin SA, Benes V, Nolan T, Pfaffl MW: Quantitative real-time RT-PCR – a perspective. J Mol Endocrinol. 2005, 34 (3): 597-601. 10.1677/jme.1.01755
Stahlberg A, Kubista M, Pfaffl M: Comparison of Reverse Transcriptases in Gene Expression Analysis. Clin Chem. 2004, 50 (9): 1678-1680. 10.1373/clinchem.2004.035469
D'Haese CA: Were the first springtails semi-aquatic? A phylogenetic approach by means of 28S rDNA and optimization alignment. P Roy Soc Lond B-Bio. 2002, 269 (1496): 1143-1151. 10.1098/rspb.2002.1981. 10.1098/rspb.2002.1981
Infante C, Matsuoka MP, Asensio E, Canavate JP, Reith M, Manchado M: Selection of housekeeping genes for gene expression studies in larvae from flatfish using real-time PCR. BMC Mol Biol. 2008, 9 (28): 12-
Scharlaken B, de Graaf DC, Goossens K, Brunain M, Peelman LJ, Jacobs FJ: Reference gene selection for insect expression studies using quantitative real-time PCR: The head of the honeybee, Apis mellifera, after a bacterial challenge. J Insect Sci. 2008, 8 (33): 10-
Stürzenbaum SR, Kille P: Control genes in quantitative molecular biological techniques: the variability of invariance. Comp Biochem Physiol B: Biochem Mol Biol. 2001, 130 (3): 281-289. 10.1016/S1096-4959(01)00440-7
Organisation for Economic Co-operation and Development: Test no 207: Earthworm, acute toxicity tests. OECD Guidelines for the testing of chemicals. 1984, OECD: Paris, France,
Bayley M, Holmstrup M: Water vapor absorption in arthropods by accumulation of myoinositol and glucose. Science. 1999, 285 (5435): 1909-1911. 10.1126/science.285.5435.1909
Dröge STJ, Paumen ML, Bleeker EAJ, Kraak MHS, Van Gestel CAM: Chronic toxicity of polycyclic aromatic compounds to the springtail Folsomia candida and the enchytraeid Enchytraeus crypticus. Environ Toxicol Chem. 2006, 25 (9): 2423-2431. 10.1897/05-628R.1
Van Gestel CAM, Koolhaas JE: Water-extractability, free ion activity, and pH explain cadmium sorption and toxicity to Folsomia candida (Collembola) in seven soil-pH combinations. Environ Toxicol Chem. 2004, 23 (8): 1822-1833. 10.1897/03-393
Bahrndorff S, Petersen SO, Loeschcke V, Overgaard J, Holmstrup M: Differences in cold and drought tolerance of high arctic and sub-arctic populations of Megaphorura arctica Tullberg 1876 (Onychiuridae: Collembola). Cryobiology. 2007, 55 (3): 315-323. 10.1016/j.cryobiol.2007.09.001
De Boer ME, De Boer TE, Nota B, Timmermans MJTN, Mariën J, Roelofs D, Van Straalen NM: Application of real-time RT-QPCR to develop new biomarkers for soil quality assessment [abstract]. Comp Biochem Physiol A Mol Integr Physiol. 2008, 151 (1): S21-S21. 10.1016/j.cbpa.2008.05.079. 10.1016/j.cbpa.2008.05.079
Bahrndorff S, Loeschcke V, Marien J, Ellers J: Dynamics of heat-induced thermal stress resistance and hsp70 expression in the springtail, Orchesella cincta. Funct Ecol. 2009, 23 (2): 233-239. 10.1111/j.1365-2435.2009.01541.x. 10.1111/j.1365-2435.2009.01541.x
Sørensen JG, Kristensen TN, Loeschcke V: The evolutionary and ecological role of heat shock proteins. Ecol Lett. 2003, 6 (11): 1025-1037. 10.1046/j.1461-0248.2003.00528.x. 10.1046/j.1461-0248.2003.00528.x
Coyle P, Philcox JC, Carey LC, Rofe AM: Metallothionein: The multipurpose protein. Cell Mol Life Sci. 2002, 59 (4): 627-647. 10.1007/s00018-002-8454-2
MultiD_Analyses_AB: GenEx Light Software, version 4.3.5. MultiD Analyses AB Year of release. 2008, http://www.multid.se
Online Statistics: An Interactive Multimedia Course of Study. http://www.onlinestatbook.com
The authors of this paper would like to thank Mikael Kubista for help with the Genex software and the data analysis and Martin Holmstrup for providing the samples for the F. candida desiccation treatment. This work was supported by grants from the Netherlands Genomics Initiative (BSIK program "Ecogenomics: Assessing the Living Soil"). JE is supported by the Netherlands Organization for Scientific Research, VIDI-grant no. 864.03.003.
MEB, TEB, DR, JE and NMvS conceived the study, set up its design and drafted the manuscript. MEB and TEB conducted the F. candida experiments and the QRT-PCR work. JE and DR performed the O. cincta treatments and RNA extractions. JM sequenced and designed the QPCR assays for the unknown genes. MJTNT partially generated material for the F. candida desiccation analysis and helped with QPCR assay design. BN helped conducting some of the treatments. All authors helped shaping and approved the final manuscript.
Muriel E de Boer, Tjalf E de Boer contributed equally to this work.