Long-term exposure to environmental concentrations of azoxystrobin

Jan 7, 2019 - Long-term exposure to environmental concentrations of azoxystrobin delays sexual development and alters reproduction in zebrafish (Danio...
0 downloads 0 Views 2MB Size
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