Pesticide Mixture Toxicity in Surface Water Extracts in Snails

Feb 22, 2016 - Lijuan Zhao , Yuxiong Huang , Adeyemi S. Adeleye , and Arturo A. Keller .... Tan , Saskia Trump , Susanne Walter-Rohde , John F. Wambau...
0 downloads 0 Views 813KB Size
Subscriber access provided by University of Otago Library

Article

Pesticide mixture toxicity in surface water extracts in snails (Lymnaea stagnalis) by an in vitro acetylcholinesterase inhibition assay and metabolomics Sara Tufi, Pim N.H. Wassenaar, Victoria Osorio, Jacob De Boer, Pim E.G. Leonards, and Marja H. Lamoree Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b04577 • Publication Date (Web): 22 Feb 2016 Downloaded from http://pubs.acs.org on February 28, 2016

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 21

Environmental Science & Technology

85x51mm (150 x 150 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

1

2

Pesticide mixture toxicity in surface water extracts in

3

snails (Lymnaea stagnalis) by an in vitro

4

acetylcholinesterase inhibition assay and metabolomics

5

6

7

Sara Tufi1, Pim N.H. Wassenaar1, Victoria Osorio2, Jacob de Boer1, Pim E.G.

8

Leonards1, Marja H. Lamoree1*

9

10

1

11

1081 HV Amsterdam, The Netherlands

12

2

Institute for Environmental Studies (IVM), VU University Amsterdam, De Boelelaan 1087,

KWR Watercycle Research Institute, Nieuwegein, The Netherlands.

13

14 15 16 17 18 19 20 21 22 23

* Corresponding author: Dr. Marja Lamoree Institute for Environmental Studies (IVM) VU University Amsterdam De Boelelaan 1085 1081 HV Amsterdam The Netherlands e-mail: [email protected] phone: + 31 (0)20 89573 Fax:

+ 31 (0)20 59 89553

1 ACS Paragon Plus Environment

Page 2 of 21

Page 3 of 21

Environmental Science & Technology

24

Abstract

25

Many chemicals in use end up in the aquatic environment. The toxicity of water samples can be tested

26

with bioassays, but a metabolomic approach has the advantage that multiple endpoints can be

27

measured simultaneously and the affected metabolic pathways can be revealed. A current challenge in

28

metabolomics is the study of mixture effects. This study aims at investigating the toxicity of an

29

environmental extract and its most abundant chemicals identified by target chemical analysis of >100

30

organic micro-pollutants and Effect-Directed Analysis (EDA) using the acetylcholinesterase (AChE)

31

bioassay and metabolomics. Surface water from an agricultural area was sampled with a large volume

32

solid phase extraction (LVSPE) device using three cartridges containing a neutral, anionic and cationic

33

sorbents able to trap several pollutants classes like pharmaceuticals, pesticides, PAHs, PCBs and

34

perfluorinated surfactants. Targeted chemical analysis and AChE bioassay were performed on the

35

cartridge extracts. The extract of the neutral sorbent cartridge contained most of the targeted

36

chemicals, mainly imidacloprid, thiacloprid and pirimicarb, and was the most potent AChE inhibitor.

37

Using an EDA approach other AChE inhibiting candidates were identified in the neutral extract, such

38

as carbendazim and esprocarb. Additionally, a metabolomics experiment on the central nervous

39

system (CNS) of the freshwater snail Lymnaea stagnalis was conducted. The snails were exposed to

40

the extract, the three most abundant chemicals individually and a mixture of these. The extract

41

disturbed more metabolic pathways than the three most abundant chemicals individually, indicating

42

the contribution of other chemicals. Most pathways perturbed by the extract exposure overlapped with

43

those related to exposure to neonicotinoids, like the polyamine metabolism involved in CNS injuries.

44

Metabolomics for the straightforward comparison between complex mixture and single compound

45

toxicity is still challenging but, compared to traditional biotesting, is a promising tool due to its

46

increased sensitivity.

47 48

Introduction

49

The number of chemicals produced and used in our modern world is constantly growing. Many of

50

those chemicals, like pesticides, household chemicals and mixtures of various chemicals in industrial

51

waste end up in the environment and can cause a negative impact on the environment1. Chemicals in

52

the aquatic ecosystem can influence the water quality, which may potentially have adverse effects on

53

drinking water quality and biodiversity2. Effect-Directed Analysis (EDA) is one of the strategies to

54

assess water quality. It is proposed to be integrated in the Water Framework Directive (WFD)3. EDA

55

can potentially identify chemicals in complex mixtures that show a toxic response in a bioassay4,5. The

56

bioassays that are often used focus on one specific toxicological endpoint, like estrogenic activity or 2 ACS Paragon Plus Environment

Environmental Science & Technology

57

photosynthetic inhibition6–8. In the EDA approach environmental extracts that show toxicity are

58

fractionated to reduce the complexity of the mixture of chemicals. Those fractions are tested with the

59

bioassay and the active fractions are subjected to chemical analyses in order to identify the chemicals

60

that cause the toxicity. Usually, chromatographic techniques such as gas chromatography (GC) or

61

liquid chromatography (LC) coupled to mass spectrometry (MS) are applied to identify these

62

chemicals9,10. A commonly used bioassay to determine a toxic effect induced by pesticide exposure is

63

the Ellman’s acetylcholinesterase (AChE) bioassay, which is linked to organophosphate (OP) and/or

64

carbamate pesticides. These classes of pesticides inhibit the enzyme AChE11 which led to the

65

accumulation of acetylcholine in the synaptic cleft, inducing an over-stimulation of the neuronal cells.

66

As a consequence, the cholinergic function in the central nervous system (CNS), the peripheral

67

nervous system and the motor components are impaired. However, pesticides like the neonicotinoids

68

cannot be detected with the AChE bioassay since they act as acetylcholine receptor agonists12.

69

Currently, there is a gap between the chemicals mode of action and the bioassay designed to test

70

chemicals toxicity. As an example, for testing neonicotinoids toxicity a bioassay has not been designed

71

yet. Therefore, even if a battery of bioassays is applied to test the quality of water samples, the toxicity

72

of some classes of compounds will never be detected. An alternative strategy to test water quality and

73

overcome the major pitfall of current biotesting strategies is metabolomics11,12. With this approach

74

multiple endpoints (metabolites) can be measured simultaneously after exposing the whole organism

75

to an environmental sample, thus enabling the study of the metabolism in a more realistic

76

environmental exposure scenario13. Many studies in environmental metabolomics have investigated

77

the effect of single compounds or a simple mixture of compounds14–18. One of the current challenges in

78

environmental metabolomics is the study of complex chemical mixture effects, like an environmental

79

sample. In this study the toxicological effects of an environmental water extract and its most abundant

80

chemicals were investigated by applying an EDA using the AChE bioassay, and metabolomics on the

81

central nervous system (CNS) of the freshwater snail Lymnaea stagnalis.

82

Surface water was sampled in an agricultural area of the Netherlands using a large solid phase

83

extraction (LVSPE) device was equipped with three cartridges placed in series containing neutral,

84

anionic and cationic sorbents. Target chemical analysis on 119 chemicals was performed on the

85

extracts of the three cartridges. Since the sampling location was known to be affected by pesticide

86

pollution, the AChE inhibition potency was tested which is commonly used to detect toxicity of

87

organophosphate and carbamate pesticides. To identify other chemical in the environmental extract, an

88

EDA study was carried out using an on-line high pressure liquid chromatography - time of flight

89

(HPLC-ToF) micro-fractionation system. Since the AChE bioassay can target only certain classes of

90

chemicals, targeted and non-targeted metabolomics was used as alternative test of contaminant 3 ACS Paragon Plus Environment

Page 4 of 21

Page 5 of 21

Environmental Science & Technology

91

mixtures to evaluate the toxic effects of the environmental extract. The effects on the snail L. stagnalis

92

CNS were studied since most pesticides are affecting the nervous system. The snails were exposed to

93

the extract of the neutral sorbent cartridge, to the three most abundant chemicals individually, and to a

94

mixture of these chemicals. The biological interpretation was facilitated by the use of biochemical

95

mapping to identify the correlation between the surface water toxicity and its most abundant

96

chemicals.

97

4 ACS Paragon Plus Environment

Environmental Science & Technology

98

Materials and Methods

99

Water sampling and extraction

100

Forty-five L of surface water was sampled at a location in the province Zuid-Holland in the

101

Netherlands (geographical coordinates of 52°00'08.5"N and 4°17'35.2"E). The sampling was carried

102

out in a ditch close to an agricultural area. Chemicals in this surface water were extracted by a LVSPE

103

device, in which three different filter cartridges were placed in series. The first cartridge contained a

104

hydrophobic polystyrene-divinylbenzene copolymer (HR-X) (neutral sorbent); the second contained a

105

weak anion exchanger (HR-XAW) and the last cartridge contained a cationic exchanger (HR-XCW).

106

The choice of these sorbents has been made in order to increase the retention capability of different

107

classes of chemicals from neutral compounds to strong acids with a pKa < 1 and strong bases with pKa

108

> 10. In the manuscript, the neutral extract refers to the extract from the neutral sorbent cartridge

109

whereas cationic and anionic extract correspond to the extracts from the cationic and anionic sorbent

110

cartridges, respectively. The chemicals in the neutral cartridge were eluted with 200 mL methanol

111

(MeOH):ethylacetate 50:50 v/v, the chemicals in the anionic cartridge with 100 mL MeOH 2%

112

ammonia 7N, and the chemicals in the cationic cartridge with 100 mL 1.7% formic acid (Fluka,

113

Steinheim, Germany). As blank, 45 L of bottled water (Gerolsteiner, Gerolstein, Germany), underwent

114

the same processes.

115

Each cartridge extract was divided into aliquots for the different experiments. For the AChE bioassay

116

2 mL of each extract were concentrated by evaporation to 0.2 mL, corresponding to a 150 times higher

117

concentration than the original water sample (environmental concentration factor, ECF). The same

118

process was used for the EDA study. For the metabolomics experiment, 84 mL of the extract were

119

dissolved in 1.2 L of copper free water suitable for the snail exposure. The final concentration in the

120

snail exposure media corresponded to 16 ECF.

121

Target chemical analysis

122

A list of 119 emerging organic micro-pollutants (see Table S1) was selected for chemical screening of

123

the surface water extracts. The selection criteria of the candidate substances were based on: (i)

124

occurrence in surface water, (ii) inclusion of a broad range of physicochemical properties and (iii)

125

biological activity. This long list was classified according to their different uses, namely:

126

pharmaceuticals and metabolites (55), herbicides and metabolites (32), insecticides (7), fungicides (2),

127

flame retardants (5), plasticizers (1), industrial chemicals (14), and artificial sweeteners (3). This target

128

chemical analysis was conducted by using an Orbitrap interfaced with a HPLC-system (Thermo

129

Electron GmbH, Bremen, Germany) (see Supporting information). 5 ACS Paragon Plus Environment

Page 6 of 21

Page 7 of 21

Environmental Science & Technology

130

Acetylcholinesterase bioassay

131

The extracts of the neutral, anionic and cationic cartridges were tested with regard to AChE inhibition

132

with the AChE bioassay. As a first step glycerol (purity ≥ 99%, Sigma-Aldrich) was added to all

133

extracts, at a final concentration of 50 mg/mL. To the cationic extracts an additional 10% of

134

ammonium hydroxide was added to buffer the pH. The extracts were dried with a gentle nitrogen flow

135

and reconstituted with potassium-phosphate buffer (0.1 M, 7.5 pH) in such an amount that the extracts

136

were concentrated with a factor of 10. In addition, three procedural blanks consisting of the same

137

solvents as the cartridge extracts were prepared to check the effect of the solvent on AChE inhibition.

138

The bioassay procedure is described in the Supporting information.

139

Effect-directed analysis (EDA)

140

The most active cartridge extract and its corresponding blank were fractionated in an EDA approach to

141

identify the causative compounds of the AChE inhibition. Before fractionation, the neutral and the

142

neutral blank extracts were dried completely with a gentle nitrogen flow and subsequently

143

reconstituted with water (H2O):acetonitrile (ACN) 90:10 v/v. These two extracts were fractionated

144

using an on-line HPLC-micro-fractionation system with a ToF-MS in parallel to identify the

145

compounds of interest. The experimental setup is described in the Supporting information. Analytical

146

standards were used to calibrate the fractionation system (Table S2).

147

L. stagnalis exposure

148

L. stagnalis snails were taken from a synchronized population and cultured at the VU University

149

Amsterdam, The Netherlands as described by Tufi et al. (2015)19. The neutral and neutral blank

150

extracts were used to expose the snails. 84 mL of the extract the blank were completely dried using a

151

Buchi syncore coupled to a BUCHI vacuum pump V-700 with vacuum controller V885 and a Buchi

152

recirculating chiller F-108 (Buchi, Flawil, Switzerland) for 3 h at 240 mbar with 300 rpm, a rack

153

temperature of 50 °C and lid temperature of 60 °C, followed by 1 h at 120 mbar. Afterwards, the

154

extracts were diluted with copper free water to a concentration factor of 16 ECF.

155

Besides the neutral extract and the neutral blank, the snails were also exposed to i) the three most

156

abundant chemicals individually at concentrations corresponding to the neutral extract (imidacloprid,

157

thiacloprid, and pirimicarb which were exceeding the normative levels in the environmental extract),

158

and ii) to a mixture of them at the same concentration levels. In addition, a control group was assessed

