Subscriber access provided by University of Sussex Library
Article
Thyroid hormone disruption by water-accommodated fractions of crude oil and sediments affected by the Hebei Spirit oil spill in zebrafish and GH3 cells Sujin Kim, Ju Hae Sohn, Sung Yong Ha, Habyeong Kang, Un Hyuk Yim, Won Joon Shim, Jong Seong Khim, Dawoon Jung, and Kyungho Choi Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b00751 • Publication Date (Web): 04 May 2016 Downloaded from http://pubs.acs.org on May 5, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 29
Environmental Science & Technology
1
Thyroid hormone disruption by water-accommodated fractions of crude oil and
2
sediments affected by the Hebei Spirit oil spill in zebrafish and GH3 cells
3 4
Sujin Kim†, Ju Hae Sohn†, Sung Yong Ha‡, Habyeong Kang†, Un Hyuk Yim‡, Won Joon Shim‡,
5
Jong Seong Khim§, Dawoon Jung†, Kyungho Choi†,*
6
†School
7
‡Oil
8
Geoje 53201, Republic of Korea
9
§School
10
of Public Health, Seoul National University, Seoul 08826, Republic of Korea
and POPs Research Group, Korea Institute of Ocean Science and Technology (KIOST),
of Earth and Environmental Sciences & Research Institute of Oceanography, Seoul
National University, Seoul 08826, Republic of Korea
11 12
* Corresponding author
13
E-mail address:
[email protected] (K.Choi).
14
(Tel) 82-2-880-2738
15
(Fax) 82-2-745-9104
16 17
1
ACS Paragon Plus Environment
Environmental Science & Technology
18
Abstract
19
A crude oil and the coastal sediments that were affected by Hebei Spirit Oil Spill (HSOS) of
20
Taean, Korea were investigated for thyroid hormone disruption potentials. Water-
21
accommodated fractions (WAFs) of Iranian Heavy crude oil, the major oil type of HSOS, and
22
the porewater or leachate of sediment samples collected along the coast line of Taean were
23
tested for thyroid disruption using developing zebrafish, and/or rat pituitary GH3 cells. Major
24
polycyclic aromatic hydrocarbons (PAHs) and their alkylated forms were also measured from
25
the test samples. In zebrafish larvae, significant decreases in whole-body thyroxine (T4) and
26
triiodothyronine (T3) levels, along with transcriptional changes of thyroid regulating genes,
27
were observed following 5 day exposure to WAFs. In GH3 cells, transcriptions of thyroid
28
regulating genes were influenced following the exposure to the sediment samples, but the
29
pattern of the regulatory change was different from those observed from the WAFs.
30
Composition of PAHs and their alkylated homologs in the WAFs could partly explain this
31
difference. Our results clearly demonstrate that WAFs of crude oil can disrupt thyroid
32
function of larval zebrafish. Sediment samples also showed thyroid disrupting potentials in
33
the GH3 cell, even several years after the oil spill. Long-term ecosystem consequences of
34
thyroid hormone disruption due to oil spill deserve further investigation.
35 36
Keywords: Hebei Spirit Oil Spill, water-accommodated fraction, endocrine disruption,
37
zebrafish, GH3, hypothalamic-pituitary-thyroid (HPT) axis
38
2
ACS Paragon Plus Environment
Page 2 of 29
Page 3 of 29
Environmental Science & Technology
39
1. Introduction
40
Oil spills along the coastal shoreline can cause serious ecological and human health problems.
41
The Hebei Spirit Oil Spill (HSOS) accident which occurred in December 2007, near Taean,
42
Korea, is one of the largest oil spills in Korean history. Many studies have assessed ecological
43
and human health damages of the HSOS accident along the Taean coastline.1-5
44
Petroleum products, e.g., crude oil, are extremely complex mixtures, and therefore may
45
contain numerous toxic components. Polycyclic aromatic hydrocarbons (PAHs) or their
46
alkylated forms are considered to be among the primary contributors to oil-related adverse
47
effects.2, 6 Most available toxicological information on HSOS to date is limited to PAH-
48
associated responses such as aryl hydrocarbon receptor (AhR) binding affinity2 or
49
cytochrome P4501A (CYP1A) induction.5, 7 In addition, sex endocrine disruption and DNA
50
damage potentials of the oil spill have been documented from the environmental samples of
51
HSOS.3
52
Growing evidences indicate that crude oil and PAHs can cause developmental toxicity in fish,
53
e.g., cardiac malformation.8-12 Exposure to weathered Alaska North Slope crude oil led to
54
changes in cardiac morphology in developing zebrafish.9,
55
syndromes were reported in zebrafish embryos following exposure to Iranian Heavy crude
56
oil.13 Underlying mechanisms of such developmental changes are not yet clear. However,
57
thyroid hormone disruption could be in part responsible for these adverse effects, as thyroid
58
hormones are important regulators of normal differentiation and development of organs
59
including heart and vascular system. 14-16 Indeed, in teleost fish, thyroid hormones such as
60
thyroxine (T4) and triiodothyronine (T3) have known to play crucial roles in growth and
61
development of an individual.17, 18
62
Several studies have indicated that water soluble components of crude oil could disrupt
63
thyroid hormone balances in fish, including turbot (Scophthalmus maximus L.) and flounders 3
ACS Paragon Plus Environment
10
Similarly, cardiotoxicity
Environmental Science & Technology
64
(Pleuronectes flesus).19-22 Moreover, similar observations were reported from avian species.
65
In seabirds, e.g., black guillmot (Cepphus grille), herring gulls (Larus argentatus), and
66
Leach’s petrels (Oceanodroma leucorhoa), a single oral dosing to petroleum products caused
67
changes in plasma T4 level.23 Daily oral ingestion of crude oil resulted in increased plasma T3
68
and hepatic CYP levels in ducks (Anas platyrhynchos) as well.24 However, mechanisms
69
underlying thyroid disruption of oil or oil related compounds have not been fully understood
70
yet.
71
In the present study, we investigated thyroid disrupting potential of a crude oil and
72
environmental samples that were contaminated by the same crude oil spill. For this purpose,
73
the water-accommodated fractions (WAFs) of the crude oil, and the pore water or leachate of
74
the sediment samples which were affected by HSOS accident, were employed. WAFs have
75
been recommended for studying short-term aquatic effects of poorly soluble mixtures such as
76
crude oil.25, 26 These samples were first tested for thyroid hormone disruption potentials using
77
a rat pituitary gland cell line (GH3). GH3 cells have been used for screening thyroid related
78
activities of chemicals.27-29 For in vivo responses, zebrafish (Danio rerio) larvae were
79
employed but only for the WAFs of the oil, because of high salinity associated with the near-
80
coastal sediment samples. Through combined approach of zebrafish larvae and GH3 cells,
81
possible effects of the oil spill on thyroid hormone balance and its regulation via
82
hypothalamus-pituitary-thyroid (HPT) axis were investigated. In addition, the test samples
83
were analyzed for major PAHs and their alkylated analogues, in order to link to the observed
84
thyroid hormone disruption. The result of this study will help better understand the impacts of
85
oil components on thyroid system and development of aquatic organisms like fish.
86
4
ACS Paragon Plus Environment
Page 4 of 29
Page 5 of 29
Environmental Science & Technology
87
2. Materials and methods
88
2.1. Sample collection and preparation
89
WAFs were prepared using Iranian Heavy crude oil following a standardized method
90
published elsewhere.25, 26 Iranian Heavy crude oil is a dominant type of the oils released from
91
the Hebei Spirit.1, 30 Briefly, 87.5 g Iranian Heavy crude oil was carefully added onto the
92
surface of a 5 L carboy filled with 3.5 L of dechlorinated water (25 g oil/L water). The
93
mixture was stirred gently for 24 h to avoid formation of oil droplets, after which the water
94
phase was siphoned out and was utilized. Prepared WAFs were sent for chemical analysis,
95
and were immediately used for the zebrafish exposure study.
96
Two types of environmental samples were collected in September 2014 from two locations,
97
i.e., Sogeunri mudflat (latitude: 36° 48'46.47", longitude: 126° 10'57.56", henceforth ‘Site A’)
98
and Sinduri beach (latitude: 36° 50'20.95", longitude: 126° 11'10.24", henceforth ‘Site B’)
99
along the coast of Taean (Figure 1). The sampling locations were chosen based on visible
100
signs of oil residuals. From Site A, sediment samples were collected, pooled, and transported
101
to the laboratory under cold condition. Upon delivery to the laboratory, the sample was added
102
with dechlorinated water (1 g/mL), and was vigorously shaked for > 12 h at 4 oC. Then the
103
sediment pore water was obtained by centrifugation at 1000g for 1 h at 4 oC. From Site B,
104
leachate from the sand was collected on site, and moved to the laboratory under cold
105
condition. Both samples were stored at 4 oC until the toxicological studies and chemical
106
measurement.
107 108
2.2. Zebrafish embryo/larval exposure
109
Fertilized embryos of zebrafish (Danio rerio) were collected from 35 L culture tanks with
110
dechlorinated water, within 4 h post-fertilization (hpf). Dechlorinated water has been used for
111
maintaining fry and adult zebrafish in our laboratory. The eggs were randomly placed into 5
ACS Paragon Plus Environment
Environmental Science & Technology
112
glass beakers containing 500 mL of 0, 20, 30, 40, or 50% WAF solutions with four replicates
113
per group. For each replicate, 240 eggs were placed. Dechlorinated water was used for
114
control, and for dilution of WAFs. The exposure concentrations of WAFs were chosen below
115
a lethal level which was determined based on a preliminary range finding test. During 120 h
116
of the exposure period, embryo and larval survival, hatchability, and malformation rate were
117
recorded daily. Water quality parameters were measured at the beginning and the conclusion
118
of the exposure (Table S1 of Supporting information). At 120 h, wet body weight (mg) of 150
119
larvae was measured for each replicate and stored at -80 oC until thyroid hormone
120
measurement. For gene transcription analysis, 20 larvae were randomly chosen for each
121
replicate and were stored at -80 oC.
122
For thyroid hormone measurement, the zebrafish larvae (n=150) per replicate were
123
homogenized in 150 µL of sample diluent, and supernatants were collected following Yu et al.
124
with some modifications (For details refer to Supporting Information).31 Whole body T4 and
125
T3 were measured using commercial enzyme-linked immunosorbent assays (ELISA) kits
126
(Uscnlif, Wuhan, China; Cat no. E0453Ge for T3; Cat no. E0452Ge for T4) following the
127
manufacturer’s instructions. All the hormones were measured above the reported detection
128
limits (0.1 ng/mL for T3 and 1.2 ng/mL for T4).
129
For measurement of gene transcripts in zebrafish, 20 whole body larvae were homogenized
130
and total RNA was isolated using Maxwell®16 LEV simply RNA purification Tissue Kits
131
(Cat.#. AS1280, Promega, Madison, WI, USA). Complementary DNAs (cDNAs) were
132
synthesized using iScriptTM cDNA synthesis kits (BioRad Hercules, CA, USA) and
133
quantitative real-time polymerase chains reaction (qRT-PCR) were carried out with
134
LightCycler® 480 SYBR Green I Mastermix (Roche Diagnostics Ltd., Lewes, UK) and
135
LightCycler 480 instrument (Roche Applied Science, Indianapolis, IN, USA). The primer
136
sequences and PCR efficiencies are listed in Table S2. The transcription level of each target 6
ACS Paragon Plus Environment
Page 6 of 29
Page 7 of 29
Environmental Science & Technology
137
gene was normalized to that of a housekeeping gene, 18s rRNA. The housekeeping gene was
138
chosen based on our preliminary study where this gene showed the greatest stability
139
compared to other candidate genes such as β-actin, and elfa, in terms of the variability of Ct
140
values. The 18s rRNA gene has been successfully used in various zebrafish exposure studies
141
elsewhere.31, 32 Target genes were quantified by using the 2-∆∆CT method.33
142 143
2.3. GH3 cell exposure
144
Transcriptional changes of thyroid regulating genes were assessed with GH3 cells. For GH3
145
cell exposure, stock solutions of WAFs, the pore water of Site A, and the leachate of Site B,
146
were prepared in dechlorinated water with the final concentration of 0.1% v/v. T3 was chosen
147
as a positive control for GH3 cell assay, and was dissolved in a medium with dimethyl
148
sulfoxide (DMSO, 0.1% v/v). GH3 cell exposure was conducted following Kim et al.29 with
149
some modifications. In brief, trypsinized GH3 cells from culture plates were seeded into 24-
150
well plates (Costar 3526) at a density of 2.0×105 cells per well, and incubated for 16 h. After
151
the growth medium was replaced with serum-free medium, cells were incubated for
152
additional 8 h. Then the cells were dosed with serially diluted WAFs (0.0125, 0.025, and
153
0.05%), Site A pore water (0.00625, 0.0125, and 0.25%), and Site B leachate (0.0125, 0.025,
154
and 0.05%). For dilution, dechlorinated water was employed. The 0.1% v/v dechlorinated
155
water was used as a control. For positive control, 0.023, 0.23, or 2.3 nM of T3 and 0.1% v/v
156
DMSO were also dosed. The cells were incubated in triplicates for each treatment for 48 h. At
157
the end of the exposure, RNA isolation, synthesis of cDNA, and qRT-PCR for the GH3 cell
158
experiment were conducted as described elsewhere.29 Sequences and efficiencies of primers
159
are listed in Table S2. Cyclophilin was used as a housekeeping gene. The exposure
160
concentrations were determined at noncytotoxic range (cell proliferation >80% of control),
161
following a WST-1 cell proliferation assay (Roche Applied Science, Mannheim, Germany; 7
ACS Paragon Plus Environment
Environmental Science & Technology
162
data not shown).
163 164
2.4. Analysis of PAHs by gas chromatography
165
Concentrations of PAHs and alkylated PAHs were measured in the 50% WAF, Site A pore
166
water, and Site B leachate, using an HP 5890 GC equipped with an HP 5972 MS (Agilent,
167
Palo Alto, CA, USA). Detection limits for target PAHs in the samples ranged from 0.09 to
168
5.67 ng/ L depending on their physico-chemical properties. The recovery efficiencies of
169
PAHs were generally comparable to those reported in our previous works.34, 35 List of target
170
PAHs and detailed information are shown in Supporting Information.
171 172
2.5. Statistical analysis
173
Statistics were carried out using IBM SPSS Statistics (version 23.0; SPSS Inc., Chicago, IL,
174
USA). For zebrafish exposure study, normality and homogeneity of variances were evaluated
175
by Shapiro-Wilk’s test and Levene’s test, respectively. If data met the normality and
176
homogeneity assumptions, the difference among treatments and control was analyzed by one-
177
way analysis of variance (ANOVA) followed by Dunnett’s test. Otherwise, a non-parametric
178
Kruskal-Wallis test combined with Mann-Whitney U test was used. Dixon’s Q test was
179
performed for identification of outliers, which were considered for exclusion before statistical
180
analysis. For GH3 cell exposure results, ANOVA test followed by Dunnett’s test was
181
performed. Logarithmic or exponential data transformations were conducted when necessary.
182
In all statistical analyses, p values less than 0.05 was considered to be statistically significant.
183
All data are shown as mean±standard error of mean (SEM).
184
8
ACS Paragon Plus Environment
Page 8 of 29
Page 9 of 29
Environmental Science & Technology
185
3. Results
186
3.1. WAF exposure in zebrafish embryos/larvae
187
3.1.1. Survival, growth, and developmental effects
188
In zebrafish embryos/larvae exposed to 20, 30, 40, or 50% WAFs, no significant effects on
189
survival, hatchability, and body weight were observed until the conclusion of exposure (120 h)
190
(Table 1). However, among the exposure groups with higher WAF concentrations,
191
malformation with severe morphological changes such as pericardial edema or yolk sac
192
edema increased slightly (Figure S1).
193 194
3.1.2. Whole-body thyroid hormones
195
Exposure to WAFs for 120 h caused significant decrease of the whole body thyroid hormones
196
in zebrafish larvae (Figure 2). T4 concentrations were significantly decreased in zebrafish
197
exposed to 40% and 50% WAFs (Figure 2A). Treatment with 50% WAF resulted in
198
significant decrease of whole-body T3 in the larvae (Figure 2B).
199 200
3.1.3. Regulation of hypothalamic-pituitary-thyroid axis genes in zebrafish larvae
201
Following exposure to WAFs, several regulating genes of the HPT axis were affected in
202
zebrafish larvae (Figure 3). The thyroid stimulating hormone beta (tshβ) and thyroid hormone
203
receptor alpha (thrα) genes were slightly up-regulated although the results were not
204
statistically significant (Figure 3A). Regulation of corticotrophin-releasing hormone (crh),
205
thyroid stimulating hormone receptor (tshr) and thyroid hormone receptor beta (thrβ) genes
206
were not altered. However, the thyroglobulin (tg), NK2 homeobox1 (nkx2.1), paired box
207
protein 8 (pax8), and hematopoietically-expressed homeobox (hhex) genes were up-regulated
208
following exposure to WAFs (Figure 3B). 50% WAFs led to up-regulation of transthyretin
209
(ttr) gene by 20-fold and deiodinase type 1 (dio1) gene by 4.8-fold, compared to the control. 9
ACS Paragon Plus Environment
Environmental Science & Technology
210
In addition, transcription of uridine diphosphate glucuronosyltransferase (ugt) gene was
211
significantly increased by exposure to 40% and 50% WAFs (Figure 3C). The transcription of
212
cyp1a gene was markedly increased in a concentration-dependent manner, and the ahr2 gene
213
was significantly up-regulated at 40% and 50% WAFs (Figure 3D).
214 215
3.2. WAF and environmental sample exposure in GH3 cells
216
In rat pituitary GH3 cells, concentration-dependent down-regulation of tshβ gene was
217
observed following exposure to WAFs, and this pattern is similar to that observed by
218
exposure to T3 (Figure 4A). In contrast, the pore water of Site A sediment (Sogeunri mudflat)
219
and the leachate from Site B (Sinduri beach) significantly up-regulated tshβ gene of GH3
220
cells. Both Sites A and B samples also up-regulated dio1 gene in GH3 cells, unlike the
221
changes observed from WAFs (Figure 4B). The transcription of dio2, thrα, and thrβ genes
222
was not significantly changed with an exception of the lowest concentration of the Site A
223
pore water where dio2 gene was up-regulated (Figure S2).
224 225
3.3. Chemical characterization of the WAF and environmental samples
226
PAHs and alkylated PAHs measured in the WAF (50%), Site A pore water, and Site B
227
leachate were different in levels and composition (Figure 5 and Table S3). Total PAHs,
228
calculated as a sum of PAHs and alkylated PAHs, were the highest in Site A pore water
229
(16,955 ng/L), followed by the 50% WAF (8,565 ng/L), and site B leachate (672 ng/L). The
230
PAHs with two benzene rings, such as naphthalene (Na) and dibenzothiophene (DBT), and
231
their alkylated homologues were generally dominant in all three test samples, perhaps
232
because of their greater water solubility. In WAFs, non-alkylated and C1-Na were dominant
233
while concentrations of DBT compounds were very low or below the detection limit. On the
234
contrary, the pore water from Site A was dominated by the alkylated DBT, whereas parent 10
ACS Paragon Plus Environment
Page 10 of 29
Page 11 of 29
Environmental Science & Technology
235
and alkyl-Na concentrations were relatively low. In the pore water sample (Site A), alkylated
236
homologues of three to four ring PAHs such as fluorine, phenanthrene, and chrysene were
237
also detected. For the Site B leachate, although total PAH concentration was low, relative
238
composition of PAHs was similar to those of Site A, e.g., with C2- and C3-DBT as dominant
239
constituents.
240 241
4. Discussion
242
4.1. Thyroid disruption by WAFs in zebrafish and GH3 cells
243
Significant decrease of thyroid hormone levels in developing zebrafish by WAF exposure
244
(Figure 2) clearly demonstrates the thyroid endocrine disruption potential of the crude oil.
245
The regulatory changes of thyroid regulating genes in GH3 cells (Figure 4) also support the
246
thyroid disrupting effects of the WAFs.
247
The significant decreases of whole body T4 and T3 contents (Figure 2) may have implications
248
in fish development, as thyroid hormones play pivotal roles in cellular differentiation and
249
neurodevelopment.18, 36, 37 It is in agreement with a previous report that observed reduced
250
plasma T4 concentrations in flounder by the water soluble fractions (WSFs) of Omani crude
251
oil.21 On the contrary, exposure to WSFs of a crude oil (BP, Wytch Farm, Dorset, UK) for 6 h
252
elevated whole body T4 without alteration in T3 in turbot (Scophthalmus maximus L.)
253
larvae.20 These reports along with our present observation suggest that oil components can
254
disrupt thyroid hormones but the pattern of disruption may vary by species and exposure
255
duration.
256
Significant down-regulation of tshβ gene in GH3 cells by WAFs and T3 exposure (Figure 4A)
257
indicates that WAFs may act in the same way as T3 in the pituitary gland, i.e., down-
258
regulation of tshβ gene expression in the pituitary cells. TSH released from pituitary gland
259
stimulates thyroid gland and stimulates synthesis of thyroid hormones.17 Interestingly, in 11
ACS Paragon Plus Environment
Environmental Science & Technology
260
zebrafish larvae, tshβ gene showed up-regulating pattern following WAFs exposure (Figure
261
3A), suggesting compensatory efforts of the zebrafish pituitary gland in response to
262
decreased thyroid hormones (Figure 2). Several other genes that stimulate synthesis of
263
thyroid hormones or development of thyroid follicle, such as tg, nkx2.1, pax8, or hhex, were
264
also up-regulated (Figure 3B), implying compensatory efforts of stimulating thyroid hormone
265
synthesis in zebrafish.
266
The reason for decreased whole body thyroid hormones in fish (Figure 2) is not clear, but
267
could be found from up-regulation of ugt gene in zebrafish which was observed after
268
exposure to 40% and 50% WAFs (Figure 3C). UGT enzyme plays a critical role in
269
inactivation and excretion of many exogenous and endogenous compounds, including T4.
270
Induction of UGT may enhance glucuronidation and facilitate elimination of T4, finally
271
leading to decreased whole body T4 contents in zebrafish. Several reports indicated that
272
induction of UGT may possibly lead to decreased T4 levels.38, 39 In thyroidectomized male
273
rats, serum thyroid hormones could be reduced through UGT induction mechanism.38 Also in
274
zebrafish, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), one of the most potent inducers of
275
CYP1A, was reported to increase hepatic UGT and result in decreased T4.40 ,41
276
Up-regulation of ttr gene (Figure 3C) may reflect a compensatory effort against low thyroid
277
hormone levels (Figure 2). TTR is an important T3 binding protein in teleost fish and
278
amphibia,42, 43 and is responsible for transportation of the hormones to the target peripheral
279
tissues. At the same time, binding to TTR can delay metabolic elimination of free thyroid
280
hormones from the circulation.
281
Up-regulation of dio1 gene in zebrafish following exposure to WAFs (Figure 3C) may be also
282
seen as compensation against the decrease in T3 level. Among two types of deiodinases, i.e.,
283
DIO1 and DIO2,44 DIO1 especially in liver plays an important role in conversion of T4 to
284
T3.45 In tilapia, up-regulation of hepatic dio1 gene transcription was reported when the fish 12
ACS Paragon Plus Environment
Page 12 of 29
Page 13 of 29
Environmental Science & Technology
285
were under hypothyroidism.46 Meanwhile, DIO2 plays a role in catalytic conversion of T4 to
286
active T3, generally in euthyroid state.45,
287
observed in GH3 cells (Figure 4B) was different from its transcriptional changes observed in
288
whole body zebrafish larvae (Figure 3C). This discrepancy may be attributed to lack of
289
feedback in the rat pituitary cell line, which may not reflect regulatory feedback efforts that
290
could be seen in vivo. In addition, differential expression of those genes by organs should be
291
considered. The dio2 gene is predominant in pituitary gland while dio1 gene is abundantly
292
expressed in liver, kidney, and thyroid.47-49
293
Significant up-regulation of cyp1a gene in the fish larvae following the exposure to WAFs
294
(Figure 3D) is in line with our expectation, as CYP1A induction has been considered as an
295
indication of PAH exposure.50, 51 The ahr2 gene was up-regulated, although the extent of
296
increase was relatively smaller. Between two orthologs of AhR in fish, namely AhR1 and
297
AhR2, ahr2 transcripts are more abundant and seem to be functionally important in
298
xenobiotic metabolism.52-54 These transcriptional changes of cyp1a and ahr2 can be
299
supported by chemical composition of the WAFs with low molecular PAHs such as Na and
300
fluorene being dominant (Figure 5).
301
Taken together, our observations suggest that WAFs of crude oil can lower thyroid hormone
302
levels in zebrafish larvae, possibly through hepatic UGT induction. However, interpretation
303
of the results from this short exposure duration (120 h) warrants caution. Because WAFs are
304
not generally recommended for studies involving longer term exposure, we could not
305
continue the exposure long enough to understand longer term consequences of the oil
306
pollution. Study design that involves long-term exposure to oil spill, and also includes time
307
series observations on gene expression, hormones, and thyroid histology, deserves further
308
consideration.
47
In our present study, dio1 gene transcription
309 13
ACS Paragon Plus Environment
Environmental Science & Technology
310
4.2. Potential ecotoxicological impact of oil spill
311
Two environmental samples collected from the locations affected by HSOS, e.g., Site A
312
sediment pore water (Sogeunri mudflat) and Site B leachate (Sinduri beach) showed clear up-
313
regulation patterns of major thyroid genes such as tshβ and dio1 gene in GH3 cells (Figure 4),
314
suggesting thyroid disrupting effects of environmental samples. However, the regulatory
315
changes in GH3 cells from these environmental samples were different from those observed
316
in WAFs. Why are the responses observed from oil spill affected area different from those
317
observed in WAFs? First, as the sampling areas are located near a city and a popular beach,
318
we cannot rule out the possibility of anthropogenic contamination in our environmental
319
samples, not directly originated from the oil spill. On the other hand, effects of weathering
320
which might lead to clear difference in composition of PAHs between WAFs and the
321
environmental samples (Figure 5) can in part explain the observed difference of thyroid
322
related gene regulations. In our present study, the WAFs of Iranian Heavy crude oil contained
323
higher concentrations of parent and alkyl-Na, while the environmental samples collected
324
from Sites A and B showed greater proportions of alkyl-DBT, and other heavier PAHs. The
325
weathering process such as evaporation, dissolution, microbial and photochemical
326
degradation, generally shift the PAH composition to the dominance of heavier and more
327
alkylated compounds.13, 55 It should be noted that our knowledge on thyroid disrupting effects
328
of low molecular weight PAHs is quite limited. Only a few studies reported thyroid disrupting
329
potentials of PAHs with more rings, such as anthrancene, phenanthrene, pyrene, and
330
benzopyrene56 or by some hydroxylated homologues.57 Further investigations on thyroid
331
disrupting effects of PAHs and their alkylated analogs are warranted.
332
The results of our study clearly show that oil spill can cause thyroid hormone disruption in
333
fish. Developing zebrafish and GH3 cells could be used to understand the mechanisms of the
334
adverse effects. Screening areas with thyroid disrupting potential and following up long-term 14
ACS Paragon Plus Environment
Page 14 of 29
Page 15 of 29
335
Environmental Science & Technology
consequences of thyroid hormone disruption along coastal ecosystem warrant further studies.
336 337
Acknowledgment
338
This study was supported by the “Oil Spill Environmental impact Assessment and
339
Environmental Restoration (PM56951)” funded by the Ministry of Oceans and Fisheries of
340
Korea. We would like to thank NeoEnBiz Co. (Bucheon, Korea) for assistance in sediment
341
sampling.
342 343
Supporting Information Available
344
Additional information, including analytical methods for hormones, genes, and PAHs, and
345
major results such as morphology of fish, gene regulation of the cell, and the levels of PAHs
346
in the samples, are provided. This information is available free of charge via the Internet at
347
http://pubs.acs.org.
348
15
ACS Paragon Plus Environment
Environmental Science & Technology
349
References
350
(1)
351
Seo, J.; Lee, S.; Park, J.; Lee, W.; Choi, Y.; Lee, K. T.; Kim, C. K.; Shim, W. J.; Naile, J. E.;
352
Giesy, J. P., Two years after the Hebei Spirit oil spill: residual crude-derived hydrocarbons
353
and potential AhR-mediated activities in coastal sediments. Environ Sci Technol 2012, 46 (3),
354
1406-1414.
355
(2)
356
Shim, W. J.; Jung, J.; Giesy, J. P.; Khim, J. S., Effect-directed analysis and mixture effects of
357
AhR-active PAHs in crude oil and coastal sediments contaminated by the Hebei Spirit oil
358
spill. Environ Pollut 2015, 199, 110-228.
359
(3)
360
Y.; Shim, W. J.; Takeda, S.; Giesy, J. P.; Choi, K., Genotoxicity and endocrine-disruption
361
potentials of sediment near an oil spill site: two years after the Hebei Spirit oil spill. Environ
362
Sci Technol 2011, 45 (17), 7481-7488.
363
(4)
364
mounted to clean the Hebei Spirit oil spill in Taean, Korea. Mar Pollut Bull 2010, 60 (1), 51-
365
57.
366
(5)
367
N. S.; Addison, R. F.; Shim, W. J., Biomarker responses in pelagic and benthic fish over 1
368
year following the Hebei Spirit oil spill (Taean, Korea). Mar Pollut Bull 2011, 62 (8), 1859-
369
1866.
370
(6)
371
Handbook of Ecotoxicology, 2nd ed.; Hoffman, D. J.; Rattner, B. A.; Burton, A. G.; Cairns, J.,
372
Eds.; CRC Press LLC: Boca Raton, FL, 2003; pp 341– 371.
373
(7)
Hong, S.; Khim, J. S.; Ryu, J.; Park, J.; Song, S. J.; Kwon, B. O.; Choi, K.; Ji, K.;
Hong, S.; Lee, S.; Choi, K.; Kim, G. B.; Ha, S. Y.; Kwon, B. O.; Ryu, J.; Yim, U. H.;
Ji, K.; Seo, J.; Liu, X.; Lee, J.; Lee, S.; Lee, W.; Park, J.; Khim, J. S.; Hong, S.; Choi,
Sim, M. S.; Jo, I. J.; Song, H. G., Acute health problems related to the operation
Jung, J.-H.; Kim, M.; Yim, U. H.; Ha, S. Y.; An, J. G.; Won, J. H.; Han, G. M.; Kim,
Albers, P. H., Petroleum and Individual Polycyclic Aromatic Hydrocarbons.
Ramachandran, S. D.; Sweezey, M. J.; Hodson, P. V.; Boudreau, M.; Courtenay, S. C.; 16
ACS Paragon Plus Environment
Page 16 of 29
Page 17 of 29
Environmental Science & Technology
374
Lee, K.; King, T.; Dixon, J. A., Influence of salinity and fish species on PAH uptake from
375
dispersed crude oil. Mar Pollut Bull 2006, 52 (10), 1182-1189.
376
(8)
377
oil on embryonic development of Pacific herring. Trans Am Fish Soc 1979, 108 (1), 70-75.
378
(9)
379
Larsen, M.; Stekoll, M. S.; Rice, S. D.; Collier, T. K., Sublethal exposure to crude oil during
380
embryonic development alters cardiac morphology and reduces aerobic capacity in adult fish.
381
Proc Natl Acad Sci 2011, 108 (17), 7086-7090.
382
(10)
383
Aryl hydrocarbon receptor-independent toxicity of weathered crude oil during fish
384
development. Environ Health Perspect 2005, 113 (12), 1755-1762.
385
(11)
386
J. E.; Short, J. W., Delayed effects on growth and marine survival of pink salmon
387
Oncorhynchus gorbuscha after exposure to crude oil during embryonic development. Mar
388
Ecol Prog Ser 2000, 208, 205-216.
389
(12)
390
Valdez oil spill. Energy Sources 2003, 25 (6), 509-517.
391
(13)
392
Incardona, J. P., Geologically distinct crude oils cause a common cardiotoxicity syndrome in
393
developing zebrafish. Chemosphere 2013, 91 (8), 1146-1155.
394
(14)
395
N Engl J Med 2001, 344 (7), 501-509.
396
(15)
397
cold acclimation in zebrafish (Danio rerio). J Exp Biol 2014, 217 (5), 718-725.
398
(16)
Smith, R. L.; Cameron, J. A., Effect of water soluble fraction of Prudhoe Bay crude
Hicken, C. E.; Linbo, T. L.; Baldwin, D. H.; Willis, M. L.; Myers, M. S.; Holland, L.;
Incardona, J. P.; Carls, M. G.; Teraoka, H.; Sloan, C. A.; Collier, T. K.; Scholz, N. L.,
Heintz, R. A.; Rice, S. D.; Wertheimer, A. C.; Bradshaw, R. F.; Thrower, F. P.; Joyce,
Short, J., Long-term effects of crude oil on developing fish: lessons from the Exxon
Jung, J.-H.; Hicken, C. E.; Boyd, D.; Anulacion, B. F.; Carls, M. G.; Shim, W. J.;
Epstein, F. H.; Klein, I.; Ojamaa, K., Thyroid hormone and the cardiovascular system.
Little, A. G.; Seebacher, F., Thyroid hormone regulates cardiac performance during
Klein, I.; Danzi, S., Thyroid disease and the heart. Circulation 2007, 116 (15), 172517
ACS Paragon Plus Environment
Environmental Science & Technology
399
1735.
400
(17)
401
Rev Fish Biol Fish 1993, 3 (4), 299-347.
402
(18)
403
Canario, A.; Sweeney, G., Thyroid hormones in growth and development of fish. Comp
404
Biochem Physiol C Toxicol Pharmacol 2001, 130 (4), 447-459.
405
(19)
406
responses of larval and juvenile turbot exposed to the water‐soluble fraction of crude oil. J
407
Fish Biol 1997, 50 (5), 953-964.
408
(20)
409
Scophthalmus maximus L., exposed to sub-lethal concentrations of petroleum hydrocarbons.
410
Fish Physiol Biochem 1997, 17 (1-6), 433-439.
411
(21)
412
respiratory and haematological parameters in flounder exposed to the water soluble fraction
413
of crude oil. J Fish Biol 1996, 49 (6), 1291-1305.
414
(22)
415
teleost fish thyroid. Environ Toxicol Chem 2004, 23 (7), 1680-1701.
416
(23)
417
caused by ingested oil. Environ Res 1981, 24 (1), 6-14.
418
(24)
419
hormones and on the mixed function oxidase system in ducks. Comp Biochem Physiol C
420
Comp Pharmacol 1990, 95 (2), 213-216.
421
(25)
422
Standardization of the preparation and quantitation of water-accommodated fractions of
423
petroleum for toxicity testing. Mar Pollut Bull 2000, 40 (11), 1007-1016.
Eales, J.; Brown, S., Measurement and regulation of thyroidal status in teleost fish.
Power, D.; Llewellyn, L.; Faustino, M.; Nowell, M.; Björnsson, B. T.; Einarsdottir, I.;
Stephens, S.; Alkindi, A.; Waring, C.; Brown, J., Corticosteroid and thyroid
Stephens, S.; Brown, J.; Frankling, S., Stress responses of larval turbot,
Alkindi, A.; Brown, J.; Waring, C.; Collins, J., Endocrine, osmoregulatory,
Brown, S. B.; Adams, B. A.; Cyr, D. G.; Eales, J. G., Contaminant effects on the
Peakall, D.; Tremblay, J.; Kinter, W.; Miller, D., Endocrine dysfunction in seabirds
Jenssen, B. M.; Ekker, M.; Zahlsen, K., Effects of ingested crude oil on thyroid
Singer, M.; Aurand, D.; Bragin, G.; Clark, J.; Coelho, G.; Sowby, M.; Tjeerdema, R.,
18
ACS Paragon Plus Environment
Page 18 of 29
Page 19 of 29
Environmental Science & Technology
424
(26)
425
Tjeerdema, R. S. Making, measuring, and using water-accommodated fractions of petroleum
426
for toxicity testing. In International Oil Spill Conference Porceedings; American Petroleum
427
Institute: Washington 2001; pp 1269-1274.
428
(27)
429
thyroid hormone function in pituitary rat GH3 cells. Mol Cell Endocrinol 2005, 244 (1), 31-
430
41.
431
(28)
432
T-screen to quantify functional potentiating, antagonistic and thyroid hormone-like activities
433
of poly halogenated aromatic hydrocarbons (PHAHs). Toxicol In Vitro 2006, 20 (4), 490-498.
434
(29)
435
triphenyl phosphate, an organophosphate flame retardant, in zebrafish (Danio rerio)
436
embryos/larvae, and in GH3 and FRTL-5 cell lines. Aquat Toxicol 2015, 160, 188-196.
437
(30)
Korea Coast Guard, KCG 2008 White Paper, 2008.
438
(31)
Yu, L.; Deng, J.; Shi, X.; Liu, C.; Yu, K.; Zhou, B., Exposure to DE-71 alters thyroid
439
hormone levels and gene transcription in the hypothalamic–pituitary–thyroid axis of zebrafish
440
larvae. Aquat Toxicol 2010, 97 (3), 226-233.
441
(32)
442
H., Effects of tris (1, 3-dichloro-2-propyl) phosphate and triphenyl phosphate on receptor-
443
associated mRNA expression in zebrafish embryos/larvae. Aquat Toxicol 2013, 128, 147-157.
444
(33)
445
time quantitative PCR and the 2-∆∆CT method. Methods 2001, 25 (4), 402-408.
446
(34)
447
Jung, J.; Giesy, J. P., Bioaccessibility of AhR-active PAHs in sediments contaminated by the
448
Hebei Spirit oil spill: Application of Tenax extraction in effect-directed analysis.
Singer, M. M.; Aurand, D. V.; Coelho, G. M.; Bragin, G. E.; Clark, J. R.; Sowby, M.;
Ghisari, M.; Bonefeld-Jorgensen, E. C., Impact of environmental chemicals on the
Schriks, M.; Vrabie, C. M.; Gutleb, A. C.; Faassen, E. J.; Rietjens, I. M.; Murk, A. J.,
Kim, S.; Jung, J.; Lee, I.; Jung, D.; Youn, H.; Choi, K., Thyroid disruption by
Liu, C.; Wang, Q.; Liang, K.; Liu, J.; Zhou, B.; Zhang, X.; Liu, H.; Giesy, J. P.; Yu,
Livak, K. J.; Schmittgen, T. D., Analysis of relative gene Expression data using real-
Hong, S.; Yim, U. H.; Ha, S. Y.; Shim, W. J.; Jeon, S.; Lee, S.; Kim, C.; Choi, K.;
19
ACS Paragon Plus Environment
Environmental Science & Technology
Page 20 of 29
449
Chemosphere 2016, 144, 706-712.
450
(35)
451
Incardona, J. P.; Linbo, T. L.; Kwon, J.-H., Differential Toxicokinetics Determines the
452
Sensitivity of Two Marine Embryonic Fish Exposed to Iranian Heavy Crude Oil. Environ Sci
453
Technol 2015, 49 (22), 13639-13648.
454
(36)
455
implications for neurodevelopmental research on the effects of thyroid hormone disruptors.
456
Toxicol Ind Health 1998, 14 (1-2), 85-101.
457
(37)
458
and retinoic acid in timing oligodendrocyte development. Development 1994, 120 (5), 1097-
459
1108.
460
(38)
461
hormone levels in rats by an extrathyroidal mechanism. Toxicol Appl Pharmacol 1992, 113
462
(1), 36-42.
463
(39)
464
primary
465
glucuronosyltranferases by methylcholanthrene, clofibrate, and dexamethasone alone and in
466
combination. Drug Metab Dispos 2000, 28 (1), 34-37.
467
(40)
468
after a single dose of 2,3,7,8-tetrachlorodibenzo-p-dioxin. Endocrinol 1977, 101 (1), 292-296.
469
(41)
470
thyroid hormones and the thyroid gland. Endocr Rev 1991, 12 (2), 135-150.
471
(42)
472
the thyroid hormone-binding protein, transthyretin. Gen Comp Endocrinol 2000, 119 (3),
473
241-255.
Jung, J.-H.; Kim, M.; Yim, U. H.; Ha, S. Y.; Shim, W. J.; Chae, Y. S.; Kim, H.;
Hauser, P.; McMillin, J. M.; Bhatara, V. S., Resistance to thyroid hormone:
Barres, B.; Lazar, M.; Raff, M., A novel role for thyroid hormone, glucocorticoids
Barter, R. A.; Klaassen, C. D., UDP-glucuronosyltransferase inducers reduce thyroid
Jemnitz, K.; Veres, Z.; Monostory, K.; Vereczkey, L., Glucuronidation of thyroxine in monolayer
cultures
of
rat
hepatocytes:
in
vitro
induction
of
UDP-
Bastomsky, C. H., Enhanced thyroxine metabolism and high uptake goiters in rats
Curran, P. G.; DeGroot, L. J., The effect of hepatic enzyme-inducing drugs on
Power, D.; Elias, N.; Richardson, S.; Mendes, J.; Soares, C.; Santos, C., Evolution of
20
ACS Paragon Plus Environment
Page 21 of 29
Environmental Science & Technology
474
(43)
475
characterization of thyroid‐hormone‐binding protein from masu salmon serum. Eur J
476
Biochem 1999, 265 (3), 944-949.
477
(44)
478
an index of chemical disruption of thyroid hormone homeostasis and thyroidal status in fish.
479
In Environmental Toxicology and Risk Assessment: Standardization of Biomerkers for
480
Endocrine Disruption and Environmental Assessment, Henshel, D. S.; Black, M. C.; Harrase,
481
M. C., Eds. American Society for Testing and Materials: West Conshohocekn, PA, 1999; Vol.
482
Eighth Volume, pp 136-153.
483
(45)
484
M. M.; David Furlow, J.; Kavlock, R.; Köhrle, J.; Opitz, R., Mechanism-based testing
485
strategy using in vitro approaches for identification of thyroid hormone disrupting chemicals.
486
Toxicol In Vitro 2013, 27 (4), 1320-1346.
487
(46)
488
Sanders, J. P.; Visser, T. J.; Kühn, E. R.; Darras, V. M., Hypothyroidism induces type I
489
iodothyronine deiodinase expression in tilapia liver. Gen Comp Endocrinol 2001, 124 (3),
490
333-342.
491
(47)
492
hormone action. J Clin Invest 2006, 116 (10), 2571-2579.
493
(48)
494
Physiol Rev 2014, 94 (2), 355-382.
495
(49)
496
Zeöld, A.; Bianco, A. C., Cellular and Molecular Basis of Deiodinase-Regulated Thyroid
497
Hormone Signaling 1. Endocrb Rev 2008, 29 (7), 898-938.
498
(50)
Yamauchi, K.; Nakajima, J. i.; Hayashi, H.; Hara, A.; Yamauchi, K., Purification and
Eales, J. G.; Brown, S. B.; Cyr, D. B.; Adams, B. A.; Finnson, K. R., Deiodination as
Murk, A. J.; Rijntjes, E.; Blaauboer, B. J.; Clewell, R.; Crofton, K. M.; Dingemans,
Van der Geyten, S.; Toguyeni, A.; Baroiller, J.-F.; Fauconneau, B.; Fostier, A.;
Bianco, A. C.; Kim, B. W., Deiodinases: implications of the local control of thyroid
Mullur, R.; Liu, Y.-Y.; Brent, G. A., Thyroid Hormone Regulation of Metabolism.
Gereben, B.; Zavacki, A. M.; Ribich, S.; Kim, B. W.; Huang, S. A.; Simonides, W. S.;
Whyte, J. J.; Jung, R.; Schmitt, C.; Tillitt, D., Ethoxyresorufin-O-deethylase (EROD) 21
ACS Paragon Plus Environment
Environmental Science & Technology
499
activity in fish as a biomarker of chemical exposure. Crit Rev Toxicol 2000, 30 (4), 347-570.
500
(51)
501
pollution. In Toxicology in Transition, Springer: 1995; pp 80-95.
502
(52)
503
Interact 2002, 141 (1-2), 131-160.
504
(53)
505
among three members of the vertebrate aryl hydrocarbon receptor family: AHR repressor,
506
AHR1, and AHR2. J Biol Chem 2002, 277 (9), 6949-6959.
507
(54)
508
two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family.
509
Proc Natl Acad Sci 1997, 94 (25), 13743-13748.
510
(55)
511
and weathered crude oil: Toxic potential and the toxicity of saturated mixtures. Environ
512
Toxicol Chem 2007, 26 (1), 24-36.
513
(56)
514
in response to various chemicals. J Environ Monit 2012, 14 (8), 2121-2126.
515
(57)
516
L., Phthalates, Polycyclic Aromatic Hydrocarbons and Perchlorate Associations with Thyroid
517
Hormones During Pregnancy: Results from a Pilot Study. Epidemiology 2008, 19 (6), S235.
Goksøyr, A., Use of cytochrome P450 1A (CYP1A) in fish as a biomarker of aquatic
Hahn, M. E., Aryl hydrocarbon receptors: diversity and evolution. Chem Biol
Karchner, S. I.; Franks, D. G.; Powell, W. H.; Hahn, M. E., Regulatory interactions
Hahn, M. E.; Karchner, S. I.; Shapiro, M. A.; Perera, S. A., Molecular evolution of
Di Toro, D. M.; McGrath, J. A.; Stubblefield, W. A., Predicting the toxicity of neat
Song, M.; Kim, Y. J.; Park, Y. K.; Ryu, J. C., Changes in thyroid peroxidase activity
Peck, J.; Robledo, C.; Neas, B.; Calafat, A.; Sjodin, A.; Blount, B.; Wild, R.; Cowan,
518
22
ACS Paragon Plus Environment
Page 22 of 29
Page 23 of 29
Environmental Science & Technology
519
520 521
Abstract art
23
ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 29
522
Table 1. Effects of WAF exposure on survival, hatchability, malformation rate and wet body
523
weight in zebrafish at 120 hpf Control
20% WAF
30% WAF
40% WAF
50% WAF
Embryo survival (%)
89.2±2.59
93.7±1.81
91.8±3.09
92.8±0.524
90.5±3.03
Hatchability (%)
98.2±0.536
98.2±0.418
97.5±0.617
97.4±0.569
98.6±0.579
Larval survival (%)
97.8±0.429
97.5±0.442
97.2±0.732
96.9±0.295
94.1±0.921
6.32±1.69
6.07±0.806
7.71±0.892
10.4±1.87
149±12.1
149±8.36
157±9.54
161±4.67
Malformation
rate 6.99±0.876
(%) Body weight (mg)
167±14.8
524
Embryo survival (%), the percentage of surviving embryos among total fertilized eggs;
525
hatchability (%), percentage of hatchling among the live embryos; larval survival (%),
526
percentage of surviving larvae among the hatched; malformation rate (%), percentage of
527
malformed individuals, including dead fish; body weight (mg), wet weight of 150 zebrafish
528
larvae per replicate. Results are shown as mean ± SEM of four replicates.
529 530 531
24
ACS Paragon Plus Environment
Page 25 of 29
Environmental Science & Technology
532 533 534
Figure 1. Location of two sampling sites near Taean coast, Korea. Site A is located in
535
Sogeunri mudflat, and Site B in Sinduri beach. (Satellite photo provided by Google Earth
536
version 7.1.5. 2015. Image providers are shown in the bottom of each map: the upper left map
537
was provided by Data SIO, NOAA, U.S. Navy, NGA, GEBCO; Data Japan Hydrographic
538
Association; Image Landsat. The other map was provided by Image © 2016 TerraMetrics;
539
Image © 2016 CNES/Astrium.)
540 541
25
ACS Paragon Plus Environment
Environmental Science & Technology
542 543 544
Figure 2. Concentration (ng/g wet weight) of (A) T4, and (B) T3 measured in whole body
545
zebrafish larvae at 120 hpf following exposure to 0, 20, 30, 40, or 50% WAFs of Iranian
546
Heavy oil. The results are shown as mean ± SEM of four replicates (n=4), and each replicate
547
includes 150 larvae. Asterisks (p*