Environmentally relevant concentrations of carbamazepine caused

Large amounts of data have. 36 demonstrated that the ..... The agreement between the. 272. PCR-based expression of genes in the liver and the microarr...
2 downloads 11 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Environmentally relevant concentrations of carbamazepine caused endocrine-disrupting effects on non-target organisms, Chinese rare minnows (Gobiocypris rarus) Saihong Yan, Miao Wang, Jinmiao Zha, Lifei Zhu, Wei Li, Qian Luo, Jing Sun, and Zijian Wang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b06476 • Publication Date (Web): 18 Dec 2017 Downloaded from http://pubs.acs.org on December 19, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35

Environmental Science & Technology

1

Environmentally relevant concentrations of carbamazepine caused

2

endocrine-disrupting effects on non-target organisms, Chinese rare

3

minnows (Gobiocypris rarus)

4

Saihong Yan,§, †,∥Miao Wang,§, † Jinmiao Zha,§, †, * Lifei Zhu, Wei Li,‡ Qian Luo,# Jing Sun, §

5

Zijian Wang§



6

§

7

Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China

8



9

Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China

10 11





State Key Laboratory of Environmental Aquatic Chemistry, Research Center for

Key Laboratory of Drinking Water Science and Technology, Research Center for

University of Chinese Academy of Sciences, Beijing 100049, China Beijing Fisheries Research Institute, Beijing 100068, China

12



13

Sciences, Chinese Academy of Sciences, Beijing 100085, China

14

#

15

China

Key Laboratory of Environmental Biotechnology, Research Center for Eco-Environmental

Shenzhen Institutes of Advanced Technology, Chinese Academy of Science, Shenzhen 518055,

ACS Paragon Plus Environment

Environmental Science & Technology

16

ABSTRACT: In the present study, Chinese rare minnows (Gobiocypris rarus) were

17

exposed to 1, 10, and 100 µg/L carbamazepine (CBZ) under flow-through conditions

18

for 28 d. A hepatic-specific custom microarray identified 111 and 71 differentially

19

expressed genes in the livers of females and males exposed to 100 µg/L CBZ,

20

respectively (ratio ≥ 2, p ≤ 0.05). The levels of five differentially expressed genes

21

associated with the hypothalamic-pituitary-gonadal (HPG) axis were quantified by

22

qPCR, and the results indicated the feasibility of screening endocrine-disrupting

23

chemicals using a custom microarray. The mRNA levels of genes related to the HPG

24

axis differed significantly in different organs of Chinese rare minnow (p30,600 µg/L)30 and severe spinal deformations with upward tail

62

bending (29,625 µg/L)31 in zebrafish. Although CBZ can cause mortality and

63

developmental impairments, its potential adverse ecological effects on aquatic

64

organisms, particularly endocrine disruption in fish, have not been documented.

65

Chinese rare minnow (Gobiocypris rarus) is widely used as a model to assess the

66

potential effects of chemicals because of its small size, ease of culture, short life cycle

67

and high spawning rate.32-34 Therefore, the purpose of this study was to investigate the

68

endocrine-disrupting effects and related mechanisms of environmentally relevant

69

concentrations of CBZ (1, 10, and 100 µg/L) on Chinese rare minnows for 28 d by

70

measuring multiple toxicological endpoints. First, the transcriptional profiles of livers

71

from CBZ-treated (100 µg/L) fish were determined using a hepatic-specific custom

72

microarray,

73

hypothalamic-pituitary-gonadal (HPG) axis were verified by RT-PCR. Second, the

74

gonadosomatic index (GSI), histopathology, and the plasma levels of vitellogenin

75

(VTG) and steroid hormones were investigated. Finally, conclusions regarding the

76

toxic effects and possible mechanism of action of CBZ in aquatic organisms were

77

drawn based on the assessed toxicological endpoints.

and

the

expression

levels

of

genes

ACS Paragon Plus Environment

associated

with

the

Page 5 of 35

Environmental Science & Technology

78

■ MATERIALS AND METHODS

79

Chemicals. CBZ (CAS: 298-46-4, purity >97%) was bought from Sigma

80

Chemical Co. (St. Louis, Missouri, USA), and all other chemicals were purchased

81

