Subscriber access provided by ECU Libraries
Ecotoxicology and Human Environmental Health
Toxicokinetics of perfluorinated alkyl acids influences their toxic potency in the zebrafish embryo (Danio rerio) Carolina Vogs, Gunnar Johanson, Markus Näslund, Sascha Wulff, Marcus Sjödin, Magnus Hellstrandh, Johan Lindberg, and Emma Elisabeth Wincent Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b07188 • Publication Date (Web): 07 Mar 2019 Downloaded from http://pubs.acs.org on March 8, 2019
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 53
Environmental Science & Technology
1
Toxicokinetics of perfluorinated alkyl acids
2
influences their toxic potency in the zebrafish
3
embryo (Danio rerio)
4
Carolina Vogs1*, Gunnar Johanson1, Markus Näslund1,2, Sascha Wulff1,2, Marcus
5
Sjödin2, Magnus Hellstrandh2, Johan Lindberg2, Emma Wincent1,2
6
1
7 8
Institute of Environmental Medicine, Karolinska Institutet, SE-171 77 Stockholm, Sweden
2
Swedish Toxicology Sciences Research Center (Swetox), 151 36 Södertälje, Sweden
9
10
11
12
ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 2 of 53
13
14
Corresponding Author*
15
Address: Institute of Environmental Medicine, Karolinska Institutet, Sweden
16
Email:
[email protected] 17 18
19
20
Abstract image
21
ACS Paragon Plus Environment
2
Page 3 of 53
Environmental Science & Technology
22
Abstract
23
Perfluorinated alkyl acids (PFAA) are highly persistent and bioaccumulative and have
24
been associated with several adverse health effects. The chemical structure mainly differs
25
in two ways: the length of the hydrophobic alkyl chain and the type of hydrophilic end
26
group. Little is known how the chemical structure affects the toxicokinetics (TK) in different
27
organisms. We studied the TK of four PFAA (PFOS, PFHxS, PFOA, PFBA) with different
28
chain lengths (4-8 carbons) and functional groups (sulfonic or carboxylic acid) in zebrafish
29
(Danio rerio) embryo. The time courses of the external (ambient water) and internal
30
concentrations were determined at three exposure concentrations from 2 up to 120 hours
31
post fertilization (hpf). Three of the four PFAA showed a biphasic uptake pattern with slow
32
uptake before hatching (around 48 hpf) and faster uptake thereafter. A two-compartment
33
TK model adequately described the biphasic uptake pattern, suggesting that the chorion
34
functions as an uptake barrier until 48 hpf. The bioconcentration factors (BCF) determined
35
at 120 hpf varied widely between PFAA with averages of approximately 4000 (PFOS),
36
200 (PFHxS), 50 (PFOA) and 0.8 (PFBA) L kg dry weight-1, suggesting that both the alkyl
ACS Paragon Plus Environment
3
Environmental Science & Technology
Page 4 of 53
37
chain length and the functional group influence the TK. The differences in toxic potency
38
were reduced by three orders of magnitude when comparing internal effect
39
concentrations instead of effective external concentrations.
40
Keywords
41
perfluorooctanesulfonic
42
perfluorobutanoic acid, toxicokinetic model, UPLC-MS/MS
acid,
perfluorohexane
sulfonate,
perfluorooctanoic
acid,
43
ACS Paragon Plus Environment
4
Page 5 of 53
44
Environmental Science & Technology
Introduction
45
Perfluorinated alkyl acid carboxylates and sulfonates (PFAA) are a group of persistent
46
organofluorine chemicals with high thermal and chemical stability, which have been
47
broadly used in commercial and industrial products (e.g. surfactants, fluorinated
48
polymers, coatings, fire-resistant foams). PFAA are extremely persistent in the
49
environment and globally distributed. PFAA contamination have been found in drinking
50
water and river water1-3, soil and biota4-5 and food products6. Perfluorooctane sulfonate
51
(PFOS), perfluorooctanoic acid (PFOA) and perfluorohexane sulfonic acid (PFHxS) are
52
among the most frequently monitored PFAA with varying environmental concentrations
53
of levels from pg L-1 to µg L-1 1, 7-10. The PFAA differ in terms of chain lengths and functional
54
groups and are thus likely to vary with respect to the fate in the environment, toxicokinetics
55
(TK), bioaccumulation and toxicity. However, it is challenging to predict these differences
56
as many of the physicochemical properties (e.g. partitioning coefficients, acid dissociation
57
constants, water solubility) are scarce or controversial and little is known about how these
58
properties affect e.g. protein binding and biotransformation.
ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 53
59
Exposure to PFAA in the general human population occurs mainly through drinking
60
water, diet and consumer products8, 11. Human whole blood concentrations measured in
61
Swedish residents ranged between 1.7 - 37 µg L-1 for PFOS, 0.5 - 12.4 µg L-1 for PFOA,
62
and 0.4 - 28.4 µg L-1 for PFHxS12. Much higher serum levels were detected in another
63
Swedish population (PFOS 24 - 1500 µg L-1, PFOA 2.4 - 92 µg L-1, PFHxS 12 – 1660 µg
64
L-1) exposed via contaminated drinking water8. Several PFAA, including PFOS, PFHxS
65
and PFOA, were detected in serum samples in over 75% of the individuals in a Swedish
66
cohort13. Relative levels of PFAA found in different autopsy tissues (i.e. liver, kidney, lung,
67
brain, bone) indicated that the distribution in the human body varies between the PFAA,
68
which is likely due to their differences in physicochemical properties14. In addition to the
69
wide prevalence of PFAA in humans, elimination from the human body is very slow as
70
reflected by mean half times of 3.4, 2.7 and 5.3 years for PFOS, PFOA and PFHxS,
71
respectively8. Still, the TK processes, and how they differ between PFAA with different
72
chain lengths and functional groups, are poorly understood.
73 74
PFAA accumulation potential in humans may be associated to human health effects, but sparse evidence still exists15,
16.
To date, experimental studies have revealed that
ACS Paragon Plus Environment
6
Page 7 of 53
Environmental Science & Technology
75
PFAA exposure may cause multiple adverse effects such as hepatotoxicity,
76
immunotoxicity, developmental toxicity, reproductive toxicity, neurotoxicity and cancer17-
77
20.
78
PFAA have been shown to interfere with various signalling pathways suggesting that
79
several adverse outcome pathways and cross-talks may play a role in toxicity21. Due to
80
emerging health concerns, large efforts to manage potential risks associated to PFAA
81
exposure are ongoing, such as the banning of longer chain-length PFAA (≥ C8) under the
82
Stockholm Convention on Persistent Organic Pollutants22. As a result, longer chain-length
83
PFAA have been substituted by shorter chain-length PFAA such as PFHxS and
84
perfluorobutanesulfonic acid (PFBA). The shorter chain-length PFAA are assumed to be
85
less bioaccumulative, eliminate faster from the human body and less harmful. However,
86
experimental data confirming these assumptions are yet lacking.
Although the mechanisms underlying these adverse effects are still largely unknown,
87
The zebrafish (Danio rerio) embryo (ZFE) is increasingly used as an alternative model
88
to fish and other higher vertebrates for toxicity testing in risk assessment. The ZFE as
89
experimental model has many advantages, e.g. a rapid ex uterus development with major
90
organs functioning already at 72 h post fertilization (hpf), optical transparency in early
ACS Paragon Plus Environment
7
Environmental Science & Technology
Page 8 of 53
91
developmental stages and cost-efficient usage. Thus, the ZFE model has been applied
92
to study toxicity mechanisms of numerous chemicals from molecular initiating events to
93
phenotypical changes. Previous studies have shown that PFAA with longer chain length
94
and with a sulfonic functional group are more potent with respect to developmental toxicity
95
in ZFE than those with shorter chain length and with a carboxylic functional group23-25.
96
However, the differences in toxic potency could not be explained by mode of action,
97
suggesting that they are due to TK differences. To date, the uptake kinetics of PFAA in
98
ZFE have only been determined for PFOS26 and the TK of other PFAA as well as the
99
relation with e.g. chain length and functional group remain open issues.
100
The aim of this study was to describe and compare the TK profiles and bioconcentration
101
factors of four PFAA in ZFE, namely the sulfonic acids PFOS (C8) and PFHxS (C6) and
102
the carboxylic acids PFOA (C8) and PFBA (C4). Our hypothesis was that the different
103
physicochemical properties of these PFAA have an impact on TK processes such as
104
uptake, biotransformation and excretion. To this end, we measured the PFAA
105
concentrations in exposure medium and in the ZFE at three exposure levels over time. In
106
addition, a two-compartment TK model was applied to fit determined internal
ACS Paragon Plus Environment
8
Page 9 of 53
Environmental Science & Technology
107
concentrations and estimate kinetic rates. Bioconcentration factors were calculated and
108
used to estimate the internal effective concentration of the PFAA.
109 110
Materials and Methods
111 112
Chemicals
113
All chemicals were purchased from Sigma-Aldrich (Stockholm, Sweden) unless stated
114
otherwise. CAS No. and physicochemical properties for the PFAA are listed in the
115
Supplement (SI Table S1). Chemical specific purity was ≥ 98% for PFOS, PFHxS and
116
PFBA and > 96% for PFOA. Stock solutions of PFOS, PFHxS and PFOA were prepared
117
in dimethyl sulfoxide (DMSO, CAS No. 67-68-5) and PFBA was prepared in double-
118
distilled water (ddH2O, laboratory reagent ELGA Purelab flex). 3-morpholinopropane-1-
119
sulfonic acid (MOPS, CAS No. 1132-61-2) and sodium hydroxide (NaOH, CAS No. 1310-
120
73-2) was added to PFBA and PFOA exposure media to prevent pH changes. HPLC
121
grade acetonitrile (ACN, CAS No. 75-05-8, VWR, Stockholm, Sweden) and ddH2O were
122
used for sample preparation. The isotope labelled internal standards (IS) sodium
ACS Paragon Plus Environment
9
Environmental Science & Technology
Page 10 of 53
123
perfluoro-1-[1,2,3,4-13C4]octanesulfonate, perfluoro-n-[1,2,3,4-13C4]octanoic acid, sodium
124
perfluoro-1-hexane[18O2]sulfonate and perfluoro-n-[1,2,3,4-13C4]butanoic acid were
125
purchased from Wellington Laboratories (Guelph, Canada). E3-Medium were prepared
126
using sodium chloride (NaCl, CAS No 764714-5), potassium chloride (KCl, CAS No.
127
7447-40-7), calcium chloride dihydrate (CaCl2*2H20, CAS No. 10035-04-8), magnesium
128
sulfate (MgSO4, Cas No. 7487-88-9) and 4-(2-Hydroxyethyl)piperazine-1-ethanesulfonic
129
acid (HEPES, CAS No. 7365-45-9). For the UPLC-MS/MS analysis, mobile phases
130
consisted of methanol (MeOH, CAS: 67-56-1, Optima LC/MS GRADE, Fisher Scientific),
131
acetonitrile (ACN, CAS No. 75-05-8, Optima LC/MS GRADE, Fisher Scientific) and
132
ammonium acetate (NH4Ac, CAS No. 631-61-8, LC/MS LiChropur).
133 134
Zebrafish and exposure
135
ZFE of the AB strain and E3 medium (5.0 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, 0.33
136
mM MgSO4, 0.2 mM HEPES) were provided by the Zebrafish Core Facility at Comparative
137
Medicine, Karolinska Institutet. Three male and female adult zebrafish were paired and
138
the eggs were collected directly after spawning and transferred into E3 medium. The ZFE
ACS Paragon Plus Environment
10
Page 11 of 53
Environmental Science & Technology
139
batch was examined under a light microscope (Nikon SMZ1270 Stereomicroscope) to
140
discard unfertilized eggs. All exposure experiments were set up according to the test
141
guideline OECD TG 23627 for fish embryo toxicity testing. Briefly, embryos were exposed
142
from 2 hpf up to 120 hpf. 30 embryos were added to 30 mL exposure medium (one
143
embryo per mL) in 50 mL glass petri dishes. The dishes, covered with lids, were kept in
144
a climate chamber at 28 ± 1 ºC under dark conditions throughout the exposures.
145 146
Concentration-effect relationships of PFAA
147
Two or three independent concentration-effect-relationship experiments were
148
conducted with each PFAA to determine the concentration causing effects in 20% of the
149
embryos (EC20). The EC20 values were subsequently used as the upper limit for the
150
highest exposure concentrations in the TK experiments. Effects on morphology were
151
assessed at 120 hpf using a Nikon SMZ25 microscope and included lethal (coagulation,
152
lack of heartbeat, non-detachment of tail, non-hatching and lack of somite formation) as
153
well as sublethal (non-inflated swim bladder, pericardial and yolk sac edemas, scoliosis)
154
apical endpoints. A ZFE with any of these morphological changes was considered to be
ACS Paragon Plus Environment
11
Environmental Science & Technology
Page 12 of 53
155
affected. Starting at 2 hpf, fertilized ZFE were exposed to a wide range of concentrations
156
of respective PFAA (Table 1). To avoid effects from altered pH, the exposure media with
157
PFOA and PFBA were buffered with 10 mM MOPS and adjusted to pH of 7.4 with 5 M
158
NaOH23. Concentration-effect relationships were obtained by fitting the sigmoidal Hill
159
model (eq. 1) to respective PFAA data set using the GraphPad software (version 5.02,
160
La Jolla, USA) 𝐶𝑛
161
𝐸 = 𝐸𝑚𝑖𝑛 +(𝐸𝑚𝑎𝑥 ― 𝐸𝑚𝑖𝑛) × (𝐸𝐶𝑛
162
E represents the percentage of ZFE affected at concentration C, Emin and Emax the
163
minimum and maximum percentage affected, n the slope of the curve and EC50 the
164
concentration causing an effect in 50% of the ZFE. The parameters resulting in the best
165
fit to the experimental data were subsequently used to calculate the exposure
166
concentration causing an effect in 20% of the ZFE (EC20) (Table 1).
50
+ 𝐶𝑛)
Eq. 1
167
Toxicokinetic experiment
168
TK experiments were performed in order to determine “internal concentration” and
169
compare “internal effective concentration” of PFAA which are apparent values based on
ACS Paragon Plus Environment
12
Page 13 of 53
Environmental Science & Technology
170
the total chemical mass measured in the whole embryo, and includes PFAA bound e.g.
171
to the surfaces of the skin and chorion in addition to PFAA in ZFE tissue. Each TK
172
experiment consisted of three different exposure scenarios performed in parallel as
173
outlined in Figure 1. First, 30 ZFE were incubated in 30 mL E3 medium without addition
174
of PFAA as a control for vitality of the ZFE batch, referred to as “biological control”.
175
Second, triplicates of 30 mL exposure medium for three different concentrations of
176
respective PFAA were prepared in E3 medium and incubated in glass petri dishes without
177
adding ZFE, referred to as “ZFE-free medium”. The ZFE-free medium was used to check
178
for unspecific loss due to e.g. adsorption to glass surfaces. Third, fertilized ZFE were
179
exposed to three concentrations of each PFAA in E3 medium, hereafter called “ZFE
180
medium”. All ZFE exposures were performed with 30 embryos and 30 mL medium. The
181
high concentration (C1) was selected based on the EC20 values and the mid (C2) and low
182
(C3) concentrations were 10 and 50 times dilutions of C1, respectively (Table 2). ZFE
183
and medium were sampled at 0 (ZFE-free medium only), 3, 6, 9, 24, 31, 48, 72, 96 and
184
120 hpf. Regarding medium samples, 300 µL aliquots were taken in three replicates from
185
each petri dish and transferred to amber glass vials (2 mL, Agilent, Stockholm, Sweden).
ACS Paragon Plus Environment
13
Environmental Science & Technology
Page 14 of 53
186
Regarding ZFE, five replicates each consisting of five embryos were collected at each of
187
the above time points. The ZFE were transferred to microtubes (1.5 mL, Eppendorf Safe-
188
Lock tubes, Sigma-Aldrich, Stockholm, Sweden) and excess medium was removed by
189
pipetting. The ZFE were immediately washed twice with 300 µL ddH2O for a short time
190
period of maximum 5 min in total to minimize back diffusion of PFAA. Subsequently, all
191
samples were stored at -80ºC until sample preparation and chemical analysis. Two
192
independent TK studies were conducted using ZFE from different fish batches.
193 194
Sample preparation
195
All samples were thawed at room temperature before sample preparation. Chemical-
196
specific internal standards (IS) were added to all samples prior to extraction to correct for
197
possible matrix effects and/or variations in extraction efficiencies (Figure 1). Defrosted
198
medium samples (300 µL) were treated by adding 700 µL of 100 nM IS dissolved in ACN
199
(100%). Prior to ZFE sample preparation, 300 µL of IS solution (100 nM IS in ACN
200
dissolved in ddH20 (v/v 70:30%)) was added. Thereafter, the ZFE were homogenized at
201
room temperature by sonication (up to 4 times for 5 seconds to prevent evaporation, MSE
ACS Paragon Plus Environment
14
Page 15 of 53
Environmental Science & Technology
202
Soniprep 150, Measuring and Scientific Equipment, London, UK), incubated at room
203
temperature for 1 h and centrifuged for 10 min at 14000 rpm. The supernatants were
204
transferred to amber glass vials and stored at -80 ºC prior to chemical analysis. Extraction
205
recovery and matrix effect analysis are described in the Supplement.
206 207
Chemical analysis
208
All chemical analyses were performed using UPLC-MS/MS. Injection volumes were 10
209
µL, 1 µL, 2 µL and 1 µL for PFOS, PFHxS, PFOA and PFBA, respectively. The PFAA
210
were separated on an Acquity™ (Waters, Milford, MA, USA) on a 50 x 2.1 mm BEH C18
211
column (Waters, Sollentuna, Sweden) with gradient elution (solvent A: 2% MeOH, 1 mM
212
NH4Ac in H2O; solvent B: 1 mM NH4Ac in ACN) using two separate gradient methods (SI
213
Table S2). The mass spectrometer was a Xevo™ TQD (Waters), which operated in
214
multiple reaction monitoring modes with negative ion mode electrospray ionization [M-H].
215
A detailed description of the compound specific settings can be found in Supplement
216
(Table S2). Quantification was achieved by external PFAA calibrators (n ≥ 5) prepared in
217
ACN:ddH20 (v/v 70:30%) with ranges of 6-800 nM for PFOS, 8-25000 nM for PFHxS, 11–
ACS Paragon Plus Environment
15
Environmental Science & Technology
Page 16 of 53
218
2800 nM for PFOA and 16–2050 nM for PFBA. Limits of quantification, defined as ten
219
times the signal-to-noise ratios, were 6 nM, 25 nM, 6 nM and 18 nM for PFOS, PFHxS,
220
PFOA and PFBA, respectively. The method error of concentrations expressed as the
221
coefficient of variance (CoV) were 13% for PFOS, 14% for PFHxS, 13% for PFOA and
222
32% for PFBA in average (SI Table S3). The large CoV for PFBA resulted from large
223
variation of measured ZFE concentrations per time point and concentration.
224 225
Data analysis
226
The results from analyses of medium were combined from two independent
227
experiments, each consisting of three replicates per concentration and time point. The
228
PFAA amount per embryo was determined by dividing the total amount of respective
229
PFAA in ZFE extract with the number of embryos used per extract. The internal amount
230
in ZFE was thereafter based on volume in µmol per liter (µM) ZFE, dry weight in µmol per
231
kg dry ZFE and total protein content in pmol per µg ZFE protein at respective ages (SI
232
Figure S1) to calculate internal concentrations. ZFE volume and dry weight at different
233
ages were taken from Brox et al.28, and total protein content was determined using the
ACS Paragon Plus Environment
16
Page 17 of 53
Environmental Science & Technology
234
Coomassie (Bradford) Protein Assay Kit according to the manufacturer’s protocol
235
(Thermos Scientific, Rockford, USA). Data from ZFE samples are presented as means of
236
ten replicates for each concentration and sampling time point, originating from two
237
independent experiments. Outliers were excluded using the Grubb’s test.
238 239
Toxicokinetic characterization
240
Net uptake rates were calculated as the derivative of internal amount [pmol h-1]. To
241
quantify the total accumulated amount of PFAA over the exposure time. The area under
242
the curve (AUC [µM h]) was determined using the trapezoidal method and a Monte Carlo
243
simulation was applied to calculate standard errors (SE) (n = 500) (Berkley Madonna,
244
Version 8.3.18, California, USA). Bioconcentration factors (BCF) were determined as 𝐶𝑖𝑛𝑡
245
𝐵𝐶𝐹 =
𝐶𝑤 ,
Eq. 2
246
where Cint and Cw are the average concentrations in ZFE and medium, respectively, at
247
120 hpf. BCF values were calculated based on both embryo volume28 (BCFvol, unitless)
ACS Paragon Plus Environment
17
Environmental Science & Technology
Page 18 of 53
248
and dry weight28 (BCFdry, L kg-1). Internal effect concentrations (IEC50 [µM]) were obtained
249
by multiplying the average BCFvol with the EC50 value:
250
𝐼𝐸𝐶50 = 𝐵𝐶𝐹𝑣𝑜𝑙 × 𝐸𝐶50
Eq. 3
251
Kinetic rates were estimated by fitting one- and two-compartment models to the ZFE
252
concentrations using Berkley Madonna. The one-compartment model included uptake
253
(kWE [h-1]) and elimination (kEW [h-1]) kinetic rates for exchange processes between ZFE
254
medium (Cw) and ZFE (Cint) neglecting the chorion as transport barrier.
255
𝑑𝐶𝑖𝑛𝑡
256
A two-compartment model was developed in order to implement uptake (kWC [h-1]) and
257
elimination (kCW [h-1]) kinetic rates for exchange processes between ZFE medium and
258
chorion (Cc) and uptake (kCE [h-1]) and elimination (kEC [h-1]) kinetic rates for exchange
259
processes between ZFE chorion and ZFE until hatching at 48 hpf. Therefore, a step
260
function was implemented to mimic ZFE with a chorion before hatching and ZFE without
261
chorion fully exposed to PFAA after hatching.
262
𝑑𝑡
𝑑𝐶𝑐 𝑑𝑡
Eq. 4
= 𝑘𝑊𝐸 × 𝐶𝑊 ― 𝑘𝐸𝑊 × 𝐶𝑖𝑛𝑡
=
𝑘𝑊𝐶 × 𝐶𝑤 ― 𝑘𝐶𝑊 × 𝐶𝑐 + 𝑘𝐸𝐶 × 𝐶𝑖𝑛𝑡 ― 𝑘𝐶𝐸 × 𝐶𝑐 𝑡 ≤ 48ℎ 0 𝑡 > 48ℎ
Eq. 5
ACS Paragon Plus Environment
18
Page 19 of 53
Environmental Science & Technology
263 𝑘𝐶𝐸 × 𝐶𝑐 ― 𝑘𝐸𝐶 × 𝐶𝑖𝑛𝑡 =𝑘 ×𝐶 ―𝑘 ×𝐶
𝑡 ≤ 48ℎ 𝑡 > 48ℎ
264
𝑑𝐶𝑖𝑛𝑡
265
To compare the fits of the one-compartment model and the two-compartment model to
266
the observed internal concentrations, we calculated the Akaike Information Criteria (AIC
267
[-])29.
𝑑𝑡
𝑊𝐸
𝑤
𝐸𝑊
𝑖𝑛𝑡
∑(𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑 ― 𝑝𝑟𝑒𝑑𝑖𝑐𝑡𝑒𝑑)2
Eq. 6
268
𝐴𝐼𝐶 = 𝑁 × 𝐿𝑛(
269
with N representing sample size and Npar representing number of parameters.
𝑁
+2 × 𝑁𝑝𝑎𝑟)
Eq. 7
270 271
Results and Discussion
272 273
PFAA concentration-effect relationships
274
We first determined concentration-effect relationships for the four PFAA as too high
275
concentrations may influence the TK processes due to embryonic toxicity. ZFE exposed
276
to a wide range of concentrations (Table 1) were followed up to 120 hpf with regards to
277
aberrant morphology. Aberrant effects were observed mainly at 96 and 120 hpf with non-
ACS Paragon Plus Environment
19
Environmental Science & Technology
Page 20 of 53
278
inflated swim bladder and scoliosis as the most prominent sublethal endpoints (SI Figure
279
S2). The toxic potency of the PFAA, determined at 120 hpf, differed by 3000-fold with
280
PFOS being most toxic followed by PFHxS and PFOA (Figure 2). PFBA showed no
281
significant toxicity even at the highest concentration (70 mM). Moreover, all embryos
282
exposed to the highest concentrations of PFOS or PFHxS were affected, whereas 100%
283
affected ZFE was not reached at any tested concentration of PFOA or PFBA. EC50 values
284
were estimated to be 3.8, 84.5 and 509 µM for PFOS, PFHxS and PFOA, respectively.
285
For PFBA, the EC50 value could not be determined due to lack of toxicity. These data
286
indicate that both the length of the fluorinated carbon chain (8, 6, or 4 carbons) and the
287
functional group (sulfonic or carboxylic acid) of PFAA influence the toxicity. While there
288
are no reports on the toxicity of PFHxS in ZFE, the EC50 values for PFOS and PFOA
289
obtained in our study are similar to those previously reported (Table 1). Furthermore, our
290
data confirm the previous notation that PFAA with a sulfonic acid group are more toxic to
291
ZFE than those with a carboxylic acid group23-25. These authors could however not explain
292
the mechanisms underlying the previously observed differences in PFAA toxicity, but
293
suggested this might be due to TK variability. Hagenaars et al.30 reported that effects
ACS Paragon Plus Environment
20
Page 21 of 53
Environmental Science & Technology
294
caused by PFOS exposure depend on the duration of exposure suggesting that a certain
295
accumulated concentration over time in ZFE is essential to cause the observed adverse
296
effects.
297 298
Time-course of PFAA concentration in exposure medium
299
To minimize confounding effects of toxicity on TK processes, the highest exposure
300
concentrations were set equal to or below the EC20 values (Table 2). Accordingly, the
301
highest concentrations differed by four orders of magnitude between the four PFAA. In
302
total, three exposure concentrations were selected for respective PFAA, encompassing
303
low to high levels compared to human exposure levels6-13.
304
To check for unspecific losses of PFAA due to e.g. adsorption to glass surface, we first
305
determined the concentration of respective PFAA over time in ZFE-free medium. All
306
concentrations of PFHxS, PFOA and PFBA remained constant over the 120 h exposure
307
while PFOS concentrations decreased within the first 3 h of exposure (SI Figure S3).
308
Measurements of ZFE-free medium showed highest variation at the lowest
ACS Paragon Plus Environment
21
Environmental Science & Technology
Page 22 of 53
309
concentrations of PFOS and PFHxS and at the highest concentrations of PFOA and
310
PFBA (SI Table S3).
311
We next determined the time course of PFAA concentration in the ZFE medium (Figure
312
3A). The PFFA concentrations in the medium remained stable throughout the 120 h
313
exposure, with minor exceptions. ZFE medium concentration time courses did not differ
314
from those observed with the ZFE-free medium. (SI Figure S3). The lack of a clear
315
decrease indicates that the amount of PFAA absorbed by the ZFE is negligible relative to
316
the total amount of PFAA in the exposure medium. In contrast to our finding, Kühnert et
317
al.31 described chemical depletion of polycyclic aromatic hydrocarbons (PAH) in ZFE
318
medium over time. While the different exposure regimen used in the PAH study may
319
underlie some of the observed differences, the most likely explanation is the differences
320
in physicochemical properties, e.g. a much higher lipophilicity (log Kow 3.3 – 5.8)
321
compared to the PFAA in our study (measured log Kow 0.5 – 2.4, SI Table S1).
322 323
Time-course of PFAA concentration in ZFE
ACS Paragon Plus Environment
22
Page 23 of 53
Environmental Science & Technology
324
Absorption in the ZFE was detected for all PFAA and concentrations (Figure 3B and
325
3C). The time courses of accumulated amount did not differ in shape when normalized to
326
volume, dry weight or total protein content at different ZFE ages, respectively (SI Figure
327
S4). PFOS, PFHxS and PFOA showed biphasic uptake patterns with an initial slow
328
increase of internal concentrations up to 31 hpf, followed by faster increase. In contrast,
329
PFBA reached nearly maximum internal concentrations already within approximately 24
330
hpf and showed a very slow increase afterwards. The average net uptake rates (pmol h-
331
1)
332
PFOA, which reflects the biphasic uptake pattern (Figure 4). The uptake rate of PFBA
333
showed high fluctuations at the different time points. To our knowledge, a biphasic uptake
334
pattern of chemicals in the ZFE has not previously been reported, the one exception being
335
PFOS26. Furthermore, there are no published reports comparing the uptake of PFAA with
336
different chain lengths and functional groups in ZFE. In the few TK studies on ZFE
337
published to date, time-profiles of internal concentrations have been described as either
338
a very fast uptake (similar to PFBA in our study) for polar compounds or a monophasic
339
uptake until steady-state for nonpolar chemicals28,
were 0.4-4.8 times higher from 31 hpf for all three concentrations of PFOS, PFHxS and
31, 32.
In these studies, the uptake
ACS Paragon Plus Environment
23
Environmental Science & Technology
Page 24 of 53
340
kinetics were proposed to depend mainly on diffusion processes and lipophilicity
341
(disregarding biotransformation). The continuous increase in internal concentration over
342
time of PFOS, PFHxS and PFOA suggests that biotransformation is lacking or of minor
343
importance. This is in line with the literature showing no evidence that PFAA are
344
biotransformed in fish or other species.
345
The initial slow uptake rate of PFOS, PFHxS and PFOA suggests that the ZFE chorion
346
functions as a transport barrier (Figure 4). The chorion with the perivitelline space, which
347
surrounds the embryo, serves as barrier for macromolecules but is permeable enough to
348
allow supply of oxygen and nutritients33. PFAA are negatively charged organic chemicals
349
at pH 7 characterized by very low pKA values (SI Table S1), which can potentially lead to
350
repulsion from the chorion. Furthermore, PFAA have a high affinity to proteins and
351
transporters that could enable chorion binding to e.g. N-linked glycoproteins34. Brox et
352
al.35 found 5 to 30 ng higher concentrations of polar compounds in ZFE with chorion
353
compared to dechorinated ZFE, concluding that the internal concentration in the ZFE with
354
chorion was overestimated until the ZFE were hatched.
ACS Paragon Plus Environment
24
Page 25 of 53
Environmental Science & Technology
355
Lowering exposure concentrations of PFAA resulted in lower internal concentrations.
356
As demonstrated by the area under the concentration-time curve (AUC, Table 2), an
357
approximately ten times concentration dilution from high to mid exposure concentration
358
led to 6, 9.5, 2.4 and 5.8 times decrease in AUC for PFOS, PFHxS, PFOA and PFBA,
359
respectively. In contrary, an approximately two times dilution from mid to low exposure
360
concentration resulted in a proportional decrease of AUC and similar net uptake rates
361
(Figure 4). This finding suggest concentration-dependent TK processes of PFAA in ZFE.
362
Our findings agree with Chen et al.36 who reported that the BCF values of long-chain
363
PFAA (PFDoA, PFTrDA, PFTeDA) in adult zebrafish were significantly higher at low than
364
at high concentrations, likely as a result of concentration-dependent TK processes.
365
Various studies36-38 emphasized that concentration-dependent TK process is relevant for
366
PFAA due to saturable protein-mediated uptake mechanisms in adult fish in addition to
367
passive diffusion. PFAA specifically bind to numerous proteins (e.g. apolipoprotein or
368
albumin in plasma, fatty acid binding proteins in liver and organic anion transporters in
369
kidney) in fish and other species.
ACS Paragon Plus Environment
25
Environmental Science & Technology
Page 26 of 53
370
PFAA internal concentrations increased over time without reaching apparent steady-
371
state for any PFAA or concentration as indicated by positive net uptake rates at 120 hpf
372
(Figure 4). However, the net uptake rate decreased between 96 hpf and 120 hpf, the
373
decrease being more pronounced at high exposure concentrations. Maximum internal
374
concentrations were similar between PFAA with the sulfonic end group (PFOS and
375
PFHxS) and between PFAA with the carboxylic end group (PFOA and PFBA) even when
376
the exposure concentrations differed by ten times. This result corresponded with
377
approximately ten times lower net uptake rates of PFOS and PFHxS compared with those
378
of PFOA and PFBA suggesting different bioconcentration potentials between PFAA with
379
different functional groups.
380 381
Bioconcentration factors (BCF)
382
Bioconcentration factors varied four orders of magnitude between the four PFAA and
383
were higher for PFAA with sulfonic acid end group than those for PFAA with carboxylic
384
acid end group (Table 2). In addition, the high exposure concentration yielded lower BCF
385
values for all four PFAA compared to mid and low exposure concentrations. Thus, the
ACS Paragon Plus Environment
26
Page 27 of 53
Environmental Science & Technology
386
BCFvol were similar at the mid and low concentrations, but 1.5, 1.5, 4.9 and 1.7 times
387
lower for PFOS, PFHxS, PFOA and PFBA, respectively, at the high concentrations.
388
These results again suggest saturable protein binding at the high concentrations.
389
Using the ZFE, we could confirm that BCF values increased with carbon chain length
390
as previously reported for adult fish39-43. In addition, our BCFdry values were in the range of
391
reported BCF values stated for different adult fish species (BCFdry 3.1 – 3200 L kg -1, Table 2).
392
The study by Hung et al.26 on PFOS is of particular interest as they also used ZFE.
393
However, we observed two orders of magnitude higher BCFvol for PFOS compared to
394
Hung et al.26. The discrepancy may be due to different exposure conditions, e.g. five
395
embryos per mL medium by Hung et al.26 versus one embryo per mL medium in our study.
396
Since steady-state concentrations in the ZFE were not reached at 120 hpf for all PFAA,
397
the BCFvol likely underestimates the “true” steady-state BCF.
398
For adult fish, it has been shown that PFAA dominantly accumulate in the liver, blood,
399
kidney and intestine, while approximately ten times lower concentrations are found in the
400
brain, ovary, gall bladder and muscle42-45. At present, it is not possible to determine the
401
tissue distribution of chemicals within the ZFE. Uhlag et al.46 demonstrated that exposure
ACS Paragon Plus Environment
27
Environmental Science & Technology
Page 28 of 53
402
of PFAA with longer carbon chain length and with sulfonic group caused larger changes
403
in locomotor responses suggesting PFAA structure-dependent accumulation in ZFE brain
404
tissue. Overcoming the knowledge gap on PFAA tissue distribution in ZFE would increase
405
the understanding of organ-specific toxicity (i.e. developmental neurotoxicity, chemical-
406
induced hepatotoxicity) and enable comparison of target specific concentrations in other
407
species.
408
Toxicokinetic Modeling
409
We fitted a one- and a two-compartment TK model to the observed internal
410
concentrations in ZFE in order to estimate uptake and elimination kinetic rates (SI Figure
411
S5, Table S4). The one-compartment model did not adequately describe the curvature of
412
the concentration-time data but rather predicted a constant increase in internal
413
concentration without reaching steady-state. Meanwhile, the two-compartment model
414
with inclusion of chorion as a transport barrier addressed the curvature and gave a better
415
fit to the experimental values. Time trends of the fitted data from each independent TK
416
experiment were consistent and differences are mainly caused by inter-variability (SI
417
Figure S6). The Akaike Information Criteria (AIC) gave further support for the two-
ACS Paragon Plus Environment
28
Page 29 of 53
Environmental Science & Technology
418
compartment model (lower AIC values), especially for the two higher PFAA exposure
419
concentrations (SI Table S4). Modelling the time courses for PFBA (all three
420
concentrations) and for lowest concentration (C3, all PFAA) generally gave less god fit
421
and less clear preference for the two-compartment model, possibly due to larger method
422
errors. Another shortcoming may be that non-linear processes such as growth and
423
concentration-dependent uptake were not included in the model.
424
The two-compartment model described a slow uptake of PFOS, PFHxS and PFOA
425
before hatching at 48 hpf and a faster uptake after hatching, in agreement with our
426
observations of internal concentrations (Figure 3C). The two-compartment model was
427
also able to match the fast increase of internal PFBA concentrations before hatching. The
428
uptake rate (kWE) differed by several orders of magnitude between PFAA (PFOS > PFHxS
429
> PFOA > PFBA), while the elimination rate (kEW) only differed one order of magnitude
430
except for PFHxS (SI Table S4). For the mid and low PFHxS exposure concentrations,
431
we estimated much lower elimination rates as a consequence of the continued stable net
432
uptake rate between 96 hpf and 120 hpf (Figure 4). Although these rate constants should
433
be interpreted with caution, it seems that the uptake rate depends more on the carbon
ACS Paragon Plus Environment
29
Environmental Science & Technology
Page 30 of 53
434
chain length and functional group than does the elimination rate, in agreement with the
435
findings by Chen et al.36. We did however not find any general pattern for the PFAA uptake
436
and elimination rates between water and chorion and between chorion and embryo.
437 438
Implication for PFAA toxicity
439
Recalculation of external effect concentrations (EC) to internal effect concentrations
440
causing 20% (IEC20) in ZFE reduced the variability in toxic potency between PFAA from
441
3000-fold to 3-fold (Figure 5, Table 2). The incorporation of BCFvol thus resulted in a shift
442
in toxicity ranking from PFOS > PFHxS > PFOA > PFBA using EC20 values to PFBA >
443
PFOS ~ PFHxS ~ PFOA using IEC20 values. Altogether, this indicates that the functional
444
group as well as the carbon chain length influence TK and therefore the toxicity. In
445
agreement, Gomis et al.47 also showed that TK differences explained most of the PFAA
446
variability between concentration-effect curves of liver enlargement in male rats.
447
Moreover, the intrinsic toxic potential differed only 3-fold between PFAA indicating that
448
the toxic mechanism of action is the same for the PFAA tested. A toxic ratio analysis using
449
our calculated IEC50 values revealed that PFOS, PFHxS and PFOA are 100 times higher
ACS Paragon Plus Environment
30
Page 31 of 53
Environmental Science & Technology
450
in toxicity than baseline toxicants48 (SI Figure S7). This result is in agreement with
451
literature showing PFAA specifically bind to receptors such as peroxisome proliferator-
452
activated receptors49 or estrogen receptors50 suggesting a specific mechanism of action.
453
Altogether, this indicates that toxicity differences in ZFE can dominantly be explained by
454
differences in TK rather than differences in toxic mechanisms of action.
455
Our findings and those by Gomis et al.47 indicate that the intrinsic toxic potency of short-
456
chain and long-chain PFAA may be similar. Long-chain PFAA are currently substituted
457
by short-chain PFAA (e.g. PFHxS, PFBA and Gen-X) based on the view that short-chain
458
PFAA are less toxic. As PFAA are environmentally persistent, short-chain and long-chain
459
PFAA will coexist in the environment. Meanwhile, the short-chain PFAA are more
460
bioavailable and the long-chain PFAA are more bioaccumulative. This leads to the issue
461
of combined exposure to several PFAA. Wen et al.45 reported that the bioconcentration
462
potential of short-chain PFAA was reduced in adult zebrafish when exposed in the
463
presence of long-chain PFAA, likely due to competition for transporters and binding sites
464
of proteins. However, combined PFAA exposure was shown to be more toxic to ZFE
465
compared to single PFAA exposures51. Furthermore, combined PFAA exposure resulted
ACS Paragon Plus Environment
31
Environmental Science & Technology
Page 32 of 53
466
in complex interactive effects, including concentration-dependent additive, synergistic
467
and antagonistic mixture effects52, suggesting that the competition for transporters and
468
binding sites of proteins seen in adult zebrafish also occurs in ZFE. The ZFE model may
469
thus be useful to compare and better understand the influence of mixed PFAA exposure
470
on TK processes and toxicity. Such understanding is essential in the study design of
471
mixed exposure experiments.
472
Acknowledgement
473
The study was supported by a grant from the Swedish Research Council Formas (2014-
474
475
1454). We are thankful to the zebrafish core facility, Karolinska Institutet.
Supporting Information
476
Descriptions of extraction recovery and matrix effects
477
Figure S1 – S7: ZFE growth data; phenotypes of PFAA exposed ZFE; ZFE-free medium
478
and ZFE-medium concentrations; internal concentrations based on different biological
ACS Paragon Plus Environment
32
Page 33 of 53
Environmental Science & Technology
479
parameters; fits of one- and two compartment model; fits of two-compartment model for
480
each individual TK experiment; toxic ratio analysis
481
Table S1 – S4: physicochemical PFAA characterization; PFAA specific settings for
482
UPLC-MS/MS; coefficient of variance; kinetic parameter estimations of one- and two-
483
compartment model
484 485 486
References
487
1.
Gyllenhammar, I.; Berger, U.; Sundström, M.; McCleaf, P.; Eurén, K.; Eriksson, S.;
488
Ahlgren, S.; Lignell, S.; Aune, M.; Kotova, N., Influence of contaminated drinking water
489
on perfluoroalkyl acid levels in human serum – A case study from Uppsala, Sweden. Envir
490
Res. 2015, 140, 673-683.
491 492
2.
Lindim, C.; Van Gils, J.; Cousins, I. T., Europe-wide estuarine export and surface
water concentrations of PFOS and PFOA. Water Res. 2016, 103, 124-132.
ACS Paragon Plus Environment
33
Environmental Science & Technology
493
3.
Page 34 of 53
Post, G. B.; Gleason, J. A.; Cooper, K. R., Key scientific issues in developing
494
drinking water guidelines for perfluoroalkyl acids: Contaminants of emerging concern.
495
PLoS Bio. 2017, 15, (12), e2002855.
496
4.
Naile, J. E.; Khim, J. S.; Wang, T.; Chen, C.; Luo, W.; Kwon, B.-O.; Park, J.; Koh,
497
C.-H.; Jones, P. D.; Lu, Y., Perfluorinated compounds in water, sediment, soil and biota
498
from estuarine and coastal areas of Korea. Environ Pollut. 2010, 158, (5), 1237-1244.
499
5.
Greaves, A. K.; Letcher, R. J.; Sonne, C.; Dietz, R.; Born, E. W., Tissue-specific
500
concentrations and patterns of perfluoroalkyl carboxylates and sulfonates in East
501
Greenland polar bears. Environ Sci Technol 2012, 46, (21), 11575-11583.
502
6.
Noorlander, C. W.; Van Leeuwen, S. P.; Te Biesebeek, J. D.; Mengelers, M. J.;
503
Zeilmaker, M. J., Levels of perfluorinated compounds in food and dietary intake of PFOS
504
and PFOA in the Netherlands. J Agricult Food Chem 2011, 59, (13), 7496-7505.
505 506
7.
Sundström, M.; Ehresman, D. J.; Bignert, A.; Butenhoff, J. L.; Olsen, G. W.; Chang,
S.-C.; Bergman, Å., A temporal trend study (1972–2008) of perfluorooctanesulfonate,
ACS Paragon Plus Environment
34
Page 35 of 53
Environmental Science & Technology
507
perfluorohexanesulfonate, and perfluorooctanoate in pooled human milk samples from
508
Stockholm, Sweden. Environ Int 2011, 37, (1), 178-183.
509
8.
Li, Y.; Fletcher, T.; Mucs, D.; Scott, K.; Lindh, C. H.; Tallving, P.; Jakobsson, K.,
510
Half-lives of PFOS, PFHxS and PFOA after end of exposure to contaminated drinking
511
water. Occup Environ Med 2017, 75, (1), 46-51.
512
9.
Banzhaf, S.; Filipovic, M.; Lewis, J.; Sparrenbom, C. J.; Barthel, R.; A review of
513
contamination of surface-, ground-, and drinking water in Sweden by perfluoroalkyl and
514
polyfluoroalkyl substances (PFASs). Ambio 2017, 46, (3), 335-346.
515
10. Jian, J. M.; Guo, Y.; Zeng, L.; Liang-Ying, L.; Lu, X.; Wang, F.; Zeng, E. Y.; Global
516
distribution of perfluorochemicals (PFCs) in potential human exposure source - A review.
517
Environ Int 2017, 108, 51-62.
518
11. Vestergren, R.; Berger, U.; Glynn, A.; Cousins, I. T., Dietary exposure to
519
perfluoroalkyl acids for the Swedish population in 1999, 2005 and 2010. Environ Int 2012,
520
49, 120-127.
ACS Paragon Plus Environment
35
Environmental Science & Technology
Page 36 of 53
521
12. Kärrman, A.; van Bavel, B.; Järnberg, U.; Hardell, L.; Lindström, G., Perfluorinated
522
chemicals in relation to other persistent organic pollutants in human blood. Chemosphere
523
2006, 64, (9), 1582-1591.
524
13. Stubleski, J.; Salihovic, S.; Lind, L.; Lind, P. M.; van Bavel, B.; Kärrman, A.,
525
Changes in serum levels of perfluoroalkyl substances during a 10-year follow-up period
526
in a large population-based cohort. Environ Int 2016, 95, 86-92.
527
14. Pérez, F.; Nadal, M.; Navarro-Ortega, A.; Fàbrega, F.; Domingo, J. L.; Barceló, D.;
528
Farré, M., Accumulation of perfluoroalkyl substances in human tissues. Environ Int 2013,
529
59, 354-362.
530
15. Steenland, K.; Fletcher, T.; Savitz, D. A., Epidemiologic evidence on the health
531
effects of perfluorooctanoic acid (PFOA). Environ health perspectives 2010, 118, (8),
532
1100.
533
16. Lei, M.; Zhang, L.; Lei, J.; Zong, L.; Li, J.; Wu, Z.; Wang, Z., Overview of emerging
534
contaminants and associated human health effects. Biomed Res Int 2015, 2015, Article
535
ID 404796, doi: 10.1155/2015/404796.
ACS Paragon Plus Environment
36
Page 37 of 53
Environmental Science & Technology
536
17. Song, M-K.; Cho, Y. S.; Jeong, S-C.; Ryu, J.B., Analysis of gene expression
537
changes in relation to hepatotoxicity induced by perfluorinated chemicals in a human
538
hepatoma cell line. J. Toxicol. Environ. Health Sci. 2016, 8, 114-127.
539
18. Das K.P.; Grey B.E.; Rosen M.B.; Wood C.R.; Tatum-Gibbs K.R.; Zehr R.D.;
540
Strynar M.J.; Lindstrom A.B.; Lau C., Developmental toxicity of perfluorononanoic acid in
541
mice. Reprod Toxicol 2015, 51, 133-44.
542
19. Wang Y.; Wang L.; Chang W.; Zhang Y.; Zhang Y; Liu W., Neurotoxic effects of
543
perfluoroalkyl acids: Neurobehavioral deficit and its molecular mechanism. Toxicol Lett.
544
2019, 305, 65-72.
545 546
20. Grandjean P.; Clapp R., Perfluorinated Alkyl Substances: Emerging Insights Into Health Risks. New Solut. 2015, 25, (2), 147-63.
547
21. Li, K.; Gao, P.; Xiang, P.; Zhang, X.; Cui, X.; Ma, L. Q., Molecular mechanisms of
548
PFOA-induced toxicity in animals and humans: implications for health risks. Environ Int
549
2017, 99, 43-54.
ACS Paragon Plus Environment
37
Environmental Science & Technology
Page 38 of 53
550
22. UNEP, Stockholm convention on persistent organic pollutants. 2010.
551
23. Hagenaars, A.; Vergauwen, L.; De Coen, W.; Knapen, D., Structure–activity
552
relationship assessment of four perfluorinated chemicals using a prolonged zebrafish
553
early life stage test. Chemosphere 2011, 82, (5), 764-772.
554
24. Zheng, X.-M.; Liu, H.-L.; Shi, W.; Wei, S.; Giesy, J. P.; Yu, H.-X., Effects of
555
perfluorinated compounds on development of zebrafish embryos. Environ Sci Pollut Res
556
Int 2012, 19, (7), 2498-2505.
557
25. Ulhaq, M.; Carlsson, G.; Örn, S.; Norrgren, L., Comparison of developmental
558
toxicity of seven perfluoroalkyl acids to zebrafish embryos. Environ Toxicol Pharmacol
559
2013, 36, (2), 423-426.
560
26. Huang, H.; Huang, C.; Wang, L.; Ye, X.; Bai, C.; Simonich, M. T.; Tanguay, R. L.;
561
Dong, Q., Toxicity, uptake kinetics and behavior assessment in zebrafish embryos
562
following exposure to perfluorooctanesulphonicacid (PFOS). Aquat Toxicol 2010, 98, (2),
563
139-147.
ACS Paragon Plus Environment
38
Page 39 of 53
564 565
Environmental Science & Technology
27. OECD, Test No. 236: Fish Embryo Acute Toxicity (FET) Test, OECD Guidelines for the Testing of Chemicals, Section 2, 2013, OECD Publishing, Paris.
566
28. Brox, S.; Seiwert, B.; Küster, E.; Reemtsma, T., Toxicokinetics of polar chemicals
567
in zebrafish embryo (Danio rerio): Influence of physicochemical properties and of
568
biological processes. Environ Sci Tech 2016, 50, (18), 10264-10272.
569 570
29. Burnham, K.P., and Anderson, D. R., Model selection and multimodel inference. 2nd ed. New York: Springer. 2002.
571
30. Hagenaars, A.; Stinckens, E.; Vergauwen, L.; Bervoets, L.; Knapen, D., PFOS
572
affects posterior swim bladder chamber inflation and swimming performance of zebrafish
573
larvae. Aquat Toxicol. 2014, 157, 225-235.
574
31. Kühnert, A.; Vogs, C.; Altenburger, R.; Küster, E., The internal concentration of
575
organic substances in fish embryos — a toxicokinetic approach. Environ Toxicol Chem.
576
2013, 32, (8), 1819-1827.
ACS Paragon Plus Environment
39
Environmental Science & Technology
Page 40 of 53
577
32. Kirla, K. T.; Groh, K. J.; Steuer, A. E.; Poetzsch, M.; Banote, R. K.; Stadnicka-
578
Michalak, J.; Eggen, R. I.; Schirmer, K.; Kraemer, T., From the cover: Zebrafish larvae
579
are insensitive to stimulation by cocaine: Importance of exposure route and toxicokinetics.
580
Toxicol Sci 2016, 154, (1), 183-193.
581
33. Rawson, D. M.; Zhang, T.; Kalicharan, D.; Jongebloed, W. L., Field emission
582
scanning electron microscopy and transmission electron microscopy studies of the
583
chorion, plasma membrane and syncytial layers of the gastrula‐stage embryo of the
584
zebrafish Brachydanio rerio: a consideration of the structural and functional relationships
585
with respect to cryoprotectant penetration. Aquacult Res 2000, 31, (3), 325-336.
586
34. Bonsignorio, D.; Perego, L.; Del Giacco, L.; Cotelli, F., Structure and
587
macromolecular composition of the zebrafish egg chorion. Zygote 1996, 4, (2), 101-108.
588
35. Brox, S.; Ritter, A. P.; Küster, E.; Reemtsma, T., Influence of the perivitelline space
589
on the quantification of internal concentrations of chemicals in eggs of zebrafish embryos
590
(Danio rerio). Aquat Toxicol 2014, 157, 134-140.
ACS Paragon Plus Environment
40
Page 41 of 53
Environmental Science & Technology
591
36. Chen, F. F.; Gong, Z. Y.; Kelly, B. C., Bioavailability and bioconcentration potential
592
of perfluoroalkyl-phosphinic and -phosphonic acids in zebrafish (Danio rerio):
593
Comparison to perfluorocarboxylates and perfluorosulfonates. Sci Total Environ 2016,
594
568, 33-41.
595
37. Khazaee, M.; Ng, C. A., Evaluating parameter availability for physiologically based
596
pharmacokinetic (PBPK) modeling of perfluorooctanoic acid (PFOA) in zebrafish. Environ
597
Sci Process Impacts 2018, 20, (1), 105-119.
598 599
38. Ng, C. A.; Hungerbühler, K., Bioconcentration of perfluorinated alkyl acids: how important is specific binding? Environ Sci Tech 2013, 47, (13), 7214-7223.
600
39. Martin, J. W.; Mabury, S. A.; Solomon, K. R.; Muir, D. C. G., Bioconcentration and
601
tissue distribution of perfluorinated acids in rainbow trout (Oncorhynchus mykiss). Environ
602
Toxicol Chem 2003, 22, (1), 196-204.
603
40. Inoue, Y.; Hashizume, N.; Yakata, N.; Murakami, H.; Suzuki, Y.; Kikushima, E.;
604
Otsuka, M., Unique Physicochemical Properties of Perfluorinated Compounds and Their
ACS Paragon Plus Environment
41
Environmental Science & Technology
Page 42 of 53
605
Bioconcentration in Common Carp Cyprinus carpio L. Arch Environ Contam Toxicol.
606
2012, 62, (4), 672-680.
607
41. Naile, J.E., Khim, J.S., Hong, S., Park, J., Kwon, B.O., Rye, J.S., Hwang, J.H.,
608
Jones, P.D., Giesy, J.P, Distributions and bioconcentration characteristics of
609
perfluorinated compounds in environmental samples collected from the west coast of
610
Korea. Chemosphere 2013, 90, (2), 387-394.
611
42. Ulhaq, M.; Sundstrom, M.; Larsson, P.; Gabrielsson, J.; Bergman, A.; Norrgren, L.;
612
Orn, S., Tissue uptake, distribution and elimination of C-14-PFOA in zebrafish (Danio
613
rerio). Aquat Toxicol 2015, 163, 148-157.
614
43. Falk, S.; Failing, K.; Georgii, S.; Brunn, H.; Stahl, T., Tissue specific uptake and
615
elimination of perfluoroalkyl acids (PFAAs) in adult rainbow trout (Oncorhynchus mykiss)
616
after dietary exposure. Chemosphere 2015, 129, 150-156.
617
45. Wen, W.; Xia, X.; Hu, D.; Zhou, D.; Wang, H.; Zhai, Y.; Lin, H., Long-chain
618
perfluoroalkyl acids (PFAAs) affect the bioconcentration and tissue distribution of short-
619
chain PFAAs in zebrafish (Danio rerio). Environ Sci Technol 2017, 51, (21), 12358-12368.
ACS Paragon Plus Environment
42
Page 43 of 53
Environmental Science & Technology
620
46. Ulhaq, M.; Örn, S.; Carlsson, G.; Morrison, D. A.; Norrgren, L., Locomotor behavior
621
in zebrafish (Danio rerio) larvae exposed to perfluoroalkyl acids. Aquat Toxicol 2013, 144,
622
332-340.
623
47. Gomis, M. I.; Vestergren, R.; Borg, D.; Cousins, I. T., Comparing the toxic potency
624
in vivo of long-chain perfluoroalkyl acids and fluorinated alternatives. Environ Int 2018,
625
113, 1-9.
626
48. Kluever, N.; Vogs, C.; Altenburger, R.; Escher, B. I.; Scholz, S., Development of a
627
general baseline toxicity QSAR model for the fish embryo acute toxicity test.
628
Chemosphere, 2016, 164, 164-173.
629
49. Zhang, W.; Zhang, Y.; Zhang, H.; Wang, J.; Cui, R.; Dai, J., Sex differences in
630
transcriptional expression of FABPs in zebrafish liver after chronic perfluorononanoic acid
631
exposure. Environ Sci Technol 2012, 46, (9), 5175-5182.
632
50. Benninghoff, A. D.; Bisson, W. H.; Koch, D. C.; Ehresman, D. J.; Kolluri, S. K.;
633
Williams, D. E., Estrogen-like activity of perfluoroalkyl acids in vivo and interaction with
634
human and rainbow trout estrogen receptors in vitro. Toxicol Sci 2010, 120, (1), 42-58.
ACS Paragon Plus Environment
43
Environmental Science & Technology
Page 44 of 53
635
51. Rainieri, S.; Conlledo, N.; Langerholc, T.; Madorran, E.; Sala, M.; Barranco, A.,
636
Toxic effects of perfluorinated compounds at human cellular level and on a model
637
vertebrate. Food Chem Toxicol. 2017, 104, 14-25.
638
52. Ding, G.; Zhang, J.; Chen, Y.; Wang, L.; Wang, M.; Xiong, D.; Sun, Y., Combined
639
effects of PFOS and PFOA on zebrafish (Danio rerio) embryos. Arch Environ Contam
640
Toxicol 2013, 64, (4), 668-675.
ACS Paragon Plus Environment
44
Page 45 of 53
641
Environmental Science & Technology
Figures
and
Tables
642 643
Figure 1: Workflow of sample collection and sample preparation used in the TK
644
experiment (5x5 illustrates in the figure that five replicates of five ZFE were sampled per
645
time point and concentration).
ACS Paragon Plus Environment
45
Environmental Science & Technology
Page 46 of 53
646
PFOS PFHxS PFOA PFBA
Affected ZFE [%]
100 80 60 40 20 0 0.1 647
1
10 100 1000 Concentration [µM]
10000
648
Figure 2: Concentration-effect relationships for the PFAA show up to four orders of
649
magnitude differences in toxicity. Points represent the percent embryos that either were
650
dead or had one or several visible adverse effect such as non-inflated swim bladder,
651
pericardial or yolk sac edemas, and scoliosis. The curves indicate the best fit to the Hill
652
equation.
ACS Paragon Plus Environment
46
Page 47 of 53
Environmental Science & Technology
B concentration in embryo (log scale)
653
0.1
0
24
48
72
96
10 1 0.1
120
0
24
48
72
96
0
24
48
72
96
1000
10 1 0
24
48
72
96
24
48
72
96
10000
10
0
24
48
72
96
0
24
48 72 96 Time [h]
120
0
24
48
72
96
120 C1=7.8 µM C2=0.8 µM C3=0.6 µM
100 50
0
24
48
72
96
120
200
0
24
48
72
96
120
300
100 10 1
C1 = 340 µM C2 = 41 µM C3 = 21 µM
400
0
120
1000
1000
100
100
1
120
50
600
Conc in ZFE [µM]
Conc in ZFE [µM] 0
100
0
120
1000
100
10
100
0.1
120
C1=0.76 µM C2=0.08 µM C3=0.04 µM
150
Conc in ZFE [µM]
Conc in ZFE [µM]
1
150
0
120
1000
10
0.1
100
Conc in ZFE [µM]
0.01
1000
C concentration in embryo (linear scale) Conc. in ZFE [µM]
Conc. in ZFE [µM]
1
Conc in ZFE [µM]
Concentration [µM] Concentration [µM] Concentration [µM] Concentration [µM]
PFBA
PFOA
PFHxS
PFOS
A concentration in ZFE medium
0
24
48 72 96 Time [h]
120
C1=4800 µM C2=550 µM C3=240 µM
200 100 0
0
24
48 72 96 Time [h]
120
654
Figure 3: Time courses of PFAA concentrations in A) ZFE medium (n=6) and in ZFE
655
depicted on B) logarithmic scale and C) linear scale (n=5-10). The exposure
656
concentrations in the legends (C1 - C3) are given as measured initial concentrations. The
657
concentrations in the ZFE [µM (=µmol L-1 ZFE)] were obtained by dividing the measured
658
amount by published data on ZFE volume at different ages (according to Brox et al.28).
ACS Paragon Plus Environment
47
Environmental Science & Technology
Page 48 of 53
The lines represent simulated internal concentrations using a two-compartment model.
660
All values are stated as mean ± SE.
120
Net uptake rate [pmol h -1/Embryo]
PFOS
120
Net uptake rate [pmol h -1/Embryo]
Net uptake rate [pmol h -1/Embryo]
Net uptake rate [pmol h -1/Embryo]
659
0.8 0.6 0.4 0.2 0.0 -0.2
0
24
48
72
96
PFOA
4
2
0 0
24
48 72 Time [h]
96
PFHxS
0.8
C1 C2 C3
0.6 0.4 0.2 0.0 -0.2
0
24
48
72
96
120
24
48 72 Time [h]
96
120
PFBA
8 6 4 2 0 -2 -4 0
661 662
Figure 4: PFAA net uptake rates per individual embryo determined as derivative of
663
internal amount at three exposure concentrations (C1 - C3).
664 665 666
ACS Paragon Plus Environment
48
Page 49 of 53
Environmental Science & Technology
PFOS PFHxS PFOA PFBA
Affected ZFE [%]
100 80 60 40 20 0 667
1
10 100 1000 10000 100000 Internal concentration [µM]
668
Figure 5: Concentration-effect relationship based on internal concentrations for the PFAA
669
show up to three fold toxicity differences. Internal effect concentrations were gained by
670
multiplying BCFvol with effective exposure concentrations.The curves indicate the best fit
671
to the Hill equation based on the calculated internal effect concentrations.
672
ACS Paragon Plus Environment
49
Environmental Science & Technology
Page 50 of 53
673 674
Table 1: External concentration range tested for the studied PFAA and the estimated
675
effective external concentration at which 20% (EC20) and 50% (EC50) of the ZFE were
676
affected at 120 hpf. The EC20 and EC50 values are given as mean ± SE.
Present Study
Literature values
External concentration range [µM] PFOS
0.30 - 57
EC20
EC50
EC50
[µM]
[µM]
[µM]
2.3 ± 0.1
3.8 ± 0.1
2.0 a 2.8 b
PFHxS
0.4 - 330
59.8 ± 1.9
84.5 ± 1.4
-
PFOA
8.8 - 1400
355.6 ± 5.4
509 ± 5.5
270 a 840b
PFBA
3500 - 12000
~ 7000
-
> 14000 a 10200 b
677
a
Hagenaars et al.23 (effect values at 120 hpf), b Uhlaq et al.25 (effect values at 144 hpf)
ACS Paragon Plus Environment
50
Page 51 of 53
Environmental Science & Technology
678
Table 2: Initial exposure concentrations (Cw), areas under the concentration-time curve
679
(AUC), bioconcentration factors (BCFvol, BCFdry) and internal effect concentrations (IEC)
680
of the four PFAA studied at three exposure concentrations (C1 - C3). The IEC20 and IEC50
681
values were calculated as BCFvol x EC20 and BCFvol x EC50, respectively. The
682
experimental values from the present study are represented as mean ± SE.
Present study
PFOS
PFHxS
PFOA
PFBA
Literature values*
Present study
Cw
AUC
BCFvol
BCFdry
BCFdry
IEC20
IEC50
[µM]
[µM h]
[-]
[L kg-1]
[L kg-1]
[µM]
[µM]
C1
0.76 ± 0.01
8124± 18
240 ± 4
2600 ± 63
877 ± 189
1415 ± 301
C2
0.08 ± 0.01
1344 ± 13
370 ± 82
3900 ± 866
15 a 3200b 1100c 720d 2700f
C3
0.04 ± 0.002 7.8 ± 0.04
638 ± 18
510 ± 92
5400 ± 968
8483 ± 23
15 ± 0.7
160 ± 8
550b 9.6c 380f
1139± 134
1608 ± 184
170b 4.0c 3.1d 20 – 30e 11f
1805 ± 645
2584 ± 932
13b
~540
-
C1 C2
0.8 ± 0.01
893 ± 0
22 ± 0.5
230 ± 5
C3
0.6 ± 0.1 340 ± 24
762 ± 6
21 ± 3.8
220 ± 40
37924 ± 140
1.4 ± 0.1
15 ± 1
C1 C2
41 ± 2
15652 ± 55
6.9 ± 0.6
73 ± 8
C3
21 ± 1
6.8 ± 0.4
72 ± 4
C1
4800 ± 1050
10436 ± 41 19383 ± 151
0.04 ± 0.01
0.5 ± 0.1
ACS Paragon Plus Environment
51
Environmental Science & Technology
C2
550 ± 31
3344 ± 31
0.07 ± 0.01
0.8 ± 0.1
C3
240 ± 14
1710 ± 18
0.10 ± 0.01
1.2 ± 0.1
Page 52 of 53
683 684 685
et al.26 – zebrafish embryo, b Chen et al.36 – adult zebrafish, c Martin et al.39 – carcass of the rainbow trout (85 % - 96 % of total mass of chemical), d Inoue et al.40 – Cyprinus carpio L., e Uhlaq et al.42 – adult zebrafish, f Naile et al.41 – fish
686
* Whole-body BCFdry values for adult fish varied by one to two orders of magnitude for
687
the same PFAA, probably due to differences in exposure systems, exposure
688
concentrations, species and gender.
a Hung
ACS Paragon Plus Environment
52
Page 53 of 53
Environmental Science & Technology
For TOC art only
ACS Paragon Plus Environment