159

without an exposure. Imidacloprid, thiacloprid, and pirimicarb were dissolved in copper free water.

160

Each exposure group contained seven snails. Each snail was separately placed in a beaker glass with

161

150 mL copper free water and fed with 250 µL Tetraphyll (133 g/L in copper free water; Tetra, Melle, 6 ACS Paragon Plus Environment

Environmental Science & Technology

162

Germany). After 2 days of conditioning, the snails were placed in a new beaker glass with 150 mL of

163

the exposure medium. After 48 hours the snails were put in an aluminium box (Sanbio, Uden, The

164

Netherlands), snap frozen in liquid nitrogen and the CNSs dissected. The CNS samples were stored in

165

0.5 mL Precellys vials (Bertin Technologies, France) containing ± 10 ceramic beads (1.4 mm,

166

zirconium oxide; MO Bio Laboratories, CA, USA) and snap frozen in liquid nitrogen. Subsequently,

167

the samples were stored at -80 °C.

168

Actual exposure concentrations

169

Aliquots of the exposure media were collected at the start of each experiment and analysed with LC-

170

MS to determine actual exposure concentrations of pirimicarb, imidacloprid and thiacloprid (see

171

Supporting information and Tables S3-5).

172

Targeted and non-targeted metabolomics

173

The CNS samples were processed individually and analyzed with the non-targeted and targeted

