Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Variable Temperature Stress in the Nematode Caenorhabditis elegans (Maupas) and Its Implications for Sensitivity to an Additional Chemical Stressor

  • Nina Cedergreen ,

    Contributed equally to this work with: Nina Cedergreen, Nils Jakob Nørhave

    Affiliation Department of Plant and Environmental Sciences, University of Copenhagen, Thorvaldsensvej 40, 1871 Frederiksberg C, Copenhagen, Denmark

    ⨯
  • Nils Jakob Nørhave ,

    Contributed equally to this work with: Nina Cedergreen, Nils Jakob Nørhave

    Affiliation Department of Plant and Environmental Sciences, University of Copenhagen, Thorvaldsensvej 40, 1871 Frederiksberg C, Copenhagen, Denmark

    ⨯
  • Claus Svendsen,

    Affiliation Centre for Ecology and Hydrology, Maclean Building, Benson Lane, Crowmarsh Gifford, Wallingford, Oxfordshire OX10 8BB, United Kingdom

    ⨯
  • David J. Spurgeon

    dasp@ceh.ac.uk

    Affiliation Centre for Ecology and Hydrology, Maclean Building, Benson Lane, Crowmarsh Gifford, Wallingford, Oxfordshire OX10 8BB, United Kingdom

    ⨯

Abstract

A wealth of studies has investigated how chemical sensitivity is affected by temperature, however, almost always under different constant rather than more realistic fluctuating regimes. Here we compared how the nematode Caenorhabditis elegans responds to copper at constant temperatures (8–24°C) and under fluctuation conditions of low (±4°C) and high (±8°C) amplitude (averages of 12, 16, 20°C and 16°C respectively). The DEBkiss model was used to interpret effects on energy budgets. Increasing constant temperature from 12–24°C reduced time to first egg, life-span and population growth rates consistent with temperature driven metabolic rate change. Responses at 8°C did not, however, accord with this pattern (including a deviation from the Temperature Size Rule), identifying a cold stress effect. High amplitude variation and low amplitude variation around a mean temperature of 12°C impacted reproduction and body size compared to nematodes kept at the matching average constant temperatures. Copper exposure affected reproduction, body size and life-span and consequently population growth. Sensitivity to copper (EC50 values), was similar at intermediate temperatures (12, 16, 20°C) and higher at 24°C and especially the innately stressful 8°C condition. Temperature variation did not increase copper sensitivity. Indeed under variable conditions including time at the stressful 8°C condition, sensitivity was reduced. DEBkiss identified increased maintenance costs and increased assimilation as possible mechanisms for cold and higher copper concentration effects. Model analysis of combined variable temperature effects, however, demonstrated no additional joint stressor response. Hence, concerns that exposure to temperature fluctuations may sensitise species to co-stressor effects seem unfounded in this case.

Citation: Cedergreen N, Nørhave NJ, Svendsen C, Spurgeon DJ (2016) Variable Temperature Stress in the Nematode Caenorhabditis elegans (Maupas) and Its Implications for Sensitivity to an Additional Chemical Stressor. PLoS ONE 11(1): e0140277. https://doi.org/10.1371/journal.pone.0140277

Editor: Jonathan H. Freedman, NIEHS/NIH, UNITED STATES

Received: June 8, 2015; Accepted: September 22, 2015; Published: January 19, 2016

Copyright: © 2016 Cedergreen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All data have been uploaded to the NERC Environmental Information Data Centre (doi:10.5285/af125e27-3b70-4f0a-81fb-a7eb10f64ef7). Cedergreen, N.; Nørhave, N. J.; Svendsen, C.; Spurgeon, D. J. (2015). Temperature and copper effects on the nematode Caenorhabditis elegans. NERC Environmental Information Data Centre. doi:10.5285/af125e27-3b70-4f0a-81fb-a7eb10f64ef7.

Funding: DJS and CS were supported by core NERC national capability funding under the Centre for Ecology and Hydrology Pollution and Ecological Risk science theme. NJN was supported by the Research School of Environmental Chemistry, Microbiology and Toxicology in Denmark and a PhD grant given jointly by the University of Copenhagen and Centre for Ecology and Hydrology National Capability funding.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Organisms living in natural environments are regularly exposed to sub-optimal conditions leading to physiological stress that can affect life history parameters, population growth rates and ultimately carrying capacities [1,2]. Studies to describe and understand the mechanisms that link environmental stress to organism biology have to date largely focussed on single continuous stressors. While such simple scenarios may be amenable to generating mechanistic understanding, such cases rarely reflect the field situation [3]. Under the variable stressor conditions that may dominate in the field, the ability of any species to survive may be determined not only by the average conditions that prevail, but also by the amplitude and frequency of change. Indeed, variable environmental conditions may in themselves be stressful. For example, remodelling of membrane lipid composition to maintain appropriate fluidity at specific temperatures [4,5,6], or adaptation of enzyme and transporter compositions to cope with specific chemical conditions [7], can be physiologically costly under rapidly fluctuation. This can result in additional costs for adaptation compared to organisms living under more stable conditions. Such effects may affect physiology directly, or may compromise the ability of organisms to cope with the presence of other stressors in multiple stressor scenarios [8,9,10].

Temperatures can vary geographically and maximum and minimum temperature ranges can fluctuate diurnally by more than 20°C and annually by up to 50°C at any given location. Hence, for resident species, exposure to variable and potentially stressful temperatures is likely to be a reality at least some of the time. The effect of these variations on species may be both direct and indirect; the latter, for example, by changing the way that species cope with a second stressor. To date studies concerning variable temperatures and other stressors have often focussed on interactions with pathogens [11]. Few have so far looked at how variable temperature changes the consequences of chemical exposure. Indeed while many studies have looked at how constant temperature affects toxicity, most frequently noting increasing severity of effect with increasing temperature [3,12,13,14], comparisons of toxicity under constant and variable temperatures are lacking. In a few cases, results from laboratory toxicity tests conducted under constant conditions have been compared to observed toxicity in field studies [15,16]. However in such studies, temperature is one of a range of variable factors considered, and any observed differences cannot be attributed to this single variable.

If there is a physiological cost of optimising biochemistry to changing temperatures, then both a direct effect of variable temperature on a trait response, and also a changed response following exposure to a second stressor may be expected. Here we report the results of a series of toxicity tests for copper with the nematode Caenorhabditis elegans conducted at different constant temperatures in the range 8 to 24°C and for variable temperatures on two amplitudes; ± 4°C around mean 12, 16 and 20°C and ± 8°C around mean 16°C. C. elegans is particularly suited to study work because it has a short-life cycle allowing easy measurement of multiple traits. This biology involves parthenogenetic, temperature independent, productions of a total broodsize of 200–300 eggs during its lifespan, with this maximum number widely attributed to sperm limitation at least under laboratory conditions [17]. Laid eggs hatch to produce larvae that progress through four stages before reproductive age is reached. After reproduction has ended, nematode keeps may live for several days. Hence traits that can be measured included growth, reproduction (brood size) and survival to jointly yield an estimate of population growth rate (PGR). Further, a version of the dynamic energy budget model (DEBkiss) can be used to provide mechanistic understanding of energy fluxes governing changes in growth and reproduction under these different conditions [18]. Four hypotheses could be specifically tested through this approach from our analyses. First, if none of the constant temperature regimes are in themselves stressful, then the observed differences in life history at different temperatures will be simply due to temperature induced changes in chemical rates, described by the Arrhenius equation. To test this hypothesis DEBkiss was used and a temperature factor was applied to the energy fluxes [18]. Second, if a significant proportion of energy is spent adjusting to variable temperature regimes, we expect a decrease in growth and reproduction under fluctuating temperature that will be greater at the large temperature variation amplitude. Third, we may expect Cu toxicity to increase at higher temperatures. This has been shown for nematodes in a previous study (Nørhave et al, 2014), and was proposed to be due either to a higher Cu uptake at high concentrations or to a more severe effect of the Reactive Oxygen Species (ROS) catalysed by Cu. In DEBkiss this mechanisms may be linked to either a reduction in energy acquisition or an increase in maintenance costs. Fourth, interactions between temperature, both constant and variable, and chemical stress follow the model of Independent Action as described by Holmstrup et al. [3]. In DEBkiss this independence will be characterised by model optimal predictions, based on model parameters derived from DEBkiss fits for individual treatments alone.