from Sigma Chemical Co. (St. Louis, Missouri, USA). Stock solutions were prepared

82

in DMSO and stored in brown bottles. The final concentration of DMSO was less

83

than 0.01%.

84

Test Fish and Culture Conditions. The brood stock of rare minnows was

85

kept in an 18-L aquaria under 10 L/h flow-through conditions with dechlorinated tap

86

water (pH 7.2-7.6, hardness of 44.0-61.0 mg/L CaCO3, and temperature of 25±1 °C)

87

and a photoperiod of 16:8 h (light:dark).32-34 The brood stock was used for the

88

evaluation of chemicals in our laboratory over 10 years. The minnows were fed newly

89

hatched brine shrimp (Artemia nauplii) twice daily and granule food (TetraMin, Tetra

90

Werke, Melle, Germany) once daily. Food residues and feces were wiped off daily

91

with a sucking pipe, and the aquariums were cleaned weekly.

92

Experimental Design. All of the fish used in the present study are offspring of

93

the same pair of brood stock. Approximately 10-month-old sexually mature fish were

94

used in the experiment. In addition, the average bodyweight and length of the female

95

fish that were tested were 0.71±0.22 g and 39.26±3.19 mm, respectively, whereas

96

those of male fish were 0.63±0.10 g and 38.94±2.97 mm, respectively. The fish were

97

randomly divided into 12 groups of 30 fish (with a sex ratio of 1:1) and were

98

acclimated for two weeks under the same conditions as the brood stock before the

99

experiment. The fish were exposed to 1, 10, and 100 µg/L CBZ for 28 d in 18-L

ACS Paragon Plus Environment

Environmental Science & Technology

100

aquaria (under 10 L/h flow-through conditions). An independent dose of 60 mL/h

101

CBZ in the stock solutions was administered by infusion from a 4-L glass aspirator,

102

and this dose was renewed every second day. The water flow and dosing rates were

103

checked twice daily. The water and CBZ were mixed and then flowed to each tank at

104

a rate of 10 L/h. Three replicates of the control and treated groups were included, and

105

the fish were sacrificed after exposure for 28 d. The body and gonad weights of the

106

males and females were recorded. Blood from fish of the same gender was collected

107

with a heparinized microcapillary, pooled per treatment group and then immediately

108

centrifuged (8000 g, 10 min, 4 °C) to obtain plasma, and the plasma levels of steroid

109

hormones and VTG were determined. In addition, the gonad and brain from each fish

110

were carefully excised, and the liver tissue from each fish was sampled and cut into

111

halves: one-half of the liver was used for microarray, and the other half was used for

112

gene expression analysis by real-time PCR. All tissues were flash-frozen in liquid

113

nitrogen and stored at -80 °C.

114

Quantitative of Carbamazepine Concentration in Water. Water samples

115

were measured using an ultra-performance liquid chromatography (UPLC) system

116

connected to a Micromass Quattro Premier XE tandem quadruple mass spectrometer

117

with an electrospray ion source (Waters Corporation, Milford, MA, USA) as

118

described by Sun et al. (2015).35 The detection limit and linearity range of the method

119

were 0.02 ng/L and 0.05-50 µg/L, respectively. Water samples were collected every

120

week, filtered through a 0.2-µm GHP membrane filter, and diluted five-fold with

121

distilled water. The measured concentrations (mean ± standard deviation; and %

ACS Paragon Plus Environment

Page 6 of 35

Page 7 of 35

Environmental Science & Technology

122

analyzed/nominal) of CBZ in the flow-through system during the exposure period

123

were 0.91±0.03 (91%), 8.82±0.64 (88%) and 82.9 ± 4.12 (83%) µg/L, respectively,

124

which were close to the nominal concentrations. Therefore, the nominal CBZ values

125

are used in the rest of the manuscript.

126

Hepatic-specific Custom Microarray. Five hundred forty-eight clones

127

corresponding to unique genes were selected for PCR amplification of the inserts

128

(Supporting Information (SI) Table S1). The 548 clones contained 33 genes related to

129

the HPG axis, 10 genes related to the hypothalamic-pituitary-thyroid (HPT) axis, 6

130

genes related to the hypothalamic-pituitary-adrenal (HPA) axis, and 500 genes related

