Subscriber access provided by EASTERN KENTUCKY UNIVERSITY
Environmental Modeling
Polybrominated diphenyl ether (PBDE) accumulation in farmed salmon evaluated using a dynamic sea-cage production model. Carla Ng, Amélie Ritscher, Konrad Hungerbuehler, and Natalie von Goetz Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b00146 • Publication Date (Web): 26 Apr 2018 Downloaded from http://pubs.acs.org on April 27, 2018
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 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 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.
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 36
Environmental Science & Technology
1
Polybrominated diphenyl ether (PBDE) accumulation in farmed salmon evaluated
2
using a dynamic sea-cage production model.
3
Carla A. Ng1*, Amélie Ritscher2, Konrad Hungerbuehler2, and Natalie von Goetz2*
4 5
1
6
O’Hara St, Pittsburgh, PA, USA 15261.
7
2
8
Zurich, Switzerland.
Department of Civil and Environmental Engineering, University of Pittsburgh, 3700
Institute for Chemical and Bioengineering, ETH Zurich, Vladimir-Prelog-Weg 1, 8093
9 10 11 12 13 14 15 16 17 18 19 20
*corresponding authors:
21
Carla A. Ng, +1 412 383 4075,
[email protected] 22
Natalie von Goetz, +41 44 632 0975,
[email protected] 23
ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 2 of 36
24
Abstract
25
Food is an important source of human exposure to hazardous chemicals. Chemical
26
concentration in a food item depends on local environmental contamination, production
27
conditions, and, for animal-derived foods, on feed. Here, we investigate these influences
28
on the accumulation of individual polybrominated diphenyl ether congeners (PBDEs) in
29
farmed Atlantic salmon (Salmo salar). We develop a dynamic model over a full sea-cage
30
salmon production cycle. To assess the influence of metabolic debromination on PBDE
31
congener profiles, in vitro measurements of debromination rates in fish liver cells were
32
extrapolated to whole-body metabolic rate constants. Model results indicate that the
33
dominant factors governing PBDE concentration in Atlantic salmon fillet are uptake via
34
contaminated feed and fish growth, while the influence of metabolic debromination is
35
minor. PBDE concentrations in fish feed depend on several factors, including the
36
geographic origin of fish feed ingredients, which are produced and traded globally.
37
Human exposure to PBDE via salmon consumption is less influenced by environmental
38
concentrations at the location of salmon farming than by environmental concentrations
39
influencing feed components. This dependence of PBDE concentrations in salmon on the
40
origin and composition of feed reveals the complexity of predicting contaminant
41
concentrations in globally traded food.
42 43 44 45 46
ACS Paragon Plus Environment
2
Page 3 of 36
Environmental Science & Technology
47
Introduction
48
Polybrominated diphenyl ethers (PBDEs) are a class of synthetic flame retardants which
49
have been extensively used to increase the fire resistance of consumer products, such as
50
foam paddings, textiles, and plastics1. PBDEs can be released into the environment
51
during production, use and disposal of treated products, leading to their detection in a
52
variety of environmental and biological samples2,3, including human tissues4–6. PBDEs
53
comprise 209 different with different numbers of bromines attached to their diphenyl
54
ether structure. There are three major commercial mixtures: Penta-, Octa- and Deca-BDE,
55
named for the degree of bromination of the dominant congener in the mixture. Deca-
56
BDE, the most widely used mixture,3 is added to polymeric materials such as
57
polystyrene, polybutylene, nylon, polypropylene and other thermoelastic plastics, and
58
contains more than 96% BDE-209, the highest-brominated congener. Fire-protected
59
polypropylene can contain up to 23% Deca-BDE by weight7.
60 61
Due to environmental and public health concerns, Penta- and Octa-BDE mixtures were
62
banned in Europe and the United States in 2004. Some PBDEs are disruptive of thyroid
63
hormone homeostasis and some are associated with neurodevelopmental and behavioral
64
effects in rodents8. In 2009, tetra- and pentabromodiphenyl ethers were listed under the
65
Stockholm Convention as persistent organic pollutants9. In the United States,
66
manufacturers agreed to voluntarily phase out Deca-BDE by the end of 2013, and in
67
August 2014, the European Chemicals Agency (ECHA) submitted a proposal to restrict
68
the use of Deca-BDE in the European Union10,11. Following these actions, PBDE
69
manufacturing moved to less developed countries such as China, India, Indonesia,
ACS Paragon Plus Environment
3
Environmental Science & Technology
Page 4 of 36
70
Thailand and Vietnam12. As products containing PBDE flame retardants have a relatively
71
long lifetime, emissions to the environment will most likely continue despite increasing
72
regulations. Due to their hydrophobicity and preferential partitioning into particulate
73
phases, large environmental reservoirs are being created in sediments, which might
74
present a long-term threat to biota13.
75 76
The dominant pathway for human exposure to PBDEs is consumption of contaminated
77
food8. To estimate human exposure to such potentially hazardous chemicals, modeling
78
approaches are often applied that rely on measured levels in the foods of interest14.
79
However, data on food origin is normally not addressed and, due to the global nature and
80
interconnectedness of the international food trade system, local environmental emissions
81
or regulations can influence human exposure to a contaminant on a much larger scale.
82
Although there are already human exposure models in place that account for local
83
emissions, chemical fate, and bioaccumulation, existing models do not account for the
84
global transport of chemicals via food trade15. Such consideration can be particularly
85
problematic for the consumption of animal products, as the food sources of the animals
86
themselves come into play. A particularly salient example is that of farmed predatory
87
fish. Both fish consumption and aquaculture have increased rapidly over the last forty
88
years. Salmon is one of the most widely consumed fish, and more than 50% of salmon
89
sold globally—mostly Atlantic salmon, Salmo salar—is farmed16.
90 91
Halogenated organic contaminants in aquaculture fish, including PBDEs, have gained
92
considerable attention during the last decade16–19. Most of this work has focused on
ACS Paragon Plus Environment
4
Page 5 of 36
Environmental Science & Technology
93
polychlorinated biphenyls and dioxins. However, a few studies have also investigated
94
PBDEs. Feed is generally considered to be the most important source for persistent
95
organic pollutants in farmed fish, with 35-59% of the total PBDE intake via feed
96
accumulating in the salmon fillet over 8-12 months trial periods20–22. Reducing PBDE
97
concentrations in feed can result in lower PBDE concentrations in fish fillet23,24. On the
98
other hand, Meng et al. (2008)25 highlighted the importance of the surrounding
99
environment, finding higher concentrations of PBDEs in seawater-farmed fish than in
100
freshwater-farmed fish fed with the same fish feed. They hypothesized this concentration
101
difference was caused by higher concentrations of PBDE in coastal waters from riverine
102
discharge. However, no direct comparison of the influence of feed versus environmental
103
concentrations in a single study has been made.
104 105
The role of metabolic debromination on PBDE congener profiles in fish also remains
106
unclear. Most studies have identified BDE-47 as the dominant PBDE congener in salmon
107
fillet. While congener profiles in commercially available fish feed can be very similar to
108
those measured in salmon fillet16,18, Isosaari et al. (2005)20 hypothesized that the high
109
concentrations of BDE-47 found in salmon fillet could be caused by metabolic
110
debromination of higher brominated congeners. In a feeding study of Atlantic salmon
111
with fish feed cleaned of persistent organic pollutants, Olli et al. (2010)23 found BDE-209
112
in considerable concentrations in the control group of Salmon fed with conventional fish
113
feed. This indicates that BDE-209 might have the potential to accumulate in salmon
114
fillets. However, Hites et al. (2004)16 and Montory and Barra (2006)18 found very low
115
concentrations of higher brominated congeners in salmon fillets, and BDE-209 was not
ACS Paragon Plus Environment
5
Environmental Science & Technology
116
detected, despite being present in the analyzed samples of fish feed. The authors
117
hypothesized the absence of BDE-209 in fillet samples to be caused either by low
118
bioavailability or by metabolic transformation of the substance to lower brominated
119
congeners. Neither the feeding study conducted by Isosaari et al. (2005)20 nor Berntssen
120
et al. (2010)24 included higher brominated congeners than BDE-183.
Page 6 of 36
121 122
Metabolic debromination of PBDEs in fish has been demonstrated by a number of
123
studies26,27. However, a recent in vitro assessment of metabolic PBDE debromination in
124
liver sub-fractions of common carp, rainbow trout and Chinook salmon revealed
125
considerable inter-species differences in metabolic debromination capacities28, and the
126
products of debromination reactions of individual PBDE congeners differed between
127
common carp and the other two species. It is therefore unclear to what extent the
128
metabolic debromination of PBDEs influences their concentrations in individual fish
129
species.
130 131
The aim of this work is twofold: first, to investigate how local environmental
132
concentrations and production conditions—including choices about feed composition and
133
origin—influence concentrations of PBDE in farmed salmon, and second, to consider the
134
influence of metabolic debromination on congener profiles in farmed fish. For these
135
purposes, we have developed a dynamic multi-chemical bioaccumulation model for a
136
complete sea-cage production period. This approach allows us to evaluate the
137
contributions of the most important sources, and to assess the importance of other factors,
138
such as metabolic debromination, to PBDE concentrations in Atlantic salmon fillets.
ACS Paragon Plus Environment
6
Page 7 of 36
Environmental Science & Technology
139
Finally, we discuss the implications of these factors for the transport of PBDEs via global
140
food trade.
141 142
Materials and Methods
143
Mass balance model framework
144
Our dynamic, one-compartment fugacity-based box model for PBDE uptake in an
145
individual salmon is based on the dynamic bioaccumulation model for PBDEs in lake
146
trout developed by Bhavsar et al. (2008)29. Tissue-specific properties and processes
147
related to the gills, fish body, and gastrointestinal (GI) tract are used to parameterize
148
PBDE uptake and loss to the whole body compartment (Figure S1). Details on the
149
process equations, variables, and their sources are provided in SI section S1. Here, we
150
briefly describe the overall structure of the multi-chemical mass balance.
151 152
The complete PBDE mass balance in the fish body includes uptake of a congener from
153
water and food, excretion via respiration and fecal egestion, and biotransformation29
154
expressions that describe the change in the fugacity of a congener in the fish body due to
155
metabolic reactions (including debromination (loss) of the modeled congener, y, and
156
inputs due to the debromination of its parent congener, x): ௗ
157
158
ௗ௧
=
ಷ ⋅ಷ, ⋅ௗಷ, ௗ௧
= ܧௐ,,௬ ⋅ ܦ,௬ ⋅ ݂ௐ,௬ + ܧ,,௬ ⋅ ܦ,௬ ⋅ ݂,௬ + ߥ௫௬ ⋅ ܦெ,௫ ⋅ ݂ி,௫
− ቀ൫ܧௐ,,௬ ⋅ ܦ,௬ ൯ + ൫ܧ,,௬ ⋅ ܦொ,௬ ൯ + ܦெ,௬ ቁ ⋅ ݂ி,௬
Equation 1
159 160
The change in mass of a particular PBDE congener y in the fish, dmy/dt, is given by the
ACS Paragon Plus Environment
7
Environmental Science & Technology
Page 8 of 36
161
product of the fish volume (VF), its fugacity capacity for congener y (ZF,y) and the change
162
in the fugacity of congener y in the fish as a function of time, dfF,y/dt. The product EW,0,y ·
163
DV,y · fW,y describes the uptake of congener y from water, while ED,0,y · DD,y · fD,y is the
164
uptake of congener y from food. EW,0,y and ED,0,y describe the initial chemical uptake
165
efficiencies from water and food, respectively (at time 0) and DV,y and DD,y are the
166
transport parameters, or D-values, for respiratory and dietary uptake (see SI section S1).
167
EW,0,y · DV,y · fF,y describes the loss of chemical via the gills and ED,0,y · DQ,y · fF,y is the loss
168
via fecal elimination. DM,x and DM,y are the D-values for the debromination of congeners x
169
and y, fF,x is the fugacity of congener x in the fish body, fF,y is the fugacity of the
170
debromination product y in the fish body, and vxy is the fraction of the parent congener x,
171
that is debrominated to the product y. The model is solved numerically using Matlab
172
(MATLAB 8.3.0.532, The MathWorks Inc., Natick, MA, 2014) to produce the
173
concentration of each modeled congener in each compartment as a function of time.
174 175
PBDE congeners
176
Specific PBDE congeners are considered in the model if concentration data in fish feed
177
were available (see PBDE concentrations in feed, below) or if it was part of the
178
debromination pathway of a congener present in feed (see PBDE biotransformation,
179
below).
180 181
The octanol-water partition coefficient (KOW) and Henry's law constant (H) are used to
182
calculate each congener's fugacity capacity, Z, for lipid, non-lipid organic matter, and
183
water. For all PBDE congeners used in the model, we used the KOW for homologue
ACS Paragon Plus Environment
8
Page 9 of 36
Environmental Science & Technology
184
structures estimated by Wania and Dugani (2003)30 (Table S1). Henry's law constants
185
were calculated according to the regression equation derived by Tittlemeier et al. (2002)31
186
for the sub-cooled liquid phase.
187 188
Sea-cage life cycle simulations
189
The production and life cycle of aquaculture salmon consists of several stages. Eggs are
190
stripped from broodstock fish, fertilized, and hatched in freshwater facilities. The young
191
fish are kept in freshwater tanks on land for up to 12 months until they undergo a so-
192
called parr-smolt transformation. During this transformation, the metabolic system of the
193
salmon adapts to saltwater conditions. Smolted fish are then transferred to sea sites and
194
further grown in large net cages in coastal waters for another 12 to 18 months until they
195
reach a slaughter weight of 4-5 kg. Due to deterioration of flesh quality, the salmon are
196
usually harvested before sexual maturation32. We consider the fish life cycle from the
197
time smolted fish are transferred to the sea site until harvest, corresponding to a model
198
run of 500 days. Two independent simulations were run to evaluate the influence of
199
production conditions and biotransformation. The first simulation was parameterized
200
according to the feeding study of Olli et al. (2010)23, hereafter referred to as model run A.
201
The second was parameterized according to the study of Berntssen et al. (2010)24,
202
hereafter referred to as model run B. These independent simulations were necessary to
203
reduce uncertainties in the evaluation of simulation performance, because the Olli et al.
204
study had higher resolution physiological and feeding data, but no PBDE concentrations
205
reported in feed, whereas the Berntssen et al. study had only average feeding data but
206
specified PBDE concentrations in feed.
ACS Paragon Plus Environment
9
Environmental Science & Technology
Page 10 of 36
207 208
Fish physiology and feeding rates
209
In the feeding study conducted by Olli et al. (2010)23, used as the basis for model run A,
210
the authors collected a variety of physiological parameters, including fish weight,
211
volume, proximate composition (including lipids, non-lipid organic matter, and water),
212
and respiration rates over the course of a complete seawater production cycle, starting
213
from July 2007 to December 2008. Physiological parameters derived from this study are
214
detailed in SI section S3.
215 216
Feeding rates for model run A were based on actual feeding data recorded for the Olli et
217
al. (2010)23 study (SI Figure S5). These data were generously provided to us (Harald
218
Breivik, Neperdo Biomarine, Norway, 2016, personal communication). The average daily
219
ration per fish was calculated for cages fed with conventional fish feed. Where there were
220
missing data for an individual cage, it was omitted from the calculation of average daily
221
intake. For days where there were no feeding data for any cage available, daily ration was
222
linearly interpolated from measurements from the previous and following days. There
223
were feeding data available for 479 out of 498 days of the study (96%, see section S3.7).
224
Variability in reported feeding rates was propagated to model results by calculating
225
PBDE concentrations based on the 25th percentile, 75th percentile, and mean daily rations.
226 227
For model run B, values for fish body weight were fitted to a power function based on the
228
values given by Berntssen et al. (2010). The fitted function compared to the
229
measurements by Berntssen et al. (2010) is shown in Figure S6. Berntssen et al. (2010)
ACS Paragon Plus Environment
10
Page 11 of 36
Environmental Science & Technology
230
do not provide any information on proximate composition of the fish used for the trial. It
231
was therefore assumed that the proximate composition of a fish at the beginning and the
232
end of model run B is equal to the composition of a fish in the model run A with the same
233
weight. Berntssen et al. (2010) only report daily average feed data from three
234
experimental periods, comprising 8 of the 12 months of the trial. We interpolate daily
235
feed intakes for model run B from these average daily feed intakes as described in SI
236
section S3.8.
237 238
For both model runs, comparison of predicted to measured concentrations required
239
conversion from whole body (predicted by the model) to fillet (provided in the feeding
240
studies) concentrations. This was done by multiplying lipid-based whole-body
241
concentrations by the ratio of fillet to whole-body lipid contents (see SI section S3.11).
242 243
PBDE concentrations in water
244
Few measurements of dissolved PBDE concentrations in coastal seawater are available.
245
Here, the concentrations of PBDEs in the dissolved fraction measured by Möller et al.
246
(2012)33 on three different cruises in the marine coastal environment of the North Sea
247
from March to July 2010 were averaged and used as model input (Table 1). These were
248
assumed to be representative of the concentrations found in the coastal region of a
249
Norwegian fjord (see SI section S4.1). For BDE-203, no measurements in sea- or
250
freshwater are available in literature. However, in the debromination pathway of Chinook
251
salmon28, this congener is an important “upstream” parent compound. In order to
252
investigate the influence of BDE-203 on the PBDE concentrations in the fish body
ACS Paragon Plus Environment
11
Environmental Science & Technology
Page 12 of 36
253
compartment, its concentration in seawater was estimated. In commercial PBDE
254
mixtures, BDE-203 is only present in the technical Octa-BDE mixture, where its
255
concentrations range from 5-35%. BDE-183 is also present only in the technical mixture,
256
where it constitutes 40% of all PBDE congeners10. It was assumed that the ratio of BDE-
257
183 to BDE-203 in the water is the same as the ratio found in the Octa-BDE, selecting
258
35% as an upper bound for BDE-203 (Table 1).
259 260
PBDE concentrations in feed
261
Olli et al. (2010)23 did not report the concentration of individual PBDE congeners in the
262
feed used in their study. Therefore, average concentrations reported by the Norwegian
263
public monitoring and mapping program for fish food were used as inputs for model run
264
A for the indicator PBDEs (BDE-183, -154, -153, -100, -99, -47, -28)34. Concentration
265
data for higher brominated PBDE congeners in fish feed are also scarce. Hites et al.
266
(2004)16 reported concentrations for sum of PBDE in fish feed purchased in Scotland,
267
British Columbia, eastern Canada, and Chile. In these feed samples, BDE-209 was
268
reported to average 15 ± 5% of the total PBDE measured. Based on this, the
269
concentration of BDE- 209 in fish feed was modelled as 15% of the averaged total PBDE
270
concentration reported by Hites et al. (2004)16. In order to maintain consistency between
271
the two data sources, the earliest values reported by Sanden (2014)34 (from the year 2004)
272
were used for the concentrations of all other PBDE congeners. No measurements of
273
BDE-203 in feed were available in literature. For model run A, we assumed that the ratio
274
of BDE-183 to BDE- 203 in the fish feed is the same as the ratio found in Octa-BDE
275
(Table 1). The concentrations of individual PBDE congeners in fish feed as reported by
ACS Paragon Plus Environment
12
Page 13 of 36
276
Environmental Science & Technology
Berntssen et al. (2010) were used directly as input for model run B (Table 1).
277 278
PBDE Biotransformation
279
The metabolic degradation of PBDE in fish is generally considered to occur via reductive
280
debromination28,35–39. Roberts et al. (2011)28 examined the nature and amounts of
281
metabolites of in vitro debromination of various BDE congeners (28, 47, 49, 99, 100,
282
153, 154, 183, 203, 208 and 209) by microsomal fractions of adult rainbow trout
283
(Oncorhynchus mykiss), adult carp (Cyprinus carpio) and juvenile Chinook salmon
284
(Oncorhynchus tschwatcha). Based on the results of individual incubations of BDE
285
congeners, the authors constructed a metabolic pathway for each fish species. These
286
pathways were markedly different between the species, with carp liver fractions having
287
significantly higher metabolic capacities, both in terms of reaction rate and
288
debromination products formed, compared to rainbow trout and Chinook salmon. The
289
debromination pathways found for rainbow trout and Chinook salmon are shown in
290
Figure S8. Based on the amount of debromination products after 24 hours incubation
291
time, Roberts et al. (2011)28 calculated in vitro metabolite formation rates. Due to the
292
lack of any further data on rate constants for metabolic reductive debromination in fish,
293
we extrapolated these product formation rates to whole-body, in vivo first-order
294
biotransformation rate constants (kM) for salmon using a four-step process: (1) intrinsic
295
hepatic in vitro clearance rates were calculated from the depletion rate of the parent
296
metabolite in Chinook salmon28; (2) this in vitro clearance was extrapolated to in vivo
297
clearance using the average measured hepatosomatic index (HSI) of the fish in the Olli et
298
al. (2010)23 study and the amount of microsomal protein in the liver estimated for trout by
ACS Paragon Plus Environment
13
Environmental Science & Technology
Page 14 of 36
299
Nichols et al. (2006)40; (3) the intrinsic hepatic in vivo clearance rate was extrapolated to
300
total hepatic clearance based on the rate of liver blood flow, fish weight, and the blood-
301
water partition coefficient of individual PBDE congeners; (4) finally, the whole body
302
biotransformation rate constant (kM) was determined using the volume of distribution of
303
the chemical between the blood and other tissues. Details of these four steps, including all
304
equations and parameters used, are given in SI section S5.
305 306
In order to investigate the influence of debromination on model performance, each
307
feeding study was run with and without biotransformation. Simulations run with
308
biotransformation are referred to as model run A1 and B1 for the Olli et al. (2010)23 and
309
Berntssen et al. (2010)24 studies, respectively. Model runs without biotransformation, in
310
which kM for all congeners is set to zero, are referred to as A2 and B2.
311 312
Sensitivity and Uncertainty
313
Many parameters used in the model are based on estimated values and therefore are
314
associated with considerable uncertainty. In order to understand the influence of this
315
uncertainty on model predictions, model sensitivity to all input parameters was
316
investigated using a factor-at-a-time approach. Each model parameter was varied by 10%
317
of its value, while all other parameters were kept at their nominal value. Normalized
318
sensitivity coefficients for each parameter (ܰܵܥ ) were calculated as follows (MacLeod
319
et al., 2002):
320
ܰܵܥ =
ఏିఏబ ିబ
⋅
బ
ఏబ
Equation 2
ACS Paragon Plus Environment
14
Page 15 of 36
Environmental Science & Technology
321
where ܲ is the varied parameter value, ܲ is the nominal value of the parameter used in
322
the model, ߠ is the nominal model output and ߠ is the model output calculated with the
323
varied parameter. As some of the congeners are debrominated metabolically or formed
324
via metabolic debromination and other congeners are solely taken up via food and water,
325
the model structure for these congeners is different. Therefore the sensitivity analysis was
326
carried out for each individual congener considered in the mass balance. The endpoint
327
evaluated was the concentration of each congener in the fish body at the end of a 500-day
328
production period, based on model run A1.
329 330
Results
331
In the following sections, the concentration dynamics and PBDE fluxes predicted by the
332
model are first illustrated for selected congeners using model run A1. Predicted
333
concentrations and congener distributions are then compared to available data23,24 for all
334
four model runs.
335
Concentration Dynamics
336
The modeled concentrations of different PBDEs in the fish body predicted by model run
337
A1 span more than three orders of magnitude. Most congeners show similar dynamics
338
with time. In Figure 1, we illustrate these concentration dynamics with the results for
339
model run A1 for congeners BDE-47, 99 and 209.
340 341
The concentration of BDE-47 predicted by the model is one order of magnitude higher
342
than BDE-99 and -209, which are similar. For these and most other congeners,
343
concentrations increase steeply for the first 100 days of production, then decrease again.
ACS Paragon Plus Environment
15
Environmental Science & Technology
Page 16 of 36
344
Around day 250 of the production period, most congeners reach a quasi-steady state and
345
their concentrations remain constant for the next 100 days. Concentrations then increase
346
again for the last 150 days of production. Differences among congeners are largely driven
347
by differences in uptake and loss fluxes, as discussed below, which in turn are strongly
348
dependent on the KOW, which affects dietary and gill uptake efficiencies as well as
349
elimination rates. The shape of the concentration curve with time is largely driven by the
350
balance of changing consumption rates with fish size, growth dilution, temperature
351
(which seasonally affects rates of consumption), and metabolic capacity. Of these, the
352
influence of size and growth is particularly strong. Small fish have relatively high
353
consumption rates, leading to rapid uptake at the beginning of the modeling period.
354
Dietary uptake (absolute, not normalized to size) increases with size. The highest growth
355
period for the modeled fish begins around day 125, and the dilution associated with
356
strong growth drives the decline between days 100 and 200, before approaching steady
357
state.
358 359
Input and Output Fluxes
360
Figure 2 shows PBDE mass fluxes into and out of the fish body compartment for each
361
day of the production period for BDE-47, -99 and -209, based on model run A1 (see SI
362
Figures S11 to S21 for remaining congeners). For all congeners except BDE-187, -149, -
363
101, -66 and -49, the dominant input flux of BDE into the fish body is dietary uptake.
364
The variability in modeled fluxes corresponds to the scatter of simulated daily feed intake
365
(see Figure S5). Only BDE-209 is also taken up via respiration. However, compared to
366
dietary intake, this mass flux is small.
ACS Paragon Plus Environment
16
Page 17 of 36
Environmental Science & Technology
367 368
The importance of different output fluxes is more variable: for most metabolically
369
inactive congeners (BDE-209, -187, -154, -101, -100, -66, -49 and -47) the dominant
370
output flux from the fish body compartment is fecal egestion. For the lowest brominated
371
(and therefore most water soluble) congener, BDE-28, loss via respiration is equal to loss
372
via fecal egestion. For all congeners that undergo metabolic debromination (BDE-203, -
373
183, -153 and -99), this process is the dominant loss flux.
374 375
Model Evaluation
376
To evaluate our model, predicted concentrations of individual PBDE congeners in the
377
simulated salmon fillet at the end of the production cycle were compared to measured
378
concentrations (Figure 3). Olli et al. (2010)23 were not able to analytically separate BDE-
379
49 and -71, therefore the cumulative concentration of these congeners is shown in Figure
380
3A. They did not measure BDE-101, -149, -187 and -203, therefore no measured values
381
for these congeners are available. Both model runs A1 and A2 underestimate the total
382
PBDE concentration in the fish fillet, but by less than a factor of 2, except for BDE-66,
383
which is underestimated by a factor of 45. In both model runs, the dominant congener in
384
the simulated fish fillet is BDE-47, and predicted concentrations correspond well with
385
measurements. The modeled concentrations of BDE-149, -183, -187 and -203, which
386
were not measured, were very low at the end of the simulated production period. Model
387
run A1 in particular underestimates the concentrations of BDE-99 and -209, though in
388
both cases predictions are within a factor of 2 of the measured values.
389
ACS Paragon Plus Environment
17
Environmental Science & Technology
Page 18 of 36
390
In model run A2, without biotransformation, the calculated concentrations of the parent
391
congeners BDE-99 and BDE-153 are within 2% and 15%, respectively, of the values
392
measured by Olli et al. (2010)23. The predicted concentrations of the product congeners
393
BDE-49, -66, -101, -149 and -187, whose concentration in feed is unknown, and for
394
which inputs are only via debromination, are consequently 0. The significantly better
395
predictions for BDE-99 and -153 by model run A2 suggest that kM derived from in vitro
396
experiments for these congeners is too high, while the substantial underestimation for
397
BDE-66 (21 times lower than the measurements for model run A1 and predicted to be
398
zero for model run A2) suggests an as yet unaccounted for source in the feed, including
399
debromination of other congeners. Predictions for these simulations are based on
400
assumed concentrations in water and feed (shown in Table 1), as Olli et al. (2010) did not
401
specifically report these concentrations. Neither the water nor the feed concentrations
402
profiles used included BDE-66. Yet substantial BDE-66 was found in salmon fillet by
403
Olli et al. (2010), indicating that BDE-66 was indeed present in water and/or feed.
404
Alternatively, it is possible that an additional pathway for formation of BDE-66 via
405
debromination of other congeners may exist, or that the efficiency of transformation from
406
BDE-99 is substantially higher than estimated in the debromination pathways used (see
407
SI figure S8). Here the simulation of the Berntssen et al. (2010) study (model runs B1 and
408
B2, Figure 3B) provides additional insight, as they specifically reported measured PBDE
409
concentrations in the feed.
410 411
In contrast with simulations A1 and A2 for the Olli et al (2010) feeding study, model runs
412
B1 and B2 generally overestimate the concentrations of BDE-47, -100, -66 and -154
ACS Paragon Plus Environment
18
Page 19 of 36
Environmental Science & Technology
413
measured by Berntssen et al. (2010)24. However, predicted concentrations are again well
414
within a factor of 2 of measured values, except for BDE-99 for model run B1. We again
415
see discrepancies driven by biotransformation. While the concentrations of parent
416
congeners BDE-99 and -153 are underestimated by model run B1, model run B2
417
overestimates their fillet concentrations. This suggests biotransformation of these
418
congeners does occur, but at a rate that is somewhat lower than those based on our in
419
vitro extrapolation method.
420 421
Sensitivity and Uncertainty
422
Based on our sensitivity analysis, model sensitivity for specific BDE congeners can be
423
categorized into two groups. The first group consists of all congeners that are not
424
metabolically active (BDE-209, BDE-100, BDE-47 and BDE-28, see 3.7). The
425
normalized sensitivity coefficients (NSC) calculated for a 10% increase in parameter
426
values for BDE-209, representative for the first group of congeners, are shown in Figure
427
S9 for model run A1. The model output is most sensitive to the congener concentration in
428
feed, the fish weight, the daily feed ration, and the KOW, in decreasing order. The feed
429
concentration is particularly critical, because the estimated concentrations of PBDE in
430
feed for model runs A1 and A2 are among the most uncertain parameters. For KOW, an
431
increase in the parameter leads to a decrease in dietary uptake efficiency and therefore in
432
the final PBDE concentration in the body. The model is relatively insensitive to
433
parameters governing uptake and elimination via the gills.
434 435
The second group is of those congeners influenced by metabolic debromination. The
ACS Paragon Plus Environment
19
Environmental Science & Technology
Page 20 of 36
436
sensitivities for BDE-99, which are representative for this group, are shown in SI Figure
437
S10. Here the most important parameters are the concentration in feed, daily feed ration,
438
and fish weight. These three all have NSC with absolute values between 0.8 and 1. The
439
next most important parameters, with NSC > ±0.5, include those associated with
440
metabolism (in vitro microsomal protein content, in vitro assay volume, in vitro
441
degradation rate, microsomal protein content, hepatosomatic index, and fractional water
442
content of blood), and lipid assimilation efficiency, all with similar sensitivity indices.
443
Similarly to the congeners of the first group, the model is relatively insensitive to
444
parameters governing uptake and elimination via the gills. However for this group of
445
congeners, in contrast with the first, the model is also relatively less sensitive to KOW.
446 447
We use three parameters to illustrate how the impacts of uncertainty can be evaluated
448
using our sensitivity analysis. It has been previously shown that high uncertainties are
449
associated with the KOW, solubility, and vapor pressures30,31,41of PBDEs. Our sensitivity
450
analysis shows that the model is insensitive to both solubility and the Henry’s law
451
coefficient (which includes both solubility and vapor pressure contributions). Thus high
452
uncertainty in their values would not change the predictions of the model. For KOW, on
453
the other hand, the model is relatively sensitive (NSCKOW = − 0.73), for the first group of
454
congeners (represented by BDE-209, see figure S9), and less so (NSCKOW = − 0.19) for
455
the second group of congeners (represented by BDE-99, see figure S10). In both cases,
456
increasing the KOW will decrease the predicted whole body concentration due to
457
decreasing dietary uptake efficiency (see Figure S7). Assuming a one order of magnitude
458
difference in KOW, we can use Equation 2 to calculate a decrease in predicted PBDE
ACS Paragon Plus Environment
20
Page 21 of 36
459
Environmental Science & Technology
concentration of a factor of 2 at most.
460 461
As our model exhibits certain non-linear characteristics, these normalized sensitivity
462
coefficients should be treated with caution. More elaborate methods, such as Monte
463
Carlo-based global sensitivity and uncertainty analysis, might be needed to fully assess
464
the model sensitivity over a wider range of parameter values.
465 466 467
Discussion
468
The concentrations of PBDEs in salmon fillet are well predicted by our model—typically
469
within a factor of 2—despite substantial uncertainties associated with the input
470
parameters. Nevertheless, differences in model performance with and without
471
debromination can provide some useful insights. For model runs A1 and A2,
472
underestimates for specific PBDE congeners could result from assuming concentrations
473
in the modeled feed that are too low compared to those used in the feed by Olli et al.
474
(2010)23. This is most likely the case for BDE-66, where the modeled final concentration
475
in fillet is more than 45 times lower than the measured values. As Sanden (2014)34, which
476
we used to parameterize feed concentrations for these model runs, did not measure BDE-
477
66, we set the model feed concentration to 0 for this congener. Thus, the only input flux
478
for BDE-66 stemmed from metabolic debromination of BDE-99. Measurements by
479
Isosaari et al. (2005)20 and Berntssen et al. (2010)24 indicate, however, that BDE-66 is
480
also found in fish oil and meal used for the production of fish feed.
481
ACS Paragon Plus Environment
21
Environmental Science & Technology
Page 22 of 36
482
A second factor that could cause underestimation of the final PBDE concentrations for
483
some congeners is that extrapolated full-body metabolic rate constants for debromination
484
are too high. This would result in an overestimation of the loss of parent congeners,
485
which could explain the low fillet concentrations of BDE-99 and BDE-153 predicted by
486
model run A1, and is further supported by the overestimation of BDE-49 (the
487
debromination product of BDE-99). The systematic underestimation of the parent
488
congeners BDE-153 and -99 by model run B1, and the fact that the concentration profile
489
for model run B2 is more similar to Berntssen et al. (2010)24, also supports this
490
hypothesis.
491 492
The model we present here supports dietary uptake as the dominant pathway for PBDE
493
into salmon compared to the uptake via gills from local environmental contamination for
494
all congeners. Previous empirical work has highlighted the importance of dietary
495
uptake16–20. Our model shows that this is due to PBDE concentrations in fish feed being
496
significantly higher compared to environmental concentrations in the seawater. As Meng
497
et al. (2008)25 suggested, however, uptake from water could be important for fish farmed
498
in waters with higher contaminant loads.
499 500
In all model runs BDE-47 is calculated to be the dominant congener in the fish fillet, in
501
good agreement with previous studies16,19,20,23,42,43. Isosaari et al. (2005)20 hypothesized
502
the dominance of BDE-47 to be caused by the debromination of higher brominated
503
congeners and the stability of BDE-47 towards further debromination reactions.
504
According to Roberts et al. (2011)28, however, BDE-47 is not part of the debromination
ACS Paragon Plus Environment
22
Page 23 of 36
Environmental Science & Technology
505
pathway of either Chinook salmon or rainbow trout (both salmonids). While it is unclear
506
to what extent the debromination pathway for Chinook salmon also holds for Atlantic
507
salmon, our model reproduced the measured fillet concentrations of BDE-47 rather well
508
when only intake via dietary uptake and water were considered. This supports the
509
hypothesis that BDE-47 in salmonids is not part of the debromination pathway, and that
510
high fillet concentrations in Atlantic salmon fillet are likely from feed.
511 512
For BDE-209, the influence of metabolic debromination on the concentration in the fish
513
fillet is difficult to assess. Previous observations of BDE-209 accumulation in Salmo
514
salar have been inconclusive. Some authors hypothesize that low concentrations of BDE-
515
209 found in salmon fillet are due to metabolic debromination to lower brominated
516
congeners16,18,43. These findings stand in contrast to measurements by Olli et al. (2010)23,
517
which indicate that uptake of BDE-209 in Atlantic salmon from feed is possible and that
518
the metabolic potential of the salmon was not sufficient to fully debrominate the BDE-
519
209. In contrast to rainbow trout and common carp, Roberts et al. (2011)28 found no
520
degradation of BDE-209 in liver fractions of Chinook salmon. To conclusively assess the
521
metabolic potential of Atlantic salmon for BDE-209, an exposure trial needs to be
522
conducted. Our model assumes that Atlantic salmon, like Chinook salmon, do not
523
debrominate BDE-209, and demonstrates that concentrations in the range of the measured
524
values by Olli et al. (2010)23 could be reached by exposure via feed and water. A recent
525
study44 found highly variable concentrations of PBDE in fish meal, fish oil and fish feed
526
of various European origins, ranging from total PBDE concentrations below the limit of
527
detection up to 2.2 µg/kg. Additionally, the congener profiles of the examined samples
ACS Paragon Plus Environment
23
Environmental Science & Technology
Page 24 of 36
528
varied markedly. While the authors found no BDE-209 in the two samples of fish oil and
529
complete fish feed, 5 out of 10 samples of fish meal contained BDE-209 concentrations
530
varying between 1 and 15 µg/kg lipid. Assuming a fish meal content of 35%23, fish feed
531
would contain between 0.07 and 0.52 µg/kg BDE-209, which is comparable to the
532
concentrations found by Hites et al. (2004)16.
533 534
Our model has the following limitations. First, it is a single-compartment model, which
535
does not account for the distribution of PBDEs to different organs, such as the liver or
536
muscle. For some higher brominated PBDE congeners, preferential distribution to liver
537
tissues has been observed for wild common dab (Limanda limanda), whiting (Merlangius
538
merlangus) and pouting (Trisopterus luscus)45. Furthermore, the preferential
539
accumulation of BDE-209 has been observed in the liver of juvenile rainbow trout
540
(Oncorhynchus mykiss)37 and Chinese sturgeon (Acipenser sinensis)46. The reason for
541
preferential distribution of these congeners to liver is currently unclear. However, it could
542
cause PBDE concentrations of higher brominated PBDE congeners in fillet to be lower
543
than calculated by the model. Modeling the fish body as a single compartment further
544
neglects the effects of differential organ growth and lipid accumulation. Berntssen et al.
545
(2005)21 found muscle fat deposition in Atlantic salmon to rapidly increase compared to
546
the rest of the body when the fish reached a weight of approximately 1 kg. Such
547
differences in relative growth rates of the lipid compartment could significantly influence
548
the uptake and distribution kinetics of PBDE in Atlantic salmon, which is primarily lipid
549
driven. A third limitation of the model is the missing link between temperature, salmon
550
growth and feed intake. It has generally been observed that salmon feed intake in
ACS Paragon Plus Environment
24
Page 25 of 36
Environmental Science & Technology
551
aquaculture conditions decreases with temperature47. However, the growth rate of the fish
552
does not necessarily change accordingly48. This can also be observed in the feeding data
553
used as input for the model: during cold water periods, feed intake decreases, while the
554
growth rate remains approximately constant. Furthermore, in warmer months of the
555
production periods, feed intake increased, while growth rates did not increase by the
556
same order of magnitude. By linking temperature, feed intake and growth of farmed
557
Atlantic salmon, the PBDE fillet concentrations could be predicted independent of a
558
specific aquaculture scenario. This could be achieved by implementing a bioenergetic
559
model for aquaculture salmon into the bioaccumulation model.
560 561
If diet is the most important uptake route for PBDEs into salmon, the final concentrations
562
reached in the fillet are strongly influenced by the ratio of growth and feeding rate. This
563
so-called feed-conversion ratio (i.e. the ratio of total amount of weight gained to amount
564
of feed fed) is a parameter closely monitored in aquaculture practices and strongly
565
depends on the salmon strain and animal husbandry. By using strains of salmon with a
566
maximum feed-conversion ratio (which promotes growth dilution) and by optimizing
567
farming practices, the concentration of PBDEs in fillet of farmed salmon can be reduced.
568
The consumption of salmon raised under such optimal production conditions would
569
therefore lead to a lower exposure to PBDEs.
570 571
The limited influences of local sources of pollution on the concentration of PBDEs in
572
farmed salmon has implications for how hazardous chemicals are transported via food
573
trade15. In the case of farmed salmon, it is not necessarily the origin of the farmed fish but
ACS Paragon Plus Environment
25
Environmental Science & Technology
Page 26 of 36
574
of the feed or feed ingredients used in its production that must be taken into account to
575
accurately trace the origins of human dietary exposure to PBDEs. Fish oils in particular
576
have been found to contain high concentrations of persistent organic pollutants, promptim
577
decontamination studies such as those by Olli et al. (2010) and Berntssen et al.
578
(2010).23,48 The model we developed in this work seeks to make explicit the competing
579
contributions of feed concentrations, water concentrations, biotransformation, and fish
580
physiology (including also fish growth and feed conversion efficiency). By making these
581
links, the model could be extended and modified to assess the influence of environmental
582
and production-specific conditions of other fish with high global trading volumes, such as
583
Tilapia or red snapper. Particular attention should be given to food species or feeds
584
produced in the contamination "hot-spots" of densely populated areas. The model could
585
similarly be used to evaluate contaminant control strategies such as replacing fish oils
586
with plant-based materials in feed or decontamination of fish oil.
587 588
The model developed here performed well in predicting congener-specific PBDE
589
concentrations in farmed salmon. The major strength of the model is its ability to
590
simultaneously treat multiple congeners in a dynamic way. This includes the tracking of
591
lower brominated congeners that accumulate from direct uptake and the debromination of
592
higher brominated congeners. Moreover, the mechanistic approach of the model allows
593
us to discern the relative contributions of different uptake routes, and notably the
594
influence of feed. This model framework could be readily adapted to other metabolizable
595
compounds and their transformation products.
596
ACS Paragon Plus Environment
26
Page 27 of 36
597
Environmental Science & Technology
Table 1. Concentrations of PBDEs in Water and Feed Used as Model Input. Congener
Water
Feed Concentration
Concentration
[µg /kg]
[mol / m3]
A1, A2
B1, B2
BDE-209
4.36×10-13
4.72
0
BDE-203
1.93×10-14
0.026
0
BDE-183
2.21×10-14
0.03
0
BDE-154
0
0.19
0.78
BDE-153
1.31×10-14
0.08
0.24
BDE-100
1.97×10-14
0.35
0.87
BDE-99
2.37×10-13
0.31
0.93
BDE-66
0
0
0.31
BDE-47
3.33×10-13
2.07
4.5
BDE-28
0
0.09
0.18
598 599
FIGURE CAPTIONS
600 601
Figure 1. Concentration of BDE-47, -99, and -209 in the fish whole-body compartment
602
over an entire sea-cage cycle for model run A1.
603 604 605
Figure 2. Daily mass intake and loss of (A) BDE-47, (B) BDE-99, and (C) BDE-209 into
606
farmed Atlantic salmon (whole body) over 500-day sea-cage growth cycle (model run
ACS Paragon Plus Environment
27
Environmental Science & Technology
607
Page 28 of 36
A1). Scatter in predicted dietary uptake reflects variability in daily feed intake.
608 609
Figure 3. Comparisons between (A) Olli et al. (2010) data and model runs A1 and A2,
610
with and without biotransformation, respectively, and (B) Berntseen et al. (2010) data
611
and model runs B1 and B2, with and without biotransformation, respectively. Error bars
612
for Olli et al. and Berntssen et al. data indicate reported variability. Error bars in model
613
run A1 indicate predictions using the 25th and 75th percentile of reported feeding rates.
614
Stars indicate that the congener was not measured in that study.
615 616
ACKNOWLEDGMENT
617
We thank Harald Breivik of Neperdo Biomarine, Norway, for kindly providing the
618
feeding rate data used in this study.
619 620
ASSOCIATED CONTENT
621
Supporting Information
622
Details on model parameterization and sensitivity analysis are available in the Supporting
623
Information (SI). This material is available free of charge via the Internet at
624
http://pubs.acs.org/.
625 626
AUTHOR INFORMATION
627
Corresponding Authors
628
Carla A. Ng, +1 412 383 4075,
[email protected] 629
Natalie von Goetz, +41 44 632 0975,
[email protected] ACS Paragon Plus Environment
28
Page 29 of 36
Environmental Science & Technology
630 631
Notes
632
The authors declare no competing financial interest.
633 634
References Cited
635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670
(1) (2)
(3) (4) (5)
(6) (7) (8) (9) (10) (11) (12)
(13)
(14)
(15)
Hites, R. A. Polybrominated Diphenyl Ethers in the Environment and in People: A Meta-Analysis of Concentrations. Environ. Sci. Technol. 2004, 38 (4), 945–956. Hale, R. C.; La Guardia, M. J.; Harvey, E.; Mainor, T. M. Potential Role of Fire Retardent-Treated Polyurethan Foam as a Source of Brominated Diphenyl Ethers to the US Environment. Chemosphere 2002, 46, 729–735. de Wit, C. A. An Overview of Brominated Flame Retardants in the Environment. Chemosphere 2002, 46 (5), 583–624. Watanabe, I.; Sakai, S. Environmental Release and Behavior of Brominated Flame Retardants. Environ. Int. 2003, 29 (6), 665–682. Alaee, M.; Arias, P.; Sjödin, A.; Bergman, Åke. An Overview of Commercially Used Brominated Flame Retardants, Their Applications, Their Use Patterns in Different Countries/Regions and Possible Modes of Release. Environ. Int. 2003, 29 (6), 683–689. Wang, Y.; Jiang, G.; Lam, P. K. S.; Li, A. Polybrominated Diphenyl Ether in the East Asian Environment: A Critical Review. Environ. Int. 2007, 33, 963–973. Centers for Disease Control and Prevention. Biomonitoring Summary Polybrominated Diphenyl Ethers; 2013. European Food Safety Authority. Scientific Opinion on Polybrominated Diphenyl Ethers (PBDEs) in Food. EFSA J. 2011, 9 (5), 1–274. UNEP. The 9 New POPs; 2010. US EPA. Action Plan Fact Sheet April 2011. 2011. ECHA. ECHA Proposes a Restriction on DecaBDE , a Brominated Flame Retardant Used in Plastics and Textiles. 2014, 69 (6), 1–4. Li, Y.; Lin, T.; Hu, L.; Feng, J.; Guo, Z. Time Trends of Polybrominated Diphenyl Ethers in East China Seas: Response to the Booming of PBDE Pollution Industry in China. Environ. Int. 2016, 92–93, 507–514. Ross, P. S.; Couillard, C. M.; Ikonomou, M. G.; Johannessen, S. C.; Lebeuf, M.; Macdonald, R. W.; Tomy, G. T. Large and Growing Environmental Reservoirs of Deca-BDE Present an Emerging Health Risk for Fish and Marine Mammals. Mar. Pollut. Bull. 2009, 58 (1), 7–10. Trudel, D.; Tlustos, C.; Von Goetz, N.; Scheringer, M.; Hungerbühler, K. PBDE Exposure from Food in Ireland: Optimising Data Exploitation in Probabilistic Exposure Modelling. J. Expo. Sci. Environ. Epidemiol. 2011, 21 (6), 565–575. Ng, C. A.; von Goetz, N. The Global Food System as a Transport Pathway for Hazardous Chemicals: The Missing Link between Emissions and Exposure. Environ. Health Perspect. 2017, 125 (1), 1–7.
ACS Paragon Plus Environment
29
Environmental Science & Technology
671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714
Page 30 of 36
(16) Hites, R. A.; Foran, J. A.; Schwager, S. J.; Knuth, B. A.; Hamilton, M. C.; Carpenter, D. O. Global Assessment of Polybrominated Diphenyl Ethers in Farmed and Wild Salmon. Environ. Sci. Technol. 2004, 38 (19), 4945–4949. (17) Jacobs, M. N.; Covaci, A.; Schepens, P. Investigation of Selected Persistent Organic Pollutants in Farmed Atlantic Salmon (Salmo Salar), Salmon Aquaculture Feed, and Fish Oil Components of the Feed. Environ. Sci. Technol. 2002, 36 (13), 2797–805. (18) Montory, M.; Barra, R. Preliminary Data on Polybrominated Diphenyl Ethers (PBDEs) in Farmed Fish Tissues (Salmo Salar) and Fish Feed in Southern Chile. Chemosphere 2006, 63 (8), 1252–1260. (19) Shaw, S. D.; Berger, M. L.; Brenner, D.; Carpenter, D. O.; Tao, L.; Hong, C. S.; Kannan, K. Polybrominated Diphenyl Ethers (PBDEs) in Farmed and Wild Salmon Marketed in the Northeastern United States. Chemosphere 2008, 71 (8), 1422–1431. (20) Isosaari, P.; Lundebye, A.; Ritchie, G.; Lie, Ø.; Kiviranta, H.; Vartiainen, T.; Lundebye, A.; Ritchie, G.; Lie, Ø.; Kiviranta, H.; et al. Dietary Accumulation Efficiencies and Biotransformation of Polybrominated Diphenyl Ethers in Farmed Atlantic Salmon (Salmo Salar). Food Addit. Contam. Part A 2005, 22 (9), 829– 837. (21) Berntssen, M. H. G.; Lundebye, A. K.; Torstensen, B. E. Reducing the Levels of Dioxins and Dioxin-like PCBs in Farmed Atlantic Salmon by Substitution of Fish Oil with Vegetable Oil in the Feed. Aquac. Nutr. 2005, 11 (3), 219–231. (22) Berntssen, M. H. G.; Valdersnes, S.; Rosenlund, G.; Torstensen, B. E.; Zeilmaker, M. J.; van Eijkeren, J. C. H. Toxicokinetics and Carry-over Model of AlphaHexabromocyclododecane (HBCD) from Feed to Consumption-Sized Atlantic Salmon (Salmo Salar). Food Addit. Contam. Part -Chem. Anal. Control Expo. Risk Assess. 2011, 28 (9), 1274–1286. (23) Olli, J. J.; Breivik, H.; Mørkøre, T.; Ruyter, B.; Johansen, J.; Reynolds, P.; Thorstad, O.; Berge, G. Removal of Persistent Organic Pollutants from Atlantic Salmon (Salmo Salar L.) Diets : Influence on Growth , Feed Utilization Efficiency and Product Quality. Aquaculture 2010, 310 (1–2), 145–155. (24) Berntssen, M. H. G.; Julshamn, K.; Lundebye, A. K. Chemical Contaminants in Aquafeeds and Atlantic Salmon (Salmo Salar) Following the Use of Traditionalversus Alternative Feed Ingredients. Chemosphere 2010, 78 (6), 637–646. (25) Meng, X.-Z.; Yu, L.; Guo, Y.; Mai, B.-X.; Zeng, E. Y. Congener-Specific Distribution of Polybrominated Diphenyl Ethers in Fish of China: Implication for Input Sources. Environ. Toxicol. Chem. SETAC 2008, 27 (1), 67–72. (26) Stapleton, H. M.; Alaee, M.; Letcher, R. J.; Baker, J. E. Debromination of the Flame Retardant Decabromodiphenyl Ether by Juvenile Carp (Cyprinus Carpio) Following Dietary Exposure. Environ. Sci. Technol. 2004, 38 (1), 112–119. (27) Tomy, G. T.; Palace, V. P.; Halldorson, T.; Braekevelt, E.; Danell, R.; Wautier, K.; Evans, B.; Brinkworth, L.; Fisk, A. T. Bioaccumulation, Biotransformation, and Biochemical Effects of Brominated Diphenyl Ethers in Juvenile Lake Trout (Salvelinus Namaycush). Environ. Sci. Technol. 2004, 38 (5), 1496–1504.
ACS Paragon Plus Environment
30
Page 31 of 36
715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760
Environmental Science & Technology
(28) Roberts, S. C.; Noyes, P. D.; Gallagher, E. P.; Stapleton, H. M. Species-Specific Differences and Structure-Activity Relationships in the Debromination of PBDE Congeners in Three Fish Species. Environ. Sci. Technol. 2011, 45 (5), 1999–2005. (29) Bhavsar, S. P.; Gandhi, N.; Gewurtz, S. B.; Tomy, G. T. Fate of PBDEs in Juvenile Lake Trout Estimated Using a Dynamic Multichemical Fish Model. Chem. Eng. 2008, 42 (10), 717010. (30) Wania, F.; Dugani, C. B. Assessing the Long-Range Transport Potential of Polybrominated Diphenyl Ethers: A Comparison of Four Multimedia Models. Environ. Toxicol. Chem. 2003, 22 (6), 1252–1261. (31) Tittlemier, S. a; Halldorson, T.; Stern, G. a; Tomy, G. T. Vapor Pressures, Aqueous Solubilities, and Henry’s Law Constants of Some Brominated Flame Retardants. Environ. Toxicol. Chem. SETAC 2002, 21 (9), 1804–1810. (32) Principles of Salmonid Culture; Pennell, W., Barton, B. A., Eds.; Elsevier, 1996. (33) Möller, A.; Xie, Z.; Caba, A.; Sturm, R.; Ebinghaus, R. Occurrence and AirSeawater Exchange of Brominated Flame Retardants and Dechlorane Plus in the North Sea. Atmos. Environ. 2012, 46, 346–353. (34) Sanden, M. Program for Overv\a Aking Av Fiskefôr \AArsrapport 2013; 2014. (35) Stapleton, H. M.; Alaee, M.; Letcher, R. J.; Baker, J. E. Debromination of the Flame Retardant Decabromodiphenyl Ether by Juvenile Carp (Cyprinus Carpio) Following Dietary Exposure. Environ. Sci. Technol. 2004, 38 (1), 112–119. (36) Stapleton, H. M.; Letcher, R. J.; Li, J.; Baker, J. E. Dietary Accumulation and Metabolism of Polybrominated Diphenyl Ethers by Juvenile Carp (Cyprinus Carpio). Env. Toxicol Chem 2004, 23 (8), 1939–1946. (37) Stapleton, H. M.; Brazil, B.; Holbrook, R. D.; Mitchelmore, C. L.; Benedict, R.; Konstantinov, A.; Potter, D. In Vivo and in Vitro Debromination of Decabromodiphenyl Ether (BDE 209) by Juvenile Rainbow Trout and Common Carp. Environ. Sci. Technol. 2006, 40 (15), 4653–4658. (38) Browne, E. P.; Stapleton, H. M.; Kelly, S. M.; Tilton, S. C.; Gallagher, E. P. In Vitro Hepatic Metabolism of 2,2’,4,4’,5-Pentabromodiphenyl Ether (BDE 99) in Chinook Salmon (Onchorhynchus Tshawytscha). Aquat. Toxicol. 2009, 92 (4), 281–287. (39) Noyes, P. D.; Kelly, S. M.; Mitchelmore, C. L.; Stapleton, H. M. Characterizing the in Vitro Hepatic Biotransformation of the Flame Retardant BDE 99 by Common Carp. Aquat. Toxicol. 2010, 97 (2), 142–150. (40) Nichols, J. W.; Schultz, I. R.; Fitzsimmons, P. N. In Vitro-in Vivo Extrapolation of Quantitative Hepatic Biotransformation Data for Fish. I. A Review of Methods, and Strategies for Incorporating Intrinsic Clearance Estimates into Chemical Kinetic Models. Aquat. Toxicol. 2006, 78 (1), 74–90. (41) Stieger, G.; Scheringer, M.; Ng, C. A.; Hungerbühler, K. Assessing the Persistence, Bioaccumulation Potential and Toxicity of Brominated Flame Retardants: Data Availability and Quality for 36 Alternative Brominated Flame Retardants. Chemosphere 2014, 116, 118–123. (42) Berntssen, M.; Valdersnes, S.; Rosenlund, G.; Torstensen, B. E.; Zeilmaker, M.; van Eijkeren, J. C. H. Toxicokinetics and Carry-over Model of α Hexabromocyclododecane (HBCD) from Feed to Consumption-Sized Atlantic Salmon (Salmo Salar). Food Addit. Contam. Part A 2011, 28 (9), 1274–1286.
ACS Paragon Plus Environment
31
Environmental Science & Technology
761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784
Page 32 of 36
(43) Airaksinen, R.; Hallikainen, A.; Rantakokko, P.; Ruokojärvi, P.; Vuorinen, P. J.; Mannio, J.; Kiviranta, H. Levels and Congener Profiles of PBDEs in Edible Baltic, Freshwater, and Farmed Fish in Finland. Environ. Sci. Technol. 2015, 49, 3851– 3859. (44) Suominen, K.; Hallikainen, A.; Ruokojärvi, P.; Airaksinen, R.; Koponen, J.; Rannikko, R.; Kiviranta, H. Occurrence of PCDD/F, PCB, PBDE, PFAS, and Organotin Compounds in Fish Meal, Fish Oil and Fish Feed. Chemosphere 2011, 85 (3), 300–306. (45) Voorspoels, S.; Covaci, A.; Schepens, P. Polybrominated Diphenyl Ethers in Marine Species from the Belgian North Sea and the Western Scheldt Estuary: Levels, Profiles, and Distribution. Environ. Sci. Technol. 2003, 37 (19), 4348– 4357. (46) Wan, Y.; Zhang, K.; Dong, Z.; Hu, J. Distribution Is a Major Factor Affecting Bioaccumulation of Decabrominated Diphenyl Ether: Chinese Sturgeon (Acipenser Sinensis) as an Example. Environ. Sci. Technol. 2013, 47 (5), 2279– 2286. (47) Alanärä, A.; Kadri, S.; Paspatis, M. Feeding Management. Food Intake Fish 2001, 332–353. (48) Lock, E. J.; Fjelldal, P. G.; Torstensen, B. E.; Bjoslashrnevik, M.; Breck, O.; Johansen, J.; Reynolds, P.; Sigholt, T.; Joerum, N.; Jakobsen, J. V.; et al. Dietary Decontaminated Fish Oil Has No Negative Impact on Fish Performance, Flesh Quality or Production-Related Diseases in Atlantic Salmon (Salmo Salar). Aquac. Nutr. 2011, 17 (3).
ACS Paragon Plus Environment
32
Page 33 of 36 4
Concentration [pg/g fish]
10
Environmental Science & Technology
103
102
101
100 0
BDE- 209 BDE-99 BDE-47 100
200Plus Environment 300 ACS Paragon Time [day]
400
500
70
Environmental Science & Technology
A: BDE-47
60
Page 34 of 36
ng/day
50 40 30 20 10 0
0
10
100
150
200
250
300
350
400
450
500
200
250
300
350
400
450
500
200
250
300
350
400
450
500
B: BDE-99
8
ng/day
50
dietary uptake waterborne uptake loss via egestion loss via gills biotrasformation (loss) bioformation (input)
6 4 2 0
0
4
50
100
150
100
150
C: BDE-209
ng/day
3
2
1
0
0
50
ACS Paragon Plus Environment
Day
1.5
2.0
1.0 0.5 0.4
1.0 0.5 0.4
0.3
0.3
0.2
0.2
0.0
*
*
* * *
0.1
Model no metabolism
* *
0.0
ACS Paragon Plus Environment
Congener
Model with metabolism
*
28
0.1
Berntssen et al. (2010)
3.0
* * *
99 10 0 10 1 14 9 15 3 15 4 18 3 18 7 20 3 20 9 To ta l
Model with metabolism Model no metabolism
66
Olli et al. (2010)
49
2.0
A
47
2.5
Environmental Science & Technology 4.0 B
28 49 47 (+ 71 ) 66 99 10 0 10 9 11 9 14 9 15 3 15 4 18 3 18 7 20 3 20 9 To ta l
Concentration [ng/g fillet]
Page 35 of 36 3.0
Congener
PBDEs: globally-sourced feed Environmental Science & Technology Page 36 of 36 or local environment? Br
Br
Br Br
Br Br
Br
O Br
Br
deca-BDE
Br
ACS Paragon Plus Environment
Role of debromination?