Pesticide Body Burden of the Crustacean Gammarus pulex as a

Jun 22, 2018 - We calculated SPEAR using the desktop application Indicate 1.0.0. ..... (3,7,47) The toxicity derived from body burden (TUmax-Int) expl...
0 downloads 0 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

Environmental Measurements Methods

Pesticide body burden of the crustacean Gammarus pulex as a measure of toxic pressure in agricultural streams Naeem Shahid, Jeremias Martin Becker, Martin Krauss, Werner Brack, and Matthias Liess Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 22 Jun 2018 Downloaded from http://pubs.acs.org on June 22, 2018

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 29

Environmental Science & Technology

1

Title

2

Pesticide body burden of the crustacean Gammarus pulex as a measure of toxic pressure in

3

agricultural streams

4

Author Affiliation 1. Author: Naeem Shahid *,†,‡,§ ([email protected], phone: 0341 235 1494, fax: +49 341

5

235 1494)

6 7

2. Author: Jeremias Martin Becker†,‡

8

3. Author: Martin Krauss∥

9

4. Author: Werner Brack‡,∥ 5. Author: Matthias Liess*,†,‡ ([email protected], phone: +49 341 235 1578, fax: +49

10

341 235 1578)

11 12



13

Permoserstraße 15, 04318 Leipzig, Germany

14



15

Germany

16

§

17

Vehari, Pakistan

18

∥UFZ,

19

Permoserstraße 15, 04318 Leipzig, Germany

20

*

UFZ, Helmholtz Centre for Environmental Research, Department of System-Ecotoxicology,

RWTH Aachen University, Institute for Environmental Research (Biology V), Aachen,

Department of Environmental Sciences, COMSATS Institute of Information Technology,

Helmholtz Centre for Environmental Research, Department Effect-Directed Analysis,

Corresponding author

21 22 23 1 ACS Paragon Plus Environment

Environmental Science & Technology

24

ABSTRACT

25

Risk assessments of toxicants in aquatic environments are typically based on the evaluation of

26

concentrations in water or sediment. However, concentrations in water are highly variable, while

27

the body burden may provide a better time-integrated measure of pesticide exposure and

28

potential effects in aquatic organisms. Here, we quantified pesticide body burdens in a dominant

29

invertebrate species from agricultural streams, Gammarus pulex, compared them pesticide

30

concentrations in water samples, and linked the pesticide contamination with observed ecological

31

effects on macroinvertebrate communities.

32

In total, 19 of 61 targeted analytes were found in the organisms, ranging from 0.037 to 93.94 ng

33

g-1 (wet weight). Neonicotinoids caused the highest toxic pressure among the pesticides detected

34

in G. pulex. Using linear solvation energy relationships (LSERs), we derived equivalent pesticide

35

concentrations in stream water based on the body burden. These equivalent concentrations

36

correlated with the concentrations in water samples collected after run-off (65% of variance

37

explained).

38

Pesticide pressure significantly affected the aquatic macroinvertebrate community structure,

39

expressed as SPEARpesticides, and caused, on average, threefold increased insecticide tolerance in

40

G. pulex as a result of adaptation. The toxic pressure derived from body burden and from water

41

samples similarly explained the change in community structure (68% and 64%). However, the

42

increased tolerance of G. pulex to pesticides was better explained by the toxicity derived from

43

body burden (70%) than by the toxicity from water samples (53%). We conclude that the internal

44

body burden of macroinvertebrates is suitable to assess the overall pesticide exposure and effects

45

in agricultural streams.

2 ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29

Environmental Science & Technology

46

Abstract Art

47 48

INTRODUCTION

49

Pesticides in agricultural streams considerably affect the diversity and species composition of

50

aquatic macroinvertebrates1, 2 and, consequently, related ecosystem functions such as leaf litter

51

degradation.3-6 Sensitive species show significant reduction even at concentrations more than 2

52

orders of magnitude below the acute median lethal concentration (LC50) that kills 50% of the

53

individuals of laboratory standard test species.2,

54

concentrations include an increased sensitivity of individuals to pesticides under multiple stress

55

conditions8 and the culmination of effects from sequential exposure.9 The same low

56

concentrations exert a considerable pressure for adaptation, which results in increased pesticide

57

tolerance in exposed non-target species, with unknown ecological consequences on their general

58

fitness.10,

59

regarding the protection from pesticide contamination are necessary. The decision base for such

60

measures is a profound quantification and localisation of pesticide exposure and effects through

61

monitoring in agricultural streams.

11

7

Reasons for these low environmental effect

To reduce effects of pesticides in the environment, local management measures

3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 29

62

One relevant pathway of pesticides to enter agricultural streams is through run-off from arable

63

land after heavy rainfall. Accordingly aquatic organisms are typically exposed to short pulses of

64

high pesticide concentrations.12 Thus, relevant measurements of pesticide exposure must include

65

run-off events; as such peak concentrations are related to the long-term effects on the

66

macroinvertebrate community much better than the low concentrations found in event-

67

independent samples.2, 13 However, measuring these short-term pesticide pulses is challenging as

68

the peak concentrations depend on the timing and intensity of both the pesticide application and

69

the subsequent run-off event, which can vary widely within a spraying season.12, 13

70

Currently, a range of passive samplers have been suggested for integrative sampling of various

71

pollutants in surface waters14-16 and have also been used to capture short-term pollution events.17,

72

18

73

linking local toxic pressure to ecological effects on macroinvertebrate communities.3,

74

However, the installation of passive samplers is costly and time-consuming. Toxic pressure was

75

defined here as a stress exerted by pesticides in streams, using Toxic Unit20 modified by applying

76

the highest Toxic Unit - TUmax as a measure for the ecological relevant toxic pressure.2 In the

77

present study, we additionally used an alternative quantification of pesticide exposure based on

78

pesticide residues in macroinvertebrates sampled from agricultural streams. Organisms

79

accumulate pesticides in lipids and proteins from the water phase according to the compound

80

properties and therefore, may serve as natural passive samplers.21 Their body burden may

81

provide a better time-integrated picture of the overall pesticide exposure than do short-term

82

measurements from the water phase. Additionally, an assessment of pesticide exposure based on

83

internal concentrations may be more representative because it is the concentration at the site of

Generally, pesticide concentrations determined by passive samplers are highly suitable for

4 ACS Paragon Plus Environment

15, 19

Page 5 of 29

Environmental Science & Technology

84

action within organisms, rather than the external concentration in water, that triggers adverse

85

effects.22

86

We used pesticide residues in the common freshwater shrimp G. pulex to predict the time

87

integrated pesticide concentrations in stream water based on equilibrium partitioning23 with a

88

linear solvation energy relationship (LSER)24 model. This model estimates the partitioning of

89

pesticides into different compartments such as storage and membrane lipids and proteins and

90

considers polar and H-bonding interactions rather than hydrophobicity only. The predicted

91

equilibrium concentrations in water (Cw,eq) were compared with the pesticide peak

92

concentrations in water samples from the same streams following rain events. Additionally, we

93

evaluated the ecological effect of the toxic pressure estimated from body burden and from water

94

samples on the macroinvertebrate community composition according to the SPEAR approach

95

(Species at Risk)2 and on the pesticide tolerance in G. pulex as a consequence of adaptation.

96

MATERIALS AND METHODS

97

Description of the study sites

98

