Subscriber access provided by University of Newcastle, Australia
Article
Feeding activity and xenobiotics modulate oxidative status in Daphnia magna: implications for ecotoxicological testing Sara Furuhagen, Birgitta Liewenborg, Magnus Breitholtz, and Elena Gorokhova Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es5044722 • Publication Date (Web): 23 Sep 2014 Downloaded from http://pubs.acs.org on September 28, 2014
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 24
Environmental Science & Technology
1
Feeding activity and xenobiotics modulate oxidative status in
2
Daphnia magna: implications for ecotoxicological testing
3
Sara Furuhagen*, Birgitta Liewenborg, Magnus Breitholtz and Elena Gorokhova
4 5
Department of Applied Environmental Science, Stockholm University, Svante Arrhenius väg 8, 106 91 Stockholm, Sweden
6
*Corresponding author:
7
Sara Furuhagen
8 9
Department of Applied Environmental Science, Stockholm University, Svante Arrhenius väg 8, 106 91 Stockholm, Sweden
10
+46 8 674 72 79
11
[email protected] 12
13
14
15
16
17
18
19
20
21
22 1 ACS Paragon Plus Environment
Environmental Science & Technology
23
ABSTRACT
24
To apply biomarkers of oxidative stress in laboratory and field settings, an understanding of
25
their responses to changes in physiological rates is important. The evidence is accumulating
26
that caloric intake can increase production of reactive oxygen species and thus affect
27
background variability of oxidative stress biomarkers in ecotoxicological testing. This study
28
aimed to delineate effects of food intake and xenobiotics on oxidative biomarkers in Daphnia
29
magna. Antioxidant capacity measured as ORAC (oxygen radical absorbance capacity) and
30
lipid peroxidation assayed as TBARS (thiobarbituric acid reactive substances) were measured.
31
Food intake was manipulated by varying food densities or by exposing the animals to
32
chemicals inhibiting feeding rate (pharmaceutical haloperidol and pesticide lindane). Feeding
33
rate proved to affect both protein, ORAC and TBARS in unexposed daphnids. However, there
34
was no significant effect of feeding rate on the protein-specific ORAC values. Both
35
substances affected individual protein and ORAC levels and changed their relationship to
36
feeding rate. Our results show that inhibition of feeding rate influenced the interpretation of
37
biomarker response and further emphasize the importance of understanding (1) baseline
38
variability in potential biomarkers due to variations in metabolic state, and (2) the
39
contribution of feeding rate on toxic response of biomarkers.
40
41
42
43
44
45 2 ACS Paragon Plus Environment
Page 2 of 24
Page 3 of 24
Environmental Science & Technology
46
INTRODUCTION
47
To counteract pro-oxidative processes in aerobic organisms, homeostasis is maintained
48
between the cellular production of reactive oxygen species (ROS) and the endogenous
49
antioxidant defense. The balance between ROS production and the antioxidant defense can be
50
affected by physiological factors, such as age and disease1 as well as by environmental
51
factors, such as hypoxia2 and pollutants3. Caloric restriction (CR) has been identified as an
52
important factor affecting the cellular production of ROS as low-calorie diets hamper ROS
53
production4. Studies on isolated rat mitochondria showed that CR results in decreased
54
substrate oxidation activity, leading to lowered mitochondrial membrane potential and
55
increased proton leakage and thus a diminished generation of ROS5. This effect appears to be
56
a general mechanism as the effects of CR have been shown for a variety of species6. As an
57
excess of ROS can be harmful to DNA, lipids and proteins7, the levels of oxidative damages
58
on these biomolecules have been found to be lower in animals given a CR diet6.
59
Chemical substances can affect cellular ROS concentrations through different mechanisms.
60
Pro-oxidative substances can increase generation of ROS and other radicals by entering redox
61
cycles8, but ROS production can be also induced by increasing metabolic rate as a response to
62
stress or through induction of enzymes involved in xenobiotic metabolism, such as CYP450
63
and NAD(P)H9. Moreover, the oxidative homeostasis can be disrupted by depletion of
64
antioxidative substances in xenobiotic metabolism7, 9. Consequently, biomarkers of oxidative
65
stress are widely used in ecotoxicology and ecology as indicators of stressful conditions10,
66
although specific causes of the observed biochemical alterations are not always identifiable
67
because of the plethora of factors that may affect the oxidative status of organisms. Therefore,
68
to facilitate biomarker application and interpretation of measured responses to stress factors,
69
we need to understand the magnitude and causes of background variability for the biomarker
70
of interest11. In particular, variability in feeding rate has a high potential to affect biomarkers 3 ACS Paragon Plus Environment
Environmental Science & Technology
71
of oxidative stress, because of the effect of caloric intake on the production of ROS. Feeding
72
rate is a sensitive and ecologically relevant end point in ecotoxicological assays as exposure
73
to a wide range of substances has been reported to cause feeding inhibition in various test
74
species12-13. In stress assessment, feeding inhibition assay proved to be about 50-fold more
75
sensitive than standardized acute ecotoxicological assays employing survival as an end
76
point14. Moreover, although effects of ad libitum feeding on many endpoints have been
77
identified as a serious methodological issue in tests with vertebrates15-17, no attempt has been
78
made so far to study the effects of feeding rate and xenobiotic exposure on oxidative
79
biomarkers in invertebrates. The objectives of this study were to (1) establish relationships
80
between feeding rate and oxidative biomarkers in Daphnia magna and (2) understand the
81
contribution of variation in feeding rate to toxic response.
82
Two model substances, haloperidol and lindane, were used to address these objectives.
83
Haloperidol is a dopamine receptor antagonist that has been found to cause feeding inhibition
84
in rats by blocking hunger perception18; and it has indeed been found to cause feeding
85
inhibition in D. magna (Furuhagen, unpubl. data). Lindane is a neurotoxic insecticide19 that
86
has been shown to inhibit feeding in daphnids by reducing the movement of filtering limbs12,
87
20-21
88
varying food levels in the absence of a toxic compound or to varying substance concentrations
89
at a constant food level. As proxies for antioxidant capacity and oxidative damage we assayed
90
oxygen radical absorbance capacity (ORAC) and lipid peroxidation levels measured as
91
thiobarbituric acid reactive substances (TBARS). Both biomarkers have been used before to
92
assess oxidative stress in microcrustaceans, including daphnids22-23. Two sets of hypotheses
93
were tested with regard to the effects of food intake and chemical exposure in this study. We
94
hypothesized that individual protein content will be positively related to feeding rate in the
95
unexposed daphnids (Hypothesis 1), due to increased protein synthetic activity at higher food
. In our study, we manipulated feeding activity in daphnids by exposing them to either
4 ACS Paragon Plus Environment
Page 4 of 24
Page 5 of 24
Environmental Science & Technology
96
intake. Moreover, the oxidative biomarkers will be positively related to the increase in protein
97
content of animals feeding at higher rates and thus having higher ROS production (Hypothesis
98
2). In addition to feeding inhibition and corresponding effects on the biomarker response, the
99
exposure to toxicants was hypothesized to have effects on the relationships between feeding
100
rate and protein allocation, and between protein content and oxidative biomarkers. In
101
particular, a negative effect on individual protein content was expected in the haloperidol-
102
exposed daphnids, whereas lindane was hypothesized to have a positive effect on the protein
103
content. The rationale for this expectation was based on the reported negative effect of
104
haloperidol on the expression of heat shock proteins (hsp)24 and the induction of hsp after
105
lindane exposure25 (Hypothesis 3). An increase in oxidative damage, independent of feeding
106
rates, was expected in exposed daphnids, whereas the antioxidant capacity could be both
107
positively and negatively affected by xenobiotics depending on the stress levels (Hypothesis
108
4).
109
110
METHODS
111
Test organism
112
The daphnids originated from a single clone (environmental pollution test strain Klon 5, the
113
State office for nature, environment, and customer protection North-Rhine Westfalia, Bonn,
114
Germany; originally from the Federal Environment Agency, Berlin, Germany). They were
115
cultured in 2 L M7 medium (OECD standard 202 and 211) at a stock density of ~20
116
individuals in each beaker and fed three times a week with a mixture of Pseudokirchneriella
117
subcapitata and Scenedesmus subspicatus.
118
5 ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 24
119
Functional response experiments
120
Using varying densities of algal suspension we manipulated daphnid feeding rate. The
121
individuals feeding at varying rates were also used to measure individual protein content,
122
ORAC (Experiment I and II), and TBARS (Experiment II) (Table 1). The feeding rate of
123
juvenile daphnids (96 h old) as a function of food density was measured in M7 medium with 5
124
daphnids in each replicate. Fluorescence of the algae was measured using Turner designs 10-
125
AU Fluorometer at the start of the test and after 24 h to assess the feeding rate of the
126
daphnids. All experiments were conducted in darkness, at 20˚C to minimize algal growth
127
during the incubation. An additional replicate without daphnids served as a control to assess
128
algal growth. Incubation for ORAC and TBARS measurements (Experiment II) was done
129
using starved animals and those incubated at the intermediate food concentration (1.5 µg C
130
mL-1) and at the near-saturation concentration (7 µg C mL-1).
131
132
Feeding inhibition tests
133
Feeding inhibition tests (haloperidol: Experiment III and lindane: Experiment V; test
134
conditions as in Experiment I) were performed according to Barata, et al.
135
modifications. In these experiments, feeding rate was measured at a constant food
136
concentration (1.5 µg C mL-1) and the animals were used for measurements of individual
137
protein content and ORAC. Exposure for TBARS and ORAC (Experiments IV and VI for
138
haloperidol and lindane, respectively) were conducted at a food concentration of 1.5 µg C mL-
139
1
140
lindane were dissolved in dimethyl sulfoxide (DMSO) for use in the bioassays. Solvent
141
control with a DMSO concentration corresponding to a volume of 0.1‰ of the total
142
incubation volume was used in all exposure tests.
26
with minor
and using the same test conditions as in Experiment II (Table 1). Both haloperidol and
6 ACS Paragon Plus Environment
Page 7 of 24
Environmental Science & Technology
143
144
Biochemical analyses
145
At the termination of each experiment, live daphnids were pooled within replicates, rinsed in
146
potassium phosphate buffer (PPB) and transferred to Eppendorf tubes. All samples were snap
147
frozen and stored at -80˚C until analysis.
148
149
Protein quantification
150
Samples for protein and ORAC measurements were homogenized in 150 µL PPB buffer (3.1
151
mM, pH 7.4). Protein content (µg ind-1) was determined by bicinchoninic acid method27 using
152
PierceTM BCA Protein Assay kit (Thermo Scientific) according to microplate procedure with
153
some modifications. Transparent microplate was used and total volume in the wells was 150
154
µL, with 10 µL homogenate, 10 µL PPB and 130 µL working solution. Absorbance was
155
measured at 540 nm using FluoStar Optima plate reader (BMG Lab Technologies, Germany).
156
157
Oxygen radical absorbance capacity (ORAC)
158
Total antioxidant capacity was assayed as ORAC according to Ou, et al.
159
modifications. Fluorescein was used as a fluorescent probe (106 nM) and 2,2- azobis(2-
160
amidinopropane) dihydrochloride (AAPH) (152.66 mM) as a source of peroxyl radicals.
161
Trolox (218 µM, Sigma-Aldrich) was used as standard. In each assay, 20 µL sample was
162
added to each well and mixed with 30 µL AAPH and 150 µL fluorescein. ORAC levels were
163
expressed as ORAC ind-1. Additionally, ORAC values were normalized to protein content
7 ACS Paragon Plus Environment
28
, with minor
Environmental Science & Technology
Page 8 of 24
164
(ORACp) before related to feeding rate because of the high contribution of proteins to ORAC
165
values29.
166
167
Lipid peroxidation
168
Two exposures were performed for unexposed and exposed daphnids for TBARS analysis
169
(Experiment II, IV, VI; Table 1) with high (A) and low (B) population density, respectively.
170
Daphnids were assayed for lipid peroxidation by a modified TBARS method for measuring
171
the aldehydic lipid peroxidation decomposition derivatives, which form fluorescent products
172
after reacting with thiobarbituric acid (TBA)30. Daphnids were homogenized in 200 µL PPB,
173
sonicated 3×10s and aliquoted. Subsamples of 125 µL tissue homogenate were treated with
174
125 µL 10% trichloroacetic acid. The sample was further mixed with 150 µL TBA (2 mM)
175
and incubated at 100˚C for 1 h. After cooling to room temperature, 220 µL of
176
butanol:pyridine (volume ratio 15:1) mixture were added, vortexed (2×10s) and centrifuged
177
for 5 minutes at 4000 g. The organic phase was used for fluorometric determination
178
(540nm/590nm) of malondialdehyde (MDA) concentration (µM MDA equivalents ind-1). The
179
MDA concentrations were expressed as TBARS ind-1. Following the common procedure for
180
biochemical end point, TBARS were normalized to protein content (TBARSp) before related
181
to feeding rate.
182
183
Data analysis and statistics
184
Feeding rate was calculated according to Båmstedt, et al.
185
unexposed daphnids was fitted to a sigmoid function32-33
186
= /1 +
31
. The functional response of
(1) 8 ACS Paragon Plus Environment
Page 9 of 24
Environmental Science & Technology
187
where F is feeding rate (µg C h-1 ind-1), Fmax is the theoretical maximum feeding rate, e is
188
Euler´s number, k is a constant, C is food concentration (µg C mL-1) and Cm is the half-
189
saturation constant.
190
To test hypotheses 1 and 2, the relationships between feeding rate and biomarkers were
191
analyzed by linear regression analysis. The effect of chemical exposure on feeding rate was
192
evaluated using Hill equation
193
= + − × / +
194
where F is feeding rate (µg C h-1 ind-1), Fmin is the minimum feeding rate (µg C h-1 ind-1) and
195
Fmax the maximum feeding rate (µg C h-1 ind-1), C is the logarithm of the substance
196
concentration (mg L-1) +1 and n is the Hill constant34. The substance concentration
197
(logarithmic value) that causes 50% feeding inhibition is represented by EC50 (mg L-1).
198
Hypotheses 3 and 4 were addressed using general linear models to test the impact of
199
substance exposure on the relationships between feeding and protein allocation and between
200
protein content and oxidative biomarkers; experimental run (A and B) was used as a co-
201
variable to account for the differences in experimental design. All values were normalized to
202
the mean value of the control group for each treatment. To facilitate statistical comparisons
203
and interpretation, the independent variables were centered by their respective mean values
204
using response interval covering both exposed and unexposed daphnids. Centering input
205
variables removes the collinearity between the main effects and the interaction predictors,
206
therefore allowing the interpretation of main effects35. The interactions between explanatory
207
variables were removed from the model when non-significant. Assumptions of normality and
208
homoscedasticity were fulfilled for all data as confirmed by residual plot analysis. Significant
209
level was ˂0.05 for all tests and all analyses were carried out in R 2.13.2.
9 ACS Paragon Plus Environment
(2)
Environmental Science & Technology
210
Projection to latent structures by means of partial least squares (PLS) was applied to visualize
211
directions in a multivariate space for maximum separation of observations (biomarkers and
212
physiological variables across an individual) belonging to different groups (treatments). The
213
amount of variation attributed to each explanatory variable was determined by regressing each
214
explanatory variable against the response variable (Treatment as a categorical variable) in the
215
absence of the other explanatory variables36 as implemented in STATISTICA 8 (StatSoft, Inc.
216
2013).
217
218
RESULTS
219
Feeding and biomarker responses in unexposed animals
220
Feeding rate increased with increasing food concentration, reaching saturation at ~6.7 µg C
221
mL-1 (Figure 1). Moreover, individual protein content was significantly positively related to
222
feeding rate (Figure 3). Also, there was a significant positive relationship between individual
223
ORAC and protein values. However, no significant relationship between ORACp and feeding
224
rate was found. Finally, individual protein content and TBARS as well as TBARSp and
225
feeding rate were positively related to each other (Table 2; Figure 3).
226
227
Feeding and biomarker responses in exposed daphnids
228
Neither haloperidol nor lindane caused mortality within the tested concentration range.
229
However, both substances had inhibitory effect on feeding rate, following a sigmoid dose-
230
response and resulting in nearly complete feeding inhibition in the highest concentrations of
231
both substances. Haloperidol EC50 was 1.62 mg L-1 (95% confidence interval: 1.37 to 1.90 mg
10 ACS Paragon Plus Environment
Page 10 of 24
Page 11 of 24
Environmental Science & Technology
232
L-1) and lindane EC50 was 1.27 mg L-1 (95% confidence interval: 1.03 to 1.55 mg L-1) (Figure
233
2).
234
Similar to the unexposed daphnids, individual protein content was significantly positively
235
related to feeding rate in exposed daphnids. Moreover, haloperidol had a significant negative
236
effect on the baseline protein content and a positive effect on the slope on the relationship
237
between protein and feeding as indicated by the significant interaction term of feeding rate
238
and treatment. By contrast, lindane had no significant effect on the individual protein content
239
(Table 3A; Figure 3).
240
Again, similar to the unexposed daphnids, ORAC was positively related to protein content in
241
both haloperidol and lindane exposed daphnids. Moreover, haloperidol had a positive effect
242
on the ORAC – protein relationship, as indicated by the significant interaction between
243
treatment and protein content. Lindane significantly decreased the baseline ORAC values as
244
indicated by the significantly different intercept (Table 3B).
245
In contrast to the unexposed daphnids, a significant positive relationship between ORACp and
246
feeding rate was observed in the animals exposed to haloperidol or lindane. Moreover, both
247
substances had significant negative effects on the ORACp baseline as indicated by the
248
significant treatment effect (Table 3C).
249
Haloperidol had no significant effect on the relationship between TBARS and protein content.
250
However, the significant interaction for lindane and protein indicates a significantly negative
251
effect on the TBARS-protein relationship due to lindane exposure (Table 3D). Neither
252
haloperidol nor lindane had any significant effect on the relationship between TBARSp and
253
feeding rate (Table 3E).
254
The PLS analysis showed a clear pattern related to the differences between the lindane-
255
exposed and unexposed daphnids revealing relevant ORACp, feeding rate and TBARSp in 11 ACS Paragon Plus Environment
Environmental Science & Technology
256
this separation. The two PLS components accounted for 69% of the total variance. No clear
257
separation between the haloperidol-exposed and unexposed animals or animals exposed to
258
low lindane concentrations was observed, with variability in both haloperidol-exposed and
259
unexposed groups being largely affected by variations in feeding and individual protein
260
content (Figure 3).
261
262
DISCUSSION
263
In Daphnia magna, biomarkers of oxidative status changed in concert with feeding rate, that
264
thus was established as a confounding factor for oxidative biomarkers in this standard test
265
species. In addition to inducing feeding inhibition, exposure to our model substances,
266
haloperidol and lindane, altered the relationship between feeding rate, protein and oxidative
267
biomarkers, emphasizing the complexity of the biomarker responses.
268
Hypothesis 1, linking protein content to feeding rate, was confirmed by the significant
269
positive relationship between feeding rate and protein content, indicating that protein
270
synthesis increases as more resources become available for protein production37. Since protein
271
content is commonly used to normalize enzyme activities and other biochemical constituents
272
used as biomarkers in ecology and ecotoxicology, this response is important for biomarker
273
applications in general.
274
Hypothesis 2, linking oxidative biomarkers to individual protein content, was also supported.
275
In both exposed and unexposed daphnids, individual ORAC values were positively related to
276
protein content, most probably due to increased production and intake of low-molecular
277
compounds and proteins with antioxidative properties. The antioxidant defense consists of
278
enzymes as well as other proteins and low-molecular-mass agents, such as ascorbic acid,
279
reduced glutathione, methionine and uric acid7. The antioxidative activity of the water soluble 12 ACS Paragon Plus Environment
Page 12 of 24
Page 13 of 24
Environmental Science & Technology
280
fraction of these compounds is measured in the ORAC assay28. The lack of correlation
281
between ORACp and feeding rate in the unexposed daphnids suggests that the allocation of
282
proteins to the antioxidant defense did not increase in response to increased caloric intake,
283
which is in agreement with previous findings38. A baseline relationship between the feeding
284
rate and alterations in antioxidant defense has thus been established, which also facilitates
285
application of ORACp as a biomarker in ecological and ecotoxicological studies11.
286
Even though EC50 values for feeding inhibition were similar between haloperidol and lindane
287
as evidenced by the overlapping confidence intervals, their effects on protein and ORAC
288
differed, which may be related to differences in their mode of action18, 20. By including protein
289
and substance exposure in the statistical models we could confirm Hypothesis 3 and show that
290
haloperidol had a positive effect on the relationship between individual protein content and
291
feeding rate. This was likely due to decreased protein synthesis with increasing haloperidol
292
concentration39. By contrast to haloperidol, lindane did not have any effect on protein content.
293
Hence, any difference in protein content in response to increasing lindane concentration was
294
solely due to the decrease in filtering activity and thus feeding rate in the exposed animals.
295
Without including feeding rate as an explanatory variable to the model, this effect would have
296
been missed and observed alterations in protein content, and in the biomarker normalized to
297
protein content, may erroneously have been attributed direct toxic mechanisms on protein
298
metabolism.
299
The Hypothesis 4 also found support as both haloperidol and lindane affected the relationship
300
between protein and ORAC. Exposure to haloperidol resulted in lower allocation of resources
301
to the antioxidant defense in relation to protein content, which is consistent with a general
302
response to starvation, whereas lindane caused a significant decrease in the baseline ORAC
303
values. Due to the effects of haloperidol and lindane on protein allocation and ORAC, there
304
were significant positive effects on the relationship between ORACp and feeding rate in the 13 ACS Paragon Plus Environment
Environmental Science & Technology
305
exposed daphnids. The altered ORAC-protein relationship and the positive relation between
306
ORACp and feeding rate indicate reduced resource allocation to antioxidant defense and/or a
307
depletion of antioxidative compounds due to the chemical exposure. One can speculate that
308
these effects may be related to detoxification as some agents contributing to the antioxidant
309
defense are also involved in xenobiotic metabolism. One of these, glutathione (GSH), is a co-
310
factor in the detoxification of ROS7 and depletion of GSH is used as a biomarker of oxidative
311
stress due its involvement in antioxidative reactions40. Both haloperidol and lindane have been
312
shown to cause GSH depletion due to its involvement in ROS detoxification and xenobiotic
313
metabolism41-43. Hence, depletion of GSH may partly explain the decrease in ORAC in
314
relation to protein content observed for exposed daphnids.
315
Contrary to Hypothesis 4, predicting that lipid peroxidation levels would increase in response
316
to xenobiotic exposure, haloperidol did not altered the relationship between protein content
317
and TBARS, thus indicating that lipid peroxidation levels were only affected by feeding rate
318
(Figure 3). These results contradict reported effects of this drug on lipid peroxidation in
319
humans that had been administrated haloperidol (10 mg day-1) for two weeks41. However, the
320
tested concentrations in our study may have been too low to induce oxidative damages during
321
the 24-h exposure and the antioxidant defense could successfully counteracted ROS effects
322
with non-detectable effects on lipid peroxidation levels. The high survival rate in this study
323
could be indicative of such responses, albeit high survivorship does not necessarily mean that
324
an oxidative stress response is absent23. Lindane significantly lowered the TBARS levels at
325
increasing protein content compared to the unexposed daphnids. Due to the positive effect of
326
lindane on hsp, and thus protein composition25, normalizing TBARS to the total protein
327
concentrations in the sample could introduce additional uncertainty and even be misleading
328
for understanding of chemical exposure effects. Assessment of normalization strategies for
14 ACS Paragon Plus Environment
Page 14 of 24
Page 15 of 24
Environmental Science & Technology
329
oxidative stress biomarkers is needed to facilitate interpretation of these responses in the field
330
and laboratory settings.
331
The observed effects of calorie intake on TBARS levels (Table 2, Figure 3) raise questions
332
about optimal testing design for studies employing oxidative stress biomarkers to evaluate
333
effects of various stressors. The ad libitum feeding regime used in many ecotoxicological tests
334
may lead to high lipid peroxidation levels in actively feeding non-stressed individuals and,
335
further, to erroneous interpretation of the results. Hence, bioassays with D. magna should be
336
performed at moderate food levels to avoid effects on biochemical markers that are solely a
337
consequence of ad libitum feeding. Moderate dietary restriction has also been shown to
338
increase the sensitivity and response to toxic substances in both rodents15 and rotifers44, and
339
our results indicate that this is also the case for oxidative biomarkers in D. magna.
340
In conclusion, our results indicate that alterations in oxidative biomarker response due to
341
xenobiotic exposure may not only be a consequence of direct interactions with oxidative
342
processes, but can also be a result of indirect response via feeding inhibition. Since protein
343
content, ORAC and TBARS were positively correlated with food intake, it is necessary to
344
account for inhibitory effects of xenobiotics on feeding rate when interpreting responses of
345
these biomarkers. Without considering confounding factors representing nutritional status and
346
metabolism when evaluating biomarker response to toxic exposure, there is a risk of making
347
erroneous conclusions about toxicity effects and mechanisms. To further increase the value of
348
biomarkers in general and oxidative biomarkers in particular, the importance of confounding
349
factors for interpretation of biomarker responses needs to be addressed routinely in biomarker
350
validation.
351
352
ACKNOWLEDGMENTS 15 ACS Paragon Plus Environment
Environmental Science & Technology
353
This study was supported by the Swedish Research Council for Environment, Agricultural
354
Science and Spatial Planning (FORMAS), Stockholm University’s strategic marine
355
environmental research program “Baltic Ecosystem Adaptive Management” and the Swedish
356
Foundation for Strategic Environmental Research (MISTRA; MistraPharma).
357
SUPPORTING INFORMATION
358
This information is available free of charge via the internet at http://pubs.acs.org.
359 360
REFERENCES
361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390
1. Yu, B. P., Aging and oxidative stress: Modulation by dietary restriction. Free Radical Biology and Medicine 1996, 21 (5), 651-668. 2. Valko, M.; Leibfritz, D.; Moncol, J.; Cronin, M. T. D.; Mazur, M.; Telser, J., Free radicals and antioxidants in normal physiological functions and human disease. International Journal of Biochemistry & Cell Biology 2007, 39 (1), 44-84. 3. Mustafa, S. A.; Al-Subiai, S. N.; Davies, S. J.; Jha, A. N., Hypoxia-induced oxidative DNA damage links with higher level biological effects including specific growth rate in common carp, Cyprinus carpio L. Ecotoxicology 2011, 20 (6), 1455-1466. 4. Speakman, J. R.; Mitchell, S. E., Caloric restriction. Molecular aspects of medicine 2011, 32 (3), 159-221. 5. Lambert, A. J.; Merry, B. J., Effect of caloric restriction on mitochondrial reactive oxygen species production and bioenergetics: reversal by insulin. American journal of physiology. Regulatory, integrative and comparative physiology 2004, 286 (1), R71-9. 6. Sohal, R. S.; Weindruch, R., Oxidative stress, caloric restriction, and aging. Science 1996, 273 (5271), 59-63. 7. Halliwell, B.; Gutteridge, J. M. C., Free radicals in biology and medicine. 4 ed.; Oxford University press: United States, 2007. 8. Lushchak, V. I., Environmentally induced oxidative stress in aquatic animals. Aquatic Toxicology 2011, 101 (1), 13-30. 9. Livingstone, D. R., Contaminant-stimulated reactive oxygen species production and oxidative damage in aquatic organisms. Marine Pollution Bulletin 2001, 42 (8), 656-666. 10. Valavanidis, A.; Vlahogianni, T.; Dassenakis, M.; Scoullos, M., Molecular biomarkers of oxidative stress in aquatic organisms in relation to toxic environmental pollutants. Ecotoxicology and Environmental Safety 2006, 64 (2), 178-189. 11. van der Oost, R.; Beyer, J.; Vermeulen, N. P. E., Fish bioaccumulation and biomarkers in environmental risk assessment: a review. Environmental Toxicology and Pharmacology 2003, 13 (2), 57-149. 12. McWilliam, R. A.; Baird, D. J., Postexposure feeding depression: A new toxicity endpoint for use in laboratory studies with Daphnia magna. Environmental Toxicology and Chemistry 2002, 21 (6), 1198-1205.
16 ACS Paragon Plus Environment
Page 16 of 24
Page 17 of 24
391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442
Environmental Science & Technology
13. Domene, X.; Natal-da-Luz, T.; Alcaniz, J. M.; Andres, P.; Sousa, J. P., Feeding inhibition in the soil collembolan Folsomia candida as an endpoint for the estimation of organic waste ecotoxicity. Environmental Toxicology and Chemistry 2007, 26 (7), 1538-1544. 14. Barata, C.; Alanon, P.; Gutierrez-Alonso, S.; Riva, M. C.; Fernandez, C.; Tarazona, J. V., A Daphnia magna feeding bioassay as a cost effective and ecological relevant sublethal toxicity test for Environmental Risk Assessment of toxic effluents. Science of the Total Environment 2008, 405 (1-3), 78-86. 15. Hart, R. W.; Keenan, K.; Turturro, A.; Abdo, K. M.; Leakey, J.; Lyncook, B., Caloric restriction and toxicity. Fundamental and Applied Toxicology 1995, 25 (2), 184-195. 16. Keenan, K. P.; Ballam, G. C.; Soper, K. A.; Laroque, P.; Coleman, J. B.; Dixit, R., Diet, caloric restriction, and the rodent bioassay. Toxicological Sciences 1999, 52 (2), 24-34. 17. Pugh, T. D.; Klopp, R. G.; Weindruch, R., Controlling caloric consumption: protocols for rodents and rhesus monkeys. Neurobiology of Aging 1999, 20 (2), 157-165. 18. Ninan, I.; Kulkarni, S. K., Dopamine receptor sensitive effect of dizocilpine on feeding behaviour. Brain Research 1998, 812 (1-2), 157-163. 19. Videla, L. A.; Barros, S. B. M.; Junqueira, V. B. C., Lindane-induced liver oxidative stress. Free Radical Biology and Medicine 1990, 9 (2), 169-179. 20. Gliwicz, M. Z.; Sieniawska, A., Filtering activity of Daphnia in low concentrations of a pesticide. Limnology and Oceanography 1986, 31 (5), 1132-1138. 21. Hartgers, E. M.; Heugens, E. H. W.; Deneer, J. W., Effect of lindane on the clearance rate of Daphnia magna. Archives of Environmental Contamination and Toxicology 1999, 36 (4), 399404. 22. Vehmaa, A.; Hogfors, H.; Gorokhova, E.; Brutemark, A.; Holmborn, T.; Engstrom-Ost, J., Projected marine climate change: effects on copepod oxidative status and reproduction. Ecology and Evolution 2013, 3 (13), 4548-4557. 23. Barata, C.; Varo, I.; Navarro, J. C.; Arun, S.; Porte, C., Antioxidant enzyme activities and lipid peroxidation in the freshwater cladoceran Daphnia magna exposed to redox cycling compounds. Comparative Biochemistry and Physiology C-Toxicology & Pharmacology 2005, 140 (2), 175-186. 24. Nakki, R.; Nickolenko, J.; Chang, J. M.; Sagar, S. M.; Sharp, F. R., Haloperidol prevents ketamine- and phencyclidine-induced HSP70 protein expression but not microglial activation. Exp. Neurol. 1996, 137 (2), 234-241. 25. Saradha, B.; Vaithinathan, S.; Mathur, P. P., Lindane alters the levels of HSP70 and clusterin in adult rat testis. Toxicology 2008, 243 (1-2), 116-123. 26. Barata, C.; Baird, D. J.; Minarro, A.; Soares, A., Do genotype responses always converge from lethal to nonlethal toxicant exposure levels? Hypothesis tested using clones of Daphnia magna straus. Environmental Toxicology and Chemistry 2000, 19 (9), 2314-2322. 27. Smith, P. K.; Krohn, R. I.; Hermanson, G. T.; Mallia, A. K.; Gartner, F. H.; Provenzano, M. D.; Fujimoto, E. K.; Goeke, N. M.; Olson, B. J.; Klenk, D. C., Measurement of protein using bicinchonic acid. Anal. Biochem. 1985, 150 (1), 76-85. 28. Ou, B. X.; Hampsch-Woodill, M.; Prior, R. L., Development and validation of an improved oxygen radical absorbance capacity assay using fluorescein as the fluorescent probe. Journal of Agricultural and Food Chemistry 2001, 49 (10), 4619-4626. 29. Cao, G. H.; Prior, R. L., Comparison of different analytical methods for assessing total antioxidant capacity of human serum. Clinical Chemistry 1998, 44 (6), 1309-1315. 30. Hodges, D. M.; DeLong, J. M.; Forney, C. F.; Prange, R. K., Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. Planta 1999, 207 (4), 604-611. 31. Båmstedt, U.; Gifford, D. J.; Irigoien, X.; Atkinson, A.; Roman, M., 8 - Feeding. In ICES Zooplankton Methodology Manual, Huntley, R. H. W. L. R. S., Ed. Academic Press: London, 2000; pp 297-399. 32. Jeschke, J. M.; Kopp, M.; Tollrian, R., Consumer-food systems: why type I functional responses are exclusive to filter feeders. Biological Reviews 2004, 79 (2), 337-349. 17 ACS Paragon Plus Environment
Environmental Science & Technology
443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480
33. Sarnelle, O.; Wilson, A. E., Type III functional response in Daphnia. Ecology 2008, 89 (6), 1723-1732. 34. Barton, H. A.; Andersen, M. E.; Allen, B. C., Dose-response characteristics of uterine responses in rats exposed to estrogen agonists. Regulatory Toxicology and Pharmacology 1998, 28 (2), 133-149. 35. Schielzeth, H., Simple means to improve the interpretability of regression coefficients. Methods in Ecology and Evolution 2010, 1 (2), 103-113. 36. Wold, S.; Sjostrom, M.; Eriksson, L., PLS-regression: a basic tool of chemometrics. Chemometrics Intell. Lab. Syst. 2001, 58 (2), 109-130. 37. Hoppe, S.; Bierhoff, H.; Cado, I.; Weber, A.; Tiebe, M.; Grummt, I.; Voit, R., AMPactivated protein kinase adapts rRNA synthesis to cellular energy supply. Proceedings of the National Academy of Sciences of the United States of America 2009, 106 (42), 17781-17786. 38. Sohal, R. S.; Ku, H. H.; Agarwal, S.; Forster, M. J.; Lal, H., Oxidative damage, mitochondrial oxidant generation and antioxidant defense during aging and in response to food restriction in the mouse. Mechanisms of Ageing and Development 1994, 74 (1-2), 121-133. 39. Barrientos, A.; Marin, C.; Miro, O.; Casademont, J.; Gomez, M.; Nunes, V.; Tolosa, E.; Urbano-Marquez, A.; Cardellach, F., Biochemical and molecular effects of chronic haloperidol administration on brain and muscle mitochondria of rats. Journal of Neuroscience Research 1998, 53 (4), 475-481. 40. Doyotte, A.; Cossu, C.; Jacquin, M. C.; Babut, M.; Vasseur, P., Antioxidant enzymes, glutathione and lipid peroxidation as relevant biomarkers of experimental or field exposure in the gills and the digestive gland of the freshwater bivalve Unio tumidus. Aquatic Toxicology 1997, 39 (2), 93-110. 41. Pai, B. N.; Janakiramaiah, N.; Gangadhar, B. N.; Ravindranath, V., Depletion of glutathione and enhanced lipid-peroxidation in the CSF of acute psychotic following haloperidol administration. Biological Psychiatry 1994, 36 (7), 489-491. 42. Tanaka, K.; Kurihara, N.; Nakajima, M., Metabolism of lindane in house-flies metabolic desaturation, dehydrogenation and dehydrochlorination, and conjugation with glutathione. Pesticide Biochemistry and Physiology 1976, 6 (4), 392-399. 43. Anguiano, O. L.; de Castro, A. C.; de D'Angelo, A. M. P., The role of glutathion conjugation in the regulation of early toad embryos' tolerance to pesticides. Comparative Biochemistry and Physiology C-Toxicology & Pharmacology 2001, 128 (1), 35-43. 44. Kailasam, M.; Kaneko, G.; Oo, A. K. S.; Ozaki, Y.; Thirunavukkarasu, A. R.; Watabe, S., Effects of calorie restriction on the expression of manganese superoxide dismutase and catalase under oxidative stress conditions in the rotifer Brachionus plicatilis. Fisheries Science 2011, 77 (3), 403-409. 45. Rousseeuw, P. J.; Ruts, I.; Tukey, J. W., The bagplot: A bivariate boxplot. Am. Stat. 1999, 53 (4), 382-387.
481 482 483
18 ACS Paragon Plus Environment
Page 18 of 24
Page 19 of 24
Environmental Science & Technology
484
485 486 487 488 489 490
Table 1. Summary of the experiments, specific conditions and measured biomarkers
Exp. I no test substance, varying food levels
Exp. II no test substance, varying food levels
Exp. III
Exp. IV
Exp. V
Exp. VI
Haloperidol
Haloperidol
Lindane
Lindane
Concentrations of the test substance (mg L-1)
0
0
0.2-3.1
0.2-3.1
0.2-1.6
0.2-1.6
Food concentration (µg C mL-1)
0-9.4
0-7
1.5
1.5
1.5
1.5
Feeding rate, Mortality
Mortality Feeding rate (B)
Feeding rate, Mortality
Mortality Feeding rate (B)
Feeding rate, Mortality
Mortality Feeding rate (B)
Protein, ORAC
Protein, ORAC, TBARS
Protein, ORAC
Protein, ORAC, TBARS
Protein, ORAC
Protein, ORAC, TBARS
Test volume (mL)
50
900 (A) 50 (B)
50
900 (A) 50 (B)
50
900 (A) 50 (B)
Individuals/replicate
5
33-35 (A) 5 (B)
5
33-35 (A) 5 (B)
5
33-35 (A) 5 (B)
5
3 (A) 18 pooled to 3 (B)
Exposure
Physiological end points Biomarkers measured
3 (A) 3 (A) 18 pooled to 5 18 pooled to 3 3 (B) (B) A – Exposure with high population density, B – Exposure with low population density Replicates/conc.
5
19 ACS Paragon Plus Environment
Environmental Science & Technology
491
492 493 494 495
Figure 1. Functional response of juvenile D. magna.fed P. subcapitata. The horizontal broken line represents the maximum feeding rate and the vertical line represents saturation level 6.7 µg C/ml, calculated using the lower 95% confidence interval for the theoretical Fmax.
496 497 498 499 500 501 502 503
20 ACS Paragon Plus Environment
Page 20 of 24
Page 21 of 24
Environmental Science & Technology
504 505 506
Figure 2. Feeding inhibition in Daphnia magna exposed to haloperidol (A) and lindane (B). The vertical lines represent EC50 values and the grey, broken lines represent upper and lower confidence interval (95%).
507 508 Table 2. Linear regression for unexposed daphnids.
Response variable
Explanatory variable
a
b
df
R2
p
Protein
FR
3.81
5.45
23
0.70
***
ORAC
Protein
0.10
0.059
23
0.71
***
ORACp
FR
-0.014
0.12
23
0.10
˃ 0.05
TBARS
Protein
1.15
-2.99
7
0.47
*
TBARSp
FR
0.47
0.51
7
0.54
*
Linear regressions (y = ax + b) for the relationships between feeding rate and biomarker responses and between biomarkers and individual protein content. Asterisks indicate level of significance for the slope (a): p ≤ 0.05 (*); p ≤ 0.01 (**); p ≤ 0.001 (***). FR – Feeding rate, ORAC – Oxygen radical absorbance capacity, ORACp - proteinspecific ORAC, TBARS – Thiobarbituric acid reactive substances, TBARSp - protein-specific TBARS
509 510 511 512 Table 3. General linear models testing effects of the xenobiotics and feeding rate on the individual protein content, ORAC,
21 ACS Paragon Plus Environment
Environmental Science & Technology
Page 22 of 24
ORACp, TBARS and TBARSp.
Haloperidol
Lindane
Estimate
SE
t
p value
Estimate
SE
t
p value
FR
0.54
0.11
4.77
***
0.15
0.017
8.77
***
Treatment
0.16
0.043
3.80
***
-0.06
0.042
-1.42
˃ 0.05
FR*Treatment
-0.27
0.13
-2.04
*
Protein
1.7
0.15
11.45
***
0.98
0.17
5.81
***
Treatment
0.023
0.043
0.54
˃ 0.05
0.27
0.056
4.89
***
Protein*Treatment
-0.60
0.22
-2.76
**
FR
0.93
0.12
7.67
***
0.38
0.15
2.55
*
Treatment
0.11
0.046
2.48
*
0.24
0.056
4.22
***
FR*Treatment
-0.98
0.14
-7.01
***
-0.44
0.17
-2.56
*
Protein
0.51
0.10
4.89
***
0.25
0.16
1.66
˃ 0.05
Treatment
0.022
0.044
0.49
˃ 0.05
0.025
0.042
0.58
˃ 0.05
0.39
0.19
2.07
*
A) Protein
B) ORAC
C) ORACp
D) TBARS
Protein*Treatment E) TBARSp FR
-0.20
0.12
-1.72
˃ 0.05
-0.20
0.12
-1.76
˃ 0.05
Treatment
-0.063
0.099
-0.64
˃ 0.05
-0.16
0.098
-1.61
˃ 0.05
Treatment group was set as a reference. Asterisks indicate level of significance: p ≤ 0.05 (*); p ≤ 0.01 (**); p ≤ 0.001 (***).FR – Feeding rate, ORAC – Oxygen radical absorbance capacity, ORACp – protein-specific ORAC, TBARS – Thiobarbituric acid reactive substances, TBARSp - protein-specific TBARS
513 514 515
22 ACS Paragon Plus Environment
Page 23 of 24
Environmental Science & Technology
516 517 518 519 520 521 522 523 524 525 526 527
Figure 3. Standardized biplot visualizing PLS model by scatter plots termed scores or loadings plots and showing how predictors form the space of the latent variables and how they are combined with the observations. Each point on the score plot represents an individual Daphnia sample projected in the bivariate space and loading scores provide the correlation between the original variables and the new component variables. The model is for Treatment as a response variable and feeding rate (FR), individual protein content (Protein), protein-specific ORAC and TBARS values (ORACp and TBARSp, respectively) as explanatory variables. The structure of the multivariate point cloud is represented by the shaded areas in the bagplot, analogous to a box-and-whiskers plot but also visualizing the spread, correlation, skewness, and tails of the data45. The small square colored according to the treatment represents the depth median (the deepest location), the dark area corresponds to 50% of the dataset, and the light area is a fence augmented by a default factor of 1.5. Sample points outside the shaded areas are outliers.
23 ACS Paragon Plus Environment
Environmental Science & Technology
145x131mm (96 x 96 DPI)
ACS Paragon Plus Environment
Page 24 of 24