Subscriber access provided by Iowa State University | Library
Ecotoxicology and Human Environmental Health
Long-term exposure to environmental concentrations of azoxystrobin delays sexual development and alters reproduction in zebrafish (Danio rerio) Fangjie Cao, Christopher J Martyniuk, Peizhuo Wu, Feng Zhao, Sen Pang, Chengju Wang, and Lihong Qiu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b05829 • Publication Date (Web): 07 Jan 2019 Downloaded from http://pubs.acs.org on January 8, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 27
Environmental Science & Technology
1
Long-term exposure to environmental concentrations of azoxystrobin delays sexual development
2
and alters reproduction in zebrafish (Danio rerio)
3
Fangjie Cao,† Christopher J. Martyniuk,‡ Peizhuo Wu,† Feng Zhao,† Sen Pang,† Chengju Wang,† and
4
Lihong Qiu†*
5
† Department
6
China
7
‡
8
University of Florida Genetics Institute, Interdisciplinary Program in Biomedical Sciences Neuroscience,
9
College of Veterinary Medicine, University of Florida, Gainesville, FL, 32611, USA
of Applied Chemistry, College of Science, China Agricultural University, Beijing, 100193,
Department of Physiological Sciences and Center for Environmental and Human Toxicology,
10
*
11
Lihong Qiu,
12
E-mail:
[email protected]; Tel: +86 (0)10 62733924
Corresponding author:
13 14 15 16 17 18 19 20 21 22 1
ACS Paragon Plus Environment
Environmental Science & Technology
23
ABSTRACT: The strobilurin fungicide azoxystrobin (AZO) can induce adverse effects in aquatic
24
organisms, but data are lacking on endpoints associated with sexual development and reproduction
25
following a chronic exposure to AZO. In this study, zebrafish embryos (F0) at 2–4 hours post-fertilization
26
(hpf) were exposed to 0.2, 2.0 and 20.0 µg/L AZO until 120 days post-fertilization (dpf). Decreased male
27
ratio and increased intersex ratio were observed by 20.0 µg/L AZO at 42 and 60 dpf, but this effect
28
disappeared at 120 dpf. AZO at 20.0 µg/L inhibited growth, retarded gonadal development, and disrupted
29
sex hormone and vitellogenin in females at 60 and 120 dpf, and in males at 42, 60 and 120 dpf. These
30
effects were associated with altered expression of cyp19a, cyp19b, hsd3b, hsd17b, vtg1 and vtg2.
31
Exposure to 2.0 µg/L AZO altered mRNA levels of these transcripts in females at 120 dpf, and in males
32
at 60 and 120 dpf. Reproduction ability was reduced by 20.0 µg/L AZO at 120 dpf. Developmental
33
defects were observed after F1 embryos from exposed parents of 20.0 µg/L were reared in AZO-free
34
water at 96 hpf. Overall, these data provide new understanding of fish sexual development and
35
reproduction following chronic exposures to AZO.
36 37 38 39 40 41 42 43 44 2
ACS Paragon Plus Environment
Page 2 of 27
Page 3 of 27
Environmental Science & Technology
45
INTRODUCTION
46
Azoxystrobin (AZO), a strobilurin fungicide, remains one of the most widely used fungicides
47
around the world (i.e. USA, Australia, UK and China),1 and is worth $1.27 billion in global sales reported
48
by Phillips McDougall in 2016.2 AZO is also one of the frequently detected fungicides in aquatic
49
environment,3-6 with mean detected concentration of 0.16 µg/L in streams in the US,7 1.43 µg/L in 20
50
streams in Lower Saxony, Germany,8 0.15 µg/L in the surface water in Demark,9 and 34 µg/L in river
51
water in Shanghai, China.4 AZO is a contaminant of concern because it induces both lethal and sub-lethal
52
effects in aquatic organisms.10-17 Previous studies report that AZO disrupts the endocrine system of
53
zebrafish.12-14 For example, a 21-day exposure to 200 µg/L AZO decreased 17β-estradiol (E2) and
54
vitellogenin (Vtg), increased testosterone (T), rewarded ovarian development, up-regulated hsd3b,
55
hsd17b and cyp17, as well as down-regulated cyp19a, cyp19b, vtg1 and vtg2 in female zebrafish;
56
moreover, male zebrafish displayed increased E2 and Vtg, decreased T, inhibited testicular development,
57
up-regulated hsd3b, hsd17b, cyp17, cyp19a, cyp19b, vtg1 and vtg2, which impaired reproductive
58
endpoints.14 Vtg and E2 were increased in zebrafish embryos following 48 and 72 exposure to 0.1 and
59
10.0 µg/L AZO, and an increase in vtg and a decrease in cyp19b were induced after 48 and 72 h exposure
60
to 1, 10 and 100 µg/L AZO.12 Similarly, a 21-day parental exposure to 200 µg/L up-regulated cyp19a,
61
cyp19b, vtg1 and vtg2 in F1 zebrafish embryos.13 However, data are lacking regarding long-term impacts
62
of AZO at environmental concentrations on sexual development and reproduction in F0 generation, and
63
developmental success in F1 offspring of fish.
64
Zebrafish (Danio rerio) serve as a good model to assess the effects of chemicals on sexual
65
development.18 Zebrafish sexual development progresses as follows: (1) initially, juveniles develop
66
undifferentiated, ovary-like gonads at approximately 10 days post-fertilization (dpf); (2) sexual 3
ACS Paragon Plus Environment
Environmental Science & Technology
67
differentiation begins at approximately 21−23 dpf, with about 50% of the ovaries beginning to undergo
68
reabsorption and transformation into testes, and ends at about 42 dpf; (3) the development and maturation
69
of ovaries and testes are almost completed at ~60 dpf.19,
70
influenced by sex steroid hormones,21 sex determination genes (foxl2, brca2, sox9b, sox9a, amh and
71
dmrt1, etc.),22, 23 and environmental factors (e.g. chemicals).18, 24, 25
20
Zebrafish sexual development can be
72
Estrogens (i.e. E2) regulate feminization by binding to the estrogen receptor,26 while androgens (11-
73
ketotestosterone and T) regulate masculinization by binding to the androgen receptor.27 Steroidogenic
74
enzymes, such as CYP11, CYP17, CYP19, HSD3B, and HSD17B, determine the synthetic rate of
75
estrogens and androgens, and can be dysregulated by chemical exposure.28, 29 Estrogens activate the
76
expression of vtg to synthesize Vtg protein, the major precursor of egg yolk. Vtg plays an important role
77
in the oocytes.28 While vtg gene is present in male fish, its expression is very low, unless it is activated
78
by xenoestrogens.30 Thus, Vtg and vtg expression are strong and widely used biomarkers for the
79
estrogenic chemicals.
80
The objectives of this study were to evaluate the effects of a long-term exposure to environmentally
81
relevant concentrations of AZO on sexual development and reproduction in F0 individuals, beginning at
82
fertilization and progressing up to 120 days post-fertilization (dpf). Growth, gonadosomatic index (GSI),
83
sex steroid hormones, Vtg, sex ratio, gonadal histology, expression levels of genes in the HPG axis, and
84
reproductive capacity in F0 generation were measured. Lastly, developmental success in F1 offspring
85
was also assessed, as our previous study demonstrated that AZO exhibited transgenerational effects in
86
zebrafish.13
87
MATERIALS AND METHODS
88
Chemicals 4
ACS Paragon Plus Environment
Page 4 of 27
Page 5 of 27
Environmental Science & Technology
89
Azoxystrobin (AZO, CAS: 131860-33-8, 98% purity) was obtained from Shenyang Chemical
90
Industry Research Institute, China. A stock solution of 10,000 mg/L AZO was prepared in analytical-
91
grade acetone containing 5% Tween 80,14 and was stored at 4 °C for further dilution. Reconstituted water
92
with a pH value of 7.5 ± 0.5 was used to make exposure solutions, which contained 1.27 mM NaHCO3,
93
0.33 mM MgSO4, 0.33 mM CaCl2, and 0.17 mM KCl.31
94
Zebrafish maintenance and egg production
95
Parental zebrafish (wild-type AB-strain; six months old) were purchased from Beijing
96
Hongdagaofeng Aquarium Department, China and were housed in a flow-through feeding equipment
97
(Esen Corp., China). Normally fertilized embryos at 2−4 hours post-fertilization (hpf) were collected for
98
further exposure. The condition for zebrafish maintenance and AZO exposure were maintained with a
99
pH value of 7.5 ± 0.5, a conductivity value of 600 ± 100 μS.cm−1, a light/dark cycle of 14:10 h, a
100
minimum dissolved oxygen concentration of up to 80% of air saturation, and a temperature of 27 ± 1 °C.
101
AZO exposure assay
102
To determine the long-term effects of AZO on sexual differentiation and reproduction in F0
103
generation, as well as developmental success in F1 offspring, zebrafish embryos (F0) at 2–4 hpf were
104
randomly assigned into groups and exposed to water control, solvent control (0.02‱ acetone and 0.001‱
105
Tween 80, v/v), 0.2, 2.0 and 20.0 µg/L AZO (equal to 0.496, 4.96 and 49.6 nM, respectively) for 120
106
days. The growth, sex ratio, gonadal development, sex hormone and Vtg, mRNA levels of genes in HPG
107
axis, and reproductive capacity in F0 generation were assessed. Additionally, development and survival
108
in F1 embryos originating from F0 parents were determined at 96 hpf. Exposure concentrations were
109
selected to span a wide range of environmentally relevant concentrations reported in the aquatic
110
environment.4, 7-9 Moreover, our previous study demonstrated that 2.0 and 20.0 µg/L AZO altered cyp11a, 5
ACS Paragon Plus Environment
Environmental Science & Technology
111
hsd3b, cyp19a, vtg1 and vtg2 in male zebrafish.14 The experimental design is presented in Figure S1 and
112
sample size is shown in Table S1. Abnormal phenotypes and mortality were recorded daily. Dead
113
individuals were removed immediately and approximately 90% of exposure solutions were renewed
114
daily. Institutional Animal Care and Use Committee (IACUC) of China Agricultural University approved
115
the experimental procedures.
116
The procedure for F0 exposure was as follows: (1) Initially, 200 fertilized embryos at 2–4 hpf were
117
maintained in glass Petri dishes of 20-cm diameter (N = 4 replicates per treatment) containing 200 mL
118
of exposure solutions. (2) After hatching (approx. 96 hpf), zebrafish larvae were transferred to 2-L
119
beakers containing 1.6-L of exposure solutions. Zebrafish larvae were fed daily starting at 5 dpf. The
120
detailed feeding procedure is provided in the SI. (3) After 21 dpf, ~180 individuals (< 10% mortality
121
occurred before 21 dpf) from each beaker were transferred to 25-L tanks containing 20-L of exposure
122
solutions until 120 dpf.
123
At 90 dpf, 10 females and 10 males were selected randomly from each tank and placed into a 10-
124
L breeding tank containing 6-L of exposure solution. Each breeding tank was equipped with a removable
125
spawning tray of stainless steel mesh at the bottom of the breeding tank so that embryos could fall through.
126
One week prior to the end of F0 exposure (about 113 dpf), females and males were separated using
127
dividers, continuously exposed to AZO until 120 dpf. At 120 dpf, water in all breeding tanks was changed
128
to reconstituted water without AZO and the spawning was stimulated at 8:00 PM in the next morning
129
when the light was turned on, followed by egg collection one hour later. F1 embryos of each breeding
130
tank were rinsed three times with chemical-free reconstituted water and were examined under a
131
stereomicroscope to calculate average egg production per female and fertilization rate. Fifty fertilized
132
and normally developed F1 embryos at 2–4 hpf were reared in AZO-free reconstituted water until 96 hpf 6
ACS Paragon Plus Environment
Page 6 of 27
Page 7 of 27
Environmental Science & Technology
133
to calculate mortality, hatching rate, and malformation rate. The images of morphological abnormity
134
were photographed using an Aigo GE-5 digital microscope (Aigo Digital Technology Co. Ltd, China).
135
Sample collections in F0 generation at 42, 60 and 120 dpf were conducted as follows: (1) Five
136
individuals of each gender from each replicate were randomly selected and euthanized with anesthetic
137
buffer (100 mg/L tricaine and 300 mg/L NaHCO3). Wet body weight and total body length were
138
measured using an electronic balance and digital caliper (Mitutoyo Corp., Japan), respectively.
139
Additionally, weights of female ovaries and male testes in each replicate were measured to calculate GSI.
140
(2) Twenty fish were randomly selected from each replicate to assess the phenotype and sex ratio using
141
histological observation. (3) Five females and males in each replicate were sampled to gonadal
142
development. (4) Blood samples of 2 fish of each gender from each replicate were collected, pooled,
143
frozen and stored at −80 °C for sex steroid hormones and Vtg measurement. (5) The brains, gonads, and
144
livers of 3 individuals of each gender per replicate were pooled, frozen and stored at −80 °C for RNA
145
extraction, respectively.
146
Sex ratio and gonadal histological examination
147
The sample preparation for gonadal histology was conducted according to a published protocol.14
148
Gonadal histology in F0 zebrafish was photographed using an Olympus microscope (Olympus, Japan).
149
Two cross-sections of each gonad (five females and five males per replicate) were selected to evaluate
150
gonadal development.32
151
Sex steroid hormones and Vtg measurement
152
Detailed information on the sex steroid hormones (E2 and T) and Vtg measurements in zebrafish can
153
be found in SI.
154
Gene expression analysis 7
ACS Paragon Plus Environment
Environmental Science & Technology
155
Total RNA was extracted using RNApre Pure Tissue Kit with DNase digestion (Tiangen Biotech,
156
Beijing, China). The quality of all RNA samples was verified using gel electrophoresis. RNA
157
concentration was determined using NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc.,
158
USA). Total RNA (1.5 μg) per sample was reverse-transcribed into cDNA using FastQuant cDNA
159
Synthesis Kit (Tiangen Biotech, Beijing, China) and stored at −20 °C for further analysis. Real-time
160
quantitative polymerase chain reaction (RT-qPCR) was performed with an ABI 7500 PCR Detection
161
System (ABI, USA) using SYBR Green PCR Master Mix Reagent Kit (Tiangen Biotech, Beijing, China).
162
No reverse transcriptase (NRT) controls were prepared to assess genomic DNA contamination and there
163
was no genomic DNA contamination. The RT-qPCR procedure was conducted as follows: 95 °C for 15
164
min and 40 cycles of 95 °C for 10 s, 60 °C for 20 s, and 72 °C for 32 s, followed by a melt curve analysis.
165
Specific primers of target genes were obtained from published sequences of zebrafish (Table S2).
166
The amplification efficiency and linear coefficient of the standard curve of all genes are provided in SI
167
(Table S2). The gene of β-actin was selected as a housekeeping gene to normalize mRNA levels of target
168
genes due to its stable expression among different tissues and individuals of zebrafish.33 The stable
169
expression of β-actin in tissue categories, genders, and ages among different treatments is shown in
170
Figure S2. Melt curves of all tested genes are presented in Figure S3-S6. Each treatment contained four
171
biological replicates run in triplicate and relative mRNA levels of target genes were analyzed using the
172
2 −∆∆Ct method.34
173
Chemical analysis
174
Exposure solutions (10 mL) in each replicate were sampled at 0, 42, 60 and 120 dpf, and transferred
175
into 50-mL centrifuge tubes for chemical analysis. Detailed information of the analytical method of AZO
176
is presented in SI. 8
ACS Paragon Plus Environment
Page 8 of 27
Page 9 of 27
177
Environmental Science & Technology
Data analysis
178
Data analyses were conducted using Graph-Pad Prism version 6.0 (GraphPad Software Inc., USA).
179
There was no significant difference between water control and solvent control for all investigated
180
endpoints and thus the solvent control was used as the control in the following analysis. One-way analysis
181
of variance (one-way ANOVA) was conducted to analyze significant differences between the control
182
and treatment groups. Clustering maps were constructed using the MultiExperimental Viewer v4.9
183
(Dana-Farber Cancer Institute, Boston, MA, USA). Values are shown as the mean ± standard deviation
184
(SD). A significant difference was noted when p < 0.05.
185
RESULTS
186
Chemical analysis
187
Chemical analysis of AZO in the water is provided in the SI. The deviations between the actual
188
concentrations and nominal concentrations of AZO at 0, 42, 60 and 120 dpf were maintained within ±
189
20% (Table S3). Therefore, the nominal concentrations were used to represent the actual concentrations
190
of AZO.35
191
Growth and GSI in F0 zebrafish
192
AZO at 0.2 and 2.0 µg/L did not affect growth or GSI in F0 female and male zebrafish following
193
120 days exposure (Table S4 and S5). However, significant reductions in wet body weight, total body
194
length, and GSI were induced in F0 females exposed to 20.0 µg/L AZO at 60 and 120 dpf compared with
195
the control group (Table S4). In F0 males, 20.0 µg/L AZO significantly decreased wet body weight, total
196
body length and GSI at 42, 60 and 120 dpf (Table S5). No significant alteration in CF was observed in
197
either gender of zebrafish at the tested doses (Table S4 and S5).
198
Sex ratio in F0 zebrafish 9
ACS Paragon Plus Environment
Environmental Science & Technology
199 200
Figure 1. Representative images of undifferentiated and intersex gonads in F0 zebrafish exposed to 20.0
201
µg/L AZO at 42 (A and B) and 60 (C and D) dpf, and the relative percentages of F0 female, male, intersex
202
and undifferentiated zebrafish (E) at 42, 60 and 120 dpf. Blue arrows indicate the undifferentiated gonads
203
in zebrafish. PO indicates perinucleolar oocytes, and Sg, Sc, St and Sz represent spermatogonia,
204
spermatocytes, spermatids and spermatozoa, respectively (400× magnification in A, B, C and D). Sex
205
ratio was calculated as the number of individuals of one phenotype × 100/the total number of individuals
206
that were selected for all endpoints at 42, 60 and 120 dpf, respectively. Values are shown as the mean ±
207
standard deviation (SD) (N = 4 replicates per treatment, n = 50 fish per replicate at 42 and 60 dpf, and n
208
= 80 at 120 dpf). Asterisks (in red) denote significant differences between the control and treatments (*p
209
< 0.05 and **p < 0.01 by Dunnett’s post-hoc test).
210
AZO at 20.0 µg/L disrupted sex ratio and sexual differentiation in F0 zebrafish at 42 and 60 dpf.
211
Undifferentiated gonads and intersex gonads were observed in zebrafish treated with 20.0 µg/L AZO at
212
42 dpf (Figure 1A and B), and intersex phenotypes were still detected at 60 dpf when compared to the
213
control group (Figure 1C and D). Exposure to 20.0 µg/L AZO significantly decreased the percentages of
214
females and males at 42 dpf, whereas the percentages of intersex and undifferentiated phenotypes were
215
significantly increased (Figure 1E). Likewise, a significant decrease in the percentages of males and an 10
ACS Paragon Plus Environment
Page 10 of 27
Page 11 of 27
Environmental Science & Technology
216
increase in intersex phenotypes were observed in 20.0 µg/L AZO group at 60 dpf (Figure 1E). However,
217
no intersex and undifferentiated phenotype in zebrafish treated with 0.2, 2.0 or 20.0 µg/L AZO at 120
218
dpf were observed, nor was the sex ratio changed in zebrafish at 120 dpf (Figure 1E).
219
Gonadal development in F0 zebrafish
220 221
Figure 2. Representative images of the ovaries in F0 female zebrafish exposed to control (A, D and G),
222
2.0 (B, E and H) and 20.0 µg/L (C, F and I) AZO at 42, 60 and 120 dpf, respectively. The oocytes
223
included perinucleolar oocytes (PO), cortical alveolar oocytes (CO), early vitellogenic oocytes (EV)
224
and late vitellogenic oocytes (LV) (100× magnification).
11
ACS Paragon Plus Environment
Environmental Science & Technology
225 226
Figure 3. Representative images of the testes in F0 male zebrafish exposed to control (A, D and G), 2.0
227
(B, E and H) and 20.0 µg/L (C, F and I) AZO at 42, 60 and 120 dpf, respectively. The spermatocytes
228
included spermatogonia (Sg), spermatocytes (Sc), spermatids (St) and spermatozoa (Sz) (200×
229
magnification).
230 231
Figure 4. The relative percentages of different stages of oocytes in F0 females (A) and spermatocytes in 12
ACS Paragon Plus Environment
Page 12 of 27
Page 13 of 27
Environmental Science & Technology
232
males (B) exposed to control, 0.2, 2.0 and 20.0 µg/L AZO at 42, 60 and 120 dpf. Values are shown as
233
the mean ± standard deviation (SD) (N = 4 replicates per treatment, n = 5 fish per gender of each
234
replicate with 10 cross-sections). Asterisks denote significant differences between the control and
235
treatments (*p < 0.05 and **p < 0.01 by Dunnett’s post-hoc test).
236
Gonadal development in F0 females and males was inhibited by 20.0 µg/L AZO exposure (Figure 2
237
and 3). In female ovaries, exposure to 20.0 µg/L AZO significantly increased the relative percentage of
238
perinucleolar oocytes (immature oocytes) at 42, 60 and 120 dpf, whereas the relative percentage of late
239
vitellogenic oocytes (mature oocytes) were significantly decreased at 60 and 120 dpf compared to the
240
control (Figure 4A). In male testes, 20.0 µg/L AZO induced a significant increase in the percentage of
241
spermatogonia at 42 dpf (Figure 4B). A significant increase in the relative percentage of spermatogonia
242
and spermatocytes (immature spermatocytes) was observed in 20.0 µg/L AZO group at 60 and 120 dpf,
243
whereas 20.0 µg/L AZO significantly reduced the relative percentage of spermatozoa (mature
244
spermatocytes) at 60 and 120 dpf (Figure 4B).
245
Sex steroid hormone and Vtg in F0 zebrafish
246
F0 females exhibited no change in E2, T and Vtg, nor the E2/T ratio following 120 days exposure
247
to 0.2 and 2.0 µg/L AZO. However, in F0 females, exposure to 20.0 µg/L AZO induced a significant
248
decrease in E2 and an increase in T at 60 and 120 dpf compared to the control (Table S6). Moreover, the
249
E2/T ratio and Vtg were significantly reduced in females exposed to 20.0 µg/L AZO at 42, 60 and 120
250
dpf (Table S6). In F0 males, E2, T and Vtg, or the E2/T ratio were not altered following 120 days
251
exposure to 0.2 and 2.0 µg/L AZO, except for a significant increase in Vtg in the 2.0 µg/L AZO group
252
at 120 dpf (Table S7). Moreover, exposure to 20.0 µg/L AZO significantly increased E2 and Vtg, E2/T
253
ratio, and decreased T at 42, 60 and 120 dpf (Table S7). 13
ACS Paragon Plus Environment
Environmental Science & Technology
254
Gene expression alteration in F0 zebrafish
255 256
Figure 5. Clustering map of expression levels of esr1, esr2b, ar and cyp19b in the brains (B), cyp11a,
257
cyp17, cyp19a, hsd3b, and hsd17b in the gonads (G), and vtg1 and vtg2 in the livers (L) of female (F)
258
and male (M) zebrafish exposed to control (0.0), 0.2, 2.0 and 20.0 µg/L AZO at 42, 60 and 120 dpf. The
259
mean values of normalized transcript levels in each treatment were used to draw the clustering map (N
260
= 4 replicates per treatment). Asterisks denote significant differences between the control and treatments
261
(*p < 0.05 and **p < 0.01 by Dunnett’s post-hoc test). All expression data are presented in SI.
262
In F0 female brains, there were no changes in the expression of esr1, esr2b, or ar following 120 days
263
exposure to 0.2, 2.0 or 20.0 µg/L AZO, whereas cyp19b were significantly down-regulated in the 20.0
264
µg/L AZO group at 60 and 120 dpf compared to the control (Figure 5 and Figure S7A, C and E). In F0
265
male brains, cyp19b was significantly up-regulated in the 20.0 µg/L AZO group at 42, 60 and 120 dpf,
266
and in the 2.0 µg/L AZO group at 60 and 120 dpf (Figure 5 and Figure S7B, D and F).
267
In F0 female ovaries, cyp19a was down-regulated in the 20.0 µg/L AZO group at 42, 60 and 120 dpf,
268
and in the 2.0 µg/L AZO group at 120 dpf compared with the control (Figure 5 and Fig. S8A, C and E).
269
Expression levels of hsd3b and hsd17b were significantly increased in the 20.0 µg/L AZO group at 42,
270
60 and 120 dpf, and in the 2.0 µg/L AZO group at 120 dpf. Likewise, cyp17 was up-regulated in the 2.0
271
and 20.0 µg/L AZO groups at 120 dpf (Figure 5 and Figure S8A, C and E). In F0 male testes, 20.0 µg/L 14
ACS Paragon Plus Environment
Page 14 of 27
Page 15 of 27
Environmental Science & Technology
272
AZO significantly up-regulated cyp17, hsd3b, hsd17b and cyp19a at 42, 60 and 120 dpf, and up-
273
regulation of these genes were detected in the 2.0 µg/L AZO group at 120 dpf (Figure 5 and Figure S8B,
274
D and F). Conversely, AZO did not alter genes related to sex differentiation (i.e. foxl2, brca2, sox9a,
275
sox9b, amh and dmrt1) in female ovaries or male testes at 42 and 60 dpf (Figure S9).
276
In F0 female livers, vtg1 and vtg2 were significantly decreased in 20.0 µg/L AZO group at 42, 60
277
and 120 dpf, and in the 2.0 µg/L AZO group at 120 dpf compared to the control (Figure 5 and Figure
278
S10A, C and E). In F0 male livers, vtg1 and vtg2 were significantly increased in the 20.0 µg/L AZO
279
group at 42, 60 and 120 dpf, and in 2.0 µg/L AZO group at 120 dpf (Figure 5 and Figure S10B, D and
280
F).
281
Reproductive capacity in F0 zebrafish and developmental success in F1 embryos
282
AZO impaired reproductive endpoints in F0 zebrafish following 120 days exposure. The average
283
number of eggs per female and average fertilization rate were significantly decreased in the 20.0 µg/L
284
AZO group (Table S7). Compared to F1 embryos from the non-exposed F0 parents, a significant increase
285
in mortality and a significant decrease in hatching rate at 96 hpf were observed in F1 offspring from the
286
exposed parents of 20.0 µg/L AZO group, even when F1 embryos were reared in AZO-free water (Table
287
S8). Moreover, F1 larvae collecting from F0 group of 20.0 µg/L AZO displayed pericardial edema, yolk
288
sac edema and spinal deformation (Figure S11), with a significant increase in malformation rate at 96
289
hpf (Table S8).
290
DISCUSSION
291
In this study, 20.0 µg/L AZO reduced the ratios of female and male zebrafish, and increased the
292
proportion of undifferentiated and intersex zebrafish at 42 dpf. Decreased male ratio and increased
293
intersex ratio were also observed in zebrafish exposed to 20.0 µg/L AZO at 60 dpf. However, no 15
ACS Paragon Plus Environment
Environmental Science & Technology
294
undifferentiated or intersex phenotype was detected in zebrafish exposed to 20.0 µg/L AZO at 120 dpf,
295
nor was sex ratio altered. These results suggest that 20.0 µg/L AZO retards sexual differentiation at 42
296
and 60 dpf, but may not impact the sex ratio at 120 dpf. Delayed sexual differentiation would induce
297
undifferentiated or the intersex phenotype, which may lead to the alteration of sex ratio in zebrafish,36, 37
298
but these impacts may be less pronounced or not present at later stages of life. For example, exposure to
299
0.15 mM methimazole, starting from 3 dpf to 33 dpf, delayed sexual differentiation and induced female-
300
biased population in zebrafish at 45 dpf, but had no effect on sex ratio after zebrafish were maintained
301
in clean water until 60 dpf.38 The pesticides prochloraz, vinclozolin, and monocrotophos are documented
302
to affect sexual development and sex ratio in zebrafish.39-42 For example, waterborne exposure to 400
303
μg/L vinclozolin, an anti-androgenic fungicide, starting from 21 to 35 dpf permanently shifted sex ratio
304
towards females (55−65%) and impaired gonadal testicular maturation in male zebrafish.39 The
305
organophosphorus pesticide monocrotophos with estrogenic activity increased the proportion of females
306
(71%) at 0.100 mg/L after exposure from fertilization to 40 dpf.40 The fungicide prochloraz, an aromatase
307
inhibitor and androgen receptor antagonist, caused male-bias (77%) population after 60 days exposure
308
to 202 μg/L.40, 42
309
Exposure to 20.0 µg/L AZO increased the percentage of immature oocytes at 42, 60 and 120 dpf,
310
and decreased the percentage of mature oocytes in female zebrafish at 60 and 120 dpf. Likewise, 20.0
311
µg/L AZO increased the percentage of immature spermatocytes and decreased the percentage of mature
312
spermatocytes in male zebrafish at 60 and 120 dpf. These results indicate that 20.0 µg/L AZO retarded
313
gonadal development in female and male zebrafish. The decreased GSI in females and males at 60 and
314
120 dpf further support this hypothesis. These data are in good agreement with that of our previous study,
315
which reported that gonadal development was delayed in female zebrafish exposed to 200 µg/L AZO, 16
ACS Paragon Plus Environment
Page 16 of 27
Page 17 of 27
Environmental Science & Technology
316
and in male zebrafish treated with 20 and 200 µg/L AZO following a 21-day exposure.14 Thus, longer
317
exposures to lower doses of AZO also appear to delay gonadal development in both females and males.
318
Here, we found that 20.0 µg/L AZO decreased E2 and E2/T ratio, and increased T in females at 60
319
and 120 dpf. Conversely, an increase in E2 and E2/T ratio, and a decrease in T was observed in males
320
exposed to 20.0 µg/L AZO at 42, 60 and 120 dpf. Our previous study reported that AZO altered E2, T
321
and Vtg levels in females after a 21 day-exposure to 200 µg/L, and in males exposed to 20 and 200
322
µg/L.14 Sex steroid hormones and genes controlling sexual determination regulate gonadal differentiation
323
and sex ratio.43 We propose that an unbalanced E2 and T is a major underlying reason for altered sexual
324
differentiation and retarded gonadal growth in zebrafish treated with 20.0 µg/L AZO as expression levels
325
of foxl2, brca2, sox9b, amh, sox9a and dmrt1 remained unchanged in females and males at 42 and 60
326
dpf. However, there are likely multiple mechanisms that are associated with gonadal differentiation and
327
sex ratios. For instance, norgestrel was reported to affect sexual differentiation and induce male-biased
328
population in zebrafish exposed to 34 and 77 ng/L starting from 20 dpf to 60 dpf. This effect was
329
accompanied by T suppression, increased dmrt1 expression, and decreased figa transcription.21
330
AZO at 20.0 µg/L down-regulated cyp19a and cyp19b, as well as up-regulated hsd3b and hsd17b
331
in female zebrafish at 42, 60 and 120 dpf, whereas up-regulation was observed for cyp19a, cyp19b, cyp17,
332
hsd3b and hsd17b in male zebrafish treated with 20.0 µg/L AZO at 42, 60 and 120 dpf. AZO at 2.0 µg/L
333
was found to down-regulate cyp19a, and up-regulate cyp17, hsd3b and hsd17b in females at 120 dpf,
334
while exposure to 2.0 µg/L AZO increased cyp19a, cyp19b, cyp17, hsd3b and hsd17b in males at 120
335
dpf. The synthesis of sex steroid hormones is controlled by steroidogenic enzymes of the HPG axis, and
336
transcriptional changes in these key steroidogenic enzymes (e.g. cyp19, cyp17, hsd3b and hsd17b) disrupt
337
the balance between E2 and T.28, 44 Taken together, we hypothesize that the altered expression of genes 17
ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 27
338
encoding key steroidogenic enzymes affects sex hormones in females exposed to 20.0 µg/L AZO at 60
339
and 120 dpf, and in males treated with 20.0 µg/L AZO at 42, 60 and 120 dpf. This hypothesis is also
340
supported by our previous study demonstrating that the dysregulation of cyp19a, cyp19b, cyp17, hsd3b
341
and hsd17b were associated with altered E2 and T in females and males following a 21-day exposure to
342
AZO.14 In another study, zebrafish larvae displayed an increase in E2, which was associated with an up-
343
regulation of cyp19b and hsd17b after 24, 48 and 72 h exposure to 1 and 100 µg/L AZO.12 However,
344
further studies should be conducted to explore if transcriptional alterations in HPG axis caused by AZO
345
actually led to changes in protein levels and enzyme activity in steroidogenic synthesis pathway.
346
We found that AZO appeared to have more profound effects in males than females and this is
347
consistent to our previous study.14 Difenoconazole, a triazole fungicide, has been shown to exert sex-
348
specific effects in adult zebrafish.45 However, the mode of action of AZO is not yet defined in fish and
349
the mechanisms of sex-dependent sensitivity are complicated. We did not specifically set out to test this
350
hypothesis so we are cautious in our interpretation about severity or sex-specific susceptibility. A novel
351
mechanism related to sex-dependent sensitivity in zebrafish would be interesting to pursue.
352
In this study, exposure to 20.0 µg/L AZO reduced Vtg protein, and down-regulated vtg1 and vtg2
353
in female zebrafish at 42, 60, and 120 dpf. Conversely, Vtg protein and vtg mRNA were increased in
354
male zebrafish treated with 20.0 µg/L AZO at 42, 60 and 120 dpf, and 2.0 µg/L AZO at 120 dpf. Vtg is
355
vital for ovarian development, oocyte maturation and yolk biosynthesis,46,
356
reduced Vtg in females is related to the delayed ovarian development. Vtg in males can be induced by
357
estrogenic or anti-androgenic chemicals, which would inhibit testicular development and maturation.48
358
Our recent study reported that 200 µg/L AZO inhibited ovarian development, decreased Vtg, and down-
359
regulated vtg1 and vtg2 in female zebrafish following 21 days exposure, whereas AZO at 20 and 200 18
ACS Paragon Plus Environment
47
and it is proposed that
Page 19 of 27
Environmental Science & Technology
360
µg/L retarded testicular development, increased Vtg, and up-regulated vtg1 and vtg2 in male zebrafish.14
361
Noteworthy was that parental exposure to 200 µg/L AZO up-regulated vtg1 and vtg2 in F1 embryos,
362
even when embryos were reared in clear water.13 Taken together, AZO appears to alter Vtg protein and
363
vtg expression differently in the sexes, and this response is rather consistent over the duration of exposure
364
(> 21d).
365
Long-term exposure to 20.0 µg/L AZO from 2–4 hpf to 120 dpf decreased egg production and
366
fertilization rate in zebrafish. We also observed reduced survival and hatching rate, morphological
367
defects, and increased malformation rate in F1 offspring from exposed parents (20.0 µg/L AZO), despite
368
that F1 embryos were reared in clean water. Reduced fecundity may not be surprising as the number of
369
mature oocytes is influenced by circulating sex steroid levels and Vtg. These, in turn, are directly related
370
to steroidogenic enzyme expression in the gonad and vtg expression in the liver.49 Thus, we propose that
371
impaired reproduction in F0 individuals is a result of delayed gonadal development caused by altered sex
372
steroids and Vtg. The exposure of F0 zebrafish to 20.0 µg/L AZO resulted in lethality, hatching delay,
373
morphological abnormality, and increased malformation rate in F1 embryos. These results were in good
374
agreement with our previous findings that impaired reproduction was induced in adult zebrafish after a
375
21-day exposure to a higher dose of 200 µg/L AZO.14 Mortality, hatching delay, morphological
376
abnormity, elevated malformation rate, and disrupted gene expression in the HPG axis were detected in
377
F1 embryos from exposed parents of 200 µg/L AZO, with or without continuous exposure to this
378
fungicide; adverse effects were also detected after F1 embryos from exposed parents of 20 µg/L AZO
379
were continuously exposed to the same dose as their parents.13 Taken together, AZO appears to affect
380
the reproduction of zebrafish with both short-term exposures at higher concentrations and long-term
381
exposures at lower, environmentally relevant concentrations. 19
ACS Paragon Plus Environment
Environmental Science & Technology
382
In conclusion, this study showed that AZO adversely influenced a number of endpoints related to
383
sexual development and reproduction in zebrafish. Transcriptional responses may be the underlying
384
molecular mechanisms in terms of disrupted sexual development and impaired reproduction in F0
385
generation exposed to 20.0 µg/L AZO, and changes in mRNA levels were detected in females and males
386
exposed to 2.0 µg/L at 120 dpf. Parental exposure to 20.0 µg/L AZO for 120 days affected survival,
387
hatching and morphology in F1 embryos that were not exposed to AZO. These data provide new
388
understanding as to the effects of AZO on sexual development and reproduction in zebrafish after long-
389
term exposure to environmentally relevant concentrations.
390
ASSOCIATED CONTENT
391
Supporting Information (SI)
392
Descriptions of feeding procedure, sex steroid hormones and Vtg measurement, and analytical
393
chemistry are presented in Page S2−S9. Experimental design (Figure S1), sample size (Table S1),
394
primer sequences (Table S2), the stability of β-actin (Figure S2), melt curves of tested genes (Figure
395
S3−S6), concentrations of AZO in exposure solution (Table S3), growth and GSI (Table S4 and S5), Sex
396
steroid hormone and Vtg (Table S6 and S7), gene expression data (Figure S7−S10 and Table S8),
397
reproductive capacity in F0 generation (Table S9), developmental success in F1 embryos (Table S10),
398
and representative pictures in F1 offspring (Figure S11) are also presented in Supporting information.
399
AUTHOR INFORMATION
400
Corresponding author
401
*E-mail:
[email protected]; Tel: +86 (0)10 62733924
402
Notes
403
The authors have no conflict of interest to declare. 20
ACS Paragon Plus Environment
Page 20 of 27
Page 21 of 27
Environmental Science & Technology
404
ACKNOWLEDGMENTS
405
This work was financially supported by the National Key R&D Program of China (Grant No.
406
2017YFD0200504).
407
REFERENCES
408
(1) PPDB, Pesticide Properties DataBase for Azoxystrobin. University of Hertfordshire.
409
https://sitem.herts.ac.uk/aeru/ppdb/en/Reports/54.htm.
410
(2) PhillipsMcDougall, AgriService: Products section-Azoxystrobin, 2016 Market. 2016.
411
(3) Rodrigues, E. T.; Lopes, I.; Pardal, M. A., Occurrence, fate and effects of azoxystrobin in aquatic
412
ecosystems: a review. Environ Int. 2013, 53, 18-28.
413
(4) Wang, Q.; Zhong, W.; Huang, S., Determination of azoxystrobin residues in surface water by HPLC
414
with solid-phase extraction. Anhui Med and Pharm. 2009, 13, (6), 611-612 (in Chinese).
415
(5) Smalling, K. L.; Orlando, J., Occurrence of Pesticides in Surface Water from Three Central
416
California Coastal Watersheds, 2008–09. US Geological Survey Data Series 600 2011, US Department
417
of the Interior, Washington, DC, 70.
418
(6) Berenzen, N.; Lentzen-Godding, A.; Probst, M.; Schulz, H.; Schulz, R.; Liess, M., A comparison
419
of predicted and measured levels of runoff-related pesticide concentrations in small lowland streams on
420
a landscape level. Chemosphere. 2005, 58, (5), 683-91.
421
(7) Battaglin, W. A.; Sandstrom, M. W.; Kuivila, K. M.; Kolpin, D. W.; Meyer, M. T., Occurrence of
422
Azoxystrobin, Propiconazole, and Selected Other Fungicides in US Streams, 2005–2006. Water Air Soil
423
Pollut. 2011, 218, (1-4), 307-322.
424
(8) Schriever, C. A.; von der Ohe, P. C.; Liess, M., Estimating pesticide runoff in small streams.
425
Chemosphere. 2007, 68, (11), 2161-71. 21
ACS Paragon Plus Environment
2018,
Environmental Science & Technology
426
(9) Jorgensen, L. F.; Kjaer, J.; Olsen, P.; Rosenbom, A. E., Leaching of azoxystrobin and its
427
degradation product R234886 from Danish agricultural field sites. Chemosphere 2012, 88, (5), 554-62.
428
(10) Kunz, J. L.; Ingersoll, C. G.; Smalling, K. L.; Elskus, A. A.; Kuivila, K. M., Chronic toxicity of
429
azoxystrobin to freshwater amphipods, midges, cladocerans, and mussels in water-only exposures.
430
Environ Toxicol Chem. 2017, 36, (9), 2308-2315.
431
(11) Han, Y.; Liu, T.; Wang, J.; Wang, J.; Zhang, C.; Zhu, L., Genotoxicity and oxidative stress induced
432
by the fungicide azoxystrobin in zebrafish (Danio rerio) livers. Pestic Biochem Physiol. 2016, 133, 13-
433
19.
434
(12) Jiang, J.; Shi, Y.; Yu, R.; Chen, L.; Zhao, X., Biological response of zebrafish after short-term
435
exposure to azoxystrobin. Chemosphere. 2018, 202, 56-64.
436
(13) Cao, F.; Li, H.; Zhao, F.; Wu, P.; Qian, L.; Huang, L.; Pang, S.; Martyniuk, C. J.; Qiu, L., Parental
437
exposure to azoxystrobin causes developmental effects and disrupts gene expression in F1 embryonic
438
zebrafish (Danio rerio). Sci Total Environ. 2018, 646, 595-605.
439
(14) Cao, F.; Zhu, L.; Li, H.; Yu, S.; Wang, C.; Qiu, L., Reproductive toxicity of azoxystrobin to adult
440
zebrafish (Danio rerio). Environ Pollut. 2016, 219, 1109-1121.
441
(15) EPA, Azoxystrobin pesticide fact sheet. http://www.epa.gov/opprd001/factsheets/azoxystr.pdf1997.
442
1997, [accessed on 18th April 2012].
443
(16) EFSA, Peer review report to the conclusion regarding the peer review of the pesticide risk
444
assessment of the active substance azoxystrobin. EFSA. 2010, 8, (4), 1542.
445
(17) Bartlett, D. W.; Clough, J. M.; Godfrey, C. R. A.; Godwin, J. R.; Hall, A. A.; Heaney, S. P.; Maund,
446
S. J., Understanding the strobilurin fungicides. Pestic. Outlook. 2001, 12, (4), 143-148.
447
18. Santos, D.; Luzio, A.; Coimbra, A. M., Zebrafish sex differentiation and gonad development: A 22
ACS Paragon Plus Environment
Page 22 of 27
Page 23 of 27
Environmental Science & Technology
448
review on the impact of environmental factors. Aquat Toxicol. 2017, 191, 141-163.
449
(19) Luzio, A.; Monteiro, S. M.; Garcia-Santos, S.; Rocha, E.; Fontainhas-Fernandes, A. A.; Coimbra,
450
A. M., Zebrafish sex differentiation and gonad development after exposure to 17alpha-ethinylestradiol,
451
fadrozole and their binary mixture: A stereological study. Aquat Toxicol. 2015, 166, 83-95.
452
(20) Wang, X. G.; Bartfai, R.; Sleptsova-Freidrich, I.; Orban, L., The timing and extent of 'juvenile
453
ovary' phase are highly variable during zebrafish testis differentiation. J Fish Biol. 2007, 70, (sa), 33-44.
454
(21) Liang, Y. Q.; Huang, G. Y.; Liu, S. S.; Zhao, J. L.; Yang, Y. Y.; Chen, X. W.; Tian, F.; Jiang, Y.
455
X.; Ying, G. G., Long-term exposure to environmentally relevant concentrations of progesterone and
456
norgestrel affects sex differentiation in zebrafish (Danio rerio). Aquat Toxicol. 2015, 160, 172-9.
457
(22) Siegfried, K. R., In search of determinants: gene expression during gonadal sex differentiation. J
458
Fish Biol. 2010, 76, (8), 1879-902.
459
(23) Shive, H. R.; West, R. R.; Embree, L. J.; Azuma, M.; Sood, R.; Liu, P.; Hickstein, D. D., brca2 in
460
zebrafish ovarian development, spermatogenesis, and tumorigenesis. PNAS. 2010, 107, (45), 19350-
461
19355.
462
(24) Shang, E. H.; Yu, R. K.; Wu, R. S. S., Hypoxia Affects Sex Differentiation and Development,
463
Leading to a Male-Dominated Population in Zebrafish (Danio rerio). Environ Sci Technol. 2006, 40,
464
3118-3122.
465
(25) Luzio, A.; Monteiro, S. M.; Rocha, E.; Fontainhas-Fernandes, A. A.; Coimbra, A. M., Development
466
and recovery of histopathological alterations in the gonads of zebrafish (Danio rerio) after single and
467
combined exposure to endocrine disruptors (17alpha-ethinylestradiol and fadrozole). Aquat Toxicol.
468
2016, 175, 90-105.
469
(26) Bondesson, M.; Hao, R.; Lin, C. Y.; Williams, C.; Gustafsson, J. A., Estrogen receptor signaling 23
ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 27
470
during vertebrate development. Biochim Biophys Acta. 2015, 1849, (2), 142-51.
471
(27) Olsson, P. E.; Berg, A. H.; von Hofsten, J.; Grahn, B.; Hellqvist, A.; Larsson, A.; Karlsson, J.;
472
Modig, C.; Borg, B.; Thomas, P., Molecular cloning and characterization of a nuclear androgen receptor
473
activated by 11-ketotestosterone. Reprod Biol Endocrinol. 2005, 3, 37.
474
(28) Ma, Y.; Han, J.; Guo, Y.; Lam, P. K.; Wu, R. S.; Giesy, J. P.; Zhang, X.; Zhou, B., Disruption of
475
endocrine function in in vitro H295R cell-based and in in vivo assay in zebrafish by 2,4-dichlorophenol.
476
Aquat Toxicol. 2012, 106-107, 173-81.
477
(29) Hilscherova, K.; Jones, P. D.; Gracia, T.; Newsted, J. L.; Zhang, X.; Sanderson, J. T.; Yu, R. M.;
478
Wu, R. S.; Giesy, J. P., Assessment of the effects of chemicals on the expression of ten steroidogenic
479
genes in the H295R cell line using real-time PCR. Toxicol Sci 2004, 81, (1), 78-89.
480
(30) Tyler, C. R.; Jobling, S.; Sumpter, J. P., Endocrine disruption in wildlife: a critical review of the
481
evidence. Crit Rev Toxicol. 1998, 28, (4), 319-61.
482
(31) ISO, Water Quality-Determination of the Acute Lethal Toxicity of Substances to a Freshwater Fish
483
[Brachydanio rerio Hamiltone-Buchanan (Teleostei, Cyprinidae)]. In: Part 3: Flow-through Method
484
1996.
485
(32) OECD, Test no 123 Guidance Document on the Diagnosis of Endocrine-related Histopathology in
486
Fish
487
ENV/JM/MONO(2010)14.
488
(33) Ji, K.; Liu, X.; Lee, S.; Kang, S.; Kho, Y.; Giesy, J. P.; Choi, K., Effects of non-steroidal anti-
489
inflammatory drugs on hormones and genes of the hypothalamic-pituitary-gonad axis, and reproduction
490
of zebrafish. J Hazard Mater. 2013, 254-255, 242-51.
491
(34) Livak, K. J.; Schmittgen, T. D., Analysis of relative gene expression data using real-time
Gonads.
Organization
for
Economic
Cooperation
24
ACS Paragon Plus Environment
and
Development
2010,
In:
Page 25 of 27
Environmental Science & Technology
492
quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25, (4), 402-8.
493
(35) OECD, OECD Guidelines for the Testing of Chemicals. In: Section 2: Effects on Biotic Systems
494
Test No. 229: Fish Short Term Reproduction Assay. Organization for Economic Cooperation and
495
Development 2012, Paris (France).
496
(36) Baumann, L.; Knorr, S.; Keiter, S.; Rehberger, K.; Volz, S.; Schiller, V.; Fenske, M.; Holbech, H.;
497
Segner, H.; Braunbeck, T., Reversibility of endocrine disruption in zebrafish (Danio rerio) after
498
discontinued exposure to the estrogen 17alpha-ethinylestradiol. Toxicol Appl Pharmacol. 2014, 278, (3),
499
230-7.
500
(37) Lin, L. L.; Janz, D. M., Effects of binary mixtures of xenoestrogens on gonadal development and
501
reproduction in zebrafish. Aquat Toxicol. 2006, 80, (4), 382-95.
502
(38) Sharma, P.; Patino, R., Regulation of gonadal sex ratios and pubertal development by the thyroid
503
endocrine system in zebrafish (Danio rerio). Gen Comp Endocrinol. 2013, 184, 111-9.
504
(39) Lor, Y.; Revak, A.; Weigand, J.; Hicks, E.; Howard, D. R.; King-Heiden, T. C., Juvenile exposure
505
to vinclozolin shifts sex ratios and impairs reproductive capacity of zebrafish. Reprod Toxicol. 2015, 58,
506
111-8.
507
(40) Zhang, X.; Gao, L.; Yang, K.; Tian, H.; Wang, W.; Ru, S., Monocrotophos pesticide modulates the
508
expression of sexual differentiation genes and causes phenotypic feminization in zebrafish (Danio rerio).
509
Comp Biochem Physiol C Toxicol Pharmacol. 2013, 157, (1), 33-40.
510
(41) Thorpe, K. L.; Pereira, M. L.; Schiffer, H.; Burkhardt-Holm, P.; Weber, K.; Wheeler, J. R., Mode
511
of sexual differentiation and its influence on the relative sensitivity of the fathead minnow and zebrafish
512
in the fish sexual development test. Aquat Toxicol. 2011, 105, (3-4), 412-20.
513
(42) Kinnberg, K.; Holbech, H.; Petersen, G. I.; Bjerregaard, P., Effects of the fungicide prochloraz on 25
ACS Paragon Plus Environment
Environmental Science & Technology
514
the sexual development of zebrafish (Danio rerio). Comp Biochem Physiol C Toxicol Pharmacol. 2007,
515
145, (2), 165-70.
516
(43) Trukhina, A. V.; Lukina, N. A.; Wackerow-Kouzova, N. D.; Smirnov, A. F., The variety of
517
vertebrate mechanisms of sex determination. Biomed Res Int 2013, 2013, 587460.
518
44. Fenske, M.; Segner, H., Aromatase modulation alters gonadal differentiation in developing
519
zebrafish (Danio rerio). Aquat Toxicol. 2004, 67, (2), 105-126.
520
(45) Mu, X.; Wang, K.; Chai, T.; Zhu, L.; Yang, Y.; Zhang, J.; Pang, S.; Wang, C.; Li, X., Sex specific
521
response in cholesterol level in zebrafish (Danio rerio) after long-term exposure of difenoconazole.
522
Environ Pollut。 2015, 197, 278-86.
523
(46) Yuan, H. X.; Xu, X.; Sima, Y. H.; Xu, S. Q., Reproductive toxicity effects of 4-nonylphenol with
524
known endocrine disrupting effects and induction of vitellogenin gene expression in silkworm, Bombyx
525
mori. Chemosphere. 2013, 93, (2), 263-8.
526
(47) Yu, L.; Liu, C.; Chen, Q.; Zhou, B., Endocrine disruption and reproduction impairment in zebrafish
527
after long-term exposure to DE-71. Environ Toxicol Chem. 2014, 33, (6), 1354-62.
528
48. Barucca, M.; Canapa, A.; Olmo, E.; Regoli, F., Analysis of vitellogenin gene induction as a valuable
529
biomarker of estrogenic exposure in various Mediterranean fish species. Environ Res 2006, 101, (1), 68-
530
73.
531
(49) Villeneuve, D. L.; Blake, L. S.; Brodin, J. D.; Cavallin, J. E.; Durhan, E. J.; Jensen, K. M.; Kahl,
532
M. D.; Makynen, E. A.; Martinovic, D.; Mueller, N. D.; Ankley, G. T., Effects of a 3beta-hydroxysteroid
533
dehydrogenase inhibitor, trilostane, on the fathead minnow reproductive axis. Toxicol Sci. 2008, 104,
534
(1), 113-23.
535 26
ACS Paragon Plus Environment
Page 26 of 27
Page 27 of 27
Environmental Science & Technology
Table of Contents Art of this manuscript
ACS Paragon Plus Environment