The investigated streams were located in central Germany, including 9 sites from intense

99

agricultural landscapes and 6 sites from either undisturbed forested sections or from streams with

100

considerable buffer strips and non-agricultural upstream reaches (Supporting Information, Figure

101

S1). The second group was considered non-agricultural streams and used as reference sites.

102

During macroinvertebrate sampling, different environmental parameters such as water level, pH,

103

water temperature, dissolved oxygen and electrical conductivity were measured (Supporting

104

Information, Table S1).

105

Sampling of Gammarus pulex and water

5 ACS Paragon Plus Environment

Environmental Science & Technology

106

Thirty individuals of G. pulex (wet weight 900mg) per site were sampled for the analysis of

107

pesticide body burden in May 2016 during the maximum pesticide application period in the

108

study area.12, 25 Organisms of different sizes were collected to avoid any bias regarding greater

109

bioaccumulation and exposure time. After collection, individuals were transferred to the lab and

110

stored at −20°C.

111

Water samples from all streams were collected using event-driven samplers (EDS)2 from May to

112

July 2016 during the maximum pesticide application period. The number of samples from each

113

site varied from 1 to 5, depending on rainfall events. Briefly, two brown glass bottles, each with

114

a capacity of 1 L, were installed at each sampling site. Both bottles were attached to a stainless

115

rod at heights of 5 and 15 cm from the regular level of the stream water that filled when the

116

water level had risen due to heavy rainfall (see Supporting Information, Figure S2 for illustration

117

of event-driven samplers). After rainfall events, bottles were collected within 24 h, kept in a cool

118

box at 4°C and transported to the laboratory. Subsequently, 1 mL aliquots were transferred into 2

119

mL autosampler vials and stored at −20°C for analysis. Results of pesticide concentrations in

120

water samples have already been published elsewhere.11 In the present study we compared water

121

concentrations (external concentrations) with body burden (internal concentrations) to measure

122

the pesticide exposure and potential effects in aquatic organisms.

123

Sample preparation and extraction

124

Target compounds were extracted from gammarids using the method of Inostroza, et al.26

125

(pulverized liquid extraction and a QuEChERS approach with an additional hexane phase).

126

Frozen gammarids were transferred to a 10 mL centrifuge tube; 1 mL of LC-MS grade

127

acetonitrile, 1 mL of LC-MS grade water and 0.5 mL of LC grade hexane were added; then the

128

sample was homogenized using an Ultra-Turrax (IKA, Staufen, Germany) homogenizer.

6 ACS Paragon Plus Environment

Page 6 of 29

Page 7 of 29

Environmental Science & Technology

129

Subsequently, 400 mg of MgSO4 and 100 mg of NaCl were added to induce phase separation

130

between water and acetonitrile. The top-layer hexane fraction was taken for the analysis of

131

pyrethroids via gas chromatography-tandem mass spectrometry (GC-MS/MS), evaporated to

132

dryness in a nitrogen stream, reconstituted in 500 µL of ethyl acetate and filtered using a PTFE

133

syringe filter (pore size 0.45 µm) into a 2 mL autosampler vial. Permethrin-D6 was added as an

134

injection standard. The acetonitrile phase was cleaned up via dispersive solid-phase extraction

135

using primary-secondary amine and transferred into a 2 mL glass vial. The sample was

136

evaporated to dryness and reconstituted in 450 µL of methanol:water (70:30) for liquid

137

chromatography-high-resolution mass spectrometry (LC-HRMS) analysis. A volume of 50 µL of

138

internal standard solution containing 40 isotope-labelled compounds (1 µg/mL) was added.

139 140

GC-MS/MS analysis

141

The pyrethroid analysis was conducted via GC-MS/MS using an Agilent GC 7890-MS/MS 7000

142

system with negative-mode chemical ionization. For separation, we used a HP-5ms UltraInert

143

column (30 m x 0.25 mm, 0.25 µm film thickness) equipped with a guard column (5 m x 0.25

144

mm). The column oven program started at 100°C, was held for 1 min, increased at 30°C/min to

145

150°C, then at 6°C/min to 186°C and at 10°C/min to 300°C, and was then held for 2 min. The

146

transfer line temperature was 300°C. Splitless injection of 1 µL was used at an injector

147

temperature of 280°C. Pyrethroids were quantified with MassHunter software (Agilent) using

148

one quantifier and up to two qualifier MS/MS transitions. A method-matched calibration was

149

prepared using water spiked at seven concentration levels ranging from 0.028 to 28 ng (g

150

gammarids wet weight)-1, for which the same sample preparation protocol was used as for the

151

samples.

7 ACS Paragon Plus Environment

Environmental Science & Technology

152 153

LC-HRMS analysis

154

Polar pesticides were analysed using LC-HRMS in positive ion mode (Ultimate 3000 LC system

155

coupled to a QExactive Plus MS equipped with a heated electrospray ionization (ESI) source, all

156

from Thermo Scientific). Separation was conducted using a methanol:water gradient (eluent A

157

and B, respectively, both containing 0.1% formic acid) on a Kinetex C18 EVO column (50 x 2.1

158

mm, 2.6 µm particle size, Phenomenex) at a flow rate of 300 µL/min. After eluting for 1 min

159

with 5% B, the fraction of B was linearly increased to 100% within 12 min and kept for 11 min.

160

Then, the flow was diverted to waste, and the column was rinsed for 2 min using

161

isopropanol:acetone 50:50 / eluent B / eluent A (85% / 10% / 5%). Finally, the column was re-

162

equilibrated to the starting conditions for 5.7 min. The heated ESI source and the transfer

163

capillary were operated at 300°C, the spray voltage was 3.8 kV, the sheath gas flow rate was 45

164

a.u., and the auxiliary gas flow rate was 1 a.u. The runs combined a full scan experiment (100-

165

1000 m/z) at a nominal resolving power of 70,000 (referenced to m/z 200) and data-independent

166

MS/MS experiments at a nominal resolving power of 35,000. For the latter, we acquired the data

167

using isolation windows of approximately 50 mass units (i.e., m/z ranges 97-147, 144-194, 191-

168

241, 238-288, 285-335, 332-382, 379-429, 426-476) and 280 mass units (i.e., m/z ranges 460-

169

740, 730-1010). Eight method-matched calibration standards corresponding to levels of 0.056 to

170

27.8 ng (g gammarids wet weight)-1 were prepared the same way as the samples using 1 mL of

171

spiked water instead of gammarids. For water analysis, the calibration standards at

172

concentrations from 1 to 2000 ng/L were prepared using the same method as that applied to the

173

samples using 1 mL of pristine stream water, 25 µL of a methanolic analyte stock solution of the

8 ACS Paragon Plus Environment

Page 8 of 29

Page 9 of 29

Environmental Science & Technology

174

appropriate concentration, 25 µL of an internal standard solution and 10 µL of a 2 M NH4-

175

formate buffer (pH 3.5).

176

For data evaluation, the software Trace Finder 3.2 (Thermo) was used, employing internal

177

calibration with the isotope-labelled compound with the closest retention time of each analyte.

178

For quantification, we used the full scan extracted ion chromatogram (7 ppm window), while for

179

confirmation, one or two diagnostic MS/MS fragments and isotope patterns were compared with

180

those of reference standards. Table S2 provides an overview of all analysed pesticides, their

181

physicochemical properties, method detection limits and recoveries for the gammarid extraction

