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