Materials and Methods

Culture conditions for test organism

C. elegans (N2 Bristol strain) initially from the C. elegans Genetics Centre, University of Minnesota) were cultivated in darkness at 20°C on nematode growth medium (NGM) agar plates at all time with an excess of Escherichia coli of uracil deficient strain OP50 [19]. At least 2 weeks before the start of all experiments new cultures were initiated at the selected test temperatures to acclimatize the worms through several generations. Cultures were maintained by transferring a chunk of agar and associated eggs, juvenile and adult worms from an existing culture to freshly prepared NGM agar plates weekly.

Set-up of toxicity tests

The temperature range selected for study were relevant to the ecological range of C. elegans [20] and physiological tolerance limits relating to stress on growth and reproductive outputs [21]. Specifically exposures were conducted at steady temperatures of 8, 12, 16, 20 and 24°C and variable temperature scenarios with daily fluctuations of ± 4°C in the ranges 8–16°C (12°C average), 12–20°C (16°C average) and 16–24°C (20°C average) and daily fluctuations of ± 8°C in the ranges 8–24°C (16°C average). Temperatures were changed with a rate of 4°C per hour to simulate changes at two points within a 24 hr period separating two stable temperature periods to provide an simulation of a daily temperature regime under condition of high diurnal variation (Fig 1). Separate toxicity tests with copper were conducted under all temperature regimes using concentration known to result in changes in measured traits [22]. Copper concentrations used were 0, 1, 3, 8, 20 and 40 mg Cu/L agar. The highest concentrations equate to levels 2–3 orders of magnitude greater than those found in surface waters and soil pore waters from uncontaminated sites, but are comparable to soil pore water concentrations in heavily polluted soils.

thumbnail
Download:
Fig 1. The temperature treatments depicted as temperature as a function of time over 24h.

Five constant temperatures (8, 12, 16, 18 and 20°C, solid lines), three temperature ranges with an amplitude of ± 4° around average temperatures of 12°C (black circles, black line), 16°C (grey circle, grey line) and 20°C (white circle, broken line) and one with an amplitude of ± 8°C from the average (black triangles, dotted line) were used. Temperatures were shifted with 4°C per hour. All treatments were run with the minimum and maximum temperature around 12 a.m. and 12 p.m., respectively.

https://doi.org/10.1371/journal.pone.0140277.g001

A stock solution of CuCl2 in demineralized water (2 g Cu/L) was made and used throughout the experiment. This copper stock solution was added to the NGM agar while liquid and mixed. Batches of NGM were produced every five days during the course of the experiment to minimize immobilization of copper in the agar, as shown for cadmium by Álvarez et al. [23]. All tests were conducted using synchronized cohorts of nematodes produced from adults selected from the main stock cultures acclimatized to each temperature scenario. To generate these cohorts, adult hermaphrodite worms were selected from the acclimatized cultures and placed on Petri dishes with copper spiked NGM and E. coli and left at the test temperatures for 4 hours (6 hours for 8°C to ensure sufficient eggs were laid). After this time, the adults were removed and the eggs laid in the 4 (6) hour period left to hatch. The offspring (the test organisms), once hatched and grown to the L4 larval stage on the copper spiked plates, were transferred to 12-well plates with copper spiked NGM and E. coli, with one individual per well. During the test the presence of an E. coli lawn was always observed at all treatment levels, indicating the provision of food in excess. There were 12 replicate worms for each treatment, except the controls of 24°C constant, 8–16°C, 8–24°C and 16–24°C which had 36 replicates, and the 8°C test in which we could not always obtain 12 worms for each treatment, hence, between 6 and 12 worms were used.

During the reproduction phase of the test, the nematodes were moved to new wells in fresh plates every day. These transfers and the associated observations allowed calculation of the time taken until first egg was laid and subsequently daily egg production. The daily transfer and counting were repeated through both the reproductive and senescence stages until death (determined by the lack of response to being probed). Fertile eggs and hatched juveniles were counted as offspring, while visibly infertile eggs were excluded from brood size counts. At regular intervals, nematode length was measured using a Nikon DS-Fi1 camera connected to a Nikon SMZ 800 stereomicroscope. Body length was measured using the program Nikon NIS Elements Imaging Software 3.2 to provide a description of growth pattern.

Copper analysis

The copper content in the water fraction of the agar was determined one and five days after production. This allowed us to distinguish dissolved copper and agar bound copper (which may not have been fully bioavailable). A study on C. elegans focussing on cadmium exposure, has investigated routes of uptake (water versus bacterial food)[24]. This worl concluded that toxic effect correlates best with the Cd concentration in the aqueous phase. As Cu was dosed into the agar as an aqueous solution, we assumed that aqueous Cu concentrations would also be the best measure of exposure in the present study and, hence, analysed this fraction from the agar. Samples of NGM from each produced batch were saved in centrifuge tubes and, after one or five days, centrifuged at 10000 G for five minutes to separate water and agar. A sample of the water fraction was removed and acidified. The analyses were performed on a graphite furnace AAS (Perkin Elmer Zeeman 5100, Waltham, MA) as described in Cedergreen et al, 2013. Blanks, spiked reference samples and a standard reference material NIST 1577c (bovine liver) were included in the analysis to confirm validity [13]. Internal Cu concentrations of nematodes could not be measured, as their small biomass prevented reliable detection.

Statistics