182

method and Table S3 contains detected pesticides at all study sites.

183 184

Calculation of chemical activity and partition coefficients

185

We assumed proteins and lipids as the relevant phases for the absorption of organic pollutants

186

and calculated estimated equilibrium concentrations (Cw,eq) in body tissues of G. pulex (Eq. 1).

‫ܥ‬௪,௘௤ =

‫ܥ‬௕,௠ (݂ௌ௅ ‫ܭ‬ௌ௅ௐ + ݂௠௟ ‫ܭ‬ெ௅ௐ + ݂௉ ‫ܭ‬௉ௐ ) × 0.2175

(1)

187

where Cb,m is the total measured concentration of contaminants in body tissues, fSL is the storage

188

lipid fraction, KSLW represents the storage lipid-water partitioning coefficient, fML is the

189

membrane lipid fraction, KMLW represents the membrane lipid-water partition coefficient, fP is the

190

protein fraction, and KPW is the protein-water partitioning coefficient. To calculate gammarid-

191

water partition coefficients, the composition of G. pulex (per dry weight) was assumed as 47%

192

protein,27 4.5% storage lipid and 1.5% membrane lipid.28 The factor 0.217529 was used for the

193

conversion of wet weight to dry weight. This assumption is well in agreement with typical values

194

from more in-depth studies on the composition of Gammarus pulex,30 although ignoring

9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 29

195

variability between sexes, seasons and individuals. The present experimental setting did not

196

allow for measuring individual lipid contents in the specimens used in the experiments.

197

Partition coefficients were calculated using linear solvation energy relationships (LSERs, Eq. 2)

198

for (1) protein (muscle)–water, (2) storage lipid–water and (3) membrane lipid–water. The

199

system parameters (v, e, s, a, b, and c) used for the respective systems were taken from the open-

200

access UFZ-LSER database.24 The compound descriptors (V, E, S, A and B) were calculated for

201

all compounds using the software ACD/Percepta (ACD labs, v2016; Supporting Information,

202

Table S4).

203

log ‫ܭ‬௜ = ‫ ܸݒ‬+ ݁‫ ܧ‬+ ‫ ܵݏ‬+ ܽ‫ ܣ‬+ ܾ‫ ܤ‬+ ܿ

(2)

204

Calculation of pesticide exposure

205

Measured and estimated equilibrium water concentrations of pesticides were converted into toxic

206

units (TU) representing the decadic logarithm of the concentration of a compound divided by its

207

LC50 in a 48 h acute toxicity test with a standard reference organism. We estimated the overall

208

pesticide-induced toxic pressure in a stream as the maximum toxic unit of all compounds (TUmax;

209

Eq. (3)2), quantified either from internal body burdens (TUmax-Int) or from external water samples

210

(TUmax-Ext). TUmax works very well for agricultural streams but may not for other ecosystems that

211

receive contamination from many different sources.31

212

஼௜

௡ TU୫ୟ୶ = Max௜ୀଵ ቂlogଵ଴ ቀ௅஼ ቁቃ ఱబ೔

(3)

213

where TUmax is the highest toxic unit of n pesticides at a given site, Ci is the freely dissolved

214

pesticide concentration (µg L−1), and LC50i is the median lethal concentration (µg L−1) of that

215

pesticide for a reference organism after a 48 h exposure (Supporting Information, Table S5). We

216

used Daphnia magna, Chironomus tentans, Chironomus riparius and Hyalella azteca as

217

reference organisms; the LC50 values were obtained from the ECOTOX database32 and the 10 ACS Paragon Plus Environment

Page 11 of 29

Environmental Science & Technology

218

Pesticide Properties Database PPDB.33 When data from several organisms were available for the

219

same compound, the most sensitive organism was used. For instance, Chironomus riparius and

220

Hyalella azteca showed much higher sensitivity to neonicotinoid insecticides as compared to

221

Daphnia magna and were therefore considered as reference organisms. In few cases, the name of

222

the reference species was not provided in the database and therefore considered unknown

223

(Supporting Information, Table S5).

224 225

Characterizing effects of pesticides on the macroinvertebrate community

226

At each sampling site, we quantified the long-term effects of pesticides on the macroinvertebrate

227

community composition using the bio-indicator SPEARpesticides2 The SPEAR index categorizes

228

aquatic macroinvertebrates into vulnerable and non-vulnerable species on the basis of four

229

ecological traits: (1) potential for pesticide exposure, (2) physiological sensitivity, (3) generation

230

time (pace of recovery), and (4) migration ability (recolonization). We calculated SPEAR using

231

the desktop application SPEAR Calculator 0.10.0.2, 34

232

Stream macroinvertebrates were sampled from May to July 2016 along with the collection of G.

233

pulex. At each site, ten subsamples were collected from different habitats along a 50 m stream

234

section. Macroinvertebrates were collected using a 25 × 25 cm kick-net with a 500 µm mesh

235

size. The content of the Surber sampler was filtered through a set of four sieves (mesh sizes: 8

236

mm, 4 mm, 2 mm, and 500 µm; Retsch GmbH, Haan, Germany) and transferred into plastic trays

237

for the identification of taxa (family level) and determination of abundance. When the family of

238

an organism could not be identified, the individuals were preserved in ethanol and identified in

239

the laboratory with a stereo-microscope (Zeiss Discovery V20; Oberkochen, Germany). The

11 ACS Paragon Plus Environment

Environmental Science & Technology

240

abundance of each family was divided by the sampling area to calculate population densities per

241

m² used for the calculation of the SPEARpesticides values.

242 243

Identification of insecticide tolerance in Gammarus pulex

244

At each study site, 170 individuals of G. pulex within a size range of 6–10 mm were collected

245

together with the individuals used for the body burden analysis, using a 25 × 25 cm kick-net with

246

a 500 µm mesh size. The individuals were caught with a pipette and transferred into aerated and

247

cooled plastic boxes filled with stream water. On the day of collection, the organisms were

248

transported to the laboratory and left to acclimate in white trays filled with a mixture of aerated

249

ADaM (Artificial Daphnia medium)35 and steam water at 15°C for 24 h.

250

Insecticide sensitivity was assessed in acute toxicity tests following the OECD guidelines for the

251

testing of chemicals and the guidelines for rapid tests for community-level risk assessment.36, 37

252

We selected the neonicotinoid insecticide clothianidin for sensitivity testing that represents the

253

most commonly applied class of insecticides in agriculture for more than a decade.38 A 500 mg L-

254

1

255

(Spiess-Urania Chemical GmbH, Germany) in distilled water while stirring for 12 h and further

256

diluted in ADaM to obtain the required test concentrations. The test organisms were exposed to

257

seven different concentrations i.e., 0, 5, 15, 45, 135, 405 and 1215 µg/L for a period of 48 h.

258

Only undamaged and active individuals without visible infection with parasites were selected for

259

pesticide exposure. Five organisms were placed in a stainless tea strainer, and four tea strainers

260

per concentration were placed in a glass beaker containing 1000 mL of the test solution. During

261

contamination, the beakers were constantly aerated and the temperature was maintained at 15°C.

262

After exposure for 48 h, the organisms were checked for immobility. Organisms were considered

stock solution of the neonicotinoid insecticide clothianidin was prepared with DANTOP®

12 ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29

Environmental Science & Technology

263