174

metabolomics approaches described by Tufi et al. (2015)19,20. The extracted metabolites were analyzed

175

by using hydrophilic interaction chromatography (HILIC) coupled to a ToFMS (Bruker Daltonics).

176

For the non-targeted platform >500 standards from the MS Metabolite Library of Standards (MSMLS;

177

IROA Technologies, MI, USA) were analyzed to determine the retention time and mass spectrum. The

178

chromatograms were processed with Bruker Daltonics Compass PathwayScreener (Bruker Daltonics)

179

to perform batch targeted analyses by screening for metabolites of the MSMLS library. Multiple

180

hypotheses testing was performed using the areas of the detected metabolites and Bruker Daltonics

181

Compass ProfileAnalysis (Bruker Daltonics). A t-test analysis was conducted on areas of the detected

182

metabolites comparing two different exposure groups and FDR correction (p-value < 0.05, t-test).

183

Biochemical networks

184

From the detected metabolites biochemical networks were built with Metamapp setting the threshold

185

at 0.7 for chemical similarity21. The formed network, based on chemical similarities and biochemical

186

interactions, was uploaded in Cytoscape (v3.1.0) in which metabolic changes between the control

187

group and different exposure groups were visualized22. Fold changes were calculated by dividing the

188

average levels of metabolites in exposed groups by the average levels of metabolites in the control

189

group. Visualization was based on fold changes and p-values. Exposure to the neutral extract was

190

compared to the neutral blank extract, and exposures to the individual chemicals and the mixture of

191

chemicals were compared to the control group. Biological interpretation of the results was performed

192

using the small molecular databases Kyoto Encyclopedia of Genes and Genomes (KEGG), the Human

193

Metabolome Database, the Small Molecule Pathway Database and the Biocyc database collection.

194 7 ACS Paragon Plus Environment

Page 8 of 21

Page 9 of 21

Environmental Science & Technology

195

Results and Discussions

196

Selection of the most active environmental extract

197

Forty chemicals were quantified in the extracts (Table S6) and the concentrations of those compounds

198

detected above the level of 10 ng/L are shown in Table 1. The highest concentrations were found in

199

the neutral extract suggesting a higher retention of this sorbent for organic compounds. High

200

concentrations of some chemicals were also found in the anionic extract indicating a potential

201

cartridge breakthrough due to the overloaded capacity of the first neutral sorbent.

202

The most abundant chemicals in the neutral extract were the carbamate pesticide pirimicarb (0.27

203

µg/l), and the two neonicotinoid pesticides imidacloprid (0.21 µg/l) and thiacloprid (0.16 µg/l) (Table

204

1). Their concentrations exceed the normative levels. Pirimicarb exceeded the maximum permissible

205

risk (MPR) concentration 4-fold, while the concentration of imidacloprid was four times higher than

206

the annual average environmental quality standard (AA-EQS), and thiacloprid exceeded the AA-EQS

207

up to 20-fold, when the concentrations quantified in the neutral, cationic and anionic extracts were