Differences between temperature treatments for brood size, lifespan and maximal length of non-Cu exposed nematodes were tested by a one-way ANOVA with a Tukey post hoc analysis using the R statistical software version 2.12.0 (http://www.r-project.org/).

The average cumulative egg production as a function of time for each treatment was described by a three parameter log-logistic sigmoid model (Eq 1) using the statistical software R version 2.12.0 within the DRC package [25]. (1) where d is the maximum number of eggs produced, e is the time when half the maximum number of eggs is produced, b is proportional to the slope of the sigmoid curve at e, y is the accumulated number of produced eggs and x is the time. From this model the estimated time to first egg (TFE) was calculated by solving x for y equal to 1/d. The DRC package returns a standard error with the estimate.

A stage based matrix population model was constructed for each treatment, based on the daily reproduction values and lifespan. The modelling was performed with the statistical software R version 2.12.0 (http://www.r-project.org/), using the popbio package [26], which is based on the work of Caswell [27]. The models were day staged, thus, making full use of the available reproduction data. A projection of the population over 500 iterations (days) were made and after checking that the population growth had stabilized, the population growth rates (PGR) were determined, with the 95% confidence intervals (CIs) determined by bootstrapping.

To define the effect of Cu on brood size, lifespan, maximum length and PGR, these responses were described with a log-logistic three parameter concentration response model (Eq 1). To assess the joint effect of a potential decrease in any of these parameter as a result of variable temperature and the effect of Cu, the model of independent action was used [28]. The model predicts that if one stressor reduces the endpoint in question to fraction of the unstressed control (Rs1), and the second stressor affects the endpoint to another fraction of the control (Rs2), then the unaffected fraction left (Rmix) when affected by both stressors, is the product of the two unaffected fractions as: (2)

The DEBkiss Model and stress interactions

The structure of the DEBkiss model is described in detail in Jager et al [18], with model application to address the effects stressor combinations in C. elegans (cadmium and fluoranthene toxicity) further elaborated by Jager et al. [29]. The model equations are described in Table 1 and the model parameters in Table 2. Basically, food is assimilated by the maximal rate of JA, which scales with the surface area of the organism at a particular given temperature [30]. The assimilated energy is divided in a proportion allocated to maintenance (JM) and growth (JV), jointly κ, while the remaining energy (1- κ) goes into reproduction (JR) when the organism has reached its length of puberty (Lp) (S1 Fig). Both the rate of assimilation (JA) and maintenance (JM), and indirectly the rate of growth, as JV = κJA-JM, and reproduction as JV = (1-κ)JA, will be affected by temperature through the general effect of temperature on chemical reactions governed by the Arrhenius equation. The Arrhenius equation states that the chemical rate k can be described by: (3) Where A is a pre-factor, Ea is the activation energy of a chemical reaction, R is the gas constant and T is the absolute temperature given in °K [31]. If the energy fluxes of the DEBkiss model are determined at a standard temperature, as for example 20°C (293.15°K), then the change of the rate of the fluxes (JA and JM) with temperature, can be estimated by multiplying them with a temperature factor. The factor (FT), given by the difference between the rate at the reference temperature kref and the ambient temperature ka, is given by: (4) Tref is the reference temperature given in °K and Ta is the ambient temperature given in °K. Ea(1/R) is a constant named the Arrhenius temperature (TA) in a DEB-context, and typically lies within 5000 and 13,000 K [32]. To describe the data from the five constant temperatures, first robust parameter estimates were obtained on the 20°C stable temperature treatment without Cu. As the experiments were run in parallel with those published by Jager et al (2014), and the 20°C control treatments were overlapping, the parameters estimated by Jager et al (2014) were used as starting values. Contrary to Jager et al (2014) we did not include an initial food limitation, as we did not have length measurements at the initial growth stages to support parametrisation. If, however, we assumed Van Bertalanffy growth started one day after hatching, growth data was very well described without the initial food limitation. When robust parameters were obtained for the 20°C control data, the 8, 12, 16 and 24°C data was included, and the temperature factor FT was multiplied with JaAm and JvM to determine the approximate size of the Arrhenius temperature (TA). Finally, all parameters were optimised jointly using a Monte Carlo Markov Chain optimisation procedure.

thumbnail
Download:
Table 1. Equations for the basic DEBkiss model.

https://doi.org/10.1371/journal.pone.0140277.t001

thumbnail
Download:
Table 2. Parameters of the DEBkiss model as used in this study.

The fixed parameters were retrieved from Jager et al (2014) and were the same for all parameter fits. The fitted parameters are given for the three scenarios: 1) All constant temperature data (minus 8°C) are described jointly using the Arrhenius temperature and 20°C as a reference, 2) The constant 16°C and the variable 16±8°C treatments are described together applying a stress factor (s) on the somatic maintenance rate and 3) The 16°C treatment is described in the presence of the three highest Cu concentrations applying a Cu-related stress factor (sCu, Table 1) on maximum assimilation rate (JaAm) and length at puberty (Lp). Parameter estimates are given ± s.d.

https://doi.org/10.1371/journal.pone.0140277.t002

To test whether the stress induced by the high amplitude fluctuations mainly affected assimilation or the cost of maintainance, growth or reproduction, a general stress parameter, s, was applied to either JaAm, JvM, yVA or yBA as described in Jager and Zimmer [34]. First parameters were optimised for the 16°C constant treatment, and then both constant and the 16 ± 8°C treatment were fitted together including the stress parameter on one of the above variables and optimised jointly using a Monte Carlo Markov Chain optimisation procedure.

The effect of Cu on growth and reproduction was tested on the 16°C constant temperature treatment using the three highest Cu-treatments: 8, 20 and 40 mg Cu L-1. The 1 and 3 mg Cu L-1 treatments were omitted as they induced an increase in reproduction, which the DEBkiss model cannot account for as it assumes a negative effect of an applied chemical. To describe the effect of Cu, a Cu-related stress parameter, sCu, was applied first to all the main DEB parameters (JaAm, JvM, yVA or yBA), none of which described both growth and reproduction really well. Then, as it was observed that the size of the nematodes at the time where they laid their first eggs decreased with increasing Cu-concentrations (which is contrary DEB-assumptions of constant size at first reproduction) [34,35], the stress factor was also applied to the parameter: length a puberty (Lp) by multiplying it with 1-sCu. To predict the joint effect of both variable temperature and Cu-stress, the Cu-model was run including the stress factor obtained for the variable stress, and the predictions were assess together with the data from the multiple stress experiment. All modelling was done using OpenModel 2.2.0.

Results

Chemical analyses

All measured Cu concentrations were within the range of 81–110% of nominal concentrations with an average of 95 ± 5% (mean ± s.d., n = 51). This indicates a high level of agreement between nominal and measured concentration and hence hereafter nominal concentrations are used for all calculations and discussions.

Are the constant temperatures in themselves stressful?

Brood sizes were similar at 12, 16 and 24°C even though increasing temperature shortened brood period, while nematodes in the 8°C treatment produced a broodsize half those of the any temperature treatment (Fig 2A). This indicates a direct chill stress effect, but no heat stress response in the temperature range considered. The 20°C temperature was ultimately optimal, with nematodes producing on average approximately 20% more offspring than in the 12, 16 and 24°C treatments (Fig 2A). Despite different growth rates, there was only a slight, non-significant, negative relationship between temperature and size at the end of their lifetimes at all constant temperatures in the 12–16°C range. Length were, however, reduced at 8°C contrary to the expectation from the Temperature Size Relationship, further supporting indication of a chill stress (Fig 2B). Both time to first egg (TFE) and lifespan decreased with increasing temperature, making population growth rate (PGR) linearly related with temperature (Fig 2C, 2D and 2E). DEBkiss well described growth and reproduction at 20°C (R2 = 0.95). Including the temperature factor, FT, and describing all growth data with the 20°C parameters optimising TA gave a poorer fit (R2 = 0.63), with the 8°C treatment clearly lying below the fitted curves (Fig 3A and 3B). Omitting the 8°C data still using the 20°C parameters increased the R2 to 0.83, supporting the indication from the single traits of broodsize and body size that the 8°C was stressful causing deviation from expectation from the Arrhenius relationship alone. All four other constant temperatures accorded well with model predictions (Fig 3A and 3B). All parameters were thereafter optimised for 12, 16, 20 and 24°C constant temperature using a Monte Carlo Markov Chain optimisation (R2 = 0.87). The resulting fits and parameter values are given in Fig 3A and 3B and Table 2.