131

to other axes. PCR was performed using 2.5 units of ExTaq (Takara, Dalian, China) in

132

a volume of 50 µL, which included 1× PCR buffer, 250 µM of each dNTP, and 0.5

133

µM of each primer (T3 and T7). The PCR amplification program consisted of three

134

steps: (1) 5 min at 94 °C, (2) 18 cycles of 40 s at 94 °C, 40 s at 58 °C, and 1 min at

135

72 °C, and (3) 5 min at 72 °C. The PCR products were precipitated with the addition

136

of ethanol, dissolved in EasyArrayTM spotting solution (CapitalBio, Beijing, China),

137

and printed in triplicate on a PolymerSlide (CapitalBio, Beijing, China) using a

138

SmartArrayTM microarrayer (CapitalBio, Beijing, China). Eight sequences derived

139

from intergenic yeast genomic regions that exhibited no significant homologies with

140

any existing teleost sequences in GenBank were used as external controls. These

141

sequences were amplified by PCR and cloned into the pSP64 Poly(A) plasmid

142

(Promega, Madison, WI, USA) to produce Poly(A)-RNA following in vitro

143

transcription, and the RNA samples were spiked with varying amounts of

ACS Paragon Plus Environment

Environmental Science & Technology

144

Poly(A)-RNA.

145

RNA samples from the group exposed to 100 µg/L CBZ were selected for

146

microarray analysis. A total of 12 arrays were used in the experiment. A total of six

147

pooled RNA samples (three from males and three from females) were obtained by

148

pooling each set of two individual RNA samples and analyzed separately. Samples

149

from the CBZ-treated fish were co-hybridized with a common control (six individual

150

RNA samples from the male or female control group were pooled) to obtain the

151

common reference design. A dye-swap replicate was used for each sample to measure

152

reproducibility and reduce systematic bias.

153

The Eberwine linear RNA amplification method was used to yield cDNA labeled

154

with a fluorescent dye (Cy5- and Cy3-dCTP), and subsequent enzymatic reactions

155

were performed as previously described. 36 Briefly, double-stranded cDNA containing

156

a T7 RNA polymerase promoter sequence was synthesized with 5 µg of total RNA

157

using the Reverse Transcription System, RNAse H, DNA polymerase I, and T4 DNA

158

polymerase according to the manufacturer’s recommended protocol (Promega,

159

Madison, WI, USA). The double-stranded cDNA was then subjected to in vitro

160

transcription reactions at 37 °C for 3 h using the T7 RiboMAX Express large-scale

161

RNA production system (Promega, Madison, WI, USA). After reverse transcription,

162

the Klenow enzyme labeling strategy was used. Labeled controls and test samples

163

were quantitatively adjusted based on the efficiency of Cy5-dCTP or Cy3-dCTP

164

(Amersham Pharmacia Biotech, Piscataway, NJ, USA) incorporation. The samples

165

were then dissolved in 80 µL of hybridization solution containing 50% formamide.

ACS Paragon Plus Environment

Page 8 of 35

Page 9 of 35

Environmental Science & Technology

166

The arrays were hybridized in a CapitalBio BioMixerTM II Hybridization Station

167

overnight at 42 °C and subjected to two consecutive washes with 0.2% SDS plus 2×

168

SSC for 5 min at 42 °C and 0.2× SSC for 5 min at room temperature.

169

The arrays were scanned with a confocal LuxScanTM scanner (CapitalBio,

170

Beijing, China), and the images were analyzed using LuxScanTM 3.0 software

171

(CapitalBio, Beijing, China). Faint spots with an intensity less than 400 units after

172

background subtraction from both channels (Cy3 and Cy5) were removed from

173

individual channel data. The data were subjected to space- and intensity-dependent

174

normalization based on a LOWESS method.37 Fold changes were calculated relative

175

to the controls. The fold differences in gene expression levels obtained after CBZ

176

exposure were calculated as the average fold changes from both the technical and

177

biological replicates. The two-tailed Student’s t-test method was used to statistically

178

analyze the data from the 12 independently tested arrays. Data from male and female

179

fish were analyzed separately. A p-value less than 0.05 combined with an average fold

180