immobile when they neither moved their bodies within 30 seconds of undisturbed observation

264

nor after probing with a rod (fanning of gills and antenna was not considered body movement).

265 266

Data analysis

267

All data analyses were conducted using the open-source application RStudio version 1.1.383 for

268

Windows39 and the basic R version 3.4.3 for Windows.40 To compare the insecticide tolerance in

269

different populations, we calculated the EC50 (concentration that immobilized 50% of the test

270

individuals) and the 95% confidence intervals from the acute toxicity tests using a generalized

271

linear model with a quasi-binomial error distribution and a logit link function. The tolerance of

272

G. pulex populations from each site was quantified as the ratio of the local EC50/mean EC50 of all

273

populations from non-contaminated streams (TUmax-Int < -3). This EC50 ratio allowed us to

274

directly compare the tolerance of different populations to neonicotinoids and thus to quantify

275

adaptation in contaminated streams.

276

Linear regression was applied to analyse the relationships between the toxic pressure calculated

277

from pesticide body burden and water concentration (TUmax-Int and TUmax-Ext, respectively) and

278

the change in community structure or the neonicotinoid tolerance. To identify the additional

279

effect of the local species diversity on the increase in neonicotinoid tolerance, we applied

280

multiple linear regression model. The normal distribution and homoscedasticity of residuals were

281

confirmed using normal-Q-Q plots, histograms, plots of residuals vs. fitted values, and statistical

282

tests for deviations from the normal distribution and from variance homogeneity. To obtain a

283

normal distribution, the EC50 ratio was ln(x) transformed before analysis. Agricultural and non-

284

agricultural streams were compared with respect to the SPEAR index, EC50 and physicochemical

13 ACS Paragon Plus Environment

Environmental Science & Technology

285

parameters using two-sample t-tests and Welch's t-tests. In cases in which the data were not

286

normally distributed, Wilcoxon's rank sum test was used.

287 288

RESULTS

289

Pesticide residues in individuals of Gammarus pulex

290

A total of 19 chemicals of 61 targeted compounds (see Supporting Information, Table S2) were

291

quantified in all organisms collected (Supporting Information, Table S3). Concentrations ranged

292

from 0.037 to 14.28 ng g-1 (wet weight) for insecticides, 0.08 to 93.94 ng g-1 for herbicides and

293

0.06 to 12.74 ng g-1 for fungicides. The fungicide prothioconazole-desthio was the most

294

frequently detected chemical across all sites. The highest concentration was found for the

295

herbicide diflufenican with 93.94 ng g-1 (Supporting Information, Table S3).

296

Generally, insecticides were detected only in organisms from agricultural streams. The only

297

exceptions were low concentrations of the pyrethroid insecticides cyfluthrin and cypermethrin in

298

organisms from a single non-agricultural stream (Non-Agri 6). Additionally, organisms from

299

non-agricultural streams showed only very low concentrations of the fungicide prothioconazole-

300

desthio (0.17-0.28 ng g-1) and the herbicides diflufenican (2.58 ng g-1) and prosulfocarb (2.18 ng

301

g-1).

302 303

Pesticides in the water samples

304

A total of 50 chemicals of 55 targeted analytes were detected in the water samples collected from

305

the streams. All the agricultural sites were contaminated with most of the analytes, including

306

insecticides, fungicides and herbicides. The fungicide prothioconazole-desthio and herbicides

307

ethofumesate and terbuthylazine were the most frequently detected chemicals in water samples

14 ACS Paragon Plus Environment

Page 14 of 29

Page 15 of 29

Environmental Science & Technology

308

across all the streams. The highest concentration was found for the herbicide ethofumesate with

309

84.70 µg L-1. In agricultural streams, the most commonly detected insecticides were clothianidin,

310

thiamethoxam and thiacloprid. The highest concentration was determined for thiamethoxam

311

(2.83 µg L-1). In most non-agricultural streams, only certain herbicides and fungicides were

312

detected. Four of these sites were slightly contaminated (3.4 to 12 ng L-1) with the neonicotinoid

313

insecticide clothianidin (Supporting Information, Table S6).

314 315

Estimation of aquatic pesticide exposure

316

The toxicity of the pesticide mixture measured in water samples after runoff events was

317

expressed in toxic units of the most toxic compound (TUmax-Ext, see materials and methods). To

318

derive the pesticide exposure alternatively from the body burden of G. pulex, we converted the

319

internal pesticide concentrations to estimated equilibrium concentrations in water and calculated

320

the maximum toxic unit (TUmax-Int) for invertebrates. The results of both methods are compared

321

in Table 1.

322

In agricultural streams, the TUmax-Int ranged from -1.3 to 0.7 with a median value of 0.4. In all of

323

the agricultural sites, neonicotinoids represented the most toxic compound. The highest toxicity

324

was caused by imidacloprid (TUmax-Int at 5 sampling sites; mean TU = 0.6), followed by

325

thiacloprid (TUmax-Int at 4 sampling sites; mean TU = -0.8). In contrast, the non-agricultural

326

streams were contaminated only to a minor extent. In those streams, the TUmax-Int ranged from -

327

6.2 to -4.3 with a median value of -6.0.

328

The TUmax-Ext in agricultural streams ranged from -3.2 to -1.1 and was always lower than the

329

TUmax-Int, with a mean difference of 1.8 TU. In non-agricultural streams, by contrast, the TUmax-

330

Ext

ranged from -3.3 to -5.0 and was always higher than the TUmax-Int (mean difference: 1.5 TU),

15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 29

331

except for one stream in which insecticides were detected in the organisms (Non-Agri 6).

332

However, TUmax-Int and TUmax-Ext were significantly correlated with one another (adjusted R² =

333

0.65, p < 0.001; Fig. 1). Both TUmax-Int and TUmax-Ext showed a clear separation of streams from

334

intense agricultural areas and from reference streams. However, within these groups the TUmax-Int

335

was not clearly associated with the TUmax-Ext, indicating a relevant amount of residual variation.

336

Table 1

337

Pesticide-induced toxic pressure in the investigated streams, expressed as the toxic unit of the most toxic compound

338

(TUmax). The toxic units were calculated from pesticide concentrations in water samples collected after runoff events

339

(TUmax-Ext); they were also derived from internal pesticide concentrations detected in G. pulex, which were converted

340

to equivalent concentrations in water (TUmax-Int). Additionally, the toxic pressure was estimated using the

341

SPEARpesticides bio-indicator based on the macroinvertebrate community composition (SPEAR), and using the

342

increase in tolerance of G. pulex (EC50) to the neonicotinoid insecticide clothianidin as a result of adaptation. Most toxic compound

TUmax

SPEAR

Site ID

EC50 (95% CI)

Ext

Int

In water

In organisms

Agri-1

-1.8

0.4

Thiamethoxam (I)

Imidacloprid (I)

22.3

263(228-302)

Agri-2

-3.2

0.0

Diflufenican (H)

Thiacloprid (I)

24.4

298(227-388)

Agri-3

-1.9

-0.6

Pirimiphos-methyl (I) Thiacloprid (I)

10.2

321(203-506)

Agri-4

-1.8

0.4

Thiamethoxam (I)

Imidacloprid (I)

17.0

178(148-213)

Agri-5

-1.9

-1.3

Clothianidin (I)

Thiacloprid (I)

32.4

191(135-270)

Agri-6

-1.1

0.5