thumbnail
Download:
Fig 2. Temperature effect on five different endpoints: final brood size (A), Final body length (B), time to first egg (C), lifespan (D) and Population Growth Rate (PGR) (E).

All data are given as a function of the mean temperature of the treatment. Data are given as mean ± s.e.m., apart from PGR which is given with 95% Confidence Intervals (CI) obtained by bootstrapping. Significantly different treatments (ANOVA followed by a Tukey post hoc test) are denoted by different letters. Constant temperature treatments are given in black symbols, treatments varying ± 4°C are given in grey symbols and treatments varying ± 8°C are given in open symbols.

https://doi.org/10.1371/journal.pone.0140277.g002

thumbnail
Download:
Fig 3. Temperature effects on traits fitted using the DEBkiss model.

The development of body length (A) and offspring production (B) as a function of time for the five constant temperature treatments are described by the DEBkiss model including a temperature function (R2 = 0.87). Fig C and D shows bodylength and offspring production as a function of time for the constant 16°C treatment and the variable 16 ± 8°C treatment described by the DEBkiss model including an 8% increase in mainainance cost for the variable treatment (R2 = 0.87). Data are given as mean ± s.e.m. Model parameters are shown in Table 2.

https://doi.org/10.1371/journal.pone.0140277.g003

Variable versus constant temperatures

For the nematodes kept under low amplitude (± 4°C) variable temperatures for the two highest average temperatures of 16 ± 4°C and 20 ± 4°C, changes observed in traits were generally positive in direction, but were not significant between the variable and constant exposure for any endpoints (Tukey: p < 0.05, Fig 2A–2D). For the low amplitude (± 4°C) variable at the low temperature range of 12 ± 4°C, both lifespan and maximal length was significantly reduced compared to the constant 12°C temperature treatment (Tukey: p < 0.05, Fig 2A and 2D), while TFE and broodsize were not affected (Tukey: p > 0.05, Fig 2B and 2C). The high amplitude treatment (± 8°C) at 16 ± 8°C, however, showed a significant reduction in maximal length (15%) and brood size (30%) compared to the 16°C constant temperature treatment. Together with a 15% longer TFE this resulted in a reduction in population growth rates of 5% in the 16 ± 8°C compared to 16°C constant treatment from 2.04 (2.02–2.07) to 1.93 (1.91–1.95)(Fig 2E).

Testing whether an increase in maintenance cost could explain the reduced growth and reproduction using the DEBkiss model, we first fitted the constant 16°C data (R2 = 0.95) and then added the dataset from the 16 ± 8°C treatment and a stress factor on the parameter: specific somatic maintenance rate (JvM). This model described the observed patterns of body size and offspring production in time well (R2 = 0.87) if JvM was increased by approximately 8% (See parameter values in Table 2, Fig 3C and 3D). Decreasing the maximum assimilation rate (JaAm) by approximately 6% or increasing the cost of growth given by yVA by 14%, however, gave equally good fits (R2 = 0.87). Hence, while increased maintenance costs could provide an explanation for the observed effects in accordance with our initial hypothesis, this mechanism could not alone be distinguished as the most important process governing the stress effects of the high variable temperature 16 ± 8°C treatment on traits.

Cu stress at different temperatures

Copper exposure had a severe effect on growth, reproduction and lifespan mainly at the two highest tested concentrations of 20 and 40 mg Cu L-1, while the two low concentrations 1 and 3 mg Cu L-1 no effect, or even, in the case of offspring production, a slight increase in trait values was seen (see Fig 4A, 4B and 4D for example data set for the 16°C treatment). No effect on any concentration on TFE was observed (Fig 4C). Integration of all measured effects to assess Cu effect on PGR indicated a reduction also at the two highest tested concentrations (Fig 4E). Reproduction, final body length, lifespan and PGR could be described with concentration-response curves for all studies under the different temperature regimes (Table 3). Consistent with our hypothesis, increasing temperature to constant 24°C resulted in a lower EC50 for all responsive traits than at 12, 16 and 20°C, consistent with previous work and expectations [22]. The largest shift in sensitivity was seen for PGR for which the EC50 reduced from 61.2 mg Cu L-1 at 12°C to 25.5 mg Cu L-1 at 20°C. Additional inclusion of a treatment at 8°C, which was not included in the study of Norhave et al. [22], identified an increased Cu sensitivity of nematodes kept at this lower temperature, most notably for broodsize (Table 3). This sensitivity was not consistent with our initial hypothesis of increased sensitivity with increasing temperature.

thumbnail
Download:
Fig 4. The effects of Cu exposure on traits responses in C. elegans in treatments with an average temperature of 16°C.

The five different endpoints: Final brood size (A), final body length (A), time to first egg (B), lifespan (D) and Population Growth Rate (PGR) (E) for the constant 16°C treatment (filled symbols), the 16 ± 4°C (grey symbols) and the 16 ± 8°C (open symbols) as a function of Cu concentrations in the agar. Data are given as mean ± s.e.m. and are described with a three parameter log-logistic concentration response model, except for TFE. The parameters are given in Table 3, together with the concentration-response parameters of the other temperature treatments.

https://doi.org/10.1371/journal.pone.0140277.g004

thumbnail
Download:
Table 3. Dose-response parameters (Eq 1) for the measured endpoints: Length increase (Maximum length minus length at hatching), lifespan, brood size and population growth rate (PGR) for the different temperature treatments.

The curves for the three treatments with an average temperature of 16°C are shown in Fig 3. Data for the 12°C and 24°C average temperature treatments are presented in S2 and S3 Figs. The parameters are given ±SE.

https://doi.org/10.1371/journal.pone.0140277.t003

Describing the effect of Cu using the DEBkiss approach, gave good fits of the responses of body length and offspring production in time by applying a Cu-related stress parameter on on maximum assimilation rate (JaAm) and length at puberty (Lp) (R2 = 0.97) for the 16°C treatment (Fig 5A and 5B). Other combinations of stress parameters could potentially also give good fits, but a Cu related stress on energy flowing into the organism is consistent with expectations driven by Cu affects on metabolic systems and ROS affect on physiological processes including the mitochondria [36]. Size at TFE at 16°C constant temperature was observed to decrease from 0.85 mm in non-exposed nematodes to 0.77, 0.70 and 0.60 mm in nematodes exposed to 8, 20 and 40 mg Cu L-1 (using the TFE-values for the 16°C constant treatment of Fig 3B and the growth predictions of Fig 5A). These result support identification of a direct impact of copper stress on this trait.

thumbnail
Download:
Fig 5. Copper exposure effects on traits fitted using the DEBkiss model.

The development of body length (A) and offspring production (B) as a function of time for the 16°C constant treatment combined with 0, 8, 20 and 40 mg Cu L-1 described by the DEBkiss model. Fits include a Cu-related stress factor (sCu, Table 1) on maximum assimilation rate (JaAm) and length at puberty (Lp) (R2 = 0.97). A mixed stressor prediction including the 8% increase in somatic maintainance obtained when fitting the constant and variable control treatments together, and the Cu-related stress factor ontained from the constant temperature treatment (A, B) is shown together with the body length (C) and offspring production (D) data from the 16 ± 8°C treatment combined with 0, 8, 20 and 40 mg Cu L-1. Jointly the mixed stressor prediction describes 75% of the variation in the data. The 1 and 3 mg Cu L-1 data are omitted from the fit, as they increased offspring production, which the DEBkiss model is not parameterised to deal with. Data are presented as mean ± s.e.m.

