Bioaccumulation, Biotransformation, and Toxicity of BDE-47, 6-OH

Jan 7, 2015 - BDE-47 induces oxidative stress, activates MAPK signaling pathway, and elevates de novo lipogenesis in the copepod Paracyclopina nana. M...
0 downloads 20 Views 1MB Size
Article Environmental Science & Bioaccumulation, Technology is published by the American Chemical biotransformation and Society. 1155 Sixteenth Street N.W., Washington, toxicity ofDC BDE-47, 20036 by RUTGERS Subscriber access provided Published by American UNIVERSITY Chemical Society. Copyright © American Chemical Society.

6-OH-BDE-47 and 6-MeO-BDE-47 Environmental in Science & Technology is published by the Americanof Chemical early life-stages Society. 1155 Sixteenth Street N.W., rerio) Washington, zebrafish (Danio

20036 by RUTGERS Subscriber access DC provided Published by American UNIVERSITY Chemical Society. Copyright © American Chemical Society.

Hongling Liu, Song Tang, Xinmei Zheng, Yuting Zhu, Zhiyuan Environmental Science & Ma, Chunsheng Liu, Marcus Technology is published by the American Chemical Hecker, David Michael Vincent Society. 1155 Sixteenth Saunders, John Giesy, StreetP. N.W., Washington, 20036 Subscriber access DC provided by RUTGERS Xiaowei Zhang, and Hongxia Yu Published by American UNIVERSITY

Chemical Society. Copyright © American Chemical Society.

Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es503833q • Publication Date (Web): 07 Jan 2015& Environmental Science Technology is published

Downloadedby from http:// Chemical the American pubs.acs.org on January 17, 2015 Society. 1155 Sixteenth Street N.W., Washington, 20036 by RUTGERS Subscriber access DC provided Published by American UNIVERSITY Chemical Society. Copyright © American Chemical Society.

Just Accepted

Environmental Science & “Just Accepted” manuscripts have been pee Technology is published online prior to technical editing, formatting for by the American Chemical Sixteenth as a fre Society provides Society. “Just 1155 Accepted” Street N.W., Washington, dissemination of scientific material as soon a DC 20036 by Subscriber access provided RUTGERS Published by American UNIVERSITY Chemical Society. Copyright © American Chemical Society.

appear in full in PDF format accompanied by a fully peer reviewed, but should not be conside Environmental ScienceObject & readers and citable by the Digital Iden Technology is published to authors. Therefore, the “Just Accepted” W by the American Chemical Society. 1155 Sixteenthis technica in the journal. After a manuscript Street N.W., Washington, Accepted” Web site and as an ASA 20036published Subscriber access DC provided by RUTGERS Published by American UNIVERSITY Chemical Society. Copyright © American Chemical Society.

changes to the manuscript text and/or graph and ethical guidelines that apply to the jou Science or consequences Environmental arising from the &use of info Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, 20036 by RUTGERS Subscriber access DC provided Published by American UNIVERSITY Chemical Society. Copyright © American Chemical Society.

PageEnvironmental 1 of 44 Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology

1

Bioaccumulation, biotransformation and toxicity of BDE-47, 6-OH-BDE-47 and

2

6-MeO-BDE-47 in early life-stages of zebrafish (Danio rerio)

3 4

Hongling Liu1#*, Song Tang2#, Xinmei Zheng1, Yuting Zhu1, Zhiyuan Ma1, Chunsheng Liu1,

5

Markus Hecker2,3, David M.V. Saunders3, John P. Giesy1,3,4,5, Xiaowei Zhang1, Hongxia Yu1*

6 7

1

8

Nanjing University, Nanjing, Jiangsu 210023, China

9

2

State Key Laboratory of Pollution Control and Resource Reuse, School of the Environment,

School of Environment and Sustainability, University of Saskatchewan, Saskatoon, SK S7N

10

5B3, Canada

11

3

Toxicology Centre, University of Saskatchewan, Saskatoon, SK S7N 5B3, Canada

12

4

Department of Veterinary Biomedical Sciences, University of Saskatchewan, Saskatoon, SK

13

S7N 5B3, Canada

14

5

15

Kong, SAR, China

16

#

These authors contributed equally to this work.

18

*

Correspondence to: Drs. Hongling Liu and Hongxia Yu, School of the Environment,

19

Nanjing University, Nanjing, Jiangsu 210023, China

20

Tel: +86-25-89680356; Fax: +86-25-89680356; Email: [email protected] (Dr. Hongling Liu)

21

and [email protected] (Dr. Hongxia Yu)

Department of Biology and Chemistry, City University of Hong Kong, Kowloon, Hong

17

22

1

ACS Paragon Plus Environment

Page 2 of 44

Page 3 of 44

Environmental Science & Technology

23

Key words: PBDEs, endocrine disruption, nuclear receptor, aryl hydrocarbon receptor,

24

estrogen receptor, androgen receptor, thyroid hormone receptor

25

2

ACS Paragon Plus Environment

Environmental Science & Technology

26

Abstract

27

2,2',4,4'-Tetrabromodiphenyl ether (BDE-47), 6-hydroxy-tetrabromodiphenyl ether

28

(6-OH-BDE-47), and 6-methoxy-tetrabromodiphenyl ether (6-MeO-BDE-47) are the most

29

detected congeners of polybrominated diphenyl ethers (PBDEs), OH-BDEs, and MeO-BDEs,

30

respectively, in aquatic organisms. Although it has been demonstrated that BDE-47 can

31

interfere with certain endocrine functions that are mediated through several nuclear hormone

32

receptors (NRs), most of these findings were from mammalian cell lines exposed in vitro. In

33

the present study, embryos and larvae of zebrafish were exposed to BDE-47, 6-OH-BDE-47,

34

and 6-MeO-BDE-47 to compare their accumulation, biotransformation, and bioconcentration

35

factors (BCF) from 4 to 120 hpf. In addition, effects on expression of genes associated with

36

eight different pathways regulated by NRs were investigated at 120 hpf. 6-MeO-BDE-47 was

37

most bioaccumulated and 6-OH-BDE-47, which was the most potent BDE, was least

38

bioaccumulated. Moreover, the amount of 6-MeO-BDE-47, but not BDE-47, transformed to

39

6-OH-BDE-47 increased in a time-dependent manner, approximately 0.01%, 0.04%, and 0.08%

40

at 48, 96, and 120 hpf, respectively. Expression of genes regulated by the aryl hydrocarbon

41

receptor (AhR), estrogen receptor (ER), and mineralocorticoid receptor (MR) was affected in

42

larvae exposed to 6-OH-BDE-47, while genes regulated by AhR, ER, and the glucocorticoid

43

receptor (GR) were altered in larvae exposed to BDE-47. The greatest effect on expression of

44

genes was observed in larvae exposed to 6-MeO-BDE-47. Specifically, 6-MeO-BDE-47

45

affected the expression of genes regulated by AhR, ER, AR, GR, and thyroid hormone

46

receptor alpha (TRα). These pathways were mostly down-regulated at 2.5 µM. Taken together,

47

these results demonstrate the importance of usage of an internal dose to assess the toxic

48

effects of PBDEs. BDE-47 and its analogs elicited distinct effects on expression of genes of

49

different hormone receptor-mediated pathways, which have expanded the knowledge of

50

different mechanisms of endocrine disrupting effects in aquatic vertebrates. Because some of 3

ACS Paragon Plus Environment

Page 4 of 44

Page 5 of 44

Environmental Science & Technology

51

these homologues are natural products assessments of risks of anthropogenic PBDE need to

52

be made against the background of concentrations from naturally occurring products. Even

53

though PBDEs are being phased out as flame retardants, the natural products remain.

4

ACS Paragon Plus Environment

Environmental Science & Technology

54

Introduction

55

Polybrominated diphenyl ethers (PBDEs) have been extensively employed as flame

56

retardants (FRs) in various consumer and commercial products for decades.1, 2 As a result of

57

the substantial production, long-term use, disposal, and recycling processes, these chemicals

58

are now frequently found in the environment.3 The persistence, bioaccumulation potential

59

and toxic potency (PBT criteria) of 2,2',4,4' -tetrabromodiphenyl ether (BDE-47), one of the

60

primary PBDEs found in the environment, has raised concern about its potential adverse

61

effects to ecosystems and human health.3-5 In addition to the synthetic BDE-47, its

62

hydroxylated (OH-) or methoxylated (MeO-) forms, 6-OH-BDE-47 and 6-MeO-BDE-47,

63

have been suggested to be natural products of marine organisms6 and have been detected in a

64

wide variety of freshwater and marine organisms including mollusks, mussels, shellfish, clam,

65

fish, seal, dolphin, and whale.4, 7-13 Moreover, it has been conclusively demonstrated that

66

MeO-BDEs, and not PBDEs, are precursors of OH-BDEs.6, 17-19 In zebrafish (Danio rerio),

67

6-MeO-BDE-47 can be transformed into 6-OH-BDE-47; however, BDE-47 cannot be

68

transformed into 6-OH-BDE-47.20 In addition, inter-conversion between 6-MeO-BDE-47 and

69

6-OH-BDE-47 has been observed during dietary exposure of Japanese medaka (Oryzias

70

latipes).17

71

To date, increasing evidence has shown that exposure to BDE-47 and its two natural

72

analogs, 6-OH-BDE-47 and 6-MeO-BDE-47, can elicit a number of adverse effects in

73

aquatic organisms including disruption of the endocrine system,21, 22 disruption of molting,23

74

developmental defects,20, 24-26 and neurobehavioral toxicity.27-29 OH- and MeO-BDEs have

75

been shown to exhibit greater toxic potencies than PBDEs for certain endpoints such as 5

ACS Paragon Plus Environment

Page 6 of 44

Page 7 of 44

Environmental Science & Technology

76

estrogenicity and androgenicity.21, 30 However, mechanisms of their toxicity are complex and

77

have not been fully resolved.31 PBDEs as well as OH- and MeO-BDEs are structural analogs

78

to thyroid hormones, T3 and T4, as well as dioxin-like chemicals such as polychlorinated

79

biphenyls (PCBs), dioxins (TCDD), and furans (PCDF).31 This raises the question of whether

80

the adverse biological outcomes resulting from exposure to PBDEs and OH- or MeO-BDEs

81

are due to their ability to simulate thyroid hormones or whether they elicit effects similar to

82

those of the above dioxin-like chemicals.

83

Nuclear receptors (NRs) are a superfamily of ligand-activated, transcription factors that

84

act globally to regulate a broad range of biological processes, including development,

85

reproduction, and metabolism.32, 33 NRs mediate signaling by ligands such as endogenous

86

hormones, lipids, and xenobiotics.34, 35 Upon binding of a ligand to the ligand binding domain

87

of several kinds of NRs, a complex array of cellular responses is initiated. Recently, several

88

in vivo and in vitro studies have investigated effects of BDE-47, OH-BDE-47 or

89

MeO-BDE-47 on certain NR mediated physiological pathways, in particular the pathways

90

involving the thyroid hormone receptor (TR), estrogen receptor (ER), androgen receptor

91

(AR), and aryl hydrocarbon receptor (AhR). For example, in adult fathead minnows

92

(Pimephales promelas), dietary exposure to BDE-47 induced transcription of TRα in the

93

brain of females, and decreased the transcription of TRβ in the brain of fish of both sexes.41

94

In porcine ovarian follicles, both BDE-47 and 6-OH-BDE-47 did not alter expression of AR

95

mRNA or associated protein, but decreased expression of ERβ mRNA and protein following

96

exposure to BDE-47 and increase both ERα and ERβ gene and protein expression following

97

exposure to 6-OH-BDE-47.42 In an AhR-responsive luciferase reporter assay, 6-OH-BDE-47 6

ACS Paragon Plus Environment

Environmental Science & Technology

98

exhibited greater potency to induce AhR activity than that of 6-MeO-PBDEs and BDE-47.43

99

However, so far, most NR studies of BDE-47, 6-OH-BDE-47, and 6-MeO-BDE-47 were

100 101

completed by use of mammalian or cellular assays. The zebrafish represents an excellent vertebrate model organism in environmental

102

toxicology studies,44-46 especially in context with investigating effects of endocrine disrupting

103

chemicals (EDCs) on reproductive and developmental systems.47-49 Moreover, developmental

104

profiling of zebrafish gene expression patterns has confirmed a high degree of conservation

105

in NR expression patterns between zebrafish and other vertebrate models.50 Therefore, in the

106

present study, zebrafish embryos and larvae were used to determine the time-course of

107

accumulation, biotransformation, and bioconcentration factors (BCFs) of BDE-47,

108

6-OH-BDE-47, and 6-MeO-BDE-47. Additionally, in order to gain a more comprehensive

109

understanding of the molecular mechanisms of the toxicity of BDE-47 and related OH- and

110

MeO-analogs on the endocrine system, their effects on expression of genes associated with

111

eight nuclear hormone receptor pathways, particularly ER, AR, AhR, TRα, peroxisome

112

proliferator-activated receptor alpha (PPARα), glucocorticoid receptor (GR),

113

mineralocorticoid receptor (MR), and pregnane x receptor (PxR), were investigated and

114

compared.

7

ACS Paragon Plus Environment

Page 8 of 44

Page 9 of 44

Environmental Science & Technology

115

Materials and methods

116

Materials and reagents

117

BDE-47 (98% purity) was purchased from Chem Service (West Chester, PA, USA).

118

6-MeO-BDE-47 and 6-OH-BDE-47 were synthesized at City University of Hong Kong, and

119

purities were more than 98% as described previously.51 13C-PCB-178, 13C-2’-OH-BDE-99

120

and 13C-BDE-139 were purchased from Cambridge Isotope Laboratories (Andover, MA,

121

USA). BDE-47, 6-OH-BDE-47, and 6-MeO-BDE-47 were dissolved in dimethyl sulfoxide

122

(DMSO, Generay Biotech, Shanghai, China) to prepare stock solutions and then diluted with

123

embryonic rearing water (60 mg/L instant ocean salt in aerated distilled water) to the desired

124

test concentrations. Concentration of DMSO in final test solutions did not exceed 0.1%.

125

RNAlater, RNA Stabilization Reagents, and RNeasy® Mini Kit were purchased from

126

QIAGEN (Hilden, Germany). Maxima® First Strand cDNA Synthesis Kits were purchased

127

from Fermentas (St Leon-Rot, Germany). SYBR® Real time PCR Master Mix Plus Kits were

128

purchased from Toyobo (Tokyo, Japan).

129

130 131

Animals and exposure experiment Adult (7-months old) AB strain zebrafish maintenance and culturing were performed as

132

previously described.20 The eggs were examined under a stereomicroscope and only normally

133

developed embryos were used for exposure experiments. Briefly, twenty embryos were

134

randomly distributed into a 25 mL glass beaker containing 20 mL of exposure solution. Fish

135

were exposed until 120 hour post-fertilization (hpf), by which time they had developed into

136

free-swimming larvae and most organs had completed development.52 The control group

137

received 0.1% DMSO (v/v) only. 100% of the exposure solutions were replaced by fresh

138

exposure solution every 48 h. For 6-OH-BDE-47, BDE-47, and 6-MeO-BDE-47 exposures, 8

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 44

139

the experiments included two parts: First, zebrafish embryos were exposed to 6-OH-BDE-47

140

(0, 0.008, 0.02, 0.05, 0.1, 0.5 µM), 6-MeO-BDE-47 (0, 0.02, 0.1, 0.5, 2.5 µM) and BDE-47

141

(0, 0.02, 0.1, 0.5, 2.5 µM) from 4 to 120 hpf to study the morphologic toxicity of compounds

142

as previously described.20 Secondly, based on the results of acute toxicity test, three

143

comparable exposure concentrations for each compound were chosen: 6-OH-BDE-47 at

144

0.008, 0.02, and 0.05 µM, and at 0.1, 0.5, and 2.5 µM for both 6-MeO-BDE-47 and BDE-47

145

from 4 to 120 hpf to study the effects on expression of 63 genes involved in eight

146

receptor-mediated pathways by q-RT-PCR. After exposure for 120 hpf larvae were

147

anesthetized with ethyl 3-aminobenzoate methanesulfonate (MS-222, Suzhou Xin Yong

148

Biological Medicine Technology Co Lt, Jiangsu, China), and were preserved in RNAlater

149

RNA Stabilization Reagents until total RNA isolation.

150 151

Bioavailability analysis and QA/QC

152

This experiment was designed to analyze bioaccumulation of the three chemicals in early

153

life-stages of zebrafish. In each treatment, 600 zebrafish embryos were exposed to

154

6-OH-BDE-47, 6-MeO-BDE-47 or BDE-47 at 300 µg/L (0.6, 0.58, and 0.62 µM,

155

respectively) from 4 to 120 hpf in glass beakers. Variation among concentrations of the three

156

compounds in both the exposure medium and the embryos/larvae was determined. 80

157

embryos/larvae and corresponding exposure solutions were collected at 12, 24, 48, 72, 96,

158

and 120 hpf. Detailed protocols for extraction, clean up, and quantification, and quality

159

assurance and quality control (QA/QC) are provided in previous studies20,

160

Supporting Information Methods and Table S1

161 9

ACS Paragon Plus Environment

53, 54

and

Page 11 of 44

162 163

Environmental Science & Technology

Quantitative RT-PCR Total RNA was isolated from zebrafish larvae using RNeasy® Mini Kit. The

164

concentration and quality of total RNA were determined in accordance with the procedures

165

described in a previous study.45 First-strand cDNA synthesis and quantitative RT-PCR were

166

performed using Maxima® First Strand cDNA Synthesis and SYBR® Realtime PCR Master

167

Mix Plus Kits.20 Quantitative RT-PCR was performed by an Applied Biosystems Stepone

168

Plus Real-time PCR System (Foster City, California, USA). The primers were either mined

169

from previous literature55 or designed using Primer 3 (http://bioinfo.ut.ee/primer3-0.4.0/).

170

Primer sequences were listed in Supporting Information Table S2. The housekeeping gene18S

171

small subunit ribosomal RNA (18S rRNA) was used as an internal control.56 The thermal

172

cycle was set at 95 oC for 2 min, followed by 40 cycles of 95 oC for 15 s and 60 oC for 1 min.

173

Melting curves were derived during RT-PCR to validate that all cDNA samples amplified

174

only a single product. Levels of expression of genes were normalized to 18S rRNA mRNA

175

contents using the 2-∆∆Ct method. Each concentration was measured in triplicate or

176

quadruplicate in a composite sample containing 20 larvae.

177 178 179

Nuclear receptor pathway analysis For genes relating to AhR and ER pathways, the Agilent Literature Search application

180

was used to construct a biological interaction network within the Cytoscape software v3.1.1

181

(Cytoscape consortium, San Diego, CA, USA).57-59 The gene networks of the other six NR

182

pathways were retrieved by either WikiPathways (http://www.wikipathways.org)60 or

183

SABioscience Gene Network Central 10

ACS Paragon Plus Environment

Environmental Science & Technology

184

(http://www.sabiosciences.com/genenetwork/genenetworkcentral.php), and integrated with

185

AhR and ER pathways as “associations” and visualized as one network by Cytoscape. Only

186

genes of interest were shown in this pathway network. The resulting network genes (nodes)

187

were colored by the Enhanced Graphics application within Cytoscape according to the

188

significant fold changes of gene expressions in the respective treatments.

189 190 191

Statistical analysis SPSS 12.0 (SPSS Inc., Chicago, IL, USA) was used for statistical analysis. A

192

Kolmogorov-Smirnov test was used to verify the normality of the data, and the homogeneity

193

of variances was analyzed by Levene’s test as previously described.20 If the data failed the

194

Kolmogorov-Smirnov test, logarithmic transformation was performed and data was checked

195

again for homogeneity of variances. A one-way analysis of variance (ANOVA) followed by

196

LSD test was used to evaluate differences between the control and exposure groups. A value

197

of p < 0.05 was considered statistically significant. To capture the likely nonlinearity in

198

concentrations in exposure water or in zebrafish embryos-larvae across different time-points,

199

generalized additive models (GAMs) were used by the ‘mgcv’ package in R software version

200

3.10 (R Core Team, Vienna, Austria). Hierarchical cluster analysis for the gene expression

201

was performed by the "complete" method in R. A heatmap of gene expression results was

202

implemented by ‘pheatmap’ package version 0.7.7 in R.

11

ACS Paragon Plus Environment

Page 12 of 44

Page 13 of 44

Environmental Science & Technology

203

Results

204

Morphologic effects of 6-OH-BDE-47, 6-MeO-BDE-47 and BDE-47

205

Among the three compounds, 6-OH-BDE-47 was the most potent to zebrafish

206

embryos/larvae (Supporting Information Figures S1. B, C, D, and F and Table S3). Exposure

207

to all concentrations caused delayed development of embryos for up to 6-8 h at 24 hpf. In

208

embryos exposed to 0.5 µM 6-OH-BDE-47, mortality significantly increased to 22.5 ± 4.08%,

209

at 48 hpf. At 48 hpf, in groups exposed to concentrations greater than 0.1 µM, the embryos

210

developed hypopigmentation. At 72 hpf, development was arrested in all embryos exposed to

211

0.5 µM 6-OH-BDE-47 (Figures S1. B and C) while development of embryos at 12-18 hpf

212

was not altered. Larvae exposed to 0.1 µM 6-OH-BDE-47 exhibited spinal curvatures (Figure

213

S1. D), decreased heartbeats, and reduced body lengths (3516 ± 250 µm in 0.1µM group and

214

4040 ± 55 µm in control). The LC50 values of 6-OH-BDE-47 for teratogenic effects were 0.28

215

µM (0.21-0.38) at 72 hpf, 0.13 µM (0.11-0.16) at 96 hpf, and 0.09 µM (0.04-0.10) at 120 hpf.

216

The most sensitive toxicological endpoint was spinal curvature (Figure S1. F), which was

217

manifested in a concentration-dependent manner with maximal effects resulting from

218

exposure to 0.08 µM (0.07-0.09) 6-OH-BDE-47 at 120 hpf.

219

There were no significant differences in developmental alterations in individuals exposed

220

to 6-MeO-BDE-47 and the control group until 96 hpf. The most sensitive toxicological

221

endpoints were concomitant spinal curvature and pericardial edema at 120 hpf at 2.5 µM

222

(Figure S1. G). Although no statistically significant differences were observed within 96 hpf

223

following exposure to BDE-47 up to 2.5 µM, concomitant spinal curvature and pericardial

224

edema occurred in embryos at 120 hpf. The NOEC of spinal curvature and pericardial edema

225

was 0.5 µM of BDE-47, although there were significant differences at 2.5 µM and the

226

proportion of affected larvae was 27.6 ± 18.3% (Figure S1 H). In addition, at 120 hpf

227

following exposure to 2.5 µM BDE-47the body lengths of larvae (3751 ± 152 µm) were 12

ACS Paragon Plus Environment

Environmental Science & Technology

228

significantly reduced compared to those in the control group (4029 ± 201µm).

229

230 231

Accumulation by zebrafish Concentrations of BDE-47, 6-OH-BDE-47, or 6-MeO-BDE-47 were below their

232

analytical method detection limits in the control group (Supporting Information Table S4). At

233

120 hpf, 100% mortality occurred following exposure to 300 µg/L 6-OH-BDE-47, therefore,

234

no data is available for this time point. GC/MS results indicated that the chemical

235

concentrations in exposure solutions decreased and the doses in embryos and larvae increased

236

in a time-dependent manner (Figure 1). Moreover, the concentrations of BDE-47 and

237

6-MeO-BDE-47 were approximately 10- to 100-fold greater than concentrations of

238

6-OH-BDE-47 in zebrafish embryos and larvae that were exposed to the same concentrations

239

of the three compounds (Figure 1). The three BDE congeners ranked as follows regarding

240

their in vivo accumulation (values from greater to lesser potential): 6-MeO-BDE-47 >

241

BDE-47 > 6-OH-BDE-47.

242

Calculated BCF values for 6-OH-BDE-47 were 4.07, 9.88, 21.7, 26.9, and 23.3 at 12 hpf,

243

24 hpf, 48 hpf, 72 hpf, and 96 hpf, respectively (Figure 2A). For 6-MeO-BDE-47, the BCF

244

values were 17.2, 42.6, 89.7, 390, 1207, and 935 at 12 hpf, 24 hpf, 48 hpf, 72 hpf, 96 hpf, and

245

120 hpf, respectively. For BDE-47, bioconcentration factor (BCF) values were 25.7, 68.0,

246

489, 750, 2489, and 2430 at 12 hpf, 24 hpf, 48 hpf, 72 hpf, 96 hpf, 120 hpf, respectively.

247

These results indicated that trends of BCF values after various durations of exposure were

248

similar for 6-MeO-BDE-47 and BDE-47, reaching a maximum at 96 hpf, and showing a

249

slight decrease at later time points.

13

ACS Paragon Plus Environment

Page 14 of 44

Page 15 of 44

Environmental Science & Technology

250 251 252

Biotransformation in zebrafish Exposure of zebrafish to 300 µg/L 6-MeO-BDE-47 resulted in increasing tissue

253

concentrations of 6-OH-BDE-47 in their tissues were quantified in a time-dependent manner,

254

with 0.02, 0.11, 0.12, and 0.25 µg/g, wm (wet mass) at 48, 72, 96, and 120 hpf, respectively

255

(Figure 2B). However, no transformation of BDE-47 into 6-OH-BDE-47 or 6-MeO-BDE

256

occurred. Also, 6-OH-BDE-47 was not transformed into either BDE-47 or 6-MeO- BDE-47.

257

Based on these transformation ratios, the amounts of biotransformed 6-OH-BDE-47 were

258

expected to be 0.04, 0.21, and 1.07 µg/g, wm after 120 hpf exposure to 0.1, 0.5, or 2.5 µM

259

6-MeO-BDE-47, respectively (Table 1). Amounts of biotransformed 6-OH-BDE-47 in larvae

260

exposed to 0.5 and 2.5 µM of 6-MeO-BDE-47 were greater than the accumulated

261

concentration of 6-OH-BDE-47 (0.21 and 1.07 µg/g, wm vs. 0.15 µg/g, wm) in larvae

262

exposed to 0.05 µM of 6-OH-BDE-47 (Table 1).

263 264

265

Transcriptional responses of NR pathways to 6-OH-BDE-47 A hierarchical cluster analysis of gene expression data showed a dendrogram which

266

highlighted five principal clusters (Figure 3A). Exposures to 2.5 µM 6-MeO-BDE-47 resulted

267

in a unique clustering of gene expression data that revealed a significantly different gene

268

expression profile from the other exposures (Figure 3A). Exposures to the same compound

269

but at different concentrations generally clustered in the same group, especially for exposures

270

to 0.008, 0.1, and 0.5 µM 6-OH-BDE-47 (Figure 3A).

271

Zebrafish embryos exposed to 6-OH-BDE-47 had significant alterations in the

272

expression of genes associated with several NR pathways. The most significant effects

273

occurred along the AhR pathway (Figures 3 and 4 and Supporting Information Table S5) with 14

ACS Paragon Plus Environment

Environmental Science & Technology

274

exposure to 0.008 µM 6-OH-BDE-47 causing a significant up-regulation in the expression of

275

ahr1a, ahr1b, ahr2, and arnt2 by 2.32-, 1.71-, 1.62-, and 1.62-fold, respectively (p < 0.05),

276

and a significant down-regulation of ahrra and cyp19b expression by 2.27- and 1.81-fold,

277

respectively (p < 0.05). Exposure to 0.02 µM 6-OH-BDE-47 significantly induced expression

278

of ahr1b and ahr2 by 1.63- and 1.67-fold, respectively, and reduced the expression of ahrra

279

and cyp1a1 by 1.96- and 1.75-fold (p < 0.05). Exposure to 0.05 µM 6-OH-BDE-47

280

significantly reduced the expression of cyp1a1, cyp1b1, arnt1la, ahrra, and cyp19b by 1.92-,

281

2.63-, 1.25-, 1.96-, and 2.23-fold, respectively (p < 0.05). In addition to the ahr receptor,

282

following exposure to 0.008 and 0.05 µM 6-OH-BDE-47 the expression of mr and er2b were

283

significantly up-regulated by 1.51- and 1.83-fold, respectively (p < 0.05).

284 285

Transcriptional responses of NR pathways to 6-MeO-BDE-47

286

Some NR-mediated pathways such as AhR, ER, AR, TR, and GR were affected by

287

exposure to 6-MeO-BDE-47 with the greatest effects occurring at the greatest concentration

288

tested, 2.5 µM (Figures 3 and 4). Exposure to lesser concentrations of 6-MeO-BDE-47, 0.1

289

µM, also induced the expression of several genes in these pathways. Specifically, arnt2, dut,

290

ugtlal, ctnnb1, pa2g4a, dap3, and rela were significantly induced by 2.08-, 1.67-, 2.23-, 1.56-,

291

1.78-, and 1.51-fold, respectively (p < 0.05). Exposure to 0.5 µM 6-MeO-BDE-47 did not

292

significantly alter the expressions of most genes, except for the down-regulation cyp1a1 and

293

er2a by 9.09- and 1.56-fold, respectively (p < 0.05). However, exposure to a greater

294

concentration of 6-MeO-BDE-47, 2.5 µM, reduced the expression of most altered genes.

295

Along the AhR pathway, ahr2 was down-regulated by 1.92-fold and associated genes such as

296

cyp1a1, cyp1b1, cyp365a, and sp1 were also down-regulated by 50-, 100-, 3.03-, and

297

1.89-fold, respectively (p < 0.05). However, both ahrra and ahrrb were significantly induced 15

ACS Paragon Plus Environment

Page 16 of 44

Page 17 of 44

Environmental Science & Technology

298

by 3.16- and 5.71-fold, respectively (p < 0.05). In the ER pathway, er2a and ccnd1 were

299

down-regulated by 1.43- and 1.85-fold, respectively, while er2b was up-regulated by

300

1.44-fold following exposure to 2.5 µM 6-MeO-BDE-47 (p < 0.05). In the AR pathway, ar,

301

ctnnb1, pa2g4b, and ncoa1 were down-regulated by 1.85-, 1.89-, 1.47-, and 1.67-fold,

302

respectively (p < 0.05). Exposure to 2.5 µM 6-MeO-BDE-47 also decreased the expression of

303

thra by 2.17-fold and TR associated genes such as ncor and fus were significantly

304

down-regulated by 2.08- and 1.85-fold (p < 0.05). Also, following exposure to 2.5 µM

305

6-MeO-BDE-47, the expression of gr and tgfb1 were down-regulated by 2.04- and 1.72-fold,

306

respectively (p < 0.05).

307 308 309

Transcriptional responses of NR pathways to BDE-47 Exposure to BDE-47 significantly affected the expression of receptors in the AhR, ER,

310

and GR pathways (Figures 3 and 4). The number of altered genes increased in a concentration

311

dependent manner. The expression of gr was significantly down-regulated by 1.28-fold

312

following exposure to 0.5 µM BDE-47 (p < 0.05). When exposed to 2.5 µM, expressions of

313

ahr1b, er2a, and er2b were significantly greater by 1.64-, 1.35-, and 1.40-fold, respectively

314

(p < 0.05) than those of the respective genes controls. Expressions of genes along AhR and

315

ER pathways, such as cyp1a1, cyp19a, cyp3a65, and ccnd1, were significantly less by factors

316

of 5.88-, 2.88-, 2.04-, and 1.52-fold, respectively, relative to that of the controls (p < 0.05).

16

ACS Paragon Plus Environment

Environmental Science & Technology

317 318

Discussion Hazard assessment of contaminants is typically based on the exposure of aquatic

319

organisms to chemical solutions for a defined exposure time and the adverse outcomes

320

observed are then correlated with the concentrations of the compounds in the ambient media.

321

Since chemicals need to be accumulated into organisms and distributed to target sites for the

322

induction of toxic effects, usage of target site effect concentrations are postulated to best

323

represent the hazards of a compound in vivo.61, 62 However, for small-bodied organisms such

324

as the zebrafish effect concentrations at the target site are difficult to determine, particularly

325

at earlier stages of development. Average body concentrations of contaminants in zebrafish

326

embryos and larvae are subject to competitive dynamic processes, which include the ability

327

of compounds to penetrate the chorion, in vivo biotransformation, distribution, and excretion.

328

Hence, internal effect concentrations need to be determined. Moreover, the ratio of the

329

internal concentration in a fish to the surrounding concentration at a steady state represents

330

the compound’s BCF, which is an important metric for regulatory assessment of chemicals.63

331

In the present study accumulation, biotransformation, varied among BCF of BDE-47,

332

6-OH-BDE-47, and 6-MeO-BDE-47 during multiple developmental stages of zebrafish. The

333

time points during which accumulation of BDE-47 and 6-MeO-BDE-47 increased

334

substantially coincide with the hatching period of zebrafish embryos (48 hpf). In the late

335

developmental periods of the larvae (after 96 hpf), the bioaccumulation of BDE-47 and

336

6-MeO-BDE-47 reached a plateau, which might be due to an increase in metabolism and/or

337

excretion activities as well as the rapid growth of larvae at this time that might dilute their

338

body concentrations over this time period. In this study, the BCF of 6-OH-BDE-47 at 96 hpf 17

ACS Paragon Plus Environment

Page 18 of 44

Page 19 of 44

Environmental Science & Technology

339

was 23.3, which is similar to previously reported results.64 in which BCF values were

340

calculated in liver of zebrafish after 96 h exposure to 100 nM 6-OH-BDE-47. At all six

341

durations of exposure, 6-OH-BDE-47 was the least accumulated into the body, though it had

342

the greatest toxic potency. It is known that a compound with greater logKow generally has

343

greater bioaccumulation potential. Values of logKow were 6.76, 7.17, and 6.59 for BDE-47,

344

6-MeO-BDE-47, and 6-OH-BDE-47, respectively.65 Thus, differences in lipophilicity were

345

considered to be the most important parameters for the different accumulation properties of

346

the test chemicals. For compounds such as 6-OH-BDE-47, the hydroxyl group by making the

347

compound more polar, result in greater excretion, and might also be an important parameter

348

for the lesser in vivo concentration. Hence, our results confirmed that in aquatic exposure

349

tests, it is not sufficient to evaluate the ecotoxicological risk of a compound based solely on

350

the exposure concentration. In addition to accumulation in the body, the results of this study

351

indicated that the amounts of 6-MeO-BDE-47, but not BDE-47, that were transformed to

352

6-OH-BDE-47 increased in a time-dependent manner, (approximately 0.01%, 0.04%, and

353

0.08% at 48, 96, and 120 hpf, respectively), which is indicative of an increasing metabolic

354

capability of zebrafish embryos/larvae with increasing age.

355

Photomicrographs demonstrated that exposures to 6-OH-BDE-47, 6-MeO-BDE-47, and

356

BDE-47 resulted in developmental abnormalities in zebrafish embryos and larvae. The

357

embryo-toxic effects of BDE47, 6-OH-BDE47 and 6-MeO-BDE47 have been investigated in

358

a previous study with zebrafish embryos exposed from 3 to 72 hpf.64 Those authors showed

359

that 6-OH-BDE-47 was the most toxic BDE-47 congener inducing a range of developmental

360

defects including pericardial edema, yolk sac deformations, lesser pigmentation, lessened 18

ACS Paragon Plus Environment

Environmental Science & Technology

361

heart rate, and delayed development at concentrations of 25-50 nM, which is consistent with

362

the morphology findings of this study (Figures S1 B, C and D).64 Furthermore, the lack of

363

toxic effects of BDE47 or 6-MeO-BDE47 at 2.5 µM until 120 hpf (Figures S1 G and H) is

364

consistent with previous findings that showed that no toxicity were observed for both BDE47

365

and 6-MeO-BDE47 in zebarfish at 72 hpf.64

366

Since disruptions of cellular molecules or processes are thought to precede adverse

367

outcomes, changes to normal molecular processes might function as sensitive biomarkers to

368

predict adverse biological outcomes.66 Moreover, alteration of NR mediated pathways has

369

been shown to be associated with adverse endocrine and developmental effects that were

370

linked with morphological deformities.20, 55 For example, previous studies have demonstrated

371

that BDE-47 can alter thyroid status and thyroid hormone-regulated gene transcription in the

372

pituitary and brain of adult fathead minnows,67 and both 6-OH-BDE-47 and 6-MeO-BDE-47

373

were shown to affect expression of TRα and TRβ genes in the TR pathway, which can result

374

in teratogenic effects such as pericardial edema, developmental retardation, and curved spine

375

in zebrafish embryos.20 Also, all BDE-47, TBBPA and BPA have been demonstrated to alter

376

expression of genes along the hypothalamus-pituitary-thyroid (HPT) axis of zebrafish larvae

377

as well as induce acute toxicity.68 Additionally, zebrafish has been used to investigate effects

378

of EDCs on the expression of genes in six NR mediated pathways.55 In this study, two

379

additional receptor pathways-AR and PxR were added, and the interactions of sixty-three

380

genes involved in eight zebrafish receptor pathways were integrated. This pathway network

381

might represent a novel tool for the examination of the molecular function of each individual

382

receptor as well as for the study of their combinatorial regulatory network within NRs. 19

ACS Paragon Plus Environment

Page 20 of 44

Page 21 of 44

383

Environmental Science & Technology

The structures of the three compounds tested are similar. The cluster dendrogram

384

showed that expression of genes in individuals exposed to various concentrations of

385

6-OH-BDE-47 clustered together. Nevertheless, clustering is a function of concentration for

386

6-MeO-BDE-47 and BDE-47. Patterns of expression of genes following exposures to 2.5 µM

387

BDE-47 and 0.5 µM 6-MeO-BDE-47 were grouped into the same cluster, which indicates

388

BDE-47 likely has fewer effects on zebrafish NR-mediated pathways than 6-MeO-BDE-47 at

389

similar waterborne exposure concentrations. Exposure to 2.5 µM of the more

390

bioaccumulative 6-MeO-BDE-47 resulted in a unique gene expression profile compared to

391

BDE-47 and 6-OH-BDE-47. In addition, 1.07 µg/g, wm of biotransformed 6-OH-BDE-47

392

were detected in larvae exposed to 2.5 µM of 6-MeO-BDE-47, which was much greater than

393

the detected amount of 6-OH-BDE-47, 0.15 µg/g, which resulted from exposure to 0.05 µM

394

6-OH-BDE-47. The greater body burden of 6-OH-BDE-47 resulting from exposure to 2.5 µM

395

of 6-MeO-BDE-47 might explain the significant and great effect on gene transcription of NR

396

pathways observed in this exposure group. The effects that occurred in the greatest exposure

397

group of 6-MeO-BDE-47, therefore, were attributed to the combined effects of

398

biotransformed 6-OH-BDE-47 as well as 6-MeO-BDE-47. The clustering of the three

399

compounds correlated well with their respective accumulation potency, indicating the great

400

importance of the usage of internal dose to assess the dose-response relationship for studies

401

of PBDEs, especially MeO-PBDEs.

402

Further analyses of endocrine pathways indicated general disruption of receptor

403

pathways by all three BDEs congeners, which correlated well with the observed teratogenic

404

effects in zebrafish. Specifically, 6-OH-BDE-47 altered the expression of AhR, ER, and MR 20

ACS Paragon Plus Environment

Environmental Science & Technology

405

receptor-mediated pathways, while AhR, ER, and GR were the primary pathways altered by

406

BDE-47. Yet, exposure to the more bioaccumulative 6-MeO-BDE-47 affected AhR, ER, AR,

407

GR, and TRα pathways. Molecular structures of OH-BDEs closely resemble those of thyroid

408

hormones (THs) and 6-OH-BDE-47 can disrupt normal thyroid homeostasis and functions as

409

either an agonist or an antagonist.31, 69 For example, expression of genes along the

410

hypothalamus-pituitary-thyroid (HPT) axis that is responsible for regulation of metabolism

411

and early life-stage development was affected by BDE-47 and its OH- or MeO- forms in

412

zebrafish embryos-larvae.20, 68 However, in the present experiment, expression of thra was

413

not significantly altered following exposure to 6-OH-BDE-47, which contradicted previous

414

findings showing thra was reduced in zebrafish exposed to 200 µg/L 6-OH-BDE-47.20 This

415

difference might be due to the almost 10 times greater concentrations used in the previous

416

study. However, exposure to 2.5 µM 6-MeO-BDE-47 significantly decreased the expression

417

of thra, nocr and fus in TRα pathway. Since significant quantities of biotransformed

418

6-OH-BDE-47 were found in vivo, 6-MeO-BDE-47 may exert adverse effects indirectly via

419

transformation into 6-OH-BDE-47, which then can bind directly to TH targeted genes by

420

mimicking THs.

421

Apart from the TRα pathway, recent studies have also focused on disruption of the AhR

422

and ER pathways by PBDEs. Cross talk between ER- and AhR-signaling pathways in fish

423

has been hypothesized previously.72, 73 The pathway analyses conducted in this study also

424

suggest interactions between these pathways as visualized in the constructed networks. All

425

three compounds altered AhR and ER pathways in zebrafish. 6-OH-BDE-47 significantly

426

increased expression of er2b in zebrafish, which is consistent with previous in vitro findings 21

ACS Paragon Plus Environment

Page 22 of 44

Page 23 of 44

Environmental Science & Technology

427

that both ERα and ERβ gene and protein expression were induced by 6-OH-BDE-47.42

428

Exposure to 2.5 µM 6-MeO-BDE-47 significantly reduced expression of er2a but induced

429

er2b indicating the compound might cause endocrine disrupting effects through interfering

430

with the ER signaling pathway.21 In addition, 6-OH-BDE-47 increased the expression of

431

several AhR receptors including ahr1a, ahr1b, and ahr2 in vivo, while 6-MeO-BDE-47 and

432

BDE-47 only affected ahr2 and ahr1a transcription, respectively. These results have

433

confirmed that OH-BDEs can induce greater dioxin-like activity than corresponding

434

MeO-BDEs and parent PBDEs in vitro.43, 74 AR and PxR pathways have been previously shown to be affected by PBDEs in vitro.21,

435 436

75-77

437

exhibited potent antiandrogenicity, with potencies ranking as follows: 6-OH-BDE-47 (IC50 =

438

0.34 µM) > BDE47 (IC50 = 3.83 µM) > 6-MeO-BDE-47 (IC50 = 41.8 µM).21, 76 However, in

439

zebrafish, both BDE-47 and 6-OH-BDE-47 did not significantly alter AR expression, which

440

is consistent with a study in porcine ovarian follicular cells, showing BDE-47 and its OH-

441

metabolites had no effect on the expression of AR mRNA and protein expression.42 The PxR,

442

a steroid and xenobiotic nuclear receptor (SXR), can be activated by BDE-47 in mice.75

443

Nevertheless, in zebrafish, BDE-47 significantly down-regulated PxR associated genes of

444

cyp3a65, hnf4a and ugtlal, but increased pou1f1. The incongruities between these results

445

could be due to differences between species and lesser concentrations used in our

446

experiments. Furthermore, PxR associated genes cyp24a1 and hnf4a were significantly

447

up-regulated by 6-OH-BDE-47 exposure, indicating 6-OH-BDE-47 might be an agonist of

448

zebrafish PxR.

In the MDA-kb2 human cell line AR receptor binding assay, all three compounds

22

ACS Paragon Plus Environment

Environmental Science & Technology

449

Reports on the effect of PBDEs on PPARα, MR, and GR are limited. PPARα plays an

450

important role in lipid homeostasis, inflammation, adipogenesis, reproduction, and

451

carcinogenesis.78 In this study, none of three compounds significantly affected the expression

452

of PPARα in zebrafish. However, treatment with a PBDE mixture, BDE-71 and BDE-47

453

caused increases in PPARγ transcript levels at day 8th in 3T3-L1 mouse embryo fibroblast

454

cells.79 GR and MR are essential for regulation of multiple physiological functions, such as

455

glucose metabolism, mineral balance, and behavior.80 In this study, exposure to 2.5 µM

456

6-MeO-BDE-47 or 0.5 µM BDE-47 caused down-regulation of GR, while 0.008 µM

457

6-OH-BDE-47 increased MR expression, which suggests that GR or MR signaling pathways

458

might be involved in the endocrine disrupting effects induced by PBDEs.

459

Altogether, the results of present study, which compared the toxicities of BDE-47

460

with its OH- and MeO- analogs in zebrafish via multiple quantitative approaches, ranging

461

from in vivo toxicity tests, bioaccumulation and biotransformation, to the molecular analysis

462

of response patterns of genes along NR pathways, highlight the importance of the usage of

463

internal dose to evaluate the toxic effects for PBDEs, and the use of early life-stages of

464

zebrafish as an efficient and reliable vertebrate model to assess toxicological effects of

465

endocrine disruptors. Our data also elucidated several molecular aspects of BDE-47,

466

6-OH-BDE-47, and 6-MeO-BDE-47 induced toxicities. Specifically, the data provided

467

valuable insights into the early interaction of these compounds with steroid hormone receptor

468

pathways, which provided novel clues for their in vivo mechanisms of subsequent endocrine

469

disruption and developmental toxicities.

470 23

ACS Paragon Plus Environment

Page 24 of 44

Page 25 of 44

471 472

Environmental Science & Technology

Acknowledgements We thank Dr. Richard A. Erickson (Upper Midwest Environmental Sciences Center, U.S.

473

Geological Survey) for providing the expertise in statistical analyses. This work was funded

474

by National Natural Science Foundation (No. 21377053 and 20977047) and Major National

475

Science and Technology Projects (No. 2012ZX07506-001 and 2012ZX07501-003-02) of

476

China. J. P. G. and M. H. were supported by the Canada Research Chair Program. J. P. G. was

477

supported by the Program of 2012 "Great Level Foreign Experts" (#GDW20123200120)

478

funded by the State Administration of Foreign Experts Affairs, China to Nanjing University,

479

and the Einstein Professor Program of the Chinese Academy of Sciences. He was also

480

supported by a Visiting Distinguished Professorship in the Department of Biology and

481

Chemistry and State Key Laboratory in Marine Pollution at City University of Hong Kong.

482 483 484

Supporting Information Further details on the analytical methods and additional tables and figures as noted in the

485

text are available in Supporting Information. This information is available free of charge via

486

the Internet at http://pubs.acs.org/.

24

ACS Paragon Plus Environment

Environmental Science & Technology

487

References

488

1. Birnbaum, L. S.; Staskal, D. F., Brominated flame retardants: cause for concern? Environ.

489

Health Persp. 2004, 112, (1), 9-17.

490

2.

491

polybrominated diphenyl ethers among workers at an electronic waste dismantling region in

492

Guangdong, China. Environ. Int. 2007, 33, (8), 1029-34.

493

3.

494

Martin, P.; Wayland, M.; Chip Weseloh, D. V.; Wilson, L., Flame retardants in eggs of four

495

gull species (Laridae) from breeding sites spanning Atlantic to Pacific Canada. Environ.

496

Pollut. 2012, 168, 1-9.

497

4.

498

Polybrominated diphenyl ethers and hydroxylated and methoxylated brominated and

499

chlorinated analogues in the plasma of fish from the Detroit River. Environ. Sci. Technol.

500

2005, 39, (15), 5612-9.

501

5.

502

naturally produced brominated phenoxyphenols and phenoxyanisoles in aquatic plants from

503

the Philippines. J. Agric. Food Chem. 2010, 58, (23), 12385-91.

504

6.

505

Polybrominated diphenyl ethers and their hydroxylated/methoxylated analogs: Environmental

506

sources, metabolic relationships, and relative toxicities. Mar. Pollut. Bull. 2011, 63, (5-12),

507

179-88.

508

7.

Qu, W. Y.; Bi, X. H.; Sheng, G. Y.; Lu, S. Y.; Fu, H.; Yuan, J.; Li, L. P., Exposure to

Chen, D.; Letcher, R. J.; Burgess, N. M.; Champoux, L.; Elliott, J. E.; Hebert, C. E.;

Valters, K.; Li, H.; Alaee, M.; D'Sa, I.; Marsh, G.; Bergman, A.; Letcher, R. J.,

Haraguchi, K.; Kotaki, Y.; Relox, J. R., Jr.; Romero, M. L.; Terada, R., Monitoring of

Wiseman, S. B.; Wan, Y.; Chang, H.; Zhang, X. W.; Hecker, M.; Jones, P. D.; Giesy, J. P.,

Rotander, A.; van Bavel, B.; Riget, F.; Auethunsson, G. A.; Polder, A.; Gabrielsen, G. W.; 25

ACS Paragon Plus Environment

Page 26 of 44

Page 27 of 44

Environmental Science & Technology

509

Vikingsson, G.; Mikkelsen, B.; Dam, M., Methoxylated polybrominated diphenyl ethers

510

(MeO-PBDEs) are major contributors to the persistent organobromine load in sub-Arctic and

511

Arctic marine mammals, 1986-2009. Sci. Total Environ. 2012, 416, 482-9.

512

8.

513

marine environment: a review. Environ. Int. 2009, 35, (3), 655-66.

514

9.

515

methoxylated PBDEs in bivalves from Beijing markets. Chemosphere 2014, 110, 97-103.

516

10. Sun, J.; Liu, J.; Liu, Y.; Jiang, G., Hydroxylated and methoxylated polybrominated

517

diphenyl ethers in mollusks from Chinese coastal areas. Chemosphere 2013, 92, (3), 322-8.

518

11. Yin, G.; Asplund, L.; Qiu, Y.; Zhou, Y.; Wang, H.; Yao, Z.; Jiang, J.; Bergman, A.,

519

Chlorinated and brominated organic pollutants in shellfish from the Yellow Sea and East

520

China Sea. Environ. Sci. Pollut. Res. Int. 2014.

521

12. Wang, H. S.; Du, J.; Ho, K. L.; Leung, H. M.; Lam, M. H. W.; Giesy, J. P.; Wong, C. K.

522

C.; Wong, M. H., Exposure of Hong Kong residents to PBDEs and their structural analogues

523

through market fish consumption. J. Hazard. Mater. 2011, 192, (1), 374-80.

524

13. McKinney, M. A.; De Guise, S.; Martineau, D.; Beland, P.; Lebeuf, M.; Letcher, R. J.,

525

Organohalogen contaminants and metabolites in beluga whale (Delphinapterus leucas) liver

526

from two Canadian populations. Environ. Toxicol. Chem. 2006, 25, (5), 1246-57.

527

14. Malmvarn, A.; Marsh, G.; Kautsky, L.; Athanasiadou, M.; Bergman, A.; Asplund, L.,

528

Hydroxylated and methoxylated brominated diphenyl ethers in the red algae Ceramium

529

tenuicorne and blue mussels from the Baltic Sea. Environ. Sci. Technol. 2005, 39, (9),

530

2990-7.

Yogui, G. T.; Sericano, J. L., Polybrominated diphenyl ether flame retardants in the U.S.

Liu, X.; Jiao, Y.; Lin, C.; Sun, K.; Zhao, Y., PBDEs, hydroxylated PBDEs and

26

ACS Paragon Plus Environment

Environmental Science & Technology

531

15. Bowden, B. F.; Towerzey, L.; Junk, P. C., A new brominated diphenyl ether from the

532

marine sponge Dysidea herbacea. Aust. J. Chem. 2000, 53, (4), 299-301.

533

16. Malmvarn, A.; Zebuhr, Y.; Kautsky, L.; Bergman, K.; Asplund, L., Hydroxylated and

534

methoxylated polybrominated diphenyl ethers and polybrominated dibenzo-p-dioxins in red

535

alga and cyanobacteria living in the Baltic Sea. Chemosphere 2008, 72, (6), 910-6.

536

17. Wan, Y.; Liu, F.; Wiseman, S.; Zhang, X.; Chang, H.; Hecker, M.; Jones, P. D.; Lam, M.

537

H.; Giesy, J. P., Interconversion of hydroxylated and methoxylated polybrominated diphenyl

538

ethers in Japanese medaka. Environ. Sci. Technol. 2010, 44, (22), 8729-35.

539

18. Wan, Y.; Wiseman, S.; Chang, H.; Zhang, X.; Jones, P. D.; Hecker, M.; Kannan, K.;

540

Tanabe, S.; Hu, J.; Lam, M. H.; Giesy, J. P., Origin of hydroxylated brominated diphenyl

541

ethers: natural compounds or man-made flame retardants? Environ. Sci. Technol. 2009, 43,

542

(19), 7536-42.

543

19. Liu, F.; Wiseman, S.; Wan, Y.; Doering, J. A.; Hecker, M.; Lam, M. H.; Giesy, J. P.,

544

Multi-species comparison of the mechanism of biotransformation of MeO-BDEs to

545

OH-BDEs in fish. Aquat. Toxicol. 2012, 114-115, 182-8.

546

20. Zheng, X.; Zhu, Y.; Liu, C.; Liu, H.; Giesy, J. P.; Hecker, M.; Lam, M. H.; Yu, H.,

547

Accumulation and biotransformation of BDE-47 by zebrafish larvae and teratogenicity and

548

expression of genes along the hypothalamus-pituitary-thyroid axis. Environ. Sci. Technol.

549

2012, 46, (23), 12943-51.

550

21. Hu, W.; Liu, H.; Sun, H.; Shen, O.; Wang, X.; Lam, M. H.; Giesy, J. P.; Zhang, X.; Yu, H.,

551

Endocrine effects of methoxylated brominated diphenyl ethers in three in vitro models. Mar.

552

Pollut. Bull. 2011, 62, (11), 2356-61. 27

ACS Paragon Plus Environment

Page 28 of 44

Page 29 of 44

Environmental Science & Technology

553

22. Chan, W. K.; Chan, K. M., Disruption of the hypothalamic-pituitary-thyroid axis in

554

zebrafish embryo-larvae following waterborne exposure to BDE-47, TBBPA and BPA. Aquat.

555

Toxicol. 2012, 108, 106-11.

556

23. Davies, R.; Zou, E., Polybrominated diphenyl ethers disrupt molting in neonatal Daphnia

557

magna. Ecotoxicology 2012, 21, (5), 1371-80.

558

24. Dong, W.; Macaulay, L. J.; Kwok, K. W. H.; Hinton, D. E.; Stapleton, H. M., Using

559

whole mount in situ hybridization to examine thyroid hormone deiodinase expression in

560

embryonic and larval zebrafish: A tool for examining OH-BDE toxicity to early life stages.

561

Aquat. Toxicol. 2013, 132, 190-9.

562

25. Usenko, C. Y.; Hopkins, D. C.; Trumble, S. J.; Bruce, E. D., Hydroxylated PBDEs induce

563

developmental arrest in zebrafish. Toxicol. Appl. Pharmacol. 2012, 262, (1), 43-51.

564

26. Kallqvist, T.; Grung, M.; Tollefsen, K. E., Chronic toxicity of

565

2,4,2',4'-tetrabromodiphenyl ether on the marine alga Skeletonema costatum and the

566

crustacean Daphnia magna. Environ. Toxicol. Chem. 2006, 25, (6), 1657-62.

567

27. Zhao, J.; Xu, T.; Yin, D. Q., Locomotor activity changes on zebrafish larvae with

568

different 2,2',4,4'-tetrabromodiphenyl ether (PBDE-47) embryonic exposure modes.

569

Chemosphere 2014, 94, 53-61.

570

28. Chen, X.; Huang, C.; Wang, X.; Chen, J.; Bai, C.; Chen, Y.; Chen, X.; Dong, Q.; Yang,

571

D., BDE-47 disrupts axonal growth and motor behavior in developing zebrafish. Aquat.

572

Toxicol. 2012, 120-121, 35-44.

573

29. Chou, C. T.; Hsiao, Y. C.; Ko, F. C.; Cheng, J. O.; Cheng, Y. M.; Chen, T. H., Chronic

574

exposure of 2,2',4,4'-tetrabromodiphenyl ether (PBDE-47) alters locomotion behavior in 28

ACS Paragon Plus Environment

Environmental Science & Technology

575

juvenile zebrafish (Danio rerio). Aquat. Toxicol. 2010, 98, (4), 388-95.

576

30. Hamers, T.; Kamstra, J. H.; Sonneveld, E.; Murk, A. J.; Visser, T. J.; Van Velzen, M. J.;

577

Brouwer, A.; Bergman, A., Biotransformation of brominated flame retardants into potentially

578

endocrine-disrupting metabolites, with special attention to 2,2',4,4'-tetrabromodiphenyl ether

579

(BDE-47). Mol. Nutr. Food Res. 2008, 52, (2), 284-98.

580

31. Ren, X. M.; Guo, L. H., Molecular toxicology of polybrominated diphenyl ethers:

581

nuclear hormone receptor mediated pathways. Environ. Sci. Process Impact. 2013, 15, (4),

582

702-8.

583

32. Castrillo, A.; Tontonoz, P., Nuclear receptors in macrophage biology: at the crossroads of

584

lipid metabolism and inflammation. Annu. Rev. Cell Dev. Biol. 2004, 20, 455-80.

585

33. Chinenov, Y.; Gupte, R.; Rogatsky, I., Nuclear receptors in inflammation control:

586

repression by GR and beyond. Mol. Cell. Endocrinol. 2013, 380, (1-2), 55-64.

587

34. Knoedler, J. R.; Denver, R. J., Kruppel-like factors are effectors of nuclear receptor

588

signaling. Gen. Comp. Endocrinol. 2014, 203, 49-59.

589

35. Omiecinski, C. J.; Vanden Heuvel, J. P.; Perdew, G. H.; Peters, J. M., Xenobiotic

590

metabolism, disposition, and regulation by receptors: from biochemical phenomenon to

591

predictors of major toxicities. Toxicol. Sci. 2011, 120 Suppl 1, S49-75.

592

36. Aranda, A.; Pascual, A., Nuclear hormone receptors and gene expression. Physiol. Rev.

593

2001, 81, (3), 1269-304.

594

37. Jennings, P.; Limonciel, A.; Felice, L.; Leonard, M. O., An overview of transcriptional

595

regulation in response to toxicological insult. Arch. Toxicol. 2013, 87, (1), 49-72.

596

38. Jacobs, M. N.; Dickins, M.; Lewis, D. F., Homology modelling of the nuclear receptors: 29

ACS Paragon Plus Environment

Page 30 of 44

Page 31 of 44

Environmental Science & Technology

597

human oestrogen receptorbeta (hERbeta), the human pregnane-X-receptor (PXR), the Ah

598

receptor (AhR) and the constitutive androstane receptor (CAR) ligand binding domains from

599

the human oestrogen receptor alpha (hERalpha) crystal structure, and the human peroxisome

600

proliferator activated receptor alpha (PPARalpha) ligand binding domain from the human

601

PPARgamma crystal structure. J Steroid Biochem. Mol. Biol. 2003, 84, (2-3), 117-32.

602

39. Prossnitz, E. R.; Arterburn, J. B.; Smith, H. O.; Oprea, T. I.; Sklar, L. A.; Hathaway, H. J.,

603

Estrogen signaling through the transmembrane G protein-coupled receptor GrPR30. Annu.

604

Rev. Physiol. 2008, 70, 165-90.

605

40. Norman, A. W.; Mizwicki, M. T.; Norman, D. P., Steroid-hormone rapid actions,

606

membrane receptors and a conformational ensemble model. Nat. Rev. Drug Discov. 2004, 3,

607

(1), 27-41.

608

41. Lema, S. C.; Dickey, J. T.; Schultz, I. R.; Swanson, P., Dietary Exposure to 2,2 ',4,4

609

'-Tetrabromodiphenyl Ether (PBDE-47) Alters Thyroid Status and Thyroid

610

Hormone-Regulated Gene Transcription in the Pituitary and Brain. Environ. Health Persp.

611

2008, 116, (12), 1694-9.

612

42. Karpeta, A.; Ptak, A.; Gregoraszczuk, E. L., Different action of

613

2,2',4,4'-tetrabromodiphenyl ether (BDE-47) and its hydroxylated metabolites on ERalpha

614

and ERbeta gene and protein expression. Toxicol. Lett. 2014, 229, (1), 250-6.

615

43. Su, G.; Xia, J.; Liu, H.; Lam, M. H.; Yu, H.; Giesy, J. P.; Zhang, X., Dioxin-like potency

616

of HO- and MeO- analogues of PBDEs' the potential risk through consumption of fish from

617

eastern China. Environ. Sci. Technol. 2012, 46, (19), 10781-8.

618

44. Stegeman, J. J.; Goldstone, J. V.; Hahn, M. E., Perspectives on Zebrafish as a Model in 30

ACS Paragon Plus Environment

Environmental Science & Technology

619

Environmental Toxicology. Fish Physiol. 2010, 29, 367-439.

620

45. Tang, S.; Cai, Q.; Chibli, H.; Allagadda, V.; Nadeau, J. L.; Mayer, G. D., Cadmium

621

sulfate and CdTe-quantum dots alter DNA repair in zebrafish (Danio rerio) liver cells. Toxicol.

622

Appl. Pharmacol. 2013, 272, (2), 443-52.

623

46. Dai, Y. J.; Jia, Y. F.; Chen, N.; Bian, W. P.; Li, Q. K.; Ma, Y. B.; Chen, Y. L.; Pei, D. S.,

624

Zebrafish as a Model System to Study Toxicology. Environ. Toxicol. Chem. 2014, 33, (1),

625

11-7.

626

47. Baker, M. E.; Hardiman, G., Transcriptional analysis of endocrine disruption using

627

zebrafish and massively parallel sequencing. J. Mol. Endocrinol. 2014, 52, (3), R241-56.

628

48. Scholz, S.; Fischer, S.; Gundel, U.; Kuster, E.; Luckenbach, T.; Voelker, D., The

629

zebrafish embryo model in environmental risk assessment--applications beyond acute toxicity

630

testing. Environ. Sci. Pollut. Res. Int. 2008, 15, (5), 394-404.

631

49. Scholz, S.; Mayer, I., Molecular biomarkers of endocrine disruption in small model fish.

632

Mol. Cell. Endocrinol. 2008, 293, (1-2), 57-70.

633

50. Bertrand, S.; Thisse, B.; Tavares, R.; Sachs, L.; Chaumot, A.; Bardet, P. L.; Escriva, H.;

634

Duffraisse, M.; Marchand, O.; Safi, R.; Thisse, C.; Laudet, V., Unexpected novel relational

635

links uncovered by extensive developmental profiling of nuclear receptor expression. Plos

636

Genet. 2007, 3, (11), 2085-100.

637

51. Marsh, G.; Stenutz, R.; Bergman, A., Synthesis of hydroxylated and methoxylated

638

polybrominated diphenyl ethers - Natural products and potential polybrominated diphenyl

639

ether metabolites Eur. J. Org. Chem. 2003, 14, 2566-76.

640

52. Amsterdam, A.; Nissen, R. M.; Sun, Z.; Swindell, E. C.; Farrington, S.; Hopkins, N., 31

ACS Paragon Plus Environment

Page 32 of 44

Page 33 of 44

Environmental Science & Technology

641

Identification of 315 genes essential for early zebrafish development. Proc. Natl. Acad. Sci. U.

642

S. A. 2004, 101, (35), 12792-7.

643

53. Wen, Q.; Liu, H. L.; Zhu, Y. T.; Zheng, X. M.; Su, G. Y.; Zhang, X. W.; Yu, H. X.; Giesy,

644

J. P.; Lam, M. H., Maternal transfer, distribution, and metabolism of BDE-47 and its related

645

hydroxylated, methoxylated analogs in zebrafish (Danio rerio). Chemosphere 2014, 120C,

646

31-6.

647

54. Wen, Q.; Liu, H. L.; Su, G. Y.; Wei, S.; Feng, J. F.; Yu, H. X., Determination of

648

Polybrominated Diphenyl Ethers and Their Derivates in Zebrafish Eggs. Chinese J. Anal.

649

Chem. 2012, 40, (11), 1698-702.

650

55. Liu, C.; Wang, Q.; Liang, K.; Liu, J.; Zhou, B.; Zhang, X.; Liu, H.; Giesy, J. P.; Yu, H.,

651

Effects of tris(1,3-dichloro-2-propyl) phosphate and triphenyl phosphate on

652

receptor-associated mRNA expression in zebrafish embryos/larvae. Aquat. Toxicol. 2013,

653

128-129, 147-57.

654

56. Tang, S.; Allagadda, V.; Chibli, H.; Nadeau, J. L.; Mayer, G. D., Comparison of

655

cytotoxicity and expression of metal regulatory genes in zebrafish (Danio rerio) liver cells

656

exposed to cadmium sulfate, zinc sulfate and quantum dots. Metallomics 2013, 5, (10),

657

1411-22.

658

57. Cline, M. S.; Smoot, M.; Cerami, E.; Kuchinsky, A.; Landys, N.; Workman, C.;

659

Christmas, R.; Avila-Campilo, I.; Creech, M.; Gross, B.; Hanspers, K.; Isserlin, R.; Kelley, R.;

660

Killcoyne, S.; Lotia, S.; Maere, S.; Morris, J.; Ono, K.; Pavlovic, V.; Pico, A. R.; Vailaya, A.;

661

Wang, P.-L.; Adler, A.; Conklin, B. R.; Hood, L.; Kuiper, M.; Sander, C.; Schmulevich, I.;

662

Schwikowski, B.; Warner, G. J.; Ideker, T.; Bader, G. D., Integration of biological networks 32

ACS Paragon Plus Environment

Environmental Science & Technology

663

and gene expression data using Cytoscape. Nat. Protoc. 2007, 2, (10), 2366-82.

664

58. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N. S.; Wang, J. T.; Ramage, D.; Amin, N.;

665

Schwikowski, B.; Ideker, T., Cytoscape: a software environment for integrated models of

666

biomolecular interaction networks. Genome Res. 2003, 13, (11), 2498-504.

667

59. Saito, R.; Smoot, M. E.; Ono, K.; Ruscheinski, J.; Wang, P. L.; Lotia, S.; Pico, A. R.;

668

Bader, G. D.; Ideker, T., A travel guide to Cytoscape plugins. Nat. Methods 2012, 9, (11),

669

1069-76.

670

60. Pico, A. R.; Kelder, T.; van Iersel, M. P.; Hanspers, K.; Conklin, B. R.; Evelo, C.,

671

WikiPathways: pathway editing for the people. PLoS Biol. 2008, 6, (7), e184.

672

61. Escher, B. I.; Hermens, J. L. M., Modes of action in ecotoxicology: Their role in body

673

burdens, species sensitivity, QSARs, and mixture effects. Environ. Sci. Technol. 2002, 36,

674

(20), 4201-17.

675

62. Escher, B. I.; Hermens, J. L. M., Internal exposure: Linking bioavailability to effects.

676

Environ. Sci. Technol. 2004, 38, (23), 455A-62A.

677

63. Adolfsson-Erici, M.; Akerman, G.; McLachlan, M. S., Internal Benchmarking Improves

678

Precision and Reduces Animal Requirements for Determination of Fish Bioconcentration

679

Factors. Environ. Sci. Technol. 2012, 46, (15), 8205-11.

680

64. van Boxtel, A. L.; Kamstra, J. H.; Cenijn, P. H.; Pieterse, B.; Wagner, J. M.; Antink, M.;

681

Krab, K.; van der Burg, B.; Marsh, G.; Brouwer, A.; Legler, J., Microarray analysis reveals a

682

mechanism of phenolic polybrominated diphenylether toxicity in zebrafish. Environ. Sci.

683

Technol. 2008, 42, (5), 1773-9.

684

65. Yu, Y.; Yang, W.; Gao, Z.; Lam, M. H. W.; Liu, X.; Wang, L.; Yu, H., RP-HPLC 33

ACS Paragon Plus Environment

Page 34 of 44

Page 35 of 44

Environmental Science & Technology

685

measurement and quantitative structure-property relationship analysis of the noctanol-water

686

partitioning coefficients of selected metabolites of polybrominated dyphenyl ethers. Environ.

687

Chem. 2008, 5, (5), 332-9.

688

66. Aardema, M. J.; MacGregor, J. T., Toxicology and genetic toxicology in the new era of

689

"toxicogenomics": impact of "-omics" technologies. Mutat. Res. 2002, 499, (1), 13-25.

690

67. Lema, S. C.; Dickey, J. T.; Schultz, I. R.; Swanson, P., Dietary exposure to

691

2,2',4,4'-tetrabromodiphenyl ether (PBDE-47) alters thyroid status and thyroid

692

hormone-regulated gene transcription in the pituitary and brain. Environ. Health Perspect.

693

2008, 116, (12), 1694-9.

694

68. Chan, W. K.; Chan, K. M., Disruption of the hypothalamic-pituitary-thyroid axis in

695

zebrafish embryo-larvae following waterborne exposure to BDE-47, TBBPA and BPA. Aquat.

696

Toxicol. 2012, 108, 106-11.

697

69. Ren, X. M.; Guo, L. H.; Gao, Y.; Zhang, B. T.; Wan, B., Hydroxylated polybrominated

698

diphenyl ethers exhibit different activities on thyroid hormone receptors depending on their

699

degree of bromination. Toxicol. Appl. Pharmacol. 2013, 268, (3), 256-63.

700

70. Zhou, H.; Wu, H.; Liao, C.; Diao, X.; Zhen, J.; Chen, L.; Xue, Q., Toxicology

701

mechanism of the persistent organic pollutants (POPs) in fish through AhR pathway. Toxicol.

702

Mech. Methods 2010, 20, (6), 279-86.

703

71. Shibata, H.; Spencer, T. E.; Onate, S. A.; Jenster, G.; Tsai, S. Y.; Tsai, M. J.; O'Malley, B.

704

W., Role of co-activators and co-repressors in the mechanism of steroid/thyroid receptor

705

action. Recent Prog. Horm. Res. 1997, 52, 141-64; discussion 164-5.

706

72. Gjernes, M. H.; Schlenk, D.; Arukwe, A., Estrogen receptor-hijacking by dioxin-like 34

ACS Paragon Plus Environment

Environmental Science & Technology

707

3,3'4,4',5-pentachlorobiphenyl (PCB126) in salmon hepatocytes involves both receptor

708

activation and receptor protein stability. Aquat. Toxicol. 2012, 124-125, 197-208.

709

73. Grans, J.; Wassmur, B.; Celander, M. C., One-way inhibiting cross-talk between

710

arylhydrocarbon receptor (AhR) and estrogen receptor (ER) signaling in primary cultures of

711

rainbow trout hepatocytes. Aquat. Toxicol. 2010, 100, (3), 263-70.

712

74. Kojima, H.; Takeuchi, S.; Uramaru, N.; Sugihara, K.; Yoshida, T.; Kitamura, S., Nuclear

713

hormone receptor activity of polybrominated diphenyl ethers and their hydroxylated and

714

methoxylated metabolites in transactivation assays using Chinese hamster ovary cells.

715

Environ. Health Perspect. 2009, 117, (8), 1210-8.

716

75. Pacyniak, E. K.; Cheng, X.; Cunningham, M. L.; Crofton, K.; Klaassen, C. D.; Guo, G.

717

L., The flame retardants, polybrominated diphenyl ethers, are pregnane X receptor activators.

718

Toxicol. Sci. 2007, 97, (1), 94-102.

719

76. Liu, H.; Hu, W.; Sun, H.; Shen, O.; Wang, X.; Lam, M. H.; Giesy, J. P.; Zhang, X.; Yu, H.,

720

In vitro profiling of endocrine disrupting potency of 2,2',4,4'-tetrabromodiphenyl ether

721

(BDE47) and related hydroxylated analogs (HO-PBDEs). Mar. Pollut. Bull. 2011, 63, (5-12),

722

287-96.

723

77. Fery, Y.; Buschauer, I.; Salzig, C.; Lang, P.; Schrenk, D., Technical pentabromodiphenyl

724

ether and hexabromocyclododecane as activators of the pregnane-X-receptor (PXR).

725

Toxicology 2009, 264, (1-2), 45-51.

726

78. Abbott, B. D., Review of the expression of peroxisome proliferator-activated receptors

727

alpha (PPAR alpha), beta (PPAR beta), and gamma (PPAR gamma) in rodent and human

728

development. Reprod. Toxicol. 2009, 27, (3-4), 246-57. 35

ACS Paragon Plus Environment

Page 36 of 44

Page 37 of 44

Environmental Science & Technology

729

79. Tung, E. W.; Boudreau, A.; Wade, M. G.; Atlas, E., Induction of adipocyte differentiation

730

by polybrominated diphenyl ethers (PBDEs) in 3T3-L1 cells. Plos One 2014, 9, (4), e94583.

731

80. Greenwood, A. K.; Butler, P. C.; White, R. B.; DeMarco, U.; Pearce, D.; Fernald, R. D.,

732

Multiple corticosteroid receptors in a teleost fish: distinct sequences, expression patterns, and

733

transcriptional activities. Endocrinology 2003, 144, (10), 4226-36.

734

36

ACS Paragon Plus Environment

Environmental Science & Technology

Page 38 of 44

Table 1 Estimated internal doses of BDE-47, 6-OH-BDE-47, and 6-MeO-BDE-47 exposures in zebrafish larvae at 120 hpf. The calculation of internal dose for each compound was based on the ratio of the exposure concentration at 300 µg/L and the measured concentration in larvae at 120 hpf. Chemical

BDE-47

6-MeO-BDE-47

6-OH-BDE-47

Nominal Concentration (µM)

0.1

0.5

2.5

0.1

0.5

2.5

0.008

0.02

0.05

BDE-47 (µg/g)

9.05

45.27

226.33













6-MeO-BDE-47 (µg/g)







18.17

90.87

454.37







6-OH-BDE-47 (µg/g)













0.024

0.061

0.15

Biotransformed 6-OH-BDE-47 (µg/g)







0.04

0.21

1.07







37

ACS Paragon Plus Environment

Page 39 of 44

Environmental Science & Technology

Figure 1 Measured concentrations in exposure medium (µg/L) and in zebrafish embryos-larvae (mg/g, wm) after exposure to 300 µg/L of BDE-47 (0.62 µM), 6-MeO-BDE-47 (0.58 µM), or 6-OH-BDE-47 (0.6 µM) across time-points (hpf). Generalized additive model (GAM) plots between concentrations (after log transformation) in exposure water or in zebrafish embryos-larvae and time (hpf) are given. Shaded areas are the 95% confidence intervals for each GAM. The F-statistics, p-values, and adjusted R2 for the specific GAMs are given in each plot, while D shows the deviance explained.

µ 00 10

0µ 10

Concentration

µ 10

1

g/L

g/L

g/L

F=87.29, p=3.89e-05, R-sq.(adj)=0.935, D=94.6%

6-MeO-BDE-47-Zfish

F=53.18, p=0.0151, R-sq.(adj)=0.973, D=98.9%

6-OH-BDE-47-Water

F=169.8, p=0.00265, R-sq.(adj)=0.992, D=99.7%

6-OH-BDE-47-Zfish

/g mg r o /L µg

0.1

6-MeO-BDE-47-Water

BDE-47-Water

F=10.13, p=0.124, R-sq.(adj)=0.873, D=94.8% F=107.7, p=7.98e-06, R-sq.(adj)=0.982, D=98.9%

BDE-47-Zfish

F=40.68, p=0.00562, R-sq.(adj)=0.953, D=97.3%

/g mg

/g mg 1 0.0

g/g m 01 0.0

0

12

24

48

72

96

Time (hpf) 38

ACS Paragon Plus Environment

120

Environmental Science & Technology

Page 40 of 44

Figure 2 A). Bioconcentration factors (BCF) calculated after exposure to 300 µg/L of BDE-47 (0.62 µM), 6-MeO-BDE-47 (0.58 µM), and 6-OH-BDE-47 (0.6 µM) across time-points (hpf). BCF was calculated based on the ratio of measured concentrations in zebrafish embryos-larvae (µg/kg, wm) and measured concentrations in exposure medium (µg/L) at a specific time (hpf). B). Measured concentrations of biotransformed 6-OH-BDE-47 in zebrafish embryos-larvae (µg/g, wm) after exposure to 300 µg 6-MeO-BDE-47/L (0.58 µM) across several time-points (hpf). A linear regression between concentrations in zebrafish embryos-larvae and time (hpf) is given (R2=0.911; p < 0.05).

A.

B. 2500

0.3

R2 = 0.911 y = -0.120 + 0.003

x

0.2 1500 6-MeO-BDE-47 6-OH-BDE-47 BDE-47

1000

Concentration (µg/g)

Bioconcentration Factors (BCFs)

2000

0.1

500

0.0

0

12

24

48

72

96

120

48

72

96

Time (hpf)

Time (hpf)

39

ACS Paragon Plus Environment

120

Page 41 of 44

Environmental Science & Technology

Figure 3 Dendrogram displaying similarities of chemicals and doses based on effects on genes in nuclear receptor pathways for BDE-47, 6-MeO-BDE-47, and 6-OH-BDE-47 in zebrafish larvae at 120 hpf. A). The dendrogram of hierarchical cluster analysis was calculated using the average gene expression values (63 genes in total) of the three or four biological replicates per exposure. Samples names are composed by the name of exposure compound followed by the exposure concentration (µM). Different colors in the dendrogram denoted five clustering groups; B). The heatmap of gene expression profiles was generated using the average gene expression values of the three or four biological replicates per exposure. The fold-changes of gene expression are given in the respective cells and genes involved in different receptor pathways are given different colors (see legend). A.

100

60 40

40

ACS Paragon Plus Environment

6 - OH - 0.02

6 - OH - 0.05 6 - OH - 0.008

6 - MeO - 0 BDE - 47 - 0

6 - OH - 0

6 - MeO - 0.1 BDE - 47 - 0.1 BDE - 47 - 0.5

0

6 - MeO - 0.5 BDE - 47 - 2.5

20 6 - MeO - 2.5

Distance

80

Environmental Science & Technology

B.

41

ACS Paragon Plus Environment

Page 42 of 44

Page 43 of 44

Environmental Science & Technology

Figure 4 Interaction network of selected genes in nuclear steroid receptor pathways of zebrafish. Nodes represent single genes, edges either protein-protein or protein-DNA interactions. Statistically significant changes (p < 0.05) in gene expression following different concentrations of treatment of BDE-47, 6-MeO-BDE-47, and 6-OH-BDE-47 at 120 hpf are given in the respective boxes (see legend).

42

ACS Paragon Plus Environment

Environmental Science & Technology

43

ACS Paragon Plus Environment

Page 44 of 44