change ≥ 2 was used to identify significantly differentially expressed genes.

181

Relative Quantification by Real-time PCR. Total RNA from different

182

tissues (brain, liver and gonad) was obtained using an SV Total RNA Isolation System

183

(Promega, Madison, WI, USA) based on the manufacturer’s recommended methods,

184

and three replicates of every treatment were included. cDNA was synthesized with

185

RNA according to Chen et al.38 Real-time PCR was performed using an ABI 7500

186

real-time quantitative PCR system (Life Technologies, Carlsbad, CA, USA). The

187

mRNA expression levels of various genes (vtg, erα, erβ1, erβ2, ar, gnrh2, gnrh1α,

ACS Paragon Plus Environment

Environmental Science & Technology

188

gnrh1β, cyp11α, cyp17, cyp19α, cyp19β, hsd3β7, and star) were detected in triplicate.

189

The β-actin gene of rare minnow was used as an endogenous control. The primer

190

sequences and accession numbers are listed in SI Table S2. The RT-PCR reaction

191

system consisted of the Brilliant II SYBR Green QPCR master mix with the forward

192

and reverse primers. The temperature program of the QPCR reaction was 95 °C for 10

193

min followed by 40 cycles of 95 °C for 30 s and 57 °C for 1 min and a final step of

194

72 °C for 30 s. All samples were analyzed in triplicate, and the average value from the

195

three replicates was used to indicate the mRNA expression levels of the genes. The

196

results were analyzed according to the delta-delta Ct method.39 A dissociation curve

197

was generated for each gene to assess the amplification of untargeted fragments. The

198

gene expression levels are presented as changes relative to the control group during

199

the same treatment period.

200 201

Gonadosomatic Index and Histological Analysis. The GSI of male and female fish from each treatment group was calculated as follows:

202

GSI (%) = (gonad weights (g)/body weights (g))×100

203

The histopathology of the gonad was analyzed according to Chen et al. (2016).38

204

The gonad samples were fixed in Bouin’s solution for 48 h and transferred to 70%,

205

80%, 90%, 95%, and 100% ethanol for dehydration. The samples were then dipped

206

into transparent reagent, incorporated into soft and hard wax and finally embedded in

207

paraffin wax. The wax was then cut into 3-4-mm-thick sections including gonad

208

samples, and the sections were finally dyed with hematoxylin and eosin (H&E).

209

Section photographs were obtained with a BX53 optical microscope (Olympus,

ACS Paragon Plus Environment

Page 10 of 35

Page 11 of 35

Environmental Science & Technology

210

Tokyo, Japan) and analyzed with the cellSens Standard imaging software (Olympus,

211

Tokyo, Japan). The numbers of mature and immature oocytes in the histological

212

sections were counted, and the number of primary oocytes in males was also recorded

213

to calculate the testis-ova rates.

214

Plasma Steroid Hormones and VTG. The plasma levels of 17β-estradiol

215

(E2) and 11-ketotestosterone (11-KT) were

detected using enzyme-linked

216

immunosorbent assay (ELISA) kits (Hufeng Biotechnology, Shanghai, China) with

217

assay sensitivities of 1.5 ng/L and 7 ng/L, respectively, according to the

218

manufacturer’s instructions. In addition, the plasma VTG level was measured using an

219

ELISA kit (Cusabio, Wuhan, Hubei, China) with an assay sensitivity of 50 ng/mL

220

according to the manufacturer’s instructions. The intra- and inter-assay coefficients of

221

variation were all less than 10%.

222

Statistical Analyses. All results are presented as the means ± standard errors

223

of the mean (SEMs), and figures were drawn with Origin Pro 8.0 (OriginLab,

224

Northampton, MA, USA). All data were analyzed for normality and homogeneity of

225

variance using the Kolmogorov-Smirnov and Levene’s tests, respectively, using SPSS

226

17.0 (SPSS, Chicago, IL, USA). Significant differences between the treatment and

227

control groups were estimated using one-way ANOVA (p < 0.05) followed by

228

Dunnett’s test for multiple comparisons.

229 230

■ RESULTS AND DISCUSSION

ACS Paragon Plus Environment

Environmental Science & Technology

231