Thiamethoxam (I)

Imidacloprid (I)

7.0

213(171-265)

Agri-7

-1.7

0.4

Thiamethoxam (I)

Imidacloprid (I)

14.7

153(121-192)

Agri-8

-1.7

0.7

Thiamethoxam (I)

Imidacloprid (I)

25.0

190(144-249)

Agri-9

-1.9

-1.3

Clothianidin (I)

Thiacloprid (I)

15.6

160(121-210)

Non-Agri-1

-3.6

-6.2

Pirimicarb (I)

Prothioconazole-desthio (F)

49.7

83(62-110)

Non-Agri-2

-4.8

-6.1

Prosulfocarb (H)

Prothioconazole-desthio (F)

37.7

94(66-132)

Non-Agri-3

-3.3

-6.0

Clothianidin (I)

Prothioconazole-desthio (F)

48.4

67(37-120)

16 ACS Paragon Plus Environment

Page 17 of 29

Environmental Science & Technology

Non-Agri-4

-3.4

-6.1

Clothianidin (I)

Prothioconazole-desthio (F)

34.4

98(85-111)

Non-Agri-5

-3.8

-4.3

Clothianidin (I)

Diflufenican (H)

37.2

93(65-130)

50.3

57(38-83)

Non-Agri-6

-5.0

-4.4

-

Cyfluthrin

(I)

343

I: insecticide; F: fungicide; H: herbicide; Agri: Agricultural sites; Non-Agri: Non-agricultural sites; Ext: external;

344

Int: Internal.

345

346 347

Fig. 1. Relationship between TUmax-Int (TUmax calculated using

348

equivalent pesticide concentrations in water derived from internal

349

pesticide concentrations of G. pulex) and TUmax-Ext (TUmax calculated

350

using pesticide concentrations measured in water). The grey area

351

corresponds to the 95% confidence interval. R² = 0.67, adjusted R² =

352

0.65, F = 26.5, residual d.f. = 13, p < 0.001.

353 354

Ecological effect assessment

17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 29

355

To compare how both estimations of pesticide exposure (TUmax-Ext and TUmax-Int) could explain

356

the observed ecological effects of pesticides on macroinvertebrate communities, we applied the

357

SPEARpesticides

358