https://doi.org/10.1371/journal.pone.0140277.g005

Interactions between variable temperature regime and Cu stress

To assess how variable temperature regimes affect sensitivity to Cu exposure, we explicitly compared responses for the three treatments with an average temperature of 16°C. Trait responses to Cu exposure in the constant 16°C and low amplitude variation 16 ± 4°C temperatures were remarkably similar for all traits, with only life-span showing difference at the lower exposure levels (longer life-span at low amplitude temperature variation). The high amplitude variable temperature (16 ± 8°C) and the 12 ± 4°C treatment showed lower trait values at 1, 3, 8 and 20 mg Cu L-1 treatments, while the 16 ± 4°C and 20 ± 4°C did not show such differences across multiple traits. Further in all treatments, difference in absolute values were not evident for the 40 mg Cu L-1 treatment (see Fig 4 for 16 ± 8°C and 16 ± 4°C, S2 Fig for 12 ± 4°C, S3 Fig for 20 ± 4°C). Applying a traditional independent action prediction for the two types of stress: variable temperature and Cu, indicated antagony for the high concentration for length, brood size and PGR (Fig 4). The evaluation of the prediction, however, was limited due to the relative few Cu concentrations producing a significant adverse effect on traits. Nonetheless, the suggested antagonism is further indicated by the EC50s for the variable temperature treatments which show higher values than any constant temperature treatment, notably for the 16 ± 8°C and also 12 ± 4°C treatments (Table 3).

The effect of variable temperature on Cu toxicity was further investigated within DEBkiss. Predicting the effect of the mixed stress of variable temperature (± 8°C) and Cu, the increase in maintenance costs obtained from the fit of the constant versus variable 16°C treatment treatments was simply added to the Cu-model to predict the comparative effect of the metal at the constant versus high amplitude variable treatment. The integration of the maintenance effect derived from high amplitude variation with the Cu effect resulted in a model fit that was able to explain 75% of the variation in the data (Fig 5C and 5D). Body length values in the 8 and 20 mg Cu L-1 treatment generally fell below the fitted line, while offspring production values at time-points in the 8 mg Cu L-1 treatment generally fell above the fitted line. These effects are consistent with the indication of antagonism between stressors that causes deviation of effects from strict independent action.

Discussion

Effects of constant temperature regime on toxicity

Constant temperature changed measured C. elegans life-cycle traits in a manner consistent with expectations from the Arrhenius relationship in the 12°C—24°C temperature range, but not for nematodes at 8°C for which DEBkiss suggested a direct stress effect. For example, maximal length of nematodes shows a slight decrease with increasing temperature in the range 12°C—24°C in a manner consistent with predictions from the Temperature Size Rule. Nematodes growing at 8°C were, however, smaller than individuals grown at higher temperatures indicating of a direct stress effect of cold exposure for this trait that is inconsistent with expectations (Fig 2B). This result questions the generality of the Temperature Size Rule for all ectotherms in a manner consistent with finding from a range of recent studies [37,38]. The linear relationship of PGR with temperature reflects the high sensitivity of this parameter to traits that were strongly positively affected by increased temperature, such at TFE, compared to traits negatively affected by temperature, such as lifespan, to which PGR in C. elegans is known to be insensitive [39].

The sensitivity to cold at 8°C in a nematode species of such a wide geographical range as C. elegans [20] observed in results from this study could be related to its native growth environment. Currently characterized populations, such as those known in France [40], demonstrate that C. elegans prefers nutrient- and microorganism-rich substrates such as rotting fruits and decomposing plant matter [41,42]. These microbial rich environments are the subject of significant warming from heat derived as a by-product of their active microbial communities. Observation of cold, but not heat (up to 24°C) stress are, thus, consistent with expectation based on the known ecology of the species.

That the warmest temperature of 24°C did not in itself prove stressful was itself quite surprising as this temperature is close to temperatures that have previously been shown to affect traits, such as reduced brood size at 25°C [22] and growth, reproduction and survival at 25.5°C [21]. Transition from optimal temperature to stressful temperature at the higher range, therefore, seems to occur within a small temperature range for the species. One practical consequence of our observation is the value of maintaining C. elegans population at temperatures in the range 12–20°C if the induction of direct temperature induced stress into nematode cohorts used for phenotypic studies is to be avoided. In this respect the implications of the practice of refrigeration storage for phenotypic measurement in later experimental studies may require investigation.

Effects of variable temperature on traits

Low temperature variations resulted in little change in trait values when compared to worms held under similar average, but constant conditions. This does not support a positive effects on traits related to stress hardening [43], but rather an absence of a variable thermal condition effect on fitness [44]. A direct stress effect of temperature variation on traits was found in two treatments, namely the high amplitude variation 16 ± 8°C and low amplitude variation at the average 12 ± 4°C. The high amplitude variation treatment reduced brood size by 30%, however, despite the magnitude of this effect, PGR was only decreased by 5% (which corresponds to a change of constant temperature from 16°C to 14.2°C based on a regression of PGR as a function of constant temperature, Fig 2E). The relatively low sensitivity of PGR to effects traits such as brood size is consistent with what is known from previous sensitivity analyses of trait relationships to PGR in C. elegans [23,45].

In was notable that the two variable treatments that showed a direct stress effect included a period where the nematodes were kept at 8°C—a constant treatment that in itself proved stressful. Our hypothesis was that any stress effects of variable temperature would be driven by increased maintenance costs accrued as individuals seek to adapt physiology to the changing temperature conditions. DEBkiss analysis confirmed the possible validity of this hypothesis, however, decreasing assimilation rates or increasing growth costs within DEBkiss could also explain the observed patterns of effect in time with almost a similar level of agreement. Hence, for this case, the DEBkiss model could not unequivocally distinguish between physiological modes of action. Further, and perhaps more fundamentally, that only some variable temperatures (those including time spent at 8°C) affected trait performance point not to a variable temperature effects, but a (maintenance) stress effect just attributable to low temperature exposure

Temperature exposures in poikilothermic organisms are known to have a range of effects on physiology. The obvious relationship is with metabolic rate, which is fundamentally behind the changes observed for some traits (notable TFE and life-span) and PGR as observed in this study. Physiological heat shock can cause the induction of a series of stress response mechanisms, such as the metallothionein, antioxidant defence mechanism and mitogen-activated protein kinase signalling, as well as the heat shock system for protein chaperoning [46]. Cold stress also has a physiological consequence for species. Rapid transfer from warmer to colder temperatures has been shown to stimulate changes in the lipid compliment of species as they attempt to maintain optimal membrane fluidity at different temperatures [4,47]. This remodelling of phospholipid content is naturally subject to time lag during the time when the require lipids are produced in the endoplasmic reticulum and Golgi. Further, during these times of rapid change, membranes may exist for a time in states outside of optimal for the currently experienced temperature conditions. While this has the potential to place organisms under stress, in this case the impacts of these changes on traits actually appear limited. It is only when species have to exist outside their physiologically normal range that stress leading to reduced trait performance results. Variable temperature exposure alone does not appear to be stressful.

Relationship between Cu toxicity and constant temperature regime