Differential Gene mRNA Expression in Response to CBZ Exposure

232

Detected Using a Microarray. The transcriptional levels in the livers of males and

233

females treated with 100 µg/L CBZ were assessed using a custom hepatic-specific

234

cDNA microarray of Chinese rare minnow (Fig. 1A and B) that contains 548 unique

235

genes, which are mainly associated with the endocrine system (Fig. 1C). Villeneuve et

236

al.40 used a cDNA microarray to assess the effects of triclocarban on the ovarian

237

transcriptome, and their results are consistent with the levels of genes assessed by

238

QPCR and ultimately indicated that exposure to 5 µg/L triclocarban for 28 d

239

significantly decreases cumulative fecundity in fathead minnows without affecting the

240

endocrine system. Kwon et al.41 evaluated the metabolism-induced toxicity of

241

THLE-2 cells using a transfected enzyme and metabolism chip. Therefore, the results

242

indicate that endocrine-disrupting chemicals can be screened using a custom

243

microarray.

244

The mRNA expression levels of 146 unique genes were significantly different in

245

male and/or female fish (SI Table S3). A total of 111 and 71 altered genes (ratio ≧ 2, p

246

< 0.05) were detected in the female and male liver, respectively. Further analysis

247

showed that 36 genes were commonly changed in both males and females after

248

treatment with 100 µg/L CBZ (Fig. 1C), and these genes included 21 genes that were

249

significantly up-regulated, 12 genes that were significantly down-regulated (ratio ≧ 2,

250

p < 0.05), and three genes that showed differential expression between males and

251

females. A KEGG signal pathway analysis indicated that the altered genes were

252

associated with the HPG axis, lipid metabolism, calcium ion pathway, and other

ACS Paragon Plus Environment

Page 12 of 35

Page 13 of 35

Environmental Science & Technology

253

pathways. In the present study, HPG axis-related genes that were differentially

254

expressed in the liver were selected to verify the microarray results. The expression

255

levels of various genes in male (erα, erβ1 and vtg) and female livers (ar, erβ1 and vtg)

256

after exposure to 100 µg/L CBZ determined by real-time PCR agreed with the levels

257

obtained in the microarray analysis (Table 1). This agreement demonstrated that the

258

rare minnow cDNA microarray can be used for screening potential EDCs. The gene

259

expression levels suggest that CBZ might affect the liver by directly binding estrogen

260

receptors (ERs) to promote the mRNA levels of vtg in female fish. However, the

261

increased levels of the vtg gene in male Chinese rare minnow might be due to the

262

endocrine-disturbing effect of CBZ.42

263

The mRNA Levels of Genes Related to the HPG Axis in Rare Minnow

264

as Determined by Real-time PCR. The mRNA expression levels of genes related

265

to the HPG axis were significantly increased in all of the treatments (p < 0.05, Fig. 2

266

and SI Table S4), and the transcriptomic responses to CBZ associated with the HPG

267

axis of rare minnow are presented in Fig. 3.

268

The mRNA expression levels of erα in the female liver were significantly

269

increased after treatment with 1 and 10 µg/L CBZ, and the hsd3β7 levels were

270

increased in the females after treatment with 1 µg/L CBZ (p < 0.05; Figs. 2 and 3). No

271

significant differences were observed in the measured genes (erα and hsd3β7) in the

272

males after all of the treatments (p < 0.05; Figs. 2 and 3). The agreement between the

273

PCR-based expression of genes in the liver and the microarray data suggested that the

274

hepatic-specific custom microarray of Chinese rare minnow can be used to screen

ACS Paragon Plus Environment

Environmental Science & Technology

275

endocrine-disrupting chemicals.

276

The different transcription responses of genes related to the HPG axis are

277

presented in Fig. 3. In the presence of external stimuli, the hypothalamus of rare

278

minnow produced GnRH that bound to GnRHR in the pituitary to produce

279

gonadotropins (FSH and LH), which acted on the gonads. GnRHs are involved in

280

regulating reproductive activity and can be affected through negative feedback

281

mechanisms of sex hormones in vertebrates.33 The results in brain tissue indicated that

282

the mRNA levels of gnrh2 in males and gnrhr1α in females were significantly

283

increased after all of the treatments (p