Quantitative proteomic analysis to understand the mechanisms of zinc

Apr 3, 2019 - Quantitative proteomic analysis to understand the mechanisms of zinc oxide nanoparticle toxicity to Daphnia pulex (Crustacea: Daphniidae...
0 downloads 0 Views 762KB Size
Subscriber access provided by UNIV AUTONOMA DE COAHUILA UADEC

Ecotoxicology and Human Environmental Health

Quantitative proteomic analysis to understand the mechanisms of zinc oxide nanoparticle toxicity to Daphnia pulex (Crustacea: Daphniidae): : Comparing with Bulk Zinc Oxide and Zinc salt Li Lin, Mingzhi Xu, Huawei Mu, Wenwen Wang, Jin Sun, Jing He, Jian-Wen Qiu, and Tiangang Luan Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 03 Apr 2019 Downloaded from http://pubs.acs.org on April 3, 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 31

Environmental Science & Technology

1

Quantitative proteomic analysis to understand the mechanisms

2

of zinc oxide nanoparticle toxicity to Daphnia pulex (Crustacea:

3

Daphniidae): Comparing with Bulk Zinc Oxide and Zinc salt

4 5

Li Lin†, Mingzhi Xu†, Huawei Mu‡, Wenwen Wang†, Jin Sun §, Jing He†, Jian-

6

Wen Qiu‖, Tiangang Luan† ,*

7 8 9 10 11 12 13 14

†State

Key Laboratory of Biocontrol/School of Life Sciences, Sun Yat-sen University,

Guangzhou 510275, China ‡School

of Life Sciences, University of Science and Technology of China, Hefei

230071, China §Department

of Ocean Science, Hong Kong University of Science and Technology,

Hong Kong, China ‖Department

of Biology, Hong Kong Baptist University, Hong Kong, P. R. China

15 16 17

* Corresponding author: Tel: +86 2084112958; E-mail address: [email protected]

18

1

ACS Paragon Plus Environment

Environmental Science & Technology

19 20

Page 2 of 31

Graphical abstract:

21

Control

+ 114

ZnO NPs

+ 115

Bulk ZnO

+ 116

ZnSO4•7H2O

+ 117

Number of proteins

Up-regulation 200 100

Down-regulation 125

113

84

0 -100 -200 -300

178 ZnO NPs

218 235 Bulk ZnO ZnSO4·7H2O

Bulk ZnO

ZnO NPs 4

29 5

8 95

36 ZnSO4·7H2O

22 23 24 25

2

ACS Paragon Plus Environment

224

Page 3 of 31

Environmental Science & Technology

26

ABSTRACT:

27

The widespread use of zinc oxide nanoparticles (ZnO NPs) has resulted in their

28

release to the environment. There has been concern about the ecotoxicity of ZnO NPs,

29

but little is known about their toxic mechanisms. In the present study, we conducted

30

acute toxicity tests to show that ZnO NPs are more toxic to the freshwater crustacean

31

Daphnia pulex compared to bulk ZnO or ZnSO4·7H2O. To provide an integrated and

32

quantitative insights into the toxicity of ZnO NPs, we conducted isobaric tags for

33

relative and absolute quantitation (iTRAQ) proteomic analysis, which detected 262,

34

331 and 360 differentially expressed proteins (DEPs) in D. pulex exposed to ZnO

35

NPs, bulk ZnO and ZnSO4·7H2O, respectively. Among the DEPs, 224 were shared

36

among the three treatments. These proteins were related to energy metabolism,

37

oxidative stress and endoplasmic reticulum stress. The three forms of Zn all caused D.

38

pulex to downregulate chitinase expression, disrupt Ca2+ homeostasis and reduce

39

expression of digestive enzymes. Nevertheless, 29 proteins were expressed only in the

40

ZnO NP treatment. In particular, histone (H3) and ribosomal proteins (L13) were

41

obviously influenced under ZnO NP treatment. However, increased expression levels

42

of h3 and l13 genes were not induced only in ZnO NP treatment, they were sensitive

43

to Zn ions under the same exposure concentration. These results indicate that the three

44

zinc substances have a similar mode of action and that released zinc ions are the main

45

contributor to ZnO NP toxicity to D. pulex under a low concentration. Further

46

investigation is needed to clarify whether a small proportion of DEPs or higher

47

bioavailability cause ZnO NPs to be more toxic compared to bulk ZnO or ionic zinc.

48 49

Keywords: Crustaceans; Ecotoxicology; Heavy metals; Proteomics; Toxicity test;

50

Nanoparticles

3

ACS Paragon Plus Environment

Environmental Science & Technology

52 53

INTRODUCTION Nanomaterials having at least one dimension between 1 and 100 nm are widely

54

used in many applications because of their unique physicochemical properties.

55

According to the Nanotechnology Consumer Products Inventory (CPI), as many as

56

1814 consumer products from 622 companies in 32 countries contain nanomaterials.1

57

Zinc oxide nanoparticles (ZnO NPs) are one of the most widely used nanoparticles

58

and are commonly found in cosmetics, paints, textiles, sensing devices, rubber, and

59

ceramics. As ZnO NPs have a broad range of applications, there has been concern

60

about their release to the environment and their risk to organisms.2, 3 ZnO NPs have

61

been shown to be toxic to a variety of organisms (e.g., bacteria, algae, plants, aquatic

62

invertebrates, terrestrial invertebrates and aquatic vertebrates),4, 5 but the majority of

63

previous studies regarding ZnO NPs have focused on their acute toxicity or the

64

oxidative stress they cause.2

65

The toxicity of ZnO NPs may be caused by the particles,6 by zinc ions released

66

from ZnO NPs,7 or their interaction.8 ZnO NPs may produce chemical radicals or

67

ROS as their surface interacts with the surrounding media.2 The particles may also

68

interact with biomacromolecules such as protein and DNA, causing dysfunction.

69

Although zinc is an essential trace element, the zinc ions released from ZnO NPs can

70

be toxic, especially to the organisms that are susceptible to zinc. However, little is

71

known about the mechanisms of ZnO NP toxicity or the relative contribution of the

72

particles and the zinc ions to the overall toxicity of ZnO NPs.

73

Branchiopod crustaceans are susceptible to environmental pollutants and are

74

widely used as toxicological model organisms. Several studies have examined the

75

toxicity of ZnO NPs to cladocerans. In acute toxicity tests using Daphnia magna,

76

Blinova et al.9 found that the 48-h EC50 of ZnO NPs varied between 1.7 and 9.0 mg

77

Zn/L, whereas Santo et al.10 reported that the 48-h EC50 of ZnO nanopowders varied

78

from 1.9 to 3.1 mg/L. Adam et al. reported that D. magna can regulate its internal zinc

79

concentration after exposure to ZnO nanoparticles and that toxicity was caused by the

80

dissolved fraction.11 In a chronic toxicity test, Adam et al.12 found that the 4

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Environmental Science & Technology

81

reproductive EC50 of ZnO NPs for D. magna was 0.112 mg Zn/L. Nevertheless, only

82

a few studies have examined the molecular mechanisms of NP toxicity. Proteomic

83

tools have been used to compare the toxic mechanisms of citrate-coated silver NPs on

84

D. magna compared to that of AgNO3 alone.13 Protein thiol contents and carbonyl

85

levels were determined, and a two-dimensional electrophoresis (2DE) analysis was

86

conducted showing that the two forms of silver had different modes of toxicity. They

87

found 17 differentially expressed proteins, among which only eight were identified by

88

MS/MS analysis. To our understanding, no proteomic study of the toxic mechanisms

89

of ZnO NPs has been conducted.

90

Daphnia pulex is a cornerstone species in freshwater food chains. Since D.

91

pulex is the first crustacean whose genome has been sequenced and annotated,14 it is

92

an ideal species for studies of toxic mechanisms. The present study thus determined

93

the acute toxicity of ZnO NPs to D. pulex and applied isobaric tags for relative and

94

absolute quantitation (iTRAQ) proteomic analysis to compare the effects of ZnO NPs,

95

bulk ZnO, and ZnSO4·7H2O on its protein expression. Based on the proteomic results,

96

the gene expression levels of h3, ll3, idh, scot and eip gene were analysed using q-

97

PCR to verify the molecular changes under the treatments.

98 99

MATERIALS AND METHODS

100

Test Reagents and Animal Species. Zinc oxide nanoparticles (nanopowder < 100

101

nm particle size, specific surface area 15-25 m2/g, 79.1-81.5% Zn, and complex-

102

ometric titration) were purchased from Sigma-Aldrich (Shanghai, China). X-Ray

103

Diffraction (XRD) analysis showed that ZnO NPs were zincite structure (a = 0.325

104

nm, c = 0.521 nm). ZnO NP stock suspensions were freshly prepared in reconstituted

105

moderately hard water (RMHW) 15 after 30 min of sonication in a water bath

106

sonicator. The RMHW (pH 7.8 ± 0.2, ionic strength 0.0051) contained NaHCO3 (96

107

mg/L), MgSO4 (60 mg/L), KCl (4 mg/L) and CaSO4·2H2O (60 mg/L).

108 109

The shape and size of ZnO NPs were characterized by using transmission electron microscopy (TEM, FEI, the Netherlands). The image demonstrates that the 5

ACS Paragon Plus Environment

Environmental Science & Technology

110

ZnO NPs used in this study were slender and spherical in shape (Figure S1), with a

111

mean size of 61 ± 12 nm (length of longest dimension ± standard deviation),

112

determined by measuring 200 particles.

113

Daphnia pulex was obtained from the School of Marine Sciences, Ningbo

114

University (Ningbo, China). The culture stock originated from a lake in Minhang

115

District (Shanghai, China) and has been maintained for more than three years in the

116

laboratory.16 We started our lab culture with three individuals. The cladocerans were

117

cultured in the RMHW with a density of 1 ind./50 mL under a 16 h light : 8 h dark

118

cycle (20 ± 1 °C) in an incubator. The medium was replenished three times per week,

119

and the animals were fed with Selenastrum capricornutum at concentrations of

120

105~106 ind./mL.

121 122

Dissolution Measurements. A 24 h dissolution experiment was performed to

123

quantify the contents of dissolved Zn released from ZnO NPs and bulk ZnO at 1

124

mg/L, 10 mg/L and 100 mg/L. At 0 (immediately after mixing the chemicals with

125

RMHW), 0.5, 1, 2, 4, 6, 12 and 24 h, samples were centrifuged at 10000 g for 15 min

126

(Thermo Fisher Scientific, Waltham, MA). The supernatant was filtered through a

127

syringe filter with a 0.22 μm pore diameter, and the Zn concentration was measured

128

using inductively coupled plasma optical emission spectrometry (ICP-OES, Thermo

129

Fisher Scientific, Waltham, MA).

130 131

Acute Toxicity Test. The acute toxicity of ZnO NPs, bulk ZnO and ZnSO4·7H2O to

132

D. pulex was determined according to the USEPA methods for measuring the acute

133

toxicity of effluents and receiving waters to freshwater and marine organisms.15

134

According to the methods, neonates were used. The effective concentration of each of

135

the three chemicals was different based on preliminary experiments. For the precise

136

calculation of EC50 values, exposure concentrations were 0.15, 0.3, 0.6, 1.2, 2.4, 4.8

137

and 9.6 mg/L for ZnO NPs, 0.1, 0.2, 0.4, 0.8, 1.6, 3.2 and 6.4 mg/L for bulk ZnO, and

138

0.25, 0.5, 1.0, 2.0, 4.0, 8.0 and 16.0 mg/L for ZnSO4·7H2O. Each concentration

139

treatment contained 4 replicates, and each replicate comprised 10 neonates (< 24 h) 6

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Environmental Science & Technology

140

and contained 50 mL of exposure solution or RMHW (control). Daphnids were

141

incubated under the above mentioned culture conditions for a 48 h exposure period

142

without feeding. At 24 h and 48 h, immobilized individuals in beakers were recorded.

143

The median effective (immobilization) concentration (EC50) was determined

144

according to Karber's method.17 The results from these toxicity tests were used to

145

determine exposure concentrations for use in proteomic analysis. Because the

146

neonates are usually more sensitive than the adults, the survival of the daphnids in the

147

following experiments could be ensured under this concentration.

148 149

Acute Exposure for Proteomics. Twenty-five daphnids (5-day-old) were transferred

150

into a beaker containing 250 mL of sublethal concentration (1/5th of the 24 h EC50) of

151

ZnO NPs, bulk ZnO, ZnSO4·7H2O or RMHW only (control). Each treatment was

152

performed in triplicate.

153

The beakers were placed in the culture conditions described above. After exposure for

154

24 h, the daphnids were collected into a 1.5-mL Eppendorf tube containing 8 M of

155

urea and were immediately frozen with liquid nitrogen. Then, the samples were stored

156

at -80°C.

157 158

Protein Extraction, Digestion, and iTRAQ Labelling. Samples were placed on ice

159

and homogenized using a plastic pestle. A mild sonication was applied using a Sonic

160

Dismembrator 300 (Thermo Fisher Scientific) to break cells and release the proteins

161

using the following sonication cycle: 1 min at a power setting of 28, followed by three

162

bursts of 15 s at a power setting of 35, each with a 30 s pause. The samples were then

163

centrifuged at 15000 g and 4°C for 15 min, and the supernatant was collected and

164

purified using a methanol and chloroform method.18 For 100 μL of sample, 400 μL of

165

methanol, 100 μL of chloroform and 300 μL of Milli-Q water were added

166

sequentially. The mixed solution was centrifuged at 14000 g for 2 min, the top layer

167

was removed, and 400 μL of methanol was added to wash the protein pellet, followed

168

by a 3-min centrifugation at 14000 g. The supernatant was removed, the protein pellet

169

was reconstituted with 8 M urea, and the protein concentration was determined using 7

ACS Paragon Plus Environment

Environmental Science & Technology

170

an RC-DC kit (Bio-Rad, Hercules, CA). Next, 100 μg of protein was reduced with 5

171

mM triscarboxyethyl phosphine hydrochloride at 60°C for 60 min and alkylated with

172

10 mM methylethanethio sulfonate at room temperature for 20 min. A 50 mM

173

triethylammonium bicarbonate (TEAB) solution was applied to dilute sample solution

174

by 8 fold. The samples were then digested using sequencing-grade trypsin (Promega,

175

Madison, WI) for 16 h at 37°C with a 1: 50 trypsin to protein mass ratio. The digested

176

samples were desalted with Sep-Pak C18 cartridges (Waters, Milford, MA) and dried

177

in a vacuum concentrator (Eppendorf, Hamburg, Germany). Samples derived from the

178

control, ZnO NP, bulk ZnO, and ZnSO4·7H2O treatments were labelled with iTRAQ

179

114, 115, 116, and 117 reagents (AB Sciex, Framingham, MA), respectively. Samples

180

from the four treatments were pooled together and dried in a vacuum concentrator.

181

Three biological replicates were used for the following analysis.

182 183

SCX Fractionation and LC-MS/MS Analysis. The procedure of SCX fractionation

184

followed that of Mu et al.19 Dried samples were dissolved in buffer A [20%

185

acetonitrile (ACN) and 10 mM KH2PO4, pH = 3.0] and fractionated using a

186

PolySULFOETHYL strong cation-exchange column (5 μm particle size, 200 × 4.6

187

mm, 200 Å pore size, PolyLC, Columbia, MD) on a Waters 2695 HPLC. The running

188

gradient was as follows: 5 min in 100% buffer A, followed by a 28 min 0-30% linear

189

gradient to buffer B (0.5 M KCl, 10 mM KH2PO4 and 20% ACN, pH = 3.0), a 5 min

190

30-100% buffer B linear gradient, 5 min in 100% buffer B, and finally 7 min in 100%

191

buffer A. The gradient flow rate was 1 mL/min. A total of 14 fractions were collected,

192

dried in a vacuum concentrator, and desalted with Sep-Pak C18 cartridges.

193

The dried fractions were dissolved in 0.1% formic acid. Each sample was

194

analysed twice with an LTQ-Orbitrap Elite coupled to an Easy-nLC (Thermo Fisher,

195

Bremen, Germany). Peptides were separated with a C18 capillary column (Michrom

196

BioResources, CA). The following running gradient was set: 5 min in 100% solution

197

A (0.1% formic acid in Milli-Q water), a 55 min 0-30% solution B (0.1% formic acid

198

in ACN) linear gradient, a 10 min 30-98% solution B gradient, 10 min in 98%

199

solution B, and 10 min in 100% solution A. Mass spectrometry scans were run from 8

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

Environmental Science & Technology

200

350 to 1600 m/z with a resolution of 60000 in the positive charge mode. The top five

201

multiple-charged ions with the highest intensities and a minimum signal threshold of

202

500.0 were selected for fragmentation using high-energy collision-induced

203

dissociation (HCD) and collision-induced dissociation (CID). Peptides were identified

204

with the CID spectra in LTQ, while iTRAQ reporters were quantified with the HCD

205

spectra in C-trap. The CID and HCD scanning approaches adopted an isolation width

206

of 2.0 m/z, a dynamic exclusion time of 30 s, and a mass window for precursor ion

207

selection of 2.0 Th. CID fragmentation was conducted under normalized collision

208

energy of 35%, activation Q of 0.25, and an activation time of 10 ms. HCD

209

fragmentation was conducted under full scan with FTMS at a resolution of 15000 in

210

centroid mode, a normalized collision energy of 45% and an activation time of 10 ms.

211 212

Mass Spectrometry Data Analysis. Mass spectrometry data were analysed as in Mu

213

et al., 2018.20 Proteome Discoverer 1.3.0.339 (Thermo Finnigan, CA) was applied to

214

convert raw data to .mgf files. Unpaired scans were filtered, and data derived from the

215

CID and HCD scans were also separated with Python scripts. Reporter groups in the

216

CID scans were replaced by those of HCD scans from the same parent ion. The

217

intensity of ions in HCD was normalized. The new CID and HCD data files were

218

submitted to Mascot version 2.3.2 (Matrix Sciences, London, UK) to search against a

219

D. pulex protein database that was downloaded from NCBI and contained 32487

220

protein sequences. The search parameters for CID files were set as: peptide charges of

221

2+, 3+, and 4+, up to 2 missed cleavage for trypsin-digested peptides, variable

222

modification: deamidated glutamine and asparagine, oxidation of methionine, fixed

223

modification: methylthio of cysteine, precursor ions of ± 5 ppm, and MS/MS ion

224

search with fragment ion tolerance of ± 0.6 Da. HCD parameters were similar to those

225

of CID except for the fragment ion tolerance, which was set to ± 20 mmu. Data with

226

ion scores lower than 27 were removed. Since short peptides might match decoy

227

sequences, peptides shorter than seven amino acids were removed. The false

228

discovery rate was set as 1%. Unlabelled peptides or peptides with an erratic ratio

229

between any two iTRAQ reporters were removed. The median ratios between 9

ACS Paragon Plus Environment

Environmental Science & Technology

230

different iTRAQ reporters (i.e., 115/114, 116/114 and 117/114) were normalized to 1.

231

Quantification was conducted with the proteins that were detected in all three

232

biological replicates and had at least four summed matched peptides. Protein

233

quantification was based on the summed intensity of matched peptides. Data from two

234

technical runs of the same biological replicate were used to determine the criteria of

235

significant differential expression. The cutoff values were defined as the log2-

236

transformed data at which 95% of all target proteins did not deviate from each other.

237

The cutoff values were 0.81 and -0.80 for 115/114, 1.01 and -0.95 for 116/114 and

238

1.00 and -0.92 for 117/114, which corresponded to 1.76- and 0.57-fold for 115/114,

239

2.01- and 0.52-fold for 116/114 and 2.00- and 0.53-fold for 117/114, respectively.

240

These fold values were used as thresholds to define significantly up- or

241

downregulated proteins between control and treatment groups.

242 243

Bioinformatics analysis. PANTHER21 was applied to classify proteins into different

244

functional groups based on their Gene Ontology (GO) annotations. Differentially

245

expressed proteins (DEPs) were assigned to several level 2 Biological Process (BP)

246

GO terms. DEPs without GO annotations were classified into “Others”. The

247

expression levels of DEPs under each BP term were compared using log2-transformed

248

data to produce a Euclidean distance similarity matrix and were clustered using the

249

centroid linkage method and visualized by Java Treeview. GO enrichment was

250

performed using the Gene Ontology Enrichment Analysis Software Toolkit

251

(GOEAST).22 Hypothetical proteins were analysed by the Basic Local Alignment

252

Search Tool (BLAST) in UniProt and searched against the UniProt Knowledgebase

253

(UniProtKB).

254 255

q-PCR analysis. Every 150 juveniles (4-6 day) of D. pulex were placed in a 500 mL

256

glass beaker containing 500 mL of exposure media and exposed for 24 h to ZnO NPs,

257

bulk ZnO or ZnSO4, at a concentration of 0.06 mg Zn/L, which was approximately

258

equal to 1/5th of the EC50 value of the ZnO NPs. The treatment was carried out in

259

triplicate. Total RNA was extracted from 100 living animals using the E.Z.N.A. 10

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

Environmental Science & Technology

260

Total® RNA Kit II (OMEGA, USA) according to the manufacturer’s protocol. The

261

samples were then treated with DNase. RNA integrity and purity were assessed by 1%

262

agarose gel electrophoresis and a spectrophotometer (NanoDrop ND-1000, USA),

263

respectively. The A260/A280 values were between 1.8 and 2.2, which indicated good

264

sample quality. Total cDNA for the real-time PCRs of each sample was generated

265

with a PrimeScript ™ RT reagent Kit (TaKaRa, Japan) according to the

266

manufacturer's instructions. The cDNA was stored at -20℃ until further use for q-

267

PCR.

268

The gapdh gene was used as a reference.23 Primer-BLAST was used to design

269

primers for h3 (gi|321463426), l13 (gi|321474005), scot (gi|321478252), idh

270

(gi|321478958) and eip (gi|321459564). These genes were selected from the DEPs

271

according to the prior proteomic results. Primers were used for qPCR, the PCR

272

products were sequenced and the results were tested using BLAST to ensure that the

273

sequence homology was above 95% (Table S1).

274

The final volume of each 10 μl qPCR system contained 5 μl of SYBR® Green

275

Real-time PCR Master Mix, 0.4 μl of each primer (10 μm), 1 μl of diluted cDNA

276

template and 3.2 μl of PCR-grade water. The experimental protocol for qPCR

277

verification consisted of a denaturation phase (10 min at 95℃), an amplification phase

278

of 40 cycles (10 s at 95℃, 20 s at 58℃ and 20 s at 72℃) and a final melting phase.

279

Relative mRNA expression levels of target genes were calculated using the Pfaffl

280

method.24 The data were analysed by independent sample t-tests to compare means

281

between treatment and control. Statistical analysis was conducted with SPSS 22.0.

282 283 284

RESULTS

285

Dissolution of ZnO NPs and bulk ZnO in exposure media. The content of

286

dissolved Zn was measured for 1, 10, and 100 mg/L ZnO NPs and bulk ZnO in

287

RMHW (Figure 1). In ZnO NP suspensions, the content of dissolved Zn increased

288

rapidly in 1 h and tended to become stable after 6 h for all three concentrations. At 24

289

h, the dissolved Zn was 0.47, 0.62 and 0.69 mg/L for the 1, 10 and 100 mg/L ZnO NP 11

ACS Paragon Plus Environment

Environmental Science & Technology

290

solutions, respectively. The amount of dissolved Zn released from the ZnO NPs was

291

in the same order of magnitude for the three concentrations. The bulk ZnO dissolved

292

in a similar pattern but had a lower dissolution rate and content of dissolved Zn. At 24

293

h, the dissolved Zn concentrations were 0.17, 0.18 and 0.23 mg/L for the 1, 10, 100

294

mg/L bulk ZnO solutions, respectively.

295 296

Acute toxicity tests. The toxicity of the ZnO NPs, bulk ZnO and ZnSO4·7H2O to D.

297

pulex was determined by acute immobilization tests (Table 1). The EC50 values of all

298

three tested chemicals were below 1 mg/L. Compared to bulk ZnO and ZnSO4·7H2O,

299

the ZnO NPs had a much lower EC50. Additional details of the toxicity tests are

300

shown in Figure S2.

301 302

Proteomic analyses. In total, 1652, 1827 and 2142 proteins were identified from the

303

three biological replicates. Of these proteins, 1371 were detected in all three

304

biological replicates. Exposure to ZnO NPs, bulk ZnO and ZnSO4·7H2O resulted in

305

262, 331 and 360 DEPs, respectively, when compared with the control (Figure 2).

306

Approximately two-thirds of DEPs were downregulated in each treatment. There was

307

a large overlap of DEPs (224 proteins) among the three treatments (Figure 3).

308

Furthermore, many these shared DEPs had identical expression trends across the three

309

treatments (Figure S3). Only 29 DEPs were unique to the ZnO NP treatment.

310

GO analyses were performed on the 224 shared DEPs (Figure 4). For cellular

311

components, 37.5% of them related to the cell, followed by organelles (21.3%) and

312

extracellular regions (16.3%). For molecular functions, 35.2% of DEPs were related

313

to catalytic activity, and 25.9% were related to binding. For biological processes,

314

37.9% were related to metabolic processes, followed by cellular processes (16.8%)

315

and localization (15.9%).

316

Pathway enrichment analysis of the 224 shared DEPs resulted in 47 GO terms

317

(Table S2), including 25 biological process GO terms (Table 2). The enriched

318

pathways included cellular amide metabolic process, ion transport, single-organism

319

transport, aminoglycan metabolic process, cellular metal ion homeostasis, organic 12

ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31

320 321

Environmental Science & Technology

acid metabolic process, tricarboxylic acid cycle and chitin metabolic process. Among the 29 DEPs which expressed only under the ZnO NP treatment, histone

322

(H3) and ribosomal proteins (L13) were obviously influenced under the ZnO NP

323

treatment. Isocitrate dehydrogenase (IDH), succinyl-CoA:3-ketoacid-coenzyme

324

ATransferase (SCOT) and ecdysone-induced protein (EIP) were changed also.

325 326

Gene expression. The expression levels of h3 and l13 were reduced in the ZnO NP

327

and ZnSO4 treatments (Figure 5), and the expression levels of scot, idh and eip were

328

reduced in the ZnSO4 treatment.

329 330

DISCUSSION

331

Published studies have demonstrated that cladocerans are sensitive to ZnO NPs.9, 10, 12

332

In a large proportion of the published toxicity studies on cladocerans, D. magna was

333

chosen as a model species. In our study, D. pulex showed more sensitivity to ZnO

334

NPs compared to the published results on D. magna.

335

The toxicity tests results showed that the ZnO NPs were more toxic than the

336

bulk ZnO and ZnSO4∙7H2O. However, the proteomic analyses were unable to

337

distinguish the toxicity effects among the ZnO NPs, bulk ZnO and ZnSO4∙7H2O.

338

When exposed to these three substances, a large overlap in differentially expressed

339

proteins was observed, and these common DEPs also had coincident expression

340

trends. By comparing the protein expression patterns in the three treatments, we found

341

that ZnO NPs, bulk ZnO and ZnSO4∙7H2O were acting through a similar mode of

342

action. The zinc ions may be the primary contributors to the toxic effect of ZnO NPs.

343

Similar to our conclusion, the importance of zinc ions to ZnO NP toxicity in D.

344

magna was demonstrated by Adam et al.25 in a transcriptomic study. These authors

345

reported that no significantly differentially transcribed gene fragments were found by

346

microarray when comparing daphnid exposure to ZnO NPs and the corresponding

347

metal salt. The bioavailability of ZnO NPs affected by their agglomeration and ion

348

dissolution. Agglomeration was not obvious under low concentration. Therefore, the 13

ACS Paragon Plus Environment

Environmental Science & Technology

349

toxicity of ZnO nanoparticles to D. magna is mainly caused by toxic metal ions under

350

the current exposure conditions. To our knowledge, the present study is the first report

351

to use proteomic analysis to compare the toxic mechanisms of ZnO NPs to

352

crustaceans with common ZnO and Zn ions. The three treatments resulted in different

353

numbers of DEPs. This might be due to the different Zn concentrations. In this study,

354

a concentration of 1/5th of the EC50 was chosen.

355

The 224 DEPs in the present study were found to be markedly associated with

356

Ca2+ homeostasis, energy metabolism, oxidative stress, endoplasmatic reticulum (ER)

357

stress and chitin metabolism. The disruption of Ca2+ homeostasis would cause

358

interference to various physiological processes. ZnO NPs were found to interfere with

359

Ca2+ homeostasis in previous studies. Wang et al.26 found that ZnO NPs increased the

360

resting Ca2+ concentration without compromising Ca2+ homeostatic mechanisms and

361

strongly inhibited store-operated Ca2+ entry in a muscarinic receptor signalling

362

pathway. Guo et al.27 found that ZnO NPs disrupted the intracellular Ca2+ homeostasis

363

in rat retinal ganglion cells by decreasing the expression and activity of plasma

364

membrane Ca2+ ATPase. In the present study, some common DEPs in the three

365

treatments were also found to be connected with Ca2+ homeostasis (Table S3). Of

366

these, Na+/K+ ATPase was downregulated. Downregulation of Na+/K+ ATPase might

367

interfere with Ca2+ homeostasis.28 In addition, some proteins such as calmodulin and

368

annexin, which are regulated by Ca2+, were also found to be differentially expressed.

369

The function of mitochondria is also affected by Ca2+ homeostasis. Some key

370

enzymes of the Krebs cycle, such as pyruvate dehydrogenase, NAD+-isocitrate

371

dehydrogenase and oxoglutarate dehydrogenase, are positively regulated by Ca2+.29

372

Cells may respond to Ca2+ signals by enhancing mitochondrial ATP production to

373

balance the increased ATP demand.30 In this study, some of the 224 common DEPs

374

were found to be enzymes of energy synthesis (Table S4). These enzymes, including

375

glyceraldehyde-3-phosphate dehydrogenase, mitochondrial malate dehydrogenase and

376

ATP synthase, were all upregulated, suggesting enhanced energy synthesis. Triboulet

377

et al.31 found that carbohydrate catabolism is one of the critical determinants of

378

sensitivity to zinc when they exposed macrophages to zinc ions and ZnO NPs. 14

ACS Paragon Plus Environment

Page 14 of 31

Page 15 of 31

379

Environmental Science & Technology

The formation of reactive oxygen species (ROS) and oxidative stress are

380

considered to be involved in the toxicity of ZnO NPs.8 Mwaanga et al.32 found that

381

ZnO NPs caused oxidative stress that is related to biochemical changes in D. magna.

382

The exposure caused decreased GST activity and increased GSH oxidation, TBARS

383

formation and MT induction. In our study, some proteins related to antioxidant

384

defence, such as vitellogenins fused with SOD, peroxiredoxin and glutathione S-

385

transferase, were upregulated in the three treatments (Table S5). Vitellogenins fused

386

with SOD can be found in daphnids33 and brine shrimp.34 Rainville et al.13 found that

387

vitellogenins fused with SOD were upregulated when D. magna was exposed to Ag

388

NPs, and the upregulation was believed to evolve in response to the oxidative stress

389

caused by Ag NPs. Vandegehuchte et al.35 found that Zn exposure resulted in the

390

upregulation of peroxiredoxin 6 and glutathione S-transferase in D. magna. The

391

upregulation of these proteins in the present study indicated the oxidative stress

392

response in D. pulex.

393

Oxidative stress, the inhibition of protein glycosylation, a reduction in disulfide

394

bond formation, calcium depletion and the accumulation of unfolded proteins in the

395

ER lumen may cause disruption of the ER function and induce an ER stress response

396

via the unfolded protein response (UPR).36 Chen et al.37 found that ZnO NPs activated

397

the ER stress-responsive pathway in human umbilical vein endothelial cells

398

(HUVECs), and the ER stress response might be used as an early and sensitive

399

endpoint for nanotoxicological studies. Yang et al.36 found that ER stress-induced

400

apoptosis is involved in hepato-toxicity induced by ZnO NPs in mice. In our study,

401

some DEPs were also found to be related to ER stress, indicating that ER stress was

402

involved in the toxicity (Table S6). The downregulation of glycosyl transferase

403

suggested the inhibition of protein glycosylation. The downregulation of RNA-

404

binding protein, rRNA methyltransferase, signal recognition particle receptor and

405

signal peptidase were also observed, indicating that the protein synthesis and

406

translocation across the ER was inhibited to reduce the stress of protein folding.

407

Protein disulfide isomerase (PDI) is an essential folding catalyst and chaperone of the

408

ER and is responsible for introducing disulfides into proteins and catalyses the 15

ACS Paragon Plus Environment

Environmental Science & Technology

409

rearrangement of incorrect disulfides.38 The upregulation of PDI suggested an

410

enhancement of protein folding in the ER.

411

The downregulation of many digestive enzymes was observed in all three

412

treatments in this study (Table S7). The expression of digestive enzymes, such as

413

dipeptidase, serine protease, chymotrypsin, and carbohydrate digestive enzymes (e.g.,

414

beta-galactosidase, alpha-amylase and endo-beta-1,4-mannanase) were significantly

415

decreased. Since the uptake of these substances might first occur in the gut of D.

416

pulex, it might cause substantial damage to the digestive system, resulting in a

417

decrease of digestive enzymes. We also found that some chitinases were

418

downregulated in all three treatments (Table S8).

419

Notably, 29 proteins were expressed only in the ZnO NP exposure test, though

420

there was a large overlap of DEPs across the three exposure conditions. Among the

421

overlapped DEPs, histone (H3) and ribosomal proteins (L13) were obviously

422

influenced under the ZnO NP treatment. Histones are the basic structural proteins of

423

chromosomes. Histones mainly contain five kinds. DNA molecules are wrapped

424

around histone octamers formed by the polymerization of H2A, H2B, H3 and H4 to

425

form the basic structure of chromatin - nucleosomes. The binding of H1 and DNA

426

between nucleosomes controls the entry and exit of DNA in order to stabilize

427

nucleosomes.39 Ribosomal Protein (Ribosomal Protein) and ribosomal RNA (rRNA)

428

constitute ribosomes, which play an important role in the protein biosynthesis. IDH,

429

SCOT and EIP were also influenced under ZnO NPs treatment. IDH and SCOT are

430

enzymes related to the tricarboxylic acid cycle, which is a key process of body energy

431

metabolism. EIP belongs to nuclear receptor superfamily as a transcription factor and

432

plays an important role in ecdysone pathway. 40

433

Similar results have been reported by other researchers. Methylation of H3 and

434

expression of some ribosomal proteins were found to be affected by gold and silver

435

nanoparticles.41, 42 From this experimental evidence, we hypothesize that

436

nanoparticles might enter cells and affect genetic processes, such as DNA/chromatin

437

structure, transcription, and translation. This might be one of the causes of the

438

significant differences in toxicity among the three treatments. From this view, the 16

ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

Environmental Science & Technology

439

effects caused by changes in proteins that were expressed only in the ZnO NP

440

treatment group may not correlate with their respective proportion in DEPs. Omics

441

studies usually draw conclusions from the proportion of differentially expressed

442

proteins or genes. This seems to be not objective and comprehensive enough. It would

443

be valuable to investigate the contribution of the changes in this small proportion of

444

DEPs to clarify the toxicity effects of ZnO NPs.

445

The expression levels of five related genes were selected for investigation among

446

the above mentioned 29 proteins. However, the results of the gene expression level

447

detection were not consistent with the results of the proteomics research. These 5

448

genes were not particularly induced in the ZnO NP treatment. Instead, they were

449

sensitive to Zn ions under the same exposure concentration. The gene expression

450

levels of h3 and l13 were reduced under ZnO NP and ZnSO4 exposure (Figure 5).

451

However, no difference was found between ZnO NP and Zn ion exposure. Expression

452

levels of scot, idh (energy metabolism-related genes) and eip (moulting-related gene)

453

mRNA showed no difference between the ZnO NP treatment and the control, but they

454

were inhibited under the ZnSO4 exposure (Figure 5). Poynton et al.43 also found that

455

Zn caused a decrease in chitinase activity in D. magna. Unlike ZnO NPs, non-

456

dissociation TiO2 NPs can promote moulting of D. magna.44 From the gene

457

expression results, it seems that Zn ions, rather than nanoparticles, are the main reason

458

for the toxic effects on moulting and genetic energy metabolism at relatively lower

459

treatment concentrations.

460

Combining the results of the acute toxicity tests and proteomic analyses, we

461

hypothesize that ZnO NPs have a higher bioavailability to D. pulex than bulk ZnO

462

and ZnSO4∙7H2O, although the three substances have a similar mode of action. The

463

relatively higher bioavailability endows the ZnO NPs with a higher toxicity, even at

464

lower concentrations, compared to bulk ZnO and ZnSO4∙7H2O. The marked

465

bioavailability of ZnO NPs was also mentioned in previously published studies.

466

Poynton et al.45 found that ZnO NPs were more toxic to the epibenthic crustacean

467

Hyalella azteca compared with zinc sulfate. However, gene expression analysis is

468

unable to distinguish between ZnO NP exposure and zinc sulfate exposure. The 17

ACS Paragon Plus Environment

Environmental Science & Technology

469

authors hypothesized that ZnO NPs provide an enhanced exposure route for Zn2+

470

uptake in H. azteca. By synthesizing 65Zn tracer, Li et al. demonstrated a distinctive

471

uptake mode of ZnO NPs compared to zinc ions.46 However, after assessing the

472

results of the present study and the previous reports, we believe that there is a need for

473

more evidence to confirm this.

474

In conclusion, proteomic research provided an integrated and quantitative view of

475

toxic mechanism at the protein level. The present study firstly revealed that the

476

pathway of ZnO NPs were the same with common ZnO and Zn ions to daphnids

477

under a low concentration by using proteomic analysis.

478 479

Implications for Toxicity Mechanisms of ZnO NPs. In the present study, we

480

investigated the toxicity of ZnO NPs, bulk ZnO and ZnSO4∙7H2O to the freshwater

481

crustacean D. pulex by proteomic analysis. The results revealed that ZnO NPs, bulk

482

ZnO and ZnSO4∙7H2O were acting through a similar mode of action and that the zinc

483

ions may primarily contribute to the toxic effect of the ZnO NPs. There were 29

484

proteins that were expressed only upon ZnO NP exposure. Though these made up a

485

small proportion of total DEPs, the contribution from the changes of these proteins are

486

valuable for further investigation of the toxic effects of ZnO NPs. Our results suggest

487

that ZnO NPs might have a higher bioavailability to D. pulex than does the bulk ZnO

488

and ionic Zn. Increased attention should be given to the difference in bioavailability

489

of ZnO NPs and Zn from other sources when investigating ZnO NP toxicity. This

490

concept could also apply to other metal-based nanoparticles with dissolving capacity.

491 492

Supporting Information

493

Materials including dose-response curves, expression patterns of the common DEPs,

494

GO enrichment analysis results, and details of the common DEPs associated with

495

different physiological processes are freely available via the Internet at

496

http://pubs.acs.org.

497 18

ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

Environmental Science & Technology

498

ACKNOWLEDGMENTS

499

This research was supported by the National Key R & D Program of China

500

(2018YFD0900803, 2018YFD0900604), National Science Foundation of China

501

(Project No. 21777197, 21625703), and the Science Foundation of Guangdong

502

Province (Project No. 2016A030313337).

503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537

REFERENCES (1) Vance, M. E.; Kuiken, T.; Vejerano, E. P.; McGinnis, S. P.; Hochella, M. F., Jr.; Rejeski, D.; Hull, M. S., Nanotechnology in the real world: Redeveloping the nanomaterial consumer products inventory. Beilstein J Nanotechnol 2015, 6, 1769-80. (2) Vale, G.; Mehennaoui, K.; Cambier, S.; Libralato, G.; Jomini, S.; Domingos, R. F., Manufactured nanoparticles in the aquatic environment-biochemical responses on freshwater organisms: A critical overview. Aquat Toxicol 2016, 170, 162-174. (3) Selck, H.; Handy, R. D.; Fernandes, T. F.; Klaine, S. J.; Petersen, E. J., Nanomaterials in the Aquatic Environment: A European Union-United States Perspective on the Status of Ecotoxicity Testing, Research Priorities, and Challenges Ahead. Environ Toxicol Chem 2016, 35, (5), 1055-1067. (4) Ma, H.; Williams, P. L.; Diamond, S. A., Ecotoxicity of manufactured ZnO nanoparticles--a review. Environ Pollut 2013, 172, 76-85. (5) Adam, N.; Schmitt, C.; De Bruyn, L.; Knapen, D.; Blust, R., Aquatic acute species sensitivity distributions of ZnO and CuO nanoparticles. Sci Total Environ 2015, 526, 233-242. (6) Poynton, H. C. L., J.; Impellitteri, C.; Smith, M. E.; Rogers, K.; Patra, M.; Hammer, K.; Allen, J.; Vulpe, C., Differential Gene Expression in Daphnia magna Suggests Distinct Modes of Action and Bioavailability for ZnO Nanoparticles and Zn Ions. Environ. Sci. Technol. 2011, 45, (2), 762-768. (7) Landa, P.; Prerostova, S.; Petrova, S.; Knirsch, V.; Vankova, R.; Vanek, T., The Transcriptomic Response of Arabidopsis thaliana to Zinc Oxide: A Comparison of the Impact of Nanoparticle, Bulk, and Ionic Zinc. Environ Sci Technol 2015, 49, (24), 14537-45. (8) Zhao, X.; Wang, S.; Wu, Y.; You, H.; Lv, L., Acute ZnO nanoparticles exposure induces developmental toxicity, oxidative stress and DNA damage in embryo-larval zebrafish. Aquat Toxicol 2013, 136-137, 49-59. (9) Blinova, I.; Ivask, A.; Heinlaan, M.; Mortimer, M.; Kahru, A., Ecotoxicity of nanoparticles of CuO and ZnO in natural water. Environ Pollut 2010, 158, (1), 41-7. (10) Santo, N.; Fascio, U.; Torres, F.; Guazzoni, N.; Tremolada, P.; Bettinetti, R.; Mantecca, P.; Bacchetta, R., Toxic effects and ultrastructural damages to Daphnia magna of two differently sized ZnO nanoparticles: Does size matter? Water Res 2014, 53, 339-350. (11) Adam, N.; Leroux, F.; Knapen, D.; Bals, S.; Blust, R., The uptake and elimination of ZnO and CuO nanoparticles in Daphnia magna under chronic exposure scenarios. Water Res 2015, 68, 249-261. (12) Adam, N.; Schmitt, C.; Galceran, J.; Companys, E.; Vakurov, A.; Wallace, R.; Knapen, D.; Blust, R., The chronic toxicity of ZnO nanoparticles and ZnCl2 to Daphnia magna and the use of different methods to assess nanoparticle aggregation and dissolution. Nanotoxicology 2014, 8, (7), 709-17. (13) Rainville, L. C.; Carolan, D.; Varela, A. C.; Doyle, H.; Sheehan, D., Proteomic evaluation of 19

ACS Paragon Plus Environment

Environmental Science & Technology

538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581

citrate-coated silver nanoparticles toxicity in Daphnia magna. Analyst 2014, 139, (7), 1678-86. (14) Colbourne, J. K.; Pfrender, M. E.; Gilbert, D.; Thomas, W. K.; Tucker, A.; Oakley, T. H.; Tokishita, S.; Aerts, A.; Arnold, G. J.; Basu, M. K.; Bauer, D. J.; Caceres, C. E.; Carmel, L.; Casola, C.; Choi, J. H.; Detter, J. C.; Dong, Q.; Dusheyko, S.; Eads, B. D.; Frohlich, T.; Geiler-Samerotte, K. A.; Gerlach, D.; Hatcher, P.; Jogdeo, S.; Krijgsveld, J.; Kriventseva, E. V.; Kultz, D.; Laforsch, C.; Lindquist, E.; Lopez, J.; Manak, J. R.; Muller, J.; Pangilinan, J.; Patwardhan, R. P.; Pitluck, S.; Pritham, E. J.; Rechtsteiner, A.; Rho, M.; Rogozin, I. B.; Sakarya, O.; Salamov, A.; Schaack, S.; Shapiro, H.; Shiga, Y.; Skalitzky, C.; Smith, Z.; Souvorov, A.; Sung, W.; Tang, Z.; Tsuchiya, D.; Tu, H.; Vos, H.; Wang, M.; Wolf, Y. I.; Yamagata, H.; Yamada, T.; Ye, Y.; Shaw, J. R.; Andrews, J.; Crease, T. J.; Tang, H.; Lucas, S. M.; Robertson, H. M.; Bork, P.; Koonin, E. V.; Zdobnov, E. M.; Grigoriev, I. V.; Lynch, M.; Boore, J. L., The ecoresponsive genome of Daphnia pulex. Science 2011, 331, (6017), 555-61. (15) EPA, U. S., Methods for Measuring the Acute Toxicity of Effluents and Receiving Waters to Freshwater and Marine Organisms; EPA-821-R-02-012;. U.S. Environmental Protection Agency, Office of Water (4303T): 1200 Pennsylvania Avenue, NW Washington, DC 2002. (16) Guo, X.; Xu, S.; Yan, X.; Zhou, W.; Dai, X.; Zou, X.; Wang, C.; Wang, D.; Zhao, Y., Cloning and expression of a Chk1 gene in Daphnia pulex during different modes of reproduction. Genes & Genomics 2015, 37, (9), 797-807. (17) Hamilton Martin A., R. R. C., and Thurton Robert V. , Trimmed Spearman-Karber method for estimating median lethal concentrations in toxicity bioassays. Environmental Science & Technology 1977, 11, (7), 714-719. (18) Mu H, S. J., Heras H, Chu K,Qiu JW. , An integrated proteomic and transcriptomic analysis of perivitelline fluid proteins in a freshwater gastropod laying aerial eggs. Journal of proteomics 2017, 155, 22-30. (19) Mu, H.; Sun, J.; Fang, L.; Luan, T.; Williams, G. A.; Cheung, S. G.; Wong, C. K.; Qiu, J. W., Genetic Basis of Differential Heat Resistance between Two Species of Congeneric Freshwater Snails: Insights from Quantitative Proteomics and Base Substitution Rate Analysis. J Proteome Res 2015, 14, (10), 4296-308. (20) Mu H, S. J., Cheung S G, Fang L, Zhou H, Luan T, Zhang H, Wong C K G, Qiu J W Comparative proteomics and codon substitution analysis reveal mechanisms of differential resistance to hypoxia in congeneric snails. Journal

of Proteomics 2018, 172, 36-48.

(21) Mi, H.; Muruganujan, A.; Casagrande, J. T.; Thomas, P. D., Large-scale gene function analysis with the PANTHER classification system. Nat Protoc 2013, 8, (8), 1551-66. (22) Zheng, Q.; Wang, X. J., GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res 2008, 36, (Web Server issue), W358-63. (23) Spanier, K. I.; Leese, F.; Mayer, C.; Colbourne, J. K.; Gilbert, D.; Pfrender, M. E.; Tollrian, R., Predator-induced defences in Daphnia pulex: Selection and evaluation of internal reference genes for gene expression studies with real-time PCR. BMC Molecular Biology 2010, 11, 50-60. (24) Pfaffl, M. W., A new mathematical model for relative quantification in real-time RT-PCR. Nucleic acids research 2001, 29, (9), 2004-2007. (25) Adam, N.; Vergauwen, L.; Blust, R.; Knapen, D., Gene transcription patterns and energy reserves in Daphnia magna show no nanoparticle specific toxicity when exposed to ZnO and CuO nanoparticles. Environ Res 2015, 138, 82-92. (26) Wang, H. J.; Growcock, A. C.; Tang, T. H.; O'Hara, J.; Huang, Y. W.; Aronstam, R. S., Zinc 20

ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625

Environmental Science & Technology

oxide nanoparticle disruption of store-operated calcium entry in a muscarinic receptor signaling pathway. Toxicol In Vitro 2010, 24, (7), 1953-61. (27) Guo, D.; Bi, H.; Wang, D.; Wu, Q., Zinc oxide nanoparticles decrease the expression and activity of plasma membrane calcium ATPase, disrupt the intracellular calcium homeostasis in rat retinal ganglion cells. Int J Biochem Cell Biol 2013, 45, (8), 1849-59. (28) DiPolo, R.; Beauge, L., Sodium/calcium exchanger: influence of metabolic regulation on ion carrier interactions. Physiol Rev 2006, 86, (1), 155-203. (29) Richard M. Denton, J. G. M., Nigel J. Edgell, Role of calcium ions in the regulation of intramitochondrial metabolism. Biochem. J 1980, 190, 107-117. (30) Laurence S. Jouaville, P. P., Carlo Bastianutto, Guy A. Rutter, Rosario Rizzuto, Regulation of mitochondrial ATP synthesis by calcium: Evidence for a long-term metabolic priming. PNAS 1999, 96, (24), 13807-13812. (31) Triboulet, S.; Aude-Garcia, C.; Armand, L.; Gerdil, A.; Diemer, H.; Proamer, F.; Collin-Faure, V.; Habert, A.; Strub, J. M.; Hanau, D.; Herlin, N.; Carriere, M.; Van Dorsselaer, A.; Rabilloud, T., Analysis of cellular responses of macrophages to zinc ions and zinc oxide nanoparticles: a combined targeted and proteomic approach. Nanoscale 2014, 6, (11), 6102-14. (32) Mwaanga, P.; Carraway, E. R.; van den Hurk, P., The induction of biochemical changes in Daphnia magna by CuO and ZnO nanoparticles. Aquat Toxicol 2014, 150, 201-9. (33) Tokishita, S.; Kato, Y.; Kobayashi, T.; Nakamura, S.; Ohta, T.; Yamagata, H., Organization and repression by juvenile hormone of a vitellogenin gene cluster in the crustacean, Daphnia magna. Biochem Biophys Res Commun 2006, 345, (1), 362-70. (34) Chen, S.; Chen, D. F.; Yang, F.; Nagasawa, H.; Yang, W. J., Characterization and processing of superoxide dismutase-fused vitellogenin in the diapause embryo formation: a special developmental pathway in the brine shrimp, Artemia parthenogenetica. Biol Reprod 2011, 85, (1), 31-41. (35) Vandegehuchte, M. B.; Vandenbrouck, T.; De Coninck, D.; De Coen, W. M.; Janssen, C. R., Gene transcription and higher-level effects of multigenerational Zn exposure in Daphnia magna. Chemosphere 2010, 80, (9), 1014-20. (36) Yang, X.; Shao, H.; Liu, W.; Gu, W.; Shu, X.; Mo, Y.; Chen, X.; Zhang, Q.; Jiang, M., Endoplasmic reticulum stress and oxidative stress are involved in ZnO nanoparticle-induced hepatotoxicity. Toxicol Lett 2015, 234, (1), 40-9. (37) Chen, R. L. H., Shi, X.; Bai, R.; Zhang, Z.; Zhao, Y.; Chang, Y.; Chen. C., Endoplasmic reticulum stress induced by zinc oxide nanoparticles is an earlier biomarker for nanotoxicological evaluation. ACS nano 2014, 8, (3), 2562-2574. (38) Wilkinson, B.; Gilbert, H. F., Protein disulfide isomerase. Biochim Biophys Acta 2004, 1699, (12), 35-44. (39) Roger, D.; Kornberg, A.; Jean O. Thomas A. Chromatin Structure: Oligomers of the Histones. Science 1974, (4139), 865-867 (40) Woodard, C. T.; Baehrecke, E. H.; Thummel, C. S. A., Molecular mechanism for the stage specificity of the Drosophila prepupal genetic response to ecdysone. Cell 1994, 79, (4), 607-615. (41) Polverino, A.; Longo, A.; Donizetti, A.; Drongitis, D.; Frucci, M.; Schiavo, L.; Carotenuto, G.; Nicolais, L.; Piscopo, M.; Vitale, E.; Fucci, L., Molecular responses of cells to 2-mercapto-1methylimidazole gold nanoparticles (AuNPs)-mmi: investigations of histone methylation changes. Journal of Nanoparticle Research 2014, 16, (7), 2516. (42) Qian, Y.; Zhang, J.; Hu, Q.; Xu, M.; Chen, Y.; Hu, G.; Zhao, M.; Liu, S., Silver nanoparticle21

ACS Paragon Plus Environment

Environmental Science & Technology

626 627 628 629 630 631 632 633 634 635 636 637 638 639

induced hemoglobin decrease involves alteration of histone 3 methylation status. Biomaterials 2015, 70, 12-22. (43) Poynton, H. C.; Varshavsky, J. R.; Chang, B.; Holman, P. S., Daphnia magna ecotoxicogenomics provides mechanistic insights into metal toxicity. Environ. Sci. Technol 2007, 41, 1044-1050. (44) Dabrunz, A.; Duester, L.; Prasse, C.; Seitz, F.; Rosenfeldt, R.; Schilde, C.; Schaumann, G. E.; Schulz, R., Biological Surface Coating and Molting Inhibition as Mechanisms of TiO2 Nanoparticle Toxicity in Daphnia magna. PLoS ONE 2011, 6, (5), 1-7. (45) Poynton, H. C.; Lazorchak, J. M.; Impellitteri, C. A.; Blalock, B.; Smith, M. E.; Struewing, K.; Unrine, J.; Roose, D., Toxicity and transcriptomic analysis in Hyalella azteca suggests increased exposure and susceptibility of epibenthic organisms to zinc oxide nanoparticles. Environ Sci Technol 2013, 47, (16), 9453-60. (46) Li, W. M.; Wang, W. X., Distinct biokinetic behavior of ZnO nanoparticles in Daphnia magna quantified by synthesizing 65Zn tracer. Water Res 2013, 47, (2), 895-902.

640 641

22

ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

642

Environmental Science & Technology

Figure Captions

643 644

Figure 1 The concentration of dissolved Zn released from ZnO NPs and bulk ZnO at

645

a 1, 10, or 100 mg/L during 24-h of incubation in RMHW. Data are shown

646

as the mean ± SD (n=3).

647 648 649 650 651 652

Figure 2 Number of differentially expressed proteins in D. pulex exposed to ZnO NPs, bulk ZnO and ZnSO4∙7H2O. Figure 3 Venn diagram showing the overlap of differentially expressed proteins of D. pulex exposed to ZnO NPs, bulk ZnO and ZnSO4∙7H2O. Figure 4 Gene ontology (GO) assignment of the shared differentially expressed proteins in D. pulex exposed to ZnO NPs, bulk ZnO and ZnSO4∙7H2O.

653

The meanings in the figure of molecular function were as follows, AC:

654

antioxidant activity (GO:0016209), NABTFA: nucleic acid binding transcription

655

factor activity (GO:0001071), PBTFA: protein binding transcription factor

656

activity (GO:0000988), RA: receptor activity (GO:0004872), TRA: translation

657

regulator activity (GO:0045182).

658

The meanings in the figure of biological process were as follows, BA: biological

659

adhesion (GO:0022610), BG: biological regulation (GO:0065007), CCO: cellular

660

component organization or biogenesis (GO:0071840), DP: developmental process

661

(GO:0032502), ISP: immune system process (GO:0002376), R: reproduction

662

(GO:0000003), RS: response to stimulus (GO:0050896).

663

Figure 5 Relative expression level of h3, l13, scot, idh and eip mRNA of D. pulex

664

exposed to ZnO NPs, bulk ZnO and ZnSO4∙7H2O at 0.06 mg Zn/L

665

(approximately equal to 1/5th of the ZnO NP EC50). Data are shown as the

666

mean ± SD (n = 3). * indicates P < 0.05, ** indicates P < 0.01

667 668 23

ACS Paragon Plus Environment

Environmental Science & Technology

669 670 671 672 673

674 675

Figure 1

24

ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

Environmental Science & Technology

676

Number of proteins

Up-regulation 200 100

125

113

84

0 -100 -200 -300

178 ZnO NPs

218 Bulk ZnO

677 678

Figure 2

679

680 681

Down-regulation

Figure 3 25

ACS Paragon Plus Environment

235 ZnSO4·7H2O

Environmental Science & Technology

682

Cellular component extracellular matrix (GO:0031012), 6.3%

membrane (GO:0016020), 6.3% macromolecular complex (GO:0032991), 12.5%

cell part (GO:0044464), 37.5%

extracellular region (GO:0005576), 16.3%

683 684

organelle (GO:0043226), 21.3%

Molecular function AC, 1.0%

NABTFA, 1.6% TRA, 2.6%

PBTFA, 0.5%

ERA, 4.1%

catalytic activity (GO:0003824), 35.2%

transporter activity (GO:0005215), 9.3% structural molecule activity (GO:0005198), 9.8%

685 686

RA, 9.8%

binding (GO:0005488), 25.9%

Biological process BA, 2.6%

ISP, 2.6%

DP, 2.2%

multicellular organismal process (GO:0032501), 3.9%

CCO, 3.9%

R, 0.4% metabolic process (GO:0008152), 37.9%

RS, 6.0%

BG, 7.8%

687 688

Page 26 of 31

cellular process (GO:0009987), 16.8%

localization (GO:0051179), 15.9%

Figure 4 26

ACS Paragon Plus Environment

Page 27 of 31

Environmental Science & Technology

689 690 691

692 693 694

Figure 5

27

ACS Paragon Plus Environment

Environmental Science & Technology

695 696 697

Table 1. Acute toxicity of ZnO NPs, bulk ZnO and ZnSO4∙7H2O to D. pulex (calculated from 4 replicates)

Treatment ZnO NPs Bulk ZnO ZnSO4·7H2O

24 h toxicity test (mg Zn/L) EC50 95% CI 0.32 0.27-0.39 0.76 0.63-0.92 0.98 0.80-1.18

48 h toxicity test (mg Zn/L)   EC50 95% CI 0.19 0.16-0.23 0.46 0.37-0.56   0.50 0.42-0.60  

698

28

ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

Environmental Science & Technology

699

Table 2. Significantly enriched biological process Gene Ontology (GO) terms for the common differentially expressed proteins of D.

700

pulex exposed to ZnO NPs, bulk ZnO and ZnSO4∙ 7H2O GO ID

Term

Level

Number of overlapped DEPs

Number of proteins P value in database

GO:0006099

tricarboxylic acid cycle

6

4

27

4.72E-03

GO:0006879

cellular iron ion homeostasis

6

3

19

2.50E-02

GO:0015988

energy coupled proton transmembrane transport, against electrochemical gradient

5

5

18

5.95E-05

GO:0015991

ATP hydrolysis coupled proton transport

5

5

18

5.95E-05

GO:0006030

chitin metabolic process

5

9

184

1.28E-03

GO:0006101

citrate metabolic process

5

4

28

5.20E-03

GO:1902600

hydrogen ion transmembrane transport

5

5

55

5.20E-03

GO:0072350

tricarboxylic acid metabolic process

5

4

29

5.80E-03

GO:0019752

carboxylic acid metabolic process

5

12

461

1.27E-02

GO:0043436

oxoacid metabolic process

5

12

468

1.42E-02

GO:0006082

organic acid metabolic process

5

12

489

1.97E-02

GO:0046916

cellular transition metal ion homeostasis

5

3

20

2.73E-02

29

ACS Paragon Plus Environment

Environmental Science & Technology

701

Page 30 of 31

Table 2 (Continued) GO ID

Term

Level

Number of overlapped DEPs

GO:0043648

dicarboxylic acid metabolic process

5

3

20

2.73E-02

GO:0006875

cellular metal ion homeostasis

5

3

25

4.89E-02

GO:0090662

ATP hydrolysis coupled transmembrane transport

4

5

18

5.95E-05

GO:0006022

aminoglycan metabolic process

4

9

196

1.84E-03

GO:0015992

proton transport

4

6

85

4.23E-03

GO:0006812

cation transport

3

11

194

6.24E-05

GO:0044765

single-organism transport

3

20

907

1.44E-03

GO:0015672

monovalent inorganic cation transport

3

8

144

1.45E-03

GO:0006818

hydrogen transport

3

6

85

4.23E-03

GO:0006811

ion transport

3

13

472

4.72E-03

GO:0006826

iron ion transport

3

3

17

1.97E-02

GO:0000041

transition metal ion transport

3

3

22

3.50E-02

GO:0043603

cellular amide metabolic process

3

8

261

4.11E-02

702 30

ACS Paragon Plus Environment

Number of proteins in database

P value

Page 31 of 31

Environmental Science & Technology

Abstract graphic 254x190mm (96 x 96 DPI)

ACS Paragon Plus Environment