208

summed.

209

The high concentrations of pesticides are consistent with the sampling location which is an agricultural

210

area with greenhouses. From other studies it is well known that pesticides can be transferred through

211

runoff waters into surface waters23. Our results showed comparable levels as found in monitoring

212

programs at this location that are reported in the Pesticides Atlas of The Netherlands24. Thiacloprid

213

was not detected in the latest available data (2013) from these monitoring programs, suggesting that

214

thiacloprid has been widely used in that area later.

8 ACS Paragon Plus Environment

Environmental Science & Technology

215 216 217 218

Page 10 of 21

Table 1. Results of the target chemical analysis, showing the most abundant chemicals identified. The columns 4-9 show the concentrations of the compounds in the extracts in µg/L (ND = not detected). Ins: Insectides, Ph: Pharmaceuticals, Herb: Herbicides, Herb M: Herbicide metabolites, Fung: Fungicides, FR: Flame Retardants, OC: Organic compounds, IC: Industrial chemicals. AA-EQS: annual average environmental quality standard, MAC-EQS: maximum allowable concentration-EQS, MPR: (.

Compound

Class

Pirimicarb

Ins

Imidacloprid

Ins

Thiacloprid

Ins

Caffeine

Ph

MCPA

Herb

Carbendazim

Fung

2.4-Dinitrophenol

Herb

4.6-dinitro-o-cresol

Herb

2.6-Dichloro benzamide Triethyl phosphate

Herb. M. FR

Mecoprop

Herb

Metalaxyl-M HMMM 4-Methyl-1H-benzotriazole Paracetamol

Fung OC IC Ph

Norm levels (µg/L)1 0.090 (MPR) 0.067 (AA) 0.2 (MAC) 0.01 (AA) 0.11 (MAC) ND 1.4 (AA) 15 (MAC) 0.6 ( AA) 0.6 (MAC) 0.001 (MPR) 9.2 (AA) 9.2 (MAC) 1000 (MPR) 440 (MPR) 18 (AA) 160 (MAC) 9.7 (MPR) 810 (MPR) ND ND

Neutral Extract (µg/L)

Neutral Blank (µg/L)

Anionic Extract (µg/L)

Anionic Blank (µg/L)

Cationic Extract (µg/L)

Cationic blank (µg/L)

0.27

ND

0.099

ND

0.001

ND

0.21

ND

0.063

ND

0.001

ND

0.16

ND

0.044

ND

ND

ND

0.052

0.001

0.012

ND

0.004

0.003

0.034

ND

0.032

ND

0.001

ND

0.031

ND

0.017

ND

ND

ND

0.017

ND

0.023

ND

ND

ND

0.017

ND

0.013

ND

ND

ND

0.015 0.015

ND 0.7

0.006 0.005

ND 0.6

0.002 0.009

ND 0.009

0.014

ND

0.011

ND

0.1

ND

0.010 0.009 0.008 0.005

ND ND ND ND

0.003 0.005 0.005 0.005

ND 0.1 ND ND

ND ND ND 0.003

ND ND ND ND

9 ACS Paragon Plus Environment

Page 11 of 21

Environmental Science & Technology

219

The neutral extract showed the inhibitoriest potential in the AChE assay (Figure 1).

220

These results are in accordance with the target chemical analysis which showed that the neutral extract

221

contained most chemicals. A surprising result was obtained for the anionic extract, which expressed

222

AChE stimulation at 15 and 60 ECF. This trend was also observed for the cationic extract, however, not

223

statistically different from the blank. In the literature AChE inhibiting effects have been described

224

extensively. However, AChE stimulating effects are seldom reported. A study conducted by López et al.

225

(2015) described an AChE activating effect of monoterpenoids, plant oils, which showed an AChE

226

stimulation effect at low concentrations (0.04 mM) but an AChE inhibiting effect at high concentrations

227

(5 mM)25. There is a possibility that the anionic and cationic extracts contain chemicals which exert such

228

effects. Another possible explanation for the observed AChE stimulation might be the presence of humic

229

substances, like humic and fulvic acids, which are major components of natural organic matter in soil and

230

water. These acids have a light yellow to dark brown colour26 and the yellow colour was indeed observed

231

in the anionic extract and might therefore indicate the presence of humic substances. Nevertheless, it has

232

been shown that humic substances, co-extracted by solid-phase extraction, can have a suppressing effect

233

on the AChE inhibition assay27. To verify the potential interference of the matrix we performed the AChE

234

bioassay with the extract and the reagents but without the enzyme. Since we did not observe any change

235

in AChE activity, we could exclude a matrix effect in the AChE bioassay.

236

Effect-Directed Analysis

237

The target chemical analysis showed that pirimicarb and the fungicide carbendazim could be responsible

238

for the observed AChE inhibition. These chemicals are mainly interacting with mitotic division, but they

239

also have AChE inhibiting properties28. In order to identify other potential candidates that contributed to

240

the AChE inhibition in the neutral extract, the neutral extract and the blank underwent an EDA

241

fractionation study coupled to the AChE bioassay. A screening approach was used to identify chemical in

242

the fractions. The most active fractions (well position 12, 13, 15, 19, 23, 36 and 37) of the neutral extract

243

are shown in Figure S1 and the identified chemicals in the active fractions are shown in Table S8. Most of

244

the identified compounds in the active fractions belong to the class of carbamate pesticides. The pesticide

245

aldicarb was found in fraction 13, pirimicarb in fraction 19, and esprocarb in fractions 36 and 37. These

246

identified pesticides were confirmed by mass accuracy, chemical formula and isotope distribution

247

matches. To provide an additional level of confidence in these compounds identification, analytical

248

standards were bought and the retention time has been compared. Even though these chemicals could

249

potentially explain the AChE inhibition of the neutral extract, the identification of the chemicals in the

250

active fractions 15 and 23 was not successful.

251

Metabolomics of the neutral surface water extract, its most abundant chemicals, and mixture 10 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 21

252

Snails were exposed to: (i) the neutral extract, (ii) the three most abundant chemicals (imidacloprid,

253

pirimicarb and thiacloprid) individually and (iii) a mixture of these chemicals at the same levels as the

254

individual concentrations. The measured exposure concentrations of the metabolomics experiments are

255

shown in Table S6. The results show that the imidacloprid concentrations were similar among the

256

exposure groups (Table S6). Thiacloprid concentration in the mixture was about 20% higher than in the

257

neutral extract and the individual exposures, whereas pirimicarb was 30% lower than in the neutral extract

258

and 20% lower than in the mixture. The neutral extract induced many alterations at the molecular level in

259

the CNS: the levels of eighteen metabolites were found to be altered. Even though some altered pathways

260

could be correlated to the three most abundant chemicals, the overall metabolic change could not be

261

explained by the simplified mixture. The environmental extract perturbed more metabolic pathways,

262

indicating therefore the contribution of other chemicals and possibly synergistic effects to the observed

263

toxicity.

264

To further investigate the results of the metabolomics experiments, the data were visualized in

265

biochemical networks to identify affected metabolite groups and pathways. The comparison between the

266

neutral extract and the neutral blank extract showed 18 statistically significant metabolites of which most

267

were down regulated (± 80%) (Figure 2). Exposure to the single pesticides, on the other hand, showed

268

fewer statistically significant differences and in the mixture exposure there were only four significantly

269

changed metabolites (Figure S2-S5).

270

Exposure to pirimicarb showed to disturb only few metabolites including methionine, 4-hydroxy-L-

271

proline, inosine and hypoxanthine (ca. 3 fold increase). Increased hypoxanthine levels were related to

272

increased oxygen radicals29. A study conducted by Wang et al. (2014) showed that pirimicarb exposure

273

caused increased oxidative stress, which might be associated to the observed increased level of

274

hypoxanthine. Furthermore, pirimicarb was related to a disturbed energy metabolism30. The interpretation

275

of the metabolite biochemical role has been based on the mammalian metabolic pathways reported in

276

KEGG. This database is increasing the number of annotated species but the metabolism of the snail

277

Lymnaea stagnalis has not been reported yet. Therefore some discrepancies might raise from the direct

278

comparison of diverse species metabolisms, even though, metabolites are highly conserved among

279

different species13.

280

The metabolite patterns in the CNS of L. stagnalis exposed to imidacloprid or thiacloprid are expected to

281

be comparable, since these chemicals have the same mode of action and therefore their effect is supposed

282

to be additive (following the concentration addition model). Therefore, it was expected that they would

283

affect the same metabolites. Thiacloprid and imidacloprid indeed showed many similar metabolic

284

alterations. They decreased the levels of the metabolites ornithine and lysine. Ornithine is involved in 11 ACS Paragon Plus Environment

Page 13 of 21

Environmental Science & Technology

285

multiple metabolic pathways, including the spermidine and spermine biosynthesis, and in arginine and

286

proline metabolism. Imidacloprid showed a significant decrease in spermidine, proline and arginine

287

levels, which indicates involvement of these pathways. A study conducted by Moser et al. (2015), who

288

exposed rats to imidacloprid, showed the same altered pattern for ornithine, proline and spermidine.

289

Furthermore, they showed that citrulline, also involved in the arginine and proline metabolism, was

290

decreased as well31. This supports the hypothesis that neonicotinoids play a role in these pathways.

291

Although imidacloprid and thiacloprid have overlapping metabolic patterns, some differences were also

292

observed. For instance, glucosamine, involved in amino sugar metabolism, was up-regulated after

293

imidacloprid exposure and downregulated after thiacloprid exposure. Glucosamine is related to chitin

294

metabolism and is for example used as building block for the exoskeletons of arthropods, like insects and

295

crustaceans, the shells of mollusks and the cell walls of fungi32. Dondero et al. (2010) showed that

296

thiacloprid exposure, in comparison to imidacloprid, resulted in severely decreased chitinase mRNA in

297

the marine mussel Mytilus galloprovincialis, which is also involved in chitin metabolism32.

298

The exposures to the three single pesticides induced changes in a higher number of metabolites compared

299

with the alteration shown after the exposure to their mixture. Indeed, only four metabolites were

300

significantly different after the mixture exposure. For example, the decrease of lysine, observed for all

301

three single compounds, was levelled out in the mixture. This effect has been already observed in other

302

studies as well30,32 and might be caused by interaction effects33. These interactions might be due to the

303

internal dose or bioavailable concentration of the chemical by its toxicokinetics (absorption, distribution,

304

metabolism and excretion) and/or to the chemical toxicodynamics through the binding of one of more

305

chemicals to a receptor through which toxicity may be mediated.

306

Moreover, it has been shown that chemicals with the same mode of action did not show additive effects

307

when mixed, suggesting that these chemicals have non-similar effects, probably due to differences in

308

compound-specific toxicodynamics14,32,34.

309

The neutral extract exposure showed many significantly changed metabolites and most of them were not

310

found after the single pesticide exposures indicating that other chemicals in this extract also play a big

311

role in the changes. Only few metabolites were also found to change significantly after the individual

312

pesticide exposures. For instance, the serotonin metabolite 5-hydroxyindoleacetic acid (5-HIAA; ca. 3

313

fold decrease) was decreased with the same fold change as after the thiacloprid exposure. The metabolites

314

guanosine, adenosine, adenine, inosine and hypoxanthine were affected by the individual compounds as

315

well. These metabolites are involved in purine nucleotide degradation, which has been linked to

316

lymphocyte deficiencies35, and lymphocyte deficiencies have been related to thiacloprid exposure32.

317

Furthermore, imidacloprid has shown to affect purine nucleotide degradation by increasing the activity of 12 ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 21

318

xanthine oxidase, converting hypoxanthine to xanthine. This reaction is responsible for the release of

319

oxygen radicals, which can cause additional damage and disturbances in tissues and organs36. Pirimicarb

320

showed to disturb the purine nucleotide metabolism by increasing hypoxanthine.

321

The neutral extract also affected spermidine and spermine biosynthesis, which was also affected by

322

imidacloprid and thiacloprid but other metabolites changed in this pathway (ornithine and spermidine

323

instead of 5’-methylthioadenosine, N(1)-acetylspermine, N-acetylputrescine and putrescine) (Figure 3).

324

It seemed that chemicals in the neutral extract were blocking the putrescine break down, resulting in an

325

increased putrescine level. Putrescine can react with S-adenosylmethioninamine to form spermidine and

326

5’-methylthioadenosine, or can be degraded to N-acetylputrescine. Spermidine, 5’-methylthioadenosine

327

and N-acetylputrescine showed decreased fold changes, supporting the hypothesis that putrescine break

328

down is inhibited. Furthermore, N(1)-acetylspermine, a break down product of spermine or a building

329

block of spermidine, was decreased as well, which might indicate a shift in the formation of spermidine

330

through a different route. Increases in putrescine levels have been related to pathological changes that

331

cause cell injury in the CNS, like severe metabolic stress, exposure to neurotoxins and seizure37. These

332

increased putrescine levels are mainly caused by an increased activity of the enzyme ornithine

333

decarboxylase (ODC), which degrades ornithine into putrescine37–39. Furthermore, S-adenosylmethionine

334

(SAM) and the enzyme S-adenosylmethionine decarboxylase which transforms SAM into S-

335

adenosylmethioninamine (dcSAM) might be depleted or inhibited, respectively, after CNS injury40.

336

Additionally, most CNS injury cases39,40 induced unaltered or only slightly changed in the levels of

337

spermine and spermidine.

338

The observed effects of the neutral extract on the spermidine and spermine biosynthesis might partially be

339

explained by the neonicotinoids imidacloprid and thiacloprid. The observed effects of the neutral extract

340

on the metabolites histidine, kynurenic acid, trigonelline and cytosine cannot be assigned to one of the

341

three most abundant chemicals and, therefore, these effects were probably caused by other chemicals in

342

the extract. For instance, decreased kynurenic acid levels have been related to pyrethroid exposure, which

343

might be present in the extract41. Furthermore, paracetamol, which was detected in the neutral extract, has

344

been related to decreased trigonelline levels as well as decreased SAM42.

345

In conclusion, the observed toxicity of the neutral extract in the metabolomics experiment could not be

346

explained entirely by the most abundant chemicals in the extract, although some metabolic patterns were

347

overlapping with the neutral extract exposure. Nevertheless, the use of metabolomics in environmental

348

extract exposure studies is very promising since it showed to be more sensitive than traditional toxicity

349

test, capable to focus on many endpoints simultaneously and able to provide information on the pathways

350

affected. The use of an aquatic model organism requiring a smaller volume of water, like Daphnia 13 ACS Paragon Plus Environment

Page 15 of 21

Environmental Science & Technology

351

magna, can be advantageous for further research. With such organisms dose response effects of water

352

extracts can be studied, even with small amounts of water extract.

353 354

Acknowledgments

355

This study was carried out within the Marie Curie Research Training Network EDA-EMERGE

356

(www.eda-emerge.eu) supported by the EU (MRTN-CT-2012-290100).

357 358

Supporting Information

359

This information is available free of charge via the Internet at http://pubs.acs.org.

360 361 362

References

363 364

(1)

Daughton, C. G. “Emerging” chemicals as pollutants in the environment: A 21st century perspective. Renew. Resour. J. 2005, 23, 6–23.

365 366

(2)

European Environmental Agency. 10 messages for 2010 Freshwater ecosystems; Copenhagen, 2010.

367 368 369 370

(3)

Brack, W.; Govender, S.; Schulze, T.; Krauss, M.; Hu, M.; Muz, M.; Hollender, J.; Schirmer, K.; Schollee, J.; Hidasi, A.; et al. EDA-EMERGE: an FP7 initial training network to equip the next generation of young scientists with the skills to address the complexity of environmental contamination with emerging pollutants. Environ. Sci. Eur. 2013, 25 (1), 18.

371 372

(4)

Brack, W. Effect-directed analysis: a promising tool for the identification of organic toxicants in complex mixtures? Anal. Bioanal. Chem. 2003, 377 (3), 397–407.

373 374 375

(5)

Weiss, J. M.; Hamers, T.; Thomas, K. V.; van der Linden, S.; Leonards, P. E. G.; Lamoree, M. H. Masking effect of anti-androgens on androgenic activity in European river sediment unveiled by effect-directed analysis. Anal. Bioanal. Chem. 2009, 394 (5), 1385–1397.

376 377 378 379

(6)

Simon, E.; van Velzen, M.; Brandsma, S. H.; Lie, E.; Løken, K.; de Boer, J.; Bytingsvik, J.; Jenssen, B. M.; Aars, J.; Hamers, T.; et al. Effect-Directed Analysis To Explore the Polar Bear Exposome: Identification of Thyroid Hormone Disrupting Compounds in Plasma. Environ. Sci. Technol. 2013, 47 (15), 8902–8912.

380 381 382 383

(7)

Booij, P.; Vethaak, A. D.; Leonards, P. E. G.; Sjollema, S. B.; Kool, J.; de Voogt, P.; Lamoree, M. H. Identification of photosynthesis inhibitors of pelagic marine algae using 96-well plate microfractionation for enhanced throughput in effect-directed analysis. Environ. Sci. Technol. 2014, 48 (14), 8003–8011. 14 ACS Paragon Plus Environment

Environmental Science & Technology

384 385 386

(8)

Fang, Y.-X.; Ying, G.-G.; Zhao, J.-L.; Chen, F.; Liu, S.; Zhang, L.-J.; Yang, B. Assessment of hormonal activities and genotoxicity of industrial effluents using in vitro bioassays combined with chemical analysis. Environ. Toxicol. Chem. 2012, 31 (6), 1273–1282.

387 388

(9)

Brack, W. Effect-Directed Analysis in Ecotoxicology. In Encyclopedia of Aquatic Ecotoxicology; Férard, J.-F., Blaise, C., Eds.; Springer Netherlands: Dordrecht, 2013; pp 387–394.

389 390 391 392

(10)

Weiss, J. M.; Simon, E.; Stroomberg, G. J.; de Boer, R.; de Boer, J.; van der Linden, S. C.; Leonards, P. E. G.; Lamoree, M. H. Identification strategy for unknown pollutants using highresolution mass spectrometry: Androgen-disrupting compounds identified through effect-directed analysis. Anal. Bioanal. Chem. 2011, 400 (9), 3141–3149.

393 394

(11)

Ellman, G. L.; Courtney, K. D.; Andres jr., V.; Featherstone, R. M. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961, 7 (2), 88–95.

395 396

(12)

Buckingham, S.; Lapied, B.; Corronc, H.; Sattelle, F. Imidacloprid actions on insect neuronal acetylcholine receptors. J. Exp. Biol. 1997, 200 (21), 2685–2692.

397 398

(13)

Bundy, J. G.; Davey, M. P.; Viant, M. R. Environmental metabolomics: a critical review and future perspectives. Metabolomics 2009, 5 (1), 3–21.

399 400 401

(14)

Baylay, A. J.; Spurgeon, D. J.; Svendsen, C.; Griffin, J. L.; Swain, S. C.; Sturzenbaum, S. R.; Jones, O. a. H. A metabolomics based test of independent action and concentration addition using the earthworm Lumbricus rubellus. Ecotoxicology 2012, 21 (5), 1436–1447.

402 403 404

(15)

Kluender, C.; Sans-Piché, F.; Riedl, J.; Altenburger, R.; Härtig, C.; Laue, G.; Schmitt-Jansen, M. A metabolomics approach to assessing phytotoxic effects on the green alga Scenedesmus vacuolatus. Metabolomics 2008, 5 (1), 59–71.

405 406 407

(16)

Ye, Y.; Wang, X.; Zhang, L.; Lu, Z.; Yan, X. Unraveling the concentration-dependent metabolic response of Pseudomonas sp. HF-1 to nicotine stress by 1H NMR-based metabolomics. Ecotoxicology 2012, 21 (5), 1314–1324.

408 409 410

(17)

Jones, O. A. H.; Swain, S., C.; Svendsen, C.; Griffin, J. L.; Sturzenbaum, S. R.; Spurgeon, D. J. Potential New Method of Mixture Effects Testing Using Metabolomics and Caenorhabditis elegans. J. Proteome Res. 2012, 11 (2), 1446–1453.

411 412 413 414

(18)

Collette, T. W.; Teng, Q.; Jensen, K. M.; Kahl, M. D.; Makynen, E. A.; Durhan, E. J.; Villeneuve, D. L.; Martinović-Weigelt, D.; Ankley, G. T.; Ekman, D. R. Impacts of an Anti-Androgen and an Androgen/Anti-Androgen Mixture on the Metabolite Profile of Male Fathead Minnow Urine. Environ. Sci. Technol. 2010, 44 (17), 6881–6886.

415 416 417

(19)

Tufi, S.; Lamoree, M. H.; De Boer, J.; Leonards, P. E. G. Cross-platform metabolic profiling: application to the aquatic model organism Lymnaea stagnalis. Anal. Bioanal. Chem. 2015, 407 (7), 1901–1912.

418 419 420

(20)

Tufi, S.; Lamoree, M.; de Boer, J.; Leonards, P. Simultaneous analysis of multiple neurotransmitters by hydrophilic interaction liquid chromatography coupled to tandem mass spectrometry. J. Chromatogr. A 2015.

15 ACS Paragon Plus Environment

Page 16 of 21

Page 17 of 21

Environmental Science & Technology

421 422 423 424

(21)

Barupal, D. K.; Haldiya, P. K.; Wohlgemuth, G.; Kind, T.; Kothari, S. L.; Pinkerton, K. E.; Fiehn, O. MetaMapp: mapping and visualizing metabolomic data by integrating information from biochemical pathways and chemical and mass spectral similarity. BMC Bioinformatics 2012, 13 (1), 99.

425 426

(22)

Lopes, C. T.; Franz, M.; Kazi, F.; Donaldson, S. L.; Morris, Q.; Bader, G. D. Cytoscape Web: an interactive web-based network browser. Bioinformatics 2010, 26 (18), 2347–2348.

427 428 429

(23)

Cerejeira, M. J.; Viana, P.; Batista, S.; Pereira, T.; Silva, E.; Valério, M. J.; Silva, A.; Ferreira, M.; Silva-Fernandes, A. M. Pesticides in Portuguese surface and ground waters. Water Res. 2003, 37 (5), 1055–1063.

430 431 432

(24)

de Snoo, G. R.; Tamis, W. L. M.; Vijver, M. G.; Musters, C.; van ’t Zelfde, M. Risk mapping of pesticides: the Dutch atlas of pesticide concentrations in surface waters: www.pesticidesatlas.nl. Commun. Agric. Appl. Biol. Sci. 2006, 71 (2 Pt A), 49–58.

433 434 435 436

(25)

López, M. D.; Campoy, F. J.; Pascual-Villalobos, M. J.; Muñoz-Delgado, E.; Vidal, C. J. Acetylcholinesterase activity of electric eel is increased or decreased by selected monoterpenoids and phenylpropanoids in a concentration-dependent manner. Chem. Biol. Interact. 2015, 229, 36– 43.

437 438

(26)

Humus Chemistry: Genesis, Composition, Reactions, Second Edition (Stevenson, F. J.). J. Chem. Educ. 1995, 72 (4), A93.

439 440

(27)

Neale, P. A.; Escher, B. I. Coextracted dissolved organic carbon has a suppressive effect on the acetylcholinesterase inhibition assay. Environ. Toxicol. Chem. SETAC 2013, 32 (7), 1526–1534.

441 442 443

(28)

Palanikumar, L.; Kumaraguru, A. K.; Ramakritinan, C. M.; Anand, M. Toxicity, biochemical and clastogenic response of chlorpyrifos and carbendazim in milkfish Chanos chanos. Int. J. Environ. Sci. Technol. 2013, 11 (3), 765–774.

444 445

(29)

Rieger, D. Relationships between energy metabolism and development of early mammalian embryos. Theriogenology 1992, 37 (1), 75–93.

446 447 448

(30)

Wang, P.; Wang, H.-P.; Xu, M.-Y.; Liang, Y.-J.; Sun, Y.-J.; Yang, L.; Li, L.; Li, W.; Wu, Y.-J. Combined subchronic toxicity of dichlorvos with malathion or pirimicarb in mice liver and serum: A metabonomic study. Food Chem. Toxicol. 2014, 70, 222–230.

449 450 451 452

(31)

Moser, V. C.; Stewart, N.; Freeborn, D. L.; Crooks, J.; MacMillan, D. K.; Hedge, J. M.; Wood, C. E.; McMahen, R. L.; Strynar, M. J.; Herr, D. W. Assessment of serum biomarkers in rats after exposure to pesticides of different chemical classes. Toxicol. Appl. Pharmacol. 2015, 282 (2), 161–174.

453 454 455

(32)

Dondero, F.; Negri, A.; Boatti, L.; Marsano, F.; Mignone, F.; Viarengo, A. Transcriptomic and proteomic effects of a neonicotinoid insecticide mixture in the marine mussel (Mytilus galloprovincialis, Lam.). Sci. Total Environ. 2010, 408 (18), 3775–3786.

456 457 458

(33)

Spurgeon, D. J.; Jones, O. A. H.; Dorne, J.-L. C. M.; Svendsen, C.; Swain, S.; Stürzenbaum, S. R. Systems toxicology approaches for understanding the joint effects of environmental chemical mixtures. Sci. Total Environ. 2010, 408 (18), 3725–3734.

16 ACS Paragon Plus Environment

Environmental Science & Technology

459 460 461

(34)

Gomez-Eyles, J. L.; Svendsen, C.; Lister, L.; Martin, H.; Hodson, M. E.; Spurgeon, D. J. Measuring and modelling mixture toxicity of imidacloprid and thiacloprid on Caenorhabditis elegans and Eisenia fetida. Ecotoxicol. Environ. Saf. 2009, 72 (1), 71–79.

462 463

(35)

Fox, I. H. Metabolic basis for disorders of purine nucleotide degradation. Metabolism. 1981, 30 (6), 616–634.

464 465

(36)

Duzguner, V.; Erdogan, S. Chronic exposure to imidacloprid induces inflammation and oxidative stress in the liver & central nervous system of rats. Pestic. Biochem. Physiol. 2012, 104 (1), 58–64.

466 467 468

(37)

Paschen, W. Polyamine metabolism in different pathological states of the brain. Mol. Chem. Neuropathol. Spons. Int. Soc. Neurochem. World Fed. Neurol. Res. Groups Neurochem. Cerebrospinal Fluid 1992, 16 (3), 241–271.

469

(38)

Seiler, N. Oxidation of Polyamines and Brain Injury. Neurochem. Res. 2000, 25 (4), 471–490.

470 471

(39)

Shohami, E.; Nates, J. L.; Glantz, L.; Trembovler, V.; Shapira, Y.; Bachrach, U. Changes in brain polyamine levels following head injury. Exp. Neurol. 1992, 117 (2), 189–195.

472 473 474

(40)

Adibhatla, R. M.; Hatcher, J. F.; Sailor, K.; Dempsey, R. J. Polyamines and central nervous system injury: spermine and spermidine decrease following transient focal cerebral ischemia in spontaneously hypertensive rats. Brain Res. 2002, 938 (1–2), 81–86.

475 476 477

(41)

Zielińska, E.; Kocki, T.; Saran, T.; Borbely, S.; Kuc, D.; Vilagi, I.; Urbańska, E. M.; Turski, W. A. Effect of pesticides on kynurenic acid production in rat brain slices. Ann. Agric. Environ. Med. AAEM 2005, 12 (2), 177–179.

478 479 480 481

(42)

Sun, J.; Schnackenberg, L. K.; Holland, R. D.; Schmitt, T. C.; Cantor, G. H.; Dragan, Y. P.; Beger, R. D. Metabonomics evaluation of urine from rats given acute and chronic doses of acetaminophen using NMR and UPLC/MS. J. Chromatogr. B 2008, 871 (2), 328–340.

17 ACS Paragon Plus Environment

Page 18 of 21

Page 19 of 21

Environmental Science & Technology

Figure 1. Results of the AChE bioassay. A) The AChE activity of the neutral extract with an IC50 of 33.39 ECF (n=3, p-values < 0.05, Tukey’s test). B) The AChE activity of the anionic and cationic extracts, not showing any AChE inhibition (n=3, p-values < 0.05, Tukey’s test). 105x39mm (300 x 300 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 2. Biochemical network based on fold changes of metabolites, comparing the CNS of L. stagnalis exposed to the neutral extract (16 ECF) and the neutral blank extract. Red and green dots indicate, respectively, a reduced or increased fold change, whereas grey dots indicate not detected metabolites. Big dots symbolize significant metabolites (n=7, p-values < 0.05, t-test). Arrows describe metabolic interactions and dashed lines chemical similarities. 203x196mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 20 of 21

Page 21 of 21

Environmental Science & Technology

Figure 3. Polyamine metabolism pathway, highlighting the catabolic and synthetic reaction and metabolites and enzymes involved (A). Metabolites of the polyamine pathway after exposure to the neutral cartridge extract (NEUTRAL) and imidacloprid (IMI) and thiacloprid (THIA) (B). The metabolites permutation is shown with a green arrow indicating an up-regulation and with a red arrow a down-regulation. The amplitude of these changes is given by the corresponding fold changes. 185x94mm (300 x 300 DPI)

ACS Paragon Plus Environment