Testing the response of C. elegans to Cu-stress at the different constant temperatures partly confirmed the study of Nørhave et al [22], who found a 1.5 to 2.5-fold decrease in Cu EC50 increasing temperatures from 11 to 24°C, and in addition found stimulations of reproduction and PGA at low Cu concentrations indicative of hormesis. However, going beyond Norhave et al. [22], the inclusion of a lower temperature treatment at 8°C (as compared to the lowest temperature treatment of 11°C used by [22]) showed an increase in sensitivity to Cu at low temperature, based on the EC50 values for the different endpoints (Table 3). Hence, an innately stressful temperature regime was shown to increase the sensitivity of C. elegans to a toxicant. This observation of an interaction between cold stress and chemical exposure has been indicated in a series of studies looking both at how cold exposure affects toxicity and also at how chemical exposure affect cold tolerance [3]. These included specific studies that have looked at these interaction for copper and low temperatures [48,49,50].

DEBkiss modelling of the data for copper exposed nematodes pointed to a number of possible effects on model parameter that could produce patterns of growth and reproduction consistent with observations. Applying a Cu-related stress parameter on maximum assimilation rate (JaAm) was one such physiological mechanism. Studies of the physiological mode of action of metals in nematodes have pointed to a mechanism that is associated with an effect of the metal on resources acquisition [51,52]. Hence, incorporation of a Cu related stress on energy assimilation is consistent with expectations of metabolic effects on metabolic systems and ROS effects on physiological processes. Copper is recognised for the particular interaction that this metal has with the mitochondria with potential effects on availability of energy through metabolism [36,53].

Effect of variable temperature regime on Cu toxicity

Our initial hypothesis was that stress related to variable temperature exposure may result in a sensitisation of the nematodes to secondary Cu exposure. However, contrary to expectation, nematodes under the high amplitude variable temperature (16 ± 8°C) actually seem better able to endure exposure to higher Cu concentration than expected by both the independent action and mixed stressor DEBkiss predictions (Figs 4 and 5). Such an antagonistic interaction points to a potential stress hardening effects resulting most likely from exposure to cold. “Stress hardening” effects by which exposure to one stressor may induce mechanisms that protect against the effect of a second stressor have been widely reported for a range of organisms and stressor combinations [9,54,55]. The increased robustness of nematodes exposed to variable conditions including cold stress was also partly supported in the 12 ± 4°C, for which EC50-values are higher than those for all other treatments except the high amplitude variable temperature for all traits except brood size (Table 3).

Use of the simple Dynamic Energy Budget model, DEBkiss, allowed a systems-based assessment of the interactions of different stressors with physiologically relevant parameters for both temperature and Cu stress effects. For the exposure to constant and variable temperatures that include a time involving exposure at 8°C, we could not identify one process that explained the decrease in growth and reproduction significantly better than other possible mechanisms. This could be attributed to the low resolution of our data with length monitored only a few points in time, or it could be that the observed temperature effects is driven by some interactions including multiple physiological parameters. Thus, though there are some chemicals and possibly also other stress factors that affect some processes very specifically [51,56], other types of stressors may affect many of the parameters included in DEBkiss simultaneously [57]. Low/variable temperature stress, potentially affecting all enzyme and lipid related processes may be one such general stressor, in contrast to the more specific effect of Cu on assimilation indicated from DEBkiss. We achieved the best combined fit of growth and reproduction by also applying the chemical stress factor to the parameter, length at puberty. This parameter is not one of the traditional parameters for application of stress in a DEB-context [30], where all organisms are expected to start reproducing at a specific size [58]. Both we, and other studies, have shown that organisms ranging from nematodes, daphnids and rotifers decrease their size at first reproduction when exposed to stress as different as starvation, cadmium and copper stress [14,59,60]. Here we do not find the DEB-assumption of constant size at the time of reproduction biologically correct in all cases and, therefore, recommend including size at puberty as one of the parameters possibly affected by stress.

It was relatively easy to combine stressors, making reliable time-dynamic mixed stress predictions using DEBkiss, as also shown in Jager et al. [29]. This identifies DEBkiss as a promising tool for research in combined stressor ecology, particularly those studies involving temperature and, therefore, the rates of physiological processes and their effects on traits over time. Challenges in the implementation of DEBkiss that are recognised include the increased amount of work involved in measuring endpoints continuously over time and the fact that measurements are sometimes in themselves stressful (for example quantifying broodsize in time for nematodes involves regular transfers between newly treated wells). As identified in this study, it may not always be possible to unequivocally identify the physiological changes on energetic parameters that underpin the changes of traits with time. However, this may be a consequence of the underlying biology of the stressor rather than a weakness of the modelling approach.

Ecological and chemical effect assessment implications

Previous studies that have developed physiological based models to understand the link between changes in resource allocation under (toxicant) stress have provided a set of expected phenotypes for animal subject to different energy allocation related processes and trait effects [61,62,63]. Direct stressor effects on traits could be related to effects on physiological parameters, with low/variable temperature effects driven by the individual stressor and the interaction showing an antagonism. In the field, environmental conditions can be subject to rapid variations that are not included in the laboratory based assessment of toxicity that dominate in the ecotoxicological literature. On the basis of the data generated here, it appears that such toxicity data at constant temperatures will provide a good assessment of effects under variable conditions in those cases where the range of the environmental variable does not include conditions beyond the physiological limits of tolerance of the species. Hence, rather than being innately stressful in its own right, temperature fluctuation is merely a vehicle that may on occasion result in exposure to a stressful condition. When temperatures remain constant and in the normal physiological range, the impacts of Cu toxicity appear to be marginal, with toxicity only affected as temperature approach or exceed physiological limits. Fluctuations, even when they may result in some stress exposure themselves, do not necessarily result in sensitisation. Indeed, as seen here a version of a “stress hardening” response may also occur [64,65]. Hence, consideration of fluctuating temperature does not in itself need a rethink of our approach to chemical risk assessment as a previously uncharacterised source of uncertainty.

Supporting Information

S1 Fig. The DEBkiss model structure.

A schematic diagram for the energy flows in DEBkiss. The equations behind the fluxes are given in Table 1.

https://doi.org/10.1371/journal.pone.0140277.s001

(TIF)

S2 Fig. The effects of Cu exposure on traits responses in C. elegans in treatments with an average temperature of 12°C.

The five different endpoints: Final brood size (A), final body length (A), time to first egg (B), lifespan (D) and Population Growth Rate (PGR) (E) for the constant 12°C treatment (filled symbols) and the 12 ± 4°C (grey symbols) as a function of Cu concentrations in the agar. Data are given as mean ± s.e.m. and are described with a three parameter log-logistic concentration response model, except for TFE. The parameters are given in Table 3, together with the concentration-response parameters of the other temperature treatments.

https://doi.org/10.1371/journal.pone.0140277.s002

(TIF)

S3 Fig. The effects of Cu exposure on traits responses in C. elegans in treatments with an average temperature of 10°C.

The five different endpoints: Final brood size (A), final body length (A), time to first egg (B), lifespan (D) and Population Growth Rate (PGR) (E) for the constant 20°C treatment (filled symbols) and the 20 ± 4°C (grey symbols) as a function of Cu concentrations in the agar. Data are given as mean ± s.e.m. and are described with a three parameter log-logistic concentration response model, except for TFE. The parameters are given in Table 3, together with the concentration-response parameters of the other temperature treatments.

https://doi.org/10.1371/journal.pone.0140277.s003

(TIF)

Acknowledgments

DJS and CSV were supported by core NERC national capability funding under the Centre for Ecology and Hydrology Pollution and Ecological Risk science theme. NN was supported by the Research School of Environmental Chemistry, Microbiology and Toxicology in Denmark and a PhD grant given jointly from the University of Copenhagen and Centre for Ecology and Hydrology National Capability funding.