(http://www.systemecology.eu/indicate/). Agricultural streams showed low to medium

359

SPEARpesticides values from 7 to 32.4, which corresponded to a range from bad (< 11) to moderate

360

(< 33) ecological status as defined by Beketov et al.41 In contrast, non-agricultural streams

361

showed significantly higher SPEARpesticides values (Wilcoxon's rank sum test, W = 54, p-value
44) ecological status (Table 1). Of all environmental parameters measured,

364

SPEARpesticides depended only on pesticide exposure and was equally well explained by the

365

toxicity calculated from body burdens (TUmax-Int, Fig. 2a) and by the toxicity in water samples

366

(TUmax-Ext, Fig. 2b).

indicator.2

Calculation

was

performed

using

the

INDICATE

tool

367

368 369

Fig. 2. Relationship between the community composition of macroinvertebrates, expressed as SPEARpesticides, and

370

water contamination, expressed as toxic units. The grey areas correspond to the 95% confidence interval. (a) The

371

calculation of TUmax-Int was based on the equivalent pesticide concentrations in water derived from internal

18 ACS Paragon Plus Environment

Page 19 of 29

Environmental Science & Technology

372

pesticide concentrations in G. pulex; R² = 0.70, adjusted R² = 0.68, F = 30.89, regression d.f. = 1, residual d.f. =

373

13, p < 0.0001. (b) The calculation of TUmax-Ext was based on pesticide concentrations in water samples collected

374

from streams after run-off events; R² = 0.67, adjusted R² = 0.64, F = 26.28, regression d.f. = 1, residual d.f. = 13, p

375

< 0.0001.

376 377

In the next step, we compared the calculated TUmax values and SPEARpesticides with a second

378

indicator of pesticide effects, the insecticide tolerance of G. pulex from agricultural and non-

379

agricultural sites. On average, populations from agricultural streams showed 3-fold higher

380

tolerance to the neonicotinoid insecticide clothianidin than those from non-agricultural control

381

streams (Welch two-sample t-test, p < 0.001, t = -6.39, d.f. = 9.67; Fig. 3).

382 383

Fig. 3. Comparison of the median effective concentrations (EC50) of G.

384

pulex collected from agricultural (red) and non-agricultural streams

385

(green) after 48 h exposure to the neonicotinoid insecticide clothianidin in

386

the laboratory. The boundaries of the central box are the 25th and 75th

19 ACS Paragon Plus Environment

Environmental Science & Technology

387

percentiles; the horizontal line is the median; the blue dot indicates the

388

mean; and whiskers of the boxplot represent the minimum and maximum

389

values. The significance level for the observations is displayed as: *** =

390

p < 0.001.

391

The observed increase in insecticide tolerance was better explained by the pesticide toxicity

392

estimated from internal body burden (TUmax-Int, adjusted R2 = 0.70, Fig. 4a) than that from

393

pesticide concentrations in water samples (TUmax-Ext, adjusted R2 = 0.53, Fig. 4b). The increase in

394

insecticide tolerance with TUmax-Int tended to be stronger when the local species diversity was

395

below average, but the difference was not significant (adjusted R² = 0.72, F = 3.28, d.f. = 1,

396

residual d.f. = 11, p = 0.097 for the interaction toxic unit:Shannon; Supporting Information,

397

Figure S4). The two indicators of pesticide effects, SPEARpesticides and the insecticide tolerance in

398

G. pulex, were closely correlated (adjusted R² = 0.67; Fig. 5).

399 400

Fig. 4. The freshwater shrimp G. pulex adapted to pesticide exposure, resulting in increased tolerance to the

401

neonicotinoid insecticide clothianidin. Pesticide exposure (in toxic units) were calculated from equivalent

402

pesticide concentrations in water derived from the internal body burden of G. pulex (a) or from pesticide

403

concentrations measured in water samples (b). The increase in tolerance was quantified as the ratio of the local

20 ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29

Environmental Science & Technology

404

EC50/mean EC50 of all populations from non-contaminated streams (TUmax-Int < -3). The grey areas correspond to

405

the 95% confidence interval. (a) R² = 0.72, adjusted R² = 0.70, F = 33.95, residual d.f. = 13, p < 0.001. (b) R² =

406

0.56, adjusted R² = 0.53, F = 16.58, residual d.f. = 13, p = 0.001.

407

408 409

Fig. 5. Relationship between the community composition of

410

macroinvertebrates, expressed as SPEARpesticides, and the increase in

411

tolerance of G. pulex to the insecticide clothianidin, expressed as ratio

412

of the local EC50/mean EC50 of all populations from non-contaminated

413

streams (TUmax-Int < -3). The grey areas correspond to the 95%

414

confidence interval. R² = 0.70, adjusted R² = 0.67, F = 29.74, d.f. = 13,

415

p < 0.001.

416 417

DISCUSSION

418

Pesticide exposure derived from water samples and body burden 21 ACS Paragon Plus Environment

Environmental Science & Technology

419

Both the pesticide exposure derived from water samples after runoff (TUmax-Ext) and from body

420

burden in G. pulex (TUmax-Int) showed that the agricultural and the non-agricultural streams

421

differed in the level of pesticide contamination (Table 1). All agricultural streams showed a

422

TUmax-Int ≥ -1.3 and a TUmax-Ext ≥ -3.2, whereas all non-agricultural streams showed a TUmax-Int ≤

423

-4.3 and a TUmax-Ext ≤ -3.3. According to Liess and von der Ohe,2 a pesticide contamination of

424

TUmax > -3 results in considerable chronic effects on macroinvertebrates, and such an effect was

425

confirmed by an increased dominance of tolerant macroinvertebrate taxa (low SPEARpesticides

426

values) and an increased insecticide tolerance of G. pulex in agricultural streams. However,

427

within the agricultural and the non-agricultural streams, no clear correlation was observed

428

between TUmax-Int and TUmax-Ext, which presumably results from the uncertainty associated with

429

both methods.

430

The neonicotinoid insecticides provided the highest toxicity (both in organisms and in the stream

431

water), confirming previous pesticide measurements from water samples in the study area.3, 10

432

Inostroza, et al.42 measured the pesticide body burden in G. pulex from the nearby Holtemme

433

River and identified neonicotinoids as the most toxic compounds. Neonicotinoids represent the

434

class of insecticides that is currently sold in the largest quantities in the study area43 and

435

worldwide.44 Although Inostroza, et al.42 reported no pesticide measurements from water

436

samples, the derived toxicity from body burden (TUmax-Int = -2.7 to -0.07) was generally lower

437

than that in our study. This difference might result from the fact that the highest pesticide

438

concentration is typically found in the small upper reaches of streams due to the dilution effects

439

in larger rivers.45, 46

440

In the present study, we found that in agricultural streams, the TUmax-Int calculated from

441

estimated equilibrium water concentrations derived from body burden was approximately one

22 ACS Paragon Plus Environment

Page 22 of 29

Page 23 of 29

Environmental Science & Technology

442

order of magnitude higher than the TUmax-Ext derived from measured water concentrations (Table

443

1). This difference was most likely due to the uncertainty of calculated partition coefficients

444

between G. pulex and the water phase but might also reflect disequilibrium at highly dynamic

445

water concentrations during a peak event. As shown previously, neonicotinoids have a tendency

446

to exhibit a higher bioaccumulation than that estimated with existing partitioning models.26 The

447

availability of effect concentrations and thus TUs based on internal concentrations in G. pulex

448

would largely remove this uncertainty and provide a better basis for the estimation of effects on

449

populations and communities. However, good correlation between body burden and community

450

effect (SPEAR) shows that even with uncertainties resulting from back calculation to water

451

concentrations, internal concentrations are a useful and robust approach in the assessment of

452

pesticide effects.

453 454

Comparison of pesticide exposure with pesticide effects

455

The effects of pesticide exposure on the macroinvertebrate community composition

456

(SPEARpesticides) were similar to the findings of previous studies in which pesticide exposure was

457

quantified from water and sediment samples.3,

458

(TUmax-Int) explained the changes in the community composition equally as effectively as the

459

toxicity derived from water samples (TUmax-Ext). These results suggest that pesticide

460

concentrations in both water samples after run-off and macroinvertebrates can be used for an

461

adequate assessment of long-term pesticide exposure in streams.

462

In contrast, the TUmax-Int explained the increase in tolerance to the insecticide clothianidin in G.

463

pulex considerably better (adj. R² = 0.70) than the TUmax-Ext (adj. R² = 0.53). This difference

464

resulted primarily from a particularly low TUmax-Ext in a single agricultural stream (Agri-2, Table

7, 47

The toxicity derived from body burden

23 ACS Paragon Plus Environment

Environmental Science & Technology

465

1), whereas TUmax-Int, SPEARpesticides and the insecticide tolerance consistently indicated a high

466

pesticide exposure in this stream. This observation suggests that assessing the pesticide exposure

467

of a stream from water samples was in our case less reliable than the assessment based on body

468

burdens or on observable effects on macroinvertebrates. Pesticide peaks in the stream water may

469

be more subject to unpredictable short-term variation caused by the erratic combination of

470

rainfall and pesticide applications than to the other exposure indicators that rely on longer-term

471

processes. However, peak concentrations in the water still provided an adequate measure of

472

pesticide exposure as the difference from the other indicators became obvious only in one of 15

473

streams (Fig. 2b, Fig. 4b).

474

Consistent with an earlier study,10 we observed that the increase in insecticide tolerance of G.

475

pulex with TUmax-Int was less strong when the local species diversity was high. However, this

476

effect was not significant, probably due to the limited number of study sites compared with

477

Becker and Liess.10 In contrast, Russo, et al.48 even reported increased sensitivity to

478

esfenvalerate in G. pulex populations collected from agricultural streams. Pesticides may weaken

479

exposed organisms,49 rendering them more sensitive to subsequent exposure unless they can

480

adapt to the toxicant. Refuge sections can serve as a source of immigrating sensitive organisms

481

and increase the local species diversity, both of which undermine the local adaptation to

482

pesticides.10, 11, 50 The higher availability of refuge sections at the sites studied by Russo et al.48

483

may therefore explain the increased sensitivity instead of adaptation, as opposed to the

484

populations investigated in Shahid et al.11 and the present study. Another study also reported

485

increased sensitivity51 to thiacloprid in G. pulex populations experiencing pre-exposure to

486

wastewater. Organisms subjected to press disturbance are more likely less tolerant to

487

superimposed pulse disturbances.52 Interestingly, with increasing duration of exposure,

24 ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29

Environmental Science & Technology

488

gammarids subjected to pure wastewater showed a slight increase in EC50 between 4 to 6 weeks.

489

It might be an indication towards evolution of co-tolerance.53

490

The increase in the insecticide tolerance of G. pulex in agricultural streams was closely

491

correlated with changes in the macroinvertebrate community composition. Both effects of

492

pesticide exposure likely represent a selection process for more tolerant organisms that acts

493

simultaneously on the population and the community levels.54 Therefore, as previously suggested

494

by Luoma,55 the local increase in pesticide tolerance of a rather tolerant species such as G.

495

pulex2 can be used as an approximation

496

community composition and vice versa.

497

We conclude that pesticide residues in macroinvertebrates are suitable to assess the overall

498

pesticide exposure and effects in agricultural streams. This method appears to be at least as

499

reliable as the assessment based on pesticide concentrations in stream water after run-off events:

500

Water concentrations underestimated the high toxicity in one stream that was indicated by body

501

burdens and by the adaptation of macroinvertebrates to insecticides at the species and

502

community levels. In contrast, the back-calculation from body burden to water concentration

503

tended to overestimate the pesticide toxicity in highly exposed streams. However, both

504

assessment methods generally provided adequate results that were confirmed by the observed

505

pesticide effects on macroinvertebrates. We suggest that the comparison between body burden

506

and synthetic passive samplers is an important area for future investigation.

of long-term effects on the macroinvertebrate

507 508

ACKNOWLEDGEMENTS

509

We thank Oliver Kaske and Klaus Seyfarth from the Department of System-Ecotoxicology,

510

Helmholtz Centre for Environmental Research GmbH – UFZ, and Stefanie Lippmann from the

25 ACS Paragon Plus Environment

Environmental Science & Technology

511

University of Koblenz-Landau for their support in the collection and identification of

512

macroinvertebrates. We also acknowledge German Academic Exchange Service (Deutscher

513

Akademischer Austauschdienst, DAAD) for financially supporting to N.S through doctoral

514

fellowship. A free academic license for JChem, InstantJChem and the Calculator Plugins was

515

kindly provided by ChemAxon (Budapest, Hungary).

516 517

AUTHOR CONTRIBUTIONS

518

N.S. collected the water samples and the invertebrates for the chemical analysis and SPEAR

519

calculation, conducted the sensitivity tests, analyzed the data and wrote the first draft of the

520

manuscript. J.M.B. contributed to the data analyses, interpretation of results and writing of the

521

manuscript. M.K. and W.B. performed the chemical analyses and contributed to the writing of

522

the manuscript. M.L. conceived the concept of research and contributed to the interpretation of

523

results and writing of the manuscript.

524

SUPPORTING INFORMATION

525

Figures showing the study map, illustration of event-driven water sampler (EDS), comparison of

526

the ecological status of streams, and effect of species diversity on the adaptation to pesticides.

527

Tables showing physicochemical properties of investigated streams, compound descriptors, list

528

of median lethal concentrations and reference organisms used in TU calculation, and lists of

529

pesticides analysed and detected in water and G. pulex (PDF).

530 531

REFERENCES

532 533

1. Beketov, M. A.; Kefford, B. J.; SchÄfer, R. B.; Liess, M., Pesticides reduce regional biodiversity of stream invertebrates. Proc. Natl. Acad. Sci. 2013, 110, (27), 11039-11043.

26 ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29

534 535 536 537 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

Environmental Science & Technology

2. Liess, M.; von der Ohe, P. C., Analyzing effects of pesticides on invertebrate communities in streams. Environ Toxicol Chem 2005, 24, (4), 954-965. 3. Münze, R.; Orlinskiy, P.; Gunold, R.; Paschke, A.; Kaske, O.; Beketov, M. A.; Hundt, M.; Bauer, C.; Schüürmann, G.; Möder, M., Pesticide impact on aquatic invertebrates identified with Chemcatcher® passive samplers and the SPEAR pesticides index. Sci. Total Environ. 2015, 537, 69-80. 4. Peters, K.; Bundschuh, M.; SchÄfer, R. B., Review on the effects of toxicants on freshwater ecosystem functions. Environ Pollut 2013, 180, 324-9. 5. Bowmer, K. H., Ecosystem effects from nutrient and pesticide pollutants: catchment care as a solution. Resources 2013, 2, (3), 439-456. 6. SchÄfer, R. B.; Caquet, T.; Siimes, K.; Mueller, R.; Lagadic, L.; Liess, M., Effects of pesticides on community structure and ecosystem functions in agricultural streams of three biogeographical regions in Europe. Sci. Total Environ. 2007, 382, (2), 272-285. 7. SchÄfer, R. B.; von der Ohe, P. C.; Rasmussen, J.; Kefford, B. J.; Beketov, M. A.; Schulz, R.; Liess, M., Thresholds for the effects of pesticides on invertebrate communities and leaf breakdown in stream ecosystems. Environ. Sci. Technol. 2012, 46, (9), 5134-5142. 8. Liess, M.; Foit, K.; Knillmann, S.; SchÄfer, R. B.; Liess, H.-D., Predicting the synergy of multiple stress effects. Sci. Rep. 2016, 6, 32965. 9. Liess, M.; Foit, K.; Becker, A.; Hassold, E.; Dolciotti, I.; Kattwinkel, M.; Duquesne, S., Culmination of low-dose pesticide effects. Environ. Sci. Technol. 2013, 47, (15), 8862-8868. 10. Becker, J. M.; Liess, M., Species diversity hinders adaptation to toxicants. Environ. Sci. Technol. 2017, 51, (17), 10195-10202. 11. Shahid, N.; Becker, J. M.; Krauss, M.; Brack, W.; Liess, M., Adaptation of Gammarus pulex to agricultural insecticide contamination in streams. Sci. Total Environ. 2018, 621, 479-485. 12. Liess, M.; Schulz, R.; Liess, M.-D.; Rother, B.; Kreuzig, R., Determination of insecticide contamination in agricultural headwater streams. Water Res. 1999, 33, (1), 239-247. 13. Liess, M.; Schulz, R., Linking insecticide contamination and population response in an agricultural stream. Environ Toxicol Chem 1999, 18, (9), 1948-1955. 14. Allan, I. J.; Booij, K.; Paschke, A.; Vrana, B.; Mills, G. A.; Greenwood, R., Field performance of seven passive sampling devices for monitoring of hydrophobic substances. Environ. Sci. Technol. 2009, 43, (14), 5383-5390. 15. SchÄfer, R. B.; Paschke, A.; Vrana, B.; Mueller, R.; Liess, M., Performance of the Chemcatcher® passive sampler when used to monitor 10 polar and semi-polar pesticides in 16 Central European streams, and comparison with two other sampling methods. Water Res. 2008, 42, (10), 2707-2717. 16. Gunold, R.; SchÄfer, R. B.; Paschke, A.; Schüürmann, G.; Liess, M., Calibration of the Chemcatcher® passive sampler for monitoring selected polar and semi-polar pesticides in surface water. Environ. Pollut. 2008, 155, (1), 52-60. 17. Fernández, D.; Vermeirssen, E. L. M.; Bandow, N.; Muñoz, K.; SchÄfer, R. B., Calibration and field application of passive sampling for episodic exposure to polar organic pesticides in streams. Environ. Pollut. 2014, 194, 196-202. 18. Greenwood, J. L.; Rosemond, A. D.; Wallace, J. B.; Cross, W. F.; Weyers, H. S., Nutrients stimulate leaf breakdown rates and detritivore biomass: bottom-up effects via heterotrophic pathways. Oecologia 2007, 151, (4), 637-49. 19. Leonard, A. W.; Hyne, R. V.; Lim, R. P.; Pablo, F.; Van den Brink, P. J., Riverine endosulfan concentrations in the Namoi River, Australia: Link to cotton field runoff and macroinvertebrate population densities. Environ Toxicol Chem 2000, 19, (6), 1540-1551. 20. Sprague, J., Measurement of pollutant toxicity to fish. II. Utilizing and applying bioassay results. Water Res. 1970, 4, (1), 3-32.

27 ACS Paragon Plus Environment

Environmental Science & Technology

581 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

21. Huckins, J. N.; Manuweera, G. K.; Petty, J. D.; Mackay, D.; Lebo, J. A., Lipid-containing semipermeable membrane devices for monitoring organic contaminants in water. Environ. Sci. Technol. 1993, 27, (12), 2489-2496. 22. Escher, B. I.; Hermens, J. L., Modes of action in ecotoxicology: their role in body burdens, species sensitivity, QSARs, and mixture effects. Environ. Sci. Technol. 2002, 36, (20), 4201-4217. 23. Inostroza, P. A.; Massei, R.; Wild, R.; Krauss, M.; Brack, W., Chemical activity and distribution of emerging pollutants: Insights from a multi-compartment analysis of a freshwater system. Environ. Pollut. 2017, 231, (Pt 1), 339-347. 24. Ulrich, N.; Endo, S.; Brown, T. N.; Watanabe, N.; Bronner, G.; Abraham, M. H.; Goss, K. U., UFZLSER database v 3.2 [Internet]. 2017. 25. Huber, A.; Bach, M.; Frede, H., Pollution of surface waters with pesticides in Germany: modeling non-point source inputs. Agric., Ecosyst. Environ. 2000, 80, (3), 191-204. 26. Inostroza, P. A.; Wicht, A.-J.; Huber, T.; Nagy, C.; Brack, W.; Krauss, M., Body burden of pesticides and wastewater-derived pollutants on freshwater invertebrates: Method development and application in the Danube River. Environ. Pollut. 2016, 214, 77-85. 27. Fredrickson, L. H.; Reid, F. A., 13.1. 1. Nutritional Values of Waterfowl Foods. 1988. 28. Ashauer, R.; Caravatti, I.; Hintermeister, A.; Escher, B. I., Bioaccumulation kinetics of organic xenobiotic pollutants in the freshwater invertebrate Gammarus pulex modeled with prediction intervals. Environ Toxicol Chem 2010, 29, (7), 1625-1636. 29. Nitchals, P., Effects of temperature on growth in Gammarus pseudolimnaeus. Project Report of Undergraduate Creative Activities and Research Experience, Program at University of Nebraska Lincoln, available online at http://biosci. unl. edu/forebeslab/Nitchals 2014. 30. Gee, J. H., Population dynamics and morphometries of Gammarus pulex L.: evidence of seasonal food limitation in a freshwater detritivore. Freshwat. Biol. 1988, 19, (3), 333-343. 31. Massei, R.; Busch, W.; Wolschke, H.; Schinkel, L.; Bitsch, M.; Schulze, T.; Krauss, M.; Brack, W., Screening of pesticide and biocide patterns as risk drivers in sediments of major European river mouths: ubiquitous or river basin-specific contamination? Environ. Sci. Technol. 2018, 52, (4), 2251-2260. 32. USEPA, ECOTOX Database. United States Environmental Protection Agency

609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625

(USEPA), Washington, DC, USA (Available from: http://cfpub.epa.gov/ecotox/quickquery.htm). 2014. 33. University of Hertfordshire, The Pesticide Properties DataBase (PPDB) Developed by the Agriculture & Environment Research Unit (AERU). University of Hertfordshire (Available from: http://sitem.herts.ac.uk/aeru/iupac/atoz.htm). 2014. 34. Department of System-Ecotoxicology SPEAR Calculator Desktop Application. Helmholtz Centre for Environmental Research, Leipzig, Germany (Available from: http://www.systemecology.eu/spear/). 0.10.0; UFZ, Leipzig, Germany, 2005. 35. Klüttgen, B.; Dülmer, U.; Engels, M.; Ratte, H., ADaM, an artificial freshwater for the culture of zooplankton. Water Res. 1994, 28, (3), 743-746. 36. Kefford, B. J., Rapid Tests for Community-Level Risk Assessments in Ecotoxicology. In Encyclopedia of Aquatic Ecotoxicology, Springer: 2013; pp 957-966. 37. OECD, Test No. 202: Daphnia sp. Acute Immobilisation Test. OECD Publishing: 2004. 38. Simon-Delso, N.; Amaral-Rogers, V.; Belzunces, L. P.; Bonmatin, J.-M.; Chagnon, M.; Downs, C.; Furlan, L.; Gibbons, D. W.; Giorio, C.; Girolami, V., Systemic insecticides (neonicotinoids and fipronil): trends, uses, mode of action and metabolites. Environmental Science and Pollution Research 2015, 22, (1), 5-34. 39. RStudio RStudio: integrated development for R, 2016.

28 ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29

626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665

Environmental Science & Technology

40. R Core Team R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2013. R Foundation for Statistical Computing. (Aavailable from http://www.rproject.org/), 2017. 41. Beketov, M.; Foit, K.; SchÄfer, R.; Schriever, C.; Sacchi, A.; Capri, E.; Biggs, J.; Wells, C.; Liess, M., SPEAR indicates pesticide effects in streams–comparative use of species-and family-level biomonitoring data. Environ. Pollut. 2009, 157, (6), 1841-1848. 42. Inostroza, P. A.; Vera-Escalona, I.; Wicht, A.-J.; Krauss, M.; Brack, W.; Norf, H., Anthropogenic Stressors Shape Genetic Structure: Insights from a Model Freshwater Population along a Land Use Gradient. Environ. Sci. Technol. 2016, 50, (20), 11346-11356. 43. BVL Inlandsabsatz und Export von Pflanzenschutzmitteln - Jahresbericht 2014; Federal Office of Consumer Protection and Food Safety Germany (BVL): 2015; p 12. 44. Jeschke, P.; Nauen, R.; Schindler, M.; Elbert, A., Overview of the status and global strategy for neonicotinoids. J. Agric. Food Chem. 2010, 59, (7), 2897-2908. 45. Szöcs, E.; Brinke, M.; Karaoglan, B.; SchÄfer, R. B., Large Scale Risks from Agricultural Pesticides in Small Streams. Environ. Sci. Technol. 2017, 51, (13), 7378-7385. 46. Schulz, R., Field studies on exposure, effects, and risk mitigation of aquatic nonpoint-source insecticide pollution. J. Environ. Qual. 2004, 33, (2), 419-448. 47. Hunt, L.; Bonetto, C.; Marrochi, N.; Scalise, A.; Fanelli, S.; Liess, M.; Lydy, M. J.; Chiu, M. C.; Resh, V. H., Species at Risk (SPEAR) index indicates effects of insecticides on stream invertebrate communities in soy production regions of the Argentine Pampas. Sci. Total Environ. 2017, 580, (Supplement C), 699709. 48. Russo, R.; Becker, J. M.; Liess, M., Sequential exposure to low levels of pesticides and temperature stress increase toxicological sensitivity of crustaceans. Sci. Total Environ. 2018, 610, 563569. 49. Ashauer, R.; Boxall, A. B.; Brown, C. D., Modeling combined effects of pulsed exposure to carbaryl and chlorpyrifos on Gammarus pulex. Environ. Sci. Technol. 2007, 41, (15), 5535-5541. 50. Gassmann, A. J.; Carrière, Y.; Tabashnik, B. E., Fitness costs of insect resistance to Bacillus thuringiensis. Annu. Rev. Entomol. 2009, 54, 147-63. 51. Zubrod, J. P.; Englert, D.; Luderwald, S.; Poganiuch, S.; Schulz, R.; Bundschuh, M., History Matters: Pre-Exposure to Wastewater Enhances Pesticide Toxicity in Invertebrates. Environ Sci Technol 2017, 51, (16), 9280-9287. 52. Parkyn, S. M.; Collier, K. J., Interaction of press and pulse disturbance on crayfish populations: flood impacts in pasture and forest streams. Hydrobiologia 2004, 527, (1), 113-124. 53. Lopes, I.; Baird, D. J.; Ribeiro, R., Genetically determined resistance to lethal levels of copper by Daphnia longispina: association with sublethal response and multiple/coresistance. Environ Toxicol Chem 2005, 24, (6), 1414-1419. 54. Grant, A., Pollution-Tolerant Species and Communities: Intriguing Toys or Invaluable Monitoring Tools? Human and Ecological Risk Assessment: An International Journal 2002, 8, (5), 955-970. 55. Luoma, S. N., Detection of trace contaminant effects in aquatic ecosystems. Journal of the Fisheries Board of Canada 1977, 34, (3), 436-439.

666

29 ACS Paragon Plus Environment