Author Contributions

Conceived and designed the experiments: NC NJN CS DJS. Performed the experiments: NC NJN. Analyzed the data: NC NJN CS DJS. Contributed reagents/materials/analysis tools: NC NJN DJS. Wrote the paper: NC DJS.

References

  1. 1. Van Straalen NM (2003) Ecotoxicology becomes stress ecology. Environmental Science & Technology 37: 324A–330A.
  2. 2. Boonstra R (2013) The ecology of stress: a marriage of disciplines. Functional Ecology 27: 7–10.
  3. 3. Holmstrup M, Bindesbol AM, Oostingh GJ, Duschl A, Scheil V, et al. (2010) Interactions between effects of environmental chemicals and natural stressors: A review. Science of the Total Environment 408: 3746–3762. pmid:19922980
  4. 4. Gracey AY, Fraser EJ, Li WZ, Fang YX, Taylor RR, et al. (2004) Coping with cold: An integrative, multitissue analysis of the transcriptome of a poikilothermic vertebrate. Proceedings of the National Academy of Sciences of the United States of America 101: 16970–16975. pmid:15550548
  5. 5. Murray P, Hayward SAL, Govan GG, Gracey AY, Cossins AR (2007) An explicit test of the phospholipid saturation hypothesis of acquired cold tolerance in Caenorhabditis elegans. Proceedings of the National Academy of Sciences of the United States of America 104: 5489–5494. pmid:17369360
  6. 6. Amarasekare P, Coutinho RM (2013) The intrinsic growth rate as a predictor of population viability under climate warming. Journal of Animal Ecology 82: 1240–1253. pmid:23926903
  7. 7. Gupta RC (2014) Biomarkers in Toxicology; Gupta RC, editors.
    • 8. Raffel TR, Halstead NT, McMahon TA, Davis AK, Rohr JR (2015) Temperature variability and moisture synergistically interact to exacerbate an epizootic disease. Proceedings of the Royal Society B-Biological Sciences 282.
    • 9. Terrell KA, Quintero RP, Murray S, Kleopfer JD, Murphy JB, et al. (2013) Cryptic impacts of temperature variability on amphibian immune function. Journal of Experimental Biology 216: 4204–4211. pmid:23948472
    • 10. Ben-Horin T, Lenihan HS, Lafferty KD (2013) Variable intertidal temperature explains why disease endangers black abalone. Ecology 94: 161–168. pmid:23600250
    • 11. Raffel TR, Romansic JM, Halstead NT, McMahon TA, Venesky MD, et al. (2013) Disease and thermal acclimation in a more variable and unpredictable climate. Nature Climate Change 3: 146–151.
    • 12. Heugens EHW, Hendriks AJ, Dekker T, Van Straalen NM, Admiraal W (2001) A review of the effects of multiple stressors on aquatic organisms and analysis of uncertainty factors for use in risk assessment. Critical Reviews in Toxicology 31: 247–284. pmid:11405441
    • 13. Cedergreen N, Norhave NJ, Nielsen K, Johansson HKL, Marcussen H, et al. (2013) Low temperatures enhance the toxicity of copper and cadmium to Enchytraeus crypticus through different mechanisms. Environmental Toxicology and Chemistry 32: 2274–2283. pmid:23661382
    • 14. Norhave NJ, Spurgeon D, Svendsen C, Cedergreen N (2012) How does growth temperature affect cadmium toxicity measured on different life history traits in the soil nematode Caenorhabditis elegans? Environmental Toxicology and Chemistry 31: 787–793. pmid:22253140
    • 15. Spurgeon DJ, Svendsen C, Lister L, Hankard PK, Kille P (2005) Earthworm responses to Cd and Cu under fluctuating environmental conditions: a comparison with results from laboratory exposures. Environmental Pollution 136: 443–452. pmid:15862398
    • 16. Jansch S, Frampton GK, Rombke J, VandenBrink PJ, ScottFordsmand JJ (2006) Effects of pesticides on soil invertebrates in model ecosystem and field studies: A review and comparison with laboratory toxicity data. Environmental Toxicology and Chemistry 25: 2490–2501. pmid:16986805
    • 17. Goranson NC, Ebersole JP, Brault S (2005) Resolving an adaptive conundrum: reproduction in Caenorhabditis elegans is not sperm-limited when food is scarce. Evolutionary Ecology Research 7: 325–333.
    • 18. Jager T, Martin BT, Zimmer EI (2013) DEBkiss or the quest for the simplest generic model of animal life history. Journal of Theoretical Biology 328: 9–18. pmid:23523873
    • 19. Brenner S (1974) Genetics of Caenorhabditis elegans. Genetics 77: 71–94. pmid:4366476
    • 20. Sivasundar A, Hey J (2003) Population genetics of Caenorhabditis elegans: The paradox of low polymorphism in a widespread species. Genetics 163: 147–157. pmid:12586703
    • 21. Klass MR (1977) Aging in nematode Caenorhabditis-elegans—major biological and environmental-factors influencing life-span. Mechanisms of Ageing and Development 6: 413–429. pmid:926867
    • 22. Norhave NJ, Spurgeon D, Svendsen C, Cedergreen N (2014) The importance of experimental time when assessing the effect of temperature on toxicity in poikilotherms. Environmental Toxicology and Chemistry 33: 1363–1371. pmid:24648200
    • 23. Alda Alvarez OA, Jager T, Kooijman SALM, Kammenga JE (2005) Responses to stress of Caenorhabditis elegans populations with different reproductive strategies. Functional Ecology 19: 656–664.
    • 24. Hoess S, Schlottmann K, Traunspurger W (2011) Toxicity of ingested cadmium to the nematode Caenorhabditis elegans. Environmental Science & Technology 45: 10219–10225.
    • 25. Ritz C, Streibig JC (2005) Bioassay analysis using R. Journal of Statistical Software 12: 1–22.
    • 26. Stubben C, Milligan B (2007) Estimating and analyzing demographic models using the popbio package in R. Journal of Statistical Software 22: 1–23.
    • 27. Caswell H (2000) Life table response experiments in ecotoxicology. In: Kammenga J Laskowski R, editors. Demography in Ecotoxicology. Chichester, UK: John Wiley.
      • 28. Cedergreen N, Svendsen C, Backhaus T (2012) Chemical Mixtures: Concepts for Predicting Toxicity. In: Jorgensen SE, editor. Encyclopedia of Environmental Management. New York, USA: Taylor & Francis. pp. 2601–2610.
        • 29. Jager T, Gudmundsdottir EM, Cedergreen N (2014) Dynamic modeling of sublethal mixture toxicity in the nematode Caenorhabditis elegans. Environmental Science & Technology 48: 7026–7033.
        • 30. Kooijman SALM (2000) Dynamic Energy Budget Models in Biological Systems. Cambridge: Cambridge University Press.
          • 31. Hill RW, Wyse G (1989) Thermal Relations. Animal Physiology 6.
          • 32. Lika K, Kearney MR, Freitas V, van der Veer HW, van der Meer J, et al. (2011) The "covariation method" for estimating the parameters of the standard Dynamic Energy Budget model I: Philosophy and approach. Journal of Sea Research 66: 270–277.
          • 33. Jager T, Albert C, Preuss TG, Ashauer R (2011) General Unified Threshold Model of Survival—a Toxicokinetic-Toxicodynamic Framework for Ecotoxicology. Environmental Science & Technology 45: 2529–2540.
          • 34. Jager T, Zimmer EI (2012) Simplified Dynamic Energy Budget model for analysing ecotoxicity data. Ecological Modelling 225: 74–81.
          • 35. Jager T, Heugens EHW, Kooijman SALM (2006) Making sense of ecotoxicological test results: Towards application of process-based models. Ecotoxicology 15: 305–314. pmid:16739032
          • 36. Bundy JG, Sidhu JK, Rana F, Spurgeon DJ, Svendsen C, et al. (2008) Systems toxicology’ approach identifies coordinated metabolic responses to copper in a terrestrial non-model invertebrate, the earthworm Lumbricus rubellus. BMC Biology 6: 1–25.
          • 37. Angilletta MJ, Dunham AE (2003) The temperature-size rule in ectotherms: Simple evolutionary explanations may not be general. American Naturalist 162: 332–342. pmid:12970841
          • 38. Aguilar-Alberola JA, Mesquita-Joanes F (2014) Breaking the temperature-size rule: Thermal effects on growth, development and fecundity of a crustacean from temporary waters. Journal of Thermal Biology 42: 15–24. pmid:24802144
          • 39. Alvarez OA, Jager T, Kooijman SALM, Kammenga JE (2005) Responses to stress of Caenorhabditis elegans populations with different reproductive strategies. Functional Ecology 19: 656–664.
          • 40. Barriere A, Felix M-A (2007) Temporal dynamics and linkage disequilibrium in natural Caenorhabditis elegans populations. Genetics 176: 999–1011. pmid:17409084
          • 41. Kiontke KC, Felix M-A, Ailion M, Rockman MV, Braendle C, et al. (2011) A phylogeny and molecular barcodes for Caenorhabditis, with numerous new species from rotting fruits. BMC Evolutionary Biology 11.
          • 42. Petersen C, Dirksen P, Prahl S, Strathmann EA, Schulenburg H (2014) The prevalence of Caenorhabditis elegans across 1.5 years in selected North German locations: the importance of substrate type, abiotic parameters, and Caenorhabditis competitors. BMC Ecology 14.
          • 43. Boardman L, Sorensen JG, Terblanche JS (2013) Physiological responses to fluctuating thermal and hydration regimes in the chill susceptible insect, Thaumatotibia leucotreta. Journal of Insect Physiology 59: 781–794. pmid:23684741
          • 44. Xing K, Hoffmann AA, Ma C-S (2014) Does Thermal Variability Experienced at the Egg Stage Influence Life History Traits across Life Cycle Stages in a Small Invertebrate? Plos One 9.
          • 45. Kammenga JE, Busschers M, Van Straalen NM, Jepson PC, Bakker J (1996) Stress induced fitness reduction is not determined by the most sensitive life-cycle trait. Functional Ecology 10: 106–111.
          • 46. Korsloot A, Van Gestel CAM, Van Straalen NM (2004) Environmental Stress and Cellular Response in Arthropods. London, UK: CRC Press.
            • 47. Hayward SAL, Manso B, Cossins AR (2014) Molecular basis of chill resistance adaptations in poikilothermic animals. Journal of Experimental Biology 217: 6–15. pmid:24353199
            • 48. Bindesbol AM, Bayley M, Damgaard C, Hedlund K, Holmstrup M (2009) Changes in membrane phospholipids as a mechanistic explanation for decreased freeze tolerance in earthworms exposed to sublethal copper concentrations. Environmental Science & Technology 43: 5495–5500.
            • 49. Bindesbol AM, Holmstrup M, Damgaard C, Bayley M (2005) Stress synergy between environmentally realistic levels of copper and frost in the earthworm Dendrobaena octaedra. Environmental Toxicology and Chemistry 24: 1462–1467. pmid:16117123
            • 50. Holmstrup M, Petersen BF, Larsen MM (1998) Combined effects of copper, desiccation, and frost on the viability of earthworm cocoons. Environmental Toxicology and Chemistry 17: 897–901.
            • 51. Swain S, Wren J, Stürzenbaum SR, Kille P, Jager T, et al. (2010) Linking toxicants mechanism of action and physiological mode of action in Caenorhabditis elegans. BMC Systems Biology 4
            • 52. Alda Alvarez OA, Jager T, Redondo EM, Kammenga JE (2006) Physiological modes of action of toxic chemicals in the nematode Acrobeloides nanus. Environmental Toxicology and Chemistry 25: 3230–3237. pmid:17220093
            • 53. Hosseini M-J, Shaki F, Ghazi-Khansari M, Pourahmad J (2014) Toxicity of copper on isolated liver mitochondria: Impairment at complexes I, II, and IV leads to increased ROS production. Cell Biochemistry and Biophysics 70: 367–381. pmid:24691927
            • 54. Shaw PJ, Tononi G, Greenspan RJ, Robinson DF (2002) Stress response genes protect against lethal effects of sleep deprivation in Drosophila. Nature 417: 287–291. pmid:12015603
            • 55. Jenkins DE, Schultz JE, Matin A (1988) Starvation-induced cross protection against heat or H2O2 challenge in Escherichia-coli. Journal of Bacteriology 170: 3910–3914. pmid:3045081
            • 56. Wren JF, Kille P, Spurgeon DJ, Swain S, Sturzenbaum SR, et al. (2011) Application of physiologically based modelling and transcriptomics to probe the systems toxicology of aldicarb for Caenorhabditis elegans (Maupas 1900). Ecotoxicology 20: 397–408. pmid:21253838
            • 57. Baas J, Jager T, Kooijman B (2010) A review of DEB theory in assessing toxic effects of mixtures. Science of the Total Environment 408: 3740–3745. pmid:19850324
            • 58. Jager T, Crommentuijn T, VanGestel CAM, Kooijman SALM (2004) Simultaneous modeling of multiple end points in life-cycle toxicity tests. Environmental Science & Technology 38: 2894–2900.
            • 59. Guinnee MA, Gardner A, Howard AE, West SA, Little TJ (2007) The causes and consequences of variation in offspring size: a case study using Daphnia. Journal of Evolutionary Biology 20: 577–587. pmid:17305824
            • 60. Stelzer CP (2001) Resource limitation and reproductive effort in a planktonic rotifer. Ecology 82: 2521–2533.
            • 61. Baas J, Jager T, Kooijman B (2010) Understanding toxicity as processes in time. Science of the Total Environment 408: 3735–3739. pmid:19969324
            • 62. Jager T, Vandenbrouck T, Baas J, De Coen WM, Kooijman S (2010) A biology-based approach for mixture toxicity of multiple endpoints over the life cycle. Ecotoxicology 19: 351–361. pmid:19771510
            • 63. Kooijman S, Bedaux JJM (1996) Analysis of toxicity tests on Daphnia survival and reproduction. Water Research 30: 1711–1723.
            • 64. Philipp EER, Wessels W, Gruber H, Strahl J, Wagner AE, et al. (2012) Gene Expression and Physiological Changes of Different Populations of the Long-Lived Bivalve Arctica islandica under Low Oxygen Conditions. Plos One 7.
            • 65. Malmendal A, Overgaard J, Bundy JG, Sorensen JG, Nielsen NC, et al. (2006) Metabolomic profiling of heat stress: hardening and recovery of homeostasis in Drosophila. American Journal of Physiology-Regulatory Integrative and Comparative Physiology 291: R205–R212.