Carry-Over Effects Across Metamorphosis of a Pesticide on Female

Carry-Over Effects Across Metamorphosis of a Pesticide on Female Lifetime Fitness Strongly Depend on Egg Hatching Phenology: A Longitudinal Study unde...
2 downloads 4 Views 787KB Size
Subscriber access provided by READING UNIV

Article

Carry-over effects across metamorphosis of a pesticide on female lifetime fitness strongly depend on egg hatching phenology: a longitudinal study under seminatural conditions Nedim Tüzün, and Robby Stoks Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04399 • Publication Date (Web): 07 Nov 2017 Downloaded from http://pubs.acs.org on November 12, 2017

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

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

Page 1 of 31

Environmental Science & Technology

1

Carry-over effects across metamorphosis of a pesticide on female lifetime fitness strongly

2

depend on egg hatching phenology: a longitudinal study under seminatural conditions

3

4

Nedim Tüzün* and Robby Stoks

5

6

Evolutionary Stress Ecology and Ecotoxicology, University of Leuven, Deberiotstraat 32, B-

7

3000 Leuven, Belgium

8

9

*

Corresponding author: [email protected]

1 ACS Paragon Plus Environment

Environmental Science & Technology

10 11

Abstract art

2 ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Environmental Science & Technology

12

Abstract

13

Current ecological risk assessment of pesticides fails to protect aquatic biodiversity. For the first

14

time, we tested two potential reasons for this failure with regard to carry-over effects across

15

metamorphosis: their dependence on hatching period, and the lack of studies quantifying adult

16

fitness under seminatural conditions. Using the damselfly Coenagrion puella sampled from six

17

populations, we designed an outdoor longitudinal one-year study starting from the egg stage. We

18

exposed the aquatic larvae to the pesticide esfenvalerate (0.11 µg/L) during the initial microcosm

19

part. Next, we monitored the lifetime fitness of the terrestrial adults in an insectary. Exposure to

20

the pesticide negatively impacted not only larval traits, but also drastically reduced lifetime

21

mating success of adult females. The impact of this post-metamorphic effect of the pesticide on

22

the population level was three times more important than the effects in the larval stage.

23

Importantly, this carry-over effect was only present in females that hatched early in the season,

24

and was not mediated by metamorphic traits (age and mass at emergence). We provide proof-of-

25

principle under seminatural conditions for two potential pitfalls that need to be considered when

26

improving risk assessment: carry-over effects on adult fitness can (i) be much more important

27

than effects during the larval stage and may not be captured by metamorphic traits, and (ii) be

28

strongly modulated by egg hatching dates.

29

30

Introduction

31

Improving ecological risk assessment of pesticides is a major challenge as it currently fails to

32

protect aquatic biodiversity.1,2 One possible reason for this failure is that pollutant effects may

33

persist after metamorphosis in the many animals with a complex life cycle.3,4 The fate of a

3 ACS Paragon Plus Environment

Environmental Science & Technology

34

pollutant-exposed population heavily depends on lifetime fitness values of its individuals, and

35

cannot be accurately judged by pre-adult traits.5 Studies of pollutant effects across

36

metamorphosis, however, rarely explicitly quantified adult fitness components (but see ref 4).

37

Instead, they focused on metamorphic traits such as age and mass at metamorphosis which do

38

not accurately predict the effects of larval stressors on adult fitness.6–8

39

Another major factor limiting risk assessment is that studies on carry-over effects of

40

pesticides are mostly limited to one population sampled at one moment in time. As a result, the

41

widespread temporal variation in egg hatching dates within populations6,9 is neglected in

42

ecotoxicology. The time constraints imposed on late-hatching larvae can, however, strongly

43

modify carry-over effects of ecological stressors.6,9–11 This is usually assumed to be a result of

44

reallocation of resources due to the accelerated growth and development imposed by time

45

constraints.

46

To advance our insights in pesticide-induced carry-over effects, we carried out a

47

longitudinal one-year study starting from the egg stage where we exposed the aquatic larvae of a

48

damselfly to a pesticide pulse and studied not only effects during the larval stage and at

49

metamorphosis, but also explicitly quantified fitness after metamorphosis in the terrestrial adult

50

stage. To specifically assess the importance of within-season temporal variation modulating

51

carry-over effects, we exposed in each of the six sampled populations two cohorts that differed

52

ca. 40 days in their natural hatching period. Given that effects of pesticides are ideally measured

53

under conditions as natural as possible,12 and to obtain realistic estimates of fitness (cfr. ref 13),

54

we integrated an outdoor microcosm part for the larval stage with a large outdoor insectary part

55

for the adult stage. As pesticide, we chose the pyrethroid insecticide esfenvalerate, one of the

56

most widely applied pyrethroids14,15 that can cause carry-over effects in damselflies.16,17 We 4 ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Environmental Science & Technology

57

applied three pulses (with three-day intervals) of esfenvalerate with a concentration of 0.11 µg/L.

58

As study species we chose the damselfly Coenagrion puella which is common in Europe,18 and

59

known to react to time constraints by accelerating growth and development rates.19

60

Materials and methods

61

Study species and populations

62

To capture regional variation in life history, we sampled damselflies from six populations,

63

situated in both urban and more natural settings in Flanders, Belgium (Appendix S1). We

64

collected eggs from each of the six populations in each of two periods, 13 – 22 June (“early

65

period”), and 23 – 31 July 2014 (“late period”), matching the peak and the end of the flight

66

season in Flanders, respectively. This resulted in two cohorts differing ca. 40 days in egg

67

hatching time. Given the average lifespan of adult C. puella is ca. 20 days20 (see also Results),

68

the two collected cohorts derived from parents that do not overlap in flight period. We collected

69

eggs from 5-10 mated females from each pond in each period. Eggs from the early and late

70

periods hatched during 4 – 12 July and 16 – 19 August 2014, respectively. To increase survival,

71

we reared the larvae in laboratory conditions for the first three weeks. Larvae were kept per

72

female together in plastic containers filled with 1.5 L of dechlorinated tap water in a

73

temperature-controlled room at 20 °C with a photoperiod of 14:10 h light:dark. Larvae were fed

74

Artemia nauplii ad libitum five days a week.

75

General experimental procedure

76

Following the three weeks period in the laboratory, the longitudinal one-year study consisted of

77

two integrated outdoor parts. In the first, outdoor microcosm part, we set up a full factorial

78

experiment consisting of 2 hatching periods (early and late) × 2 pesticide treatments (control and 5 ACS Paragon Plus Environment

Environmental Science & Technology

79

0.2 µg/L esfenvalerate). This part ran from summer 2014 until spring 2015 (with a short indoor

80

part to avoid winter freezing; see below). The pesticide application was carried out in May 2015,

81

which is within the main application period of agricultural pesticides in the study region.21

82

During this period all animals were still in the larval stage. As in nature, all larvae were

83

simultaneously exposed to the pesticide irrespective of their egg hatching date the previous

84

summer. In the second, outdoor insectary part, adults that emerged from the microcosms were

85

placed in a single outdoor insectary where we monitored adult life-history traits and lifetime

86

mating success. This part ran from late spring to summer 2015.

87

Outdoor microcosm part

88

Each of the four hatching period × pesticide treatment combinations was replicated in twelve

89

microcosms, giving a total of 48 microcosms. The set of twelve replicated microcosms per

90

treatment combination consisted each time of two microcosms of each of the six populations.

91

Hence, each microcosm contained larvae from a single population. Microcrosms were 10 L

92

polypropylene containers filled with 3 L dechlorinated tap water and 3 L of water from an

93

adjacent pond (filtered through a 500 µm mesh). Microcosms were placed at an outdoor

94

experimental area in Heverlee, Belgium.

95

To allow the development of Daphnia and protozoan populations (food resource for the

96

small damselfly larvae), the microcosms were set up ca. two weeks prior to the introduction of

97

the larvae. We inoculated each container with ca. 300 D. magna. In addition, we added grass to

98

stimulate growth of protozoa. We covered each container with a net to prevent predators entering

99

and adult damselflies escaping.

6 ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

100

Environmental Science & Technology

We introduced 30 C. puella larvae in each container (in total 1,440 larvae), corresponding

101

to typical field densities of coenagrionid damselfly larvae.22 We introduced the larvae of each

102

period when ca. three weeks old (26-27 July 2014 for the early period and 6-7 September 2014

103

for the late period), thereby maintaining the ca. 40 day difference in hatching dates of the two

104

cohorts. Throughout the outdoor microcosm experiment, we provided larvae with ad libitum

105

food by adding weekly ca. 300 D. magna obtained from outdoor stock tanks. We reduced the

106

feeding frequency to biweekly during periods when temperatures dropped below 10 °C, as larval

107

growth of the study species ceases at these temperatures.23 To prevent freezing, we moved the

108

microcosms to an unheated indoor facility for 45 days during 27 December 2014 – 12 February

109

2015, when outdoor temperatures dropped below 0 °C. Daily mean water temperatures when

110

placed inside were ca. 4.3 °C (see also Fig. S2 in Appendix S2). Afterwards, microcosms were

111

placed again outside.

112

Starting 4 May 2015, we applied three pesticide pulses to half of the microcosms, with

113

three days between pulses. This simulates spring pesticide applications and the associated run-off

114

in nearby surface waters.24 The nominal esfenvalerate concentration of 0.2 µg/L induces

115

mortality and reduces growth rate in another Coenagrion damselfly.16 Although this

116

concentration exceeds the predicted environmental concentrations (e.g. 0.06 µg/L for a realistic

117

scenario25), it falls within the range of esfenvalerate concentrations detected in natural water

118

bodies, e.g. in Denmark (up to 0.66 µg/L, ref 26) and the USA (up to 0.8 µg/L, refs 27–29).

119

Moreover, regulatory surface water models, even when simulating most realistic scenarios,

120

strongly underpredict measured field concentrations; especially for hydrophobic insecticides

121

such as esfenvalerate.30 Importantly, our main aim was to test for a proof-of-principle for carry-

122

over effects of pesticides and their dependence on the hatching period, rather than fine tuning the

7 ACS Paragon Plus Environment

Environmental Science & Technology

123

specific risk assessment of esfenvalerate. We initially prepared a 2 × 105 µg/L stock solution by

124

dissolving 1 mg esfenvalerate powder (purity >99%, Sigma-Aldrich) in 5 mL absolute ethanol.

125

This stock solution was further diluted with filtered pond water (mesh size: 500 µm) to obtain a

126

spraying solution with a concentration of 24 µg/L. Fifty millilitre of this spraying solution was

127

gently poured over the surface of the microcosms to obtain the nominal esfenvalerate

128

concentrations of 0.2 µg/L. For the control treatment, we added 50 mL of ethanol dissolved in

129

filtered pond water with a concentration of 24 µL/L ethanol, i.e. the ethanol concentration of the

130

esfenvalerate treatment, to the microcosms. Growth and behavior of damselfly larvae are not

131

affected by ethanol concentrations up to 5 × 103 µL/L (Lizanne Janssens, unpublished data). The

132

esfenvalerate concentration in the containers 20 min after application was 0.11 µg/L, whereas

133

after three days, just before the next pulse, the concentration was below the detection limit (
3 standard deviations from the

206

mean), and was excluded from all analyses. For the analysis of post-exposure survival, we

207

identified one microcosm from the late period - control treatment combination as an outlier (>3

208

standard deviations from the mean) and excluded this microcosm from all analyses. To

209

strengthen the interpretation of key results (see Results Lifetime mating success), we additionally

210

calculated effect sizes estimated as Hedges’ d with 95% confidence intervals.38

211

Path analyses

11 ACS Paragon Plus Environment

Environmental Science & Technology

212

We used structural equation modeling (SEM39) to test whether any potential effect of the

213

pesticide treatment on LMS was operating directly, or was indirectly mediated via the

214

metamorphic traits (age and mass at emergence). We specifically applied the path analysis

215

approach, which is a SEM method that deals explicitly with observed variables. Our data

216

included a non-normally distributed response variable (i.e. Poisson-distributed LMS), and was

217

hierarchically structured (i.e. nested random effects). Therefore, we applied the piecewise SEM

218

(or ‘generalized multilevel path analysis’) approach, which allows translating a path diagram into

219

a set of GLMMs.40,41

220

To specifically evaluate the role of direct and indirect effects of the pesticide, we built

221

three candidate models: (1) the pesticide affects fitness only directly; (2) the pesticide affects

222

fitness only indirectly via the metamorphic traits; and (3) the pesticide affects fitness both

223

directly and indirectly (Fig. 1). As the pesticide effect on LMS was modulated by period (see

224

results of GLMMs), we ran separate path models for early and late females with the pesticide

225

treatment as the predictor variable. All models included paths going from age and mass at

226

emergence to LMS. Further, to correct for potential density-mediated effects, we included paths

227

going from larval density per microcosm to age and mass at emergence (results not shown). As

228

the effects on LMS remained when the weather variables were removed from the GLMMs (data

229

not shown), we did not include them in the path analysis. We included correlated errors between

230

age and mass at emergence to capture the potential trade-off between both variables. The error

231

structure and random effects of the GLMMs used to build the path models were identical to the

232

ones used in the mean analyses (see above). For all path models, the categorical variable

233

‘pesticide treatment’ was binary coded: control = 0, pesticide = 1. Hence, a negative path going

12 ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31

Environmental Science & Technology

234

from pesticide treatment to LMS, for example, should be interpreted as decreased LMS for

235

pesticide-exposed females.

236

We selected the most parsimonious among the three candidate path models using the

237

AICc (Akaike Information Criterion, corrected for small sample size) implemented for path

238

models.42 For this, we calculated and compared ∆AICc scores and relative support for each

239

model (AICc weights), where lower AICc scores and higher AICc weights indicate better

240

models.43 The overall fit of path models was evaluated with Shipley’s test of d-separation,40

241

which tests for missing paths in the model. Using the combined significance of these unrealized

242

paths, a Fisher’s C statistic is calculated and compared to a χ2 distribution to reject (P < 0.05) or

243

accept (P > 0.05) the model.

244

All analyses were performed using R version 3.3.2.44 We used the ‘lme4’45 and ‘nlme’46

245

packages for GLMMs, the ‘car’ package47 for testing fixed effects with Wald- χ2, the

246

‘piecewiseSEM’ package41 for conducting path analyses, the ‘MuMIn’ package48 to calculate

247

AICc weights, and the ‘effsize’ package49 to estimate Hedges’ d. We report standardized

248

coefficients for the path models. Mass at emergence and lifespan were log-transformed to meet

249

model assumptions.

250

Results and discussion

251

Below we report effects of pesticide exposure and hatching period on life-history traits. For full

252

details of the statistical analyses, including results for covariates, see Appendix S4.

253

Larval life history and metamorphic traits

13 ACS Paragon Plus Environment

Environmental Science & Technology

254

Larval survival prior to pesticide exposure was ca. 79% and independent of hatching period.

255

During the pesticide exposure period, survival was reduced from ca. 89% (of the 79% pre-

256

exposure survival) in the control group to ca. 81% in the pesticide-exposed group (χ12 = 11.60, P

257

= 0.001, Fig. 2a). Hatching period (χ12 = 2.87, P = 0.090) and its interaction with pesticide

258

treatment (χ12 = 0.43, P = 0.514) did not influence survival during the pesticide exposure period.

259

Pesticide exposure reduced larval growth (χ12 = 21.15, P < 0.001, Fig. 2b). Larvae from the late

260

period (‘late larvae’) had faster growth rates during the exposure period than larvae from the

261

early period (‘early larvae’) (χ12 = 6.73, P = 0.009, Fig. 2b), and this did not depend on the

262

pesticide treatment (pesticide × period: χ12 = 0.11, P = 0.739). Pesticide exposure negatively

263

affected post-exposure survival (χ12 = 8.37, P = 0.004, Fig. 2c). Likewise, late larvae had a lower

264

survival during the post-exposure period than early larvae (χ12 = 12.22, P < 0.001, Fig. 2c). The

265

interaction of pesticide treatment and hatching period did not influence post-exposure survival

266

(χ12 = 1.55, P = 0.212).

267

Age at emergence strongly differed between individuals of the two hatching periods (χ12

268

= 172.25, P < 0.001), with late females developing ca. 28 days faster compared early females

269

(Fig. 3a). Pesticide exposure (χ12 = 1.17, P = 0.279) and its interaction with hatching period (χ12 =

270

0.34, P = 0.560) did not influence age at emergence. Damselflies from the pesticide treatment

271

and late hatching group emerged with a slightly lower mass compared to the control and early

272

group, respectively, yet these effects were marginally non-significant (pesticide treatment: χ12 =

273

3.43, P = 0.064; hatching period: χ12 = 3.12, P = 0.077, Fig. 3b). The interaction of pesticide

274

treatment and hatching period did not influence mass at emergence (χ12 = 0.27, P = 0.600).

275 276

As expected, the chosen dose of esfenvalerate reduced survival both during and after the exposure period. This matches studies at similar esfenvalerate concentrations in aquatic 14 ACS Paragon Plus Environment

Page 14 of 31

Page 15 of 31

Environmental Science & Technology

277

insects,14,50–52 including damselflies.16,17 Likewise, the pesticide reduced the growth rate,

278

possibly due to energy re-allocation to detoxification.52

279

As predicted by life history theory on time constraints,53,54 late-hatched larvae developed

280

and grew faster than early-hatched larvae. There is ample empirical support for this pattern from

281

a wide range of taxa,55 particularly from studies with anurans11,56 and damselflies.57,58 This

282

accelerated life history decreased post-exposure survival. Similarly, survival until emergence

283

was lower in late compared to early larvae in the damselfly Lestes sponsa,59 and may reflect the

284

trade-off between growth and survival.60

285

Lifetime mating success

286

Of the adults released in the insectary, 68 of the 248 females mated at least once. Pesticide

287

exposure and hatching period did not affect lifespan (see Appendix S5). While the main effects

288

of pesticide (χ12 = 0.15, P = 0.707) and period (χ12 = 0.05, P = 0.819) were not significant, female

289

LMS was modulated by a significant pesticide × period interaction (χ12 = 4.38, P = 0.036, Fig.

290

3a,b): previous pesticide exposure reduced LMS of early females with ca. 60% (Hedges’ d = -

291

0.391, 95% CI: [-0.725, -0.058]), while, if anything, the pesticide seemed to have a positive

292

effect on LMS of late females (Hedges’ d = 0.235, 95% CI: [-0.158, 0.627]). This drastic

293

negative pesticide effect on female LMS was not mediated via the metamorphic traits (the ‘direct

294

paths only’ model was the best supported model, Table S5 in Appendix S6). Indeed, the model

295

for early females revealed that the pesticide exposure had a direct negative effect on LMS (as

296

indicated by the significant negative path coefficient, Fig. 3c, Table S6). Further, increasing mass

297

and age at emergence resulted in a higher LMS. The model for late females confirmed the

298

pesticide had no significant effect on their LMS (Fig. 3d, Table S6).

15 ACS Paragon Plus Environment

Environmental Science & Technology

299

A key finding was that the negative effects of larval pesticide exposure carried over to the

300

adult stage and decreased LMS in early females. This result adds to the mounting evidence of

301

carry-over effects and their impact on fitness.61,62 This is in line with effects of esfenvalerate

302

often being detected long after the exposure period ends,50,63–65 and being able to bridge

303

metamorphosis16,66,67 (but see ref 63). The path analysis indicated that this effect was not

304

mediated via pesticide effects on mass or age at emergence. This is in contrast with the

305

assumption that reduction in fitness due to exposure to a stressor should be operating by

306

affecting these two key metamorphic traits.53,54 Although the mechanisms underlying latent

307

effects of pesticides are not fully understood, pesticide-mediated reductions in adult energy

308

reserves50,65 may be responsible. We hypothesize this ‘hidden’ carry-over effect on female

309

fitness is mediated by negative effects of esfenvalerate on fat and flight muscle mass, two flight-

310

related traits that are impaired by this pesticide in damselflies,16,67 thereby reducing the flight

311

performance. Similarly, studies documenting negative effects of contaminants on mating

312

success, for example in fish,68 birds,69 bed bugs,70 and moths71 suggest impaired reproductive

313

behavior as the underlying mechanism. Compared to early females, the LMS of late females did

314

not seem to be as strongly influenced by the pesticide; possibly because more vulnerable late

315

larvae were already eliminated during metamorphosis (i.e. lower post-exposure survival in the

316

late period). The apparent positive effect of the pesticide exposure on late period animals (Fig.

317

4b) may be a hormetic response, not uncommon in ecotoxicology,72 operating via a mechanism

318

not identified in the present study.

319

General implications for risk assessment of pesticides

320

Using a powerful longitudinal design that integrated an exposure part in microcosms with the

321

monitoring of lifetime fitness in a large outdoor insectary, our study provides the first data on 16 ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

Environmental Science & Technology

322

carry-over effects of a pesticide across metamorphosis on adult lifetime fitness under seminatural

323

conditions. Although the here used concentration of esfenvalerate (0.11 µg/L) is ca. two times

324

higher than the predicted environmental concentration based on worst-case scenarios25,

325

regulatory surface water models simulating realistic scenarios strongly underpredict actual

326

measured field concentrations; especially for hydrophobic insecticides such as esfenvalerate.30

327

This is in agreement with several studies reporting substantially higher (up to 0.8 µg/L)

328

esfenvalerate concentrations in natural water bodies.26–29 Nevertheless, our aim was to test under

329

seminatural conditions for mechanisms possibly responsible for the failure of current risk

330

assessment of pesticides in general (detailed below), rather than fine tuning regulatory limits for

331

esfenvalerate.

332

Aquatic habitats are usually dominated by insects with complex life cycles, i.e. with

333

aquatic larvae metamorphosing into terrestrial adults. Although current water quality standards

334

acknowledge latent effects of pesticide exposure,73 they do not explicitly take into account carry-

335

over effects across metamorphosis, hence assume that protection of aquatic stages translate to the

336

protection of adults. We provide a proof-of-principle that this assumption is not met, and suggest

337

that the current guidelines would have erroneously concluded the pesticide effect to be weak if

338

looking only at responses at the larval stage (for metals 4,74). Notably, a simple simulation based

339

on our results and the general relationships between fitness components and the number of

340

offspring that reaches adulthood (based on studies in a natural C. puella population 33,36) revealed

341

that when including the carry-over effects, the pesticide treatment is expected to reduce the

342

number of offspring reaching maturity by more than three times (see Appendix S7 for details of

343

the simulation). Importantly, this simulation shows that the pesticide-induced mortality in the

344

larval stage accounted only for ca. 24% of the total impact on the population in the next

17 ACS Paragon Plus Environment

Environmental Science & Technology

345

generation, while the rest of the impact was captured by variation in LMS (Appendix S7). These

346

delayed effects after metamorphosis could not be predicted by effects on the two key

347

metamorphic traits (age and mass at metamorphosis), thereby violating the assumption of life

348

history theory that metamorphic traits should predict the effects of stressors on adult fitness.75

349

This further complements the few studies reporting that post-metamorphic traits are weak

350

predictors of carry-over effects of ecological stressors on adult fitness.6–8 This argues against

351

using mass and age at metamorphosis as proxies for effects on adult fitness in ecotoxicological

352

studies,76–78 and underscores the importance of identifying the underlying mechanisms of carry-

353

over effects.62,79

354

We further identified an overlooked, yet important temporal component associated with

355

time constraints strongly determining the impact of pesticide exposure. Indeed, the larval

356

pesticide exposure affected lifetime fitness only in early-hatched females, and not in the more

357

time-stressed late-hatched females. This large impact of hatching phenology is striking because

358

pesticide exposure occurred ca. 10 months after egg hatching, thereby simulating the general

359

scenario where aquatic insects that hatched in summer are exposed the next spring to pesticide

360

pulses. Given the widespread variation in hatching periods in natural populations and the

361

associated differences in time constraints,6,9 this is expected to be a general factor that may

362

critically affect the consequences of pollutant exposure on adult fitness. Ecotoxicology tests on

363

field-collected animals (as frequently done in mesocosm tests) typically ignore the hatching date.

364

Hence, the here identified mechanisms may introduce considerable noise when animals differing

365

in hatching period are being tested, and may explain differences between toxicological studies;

366

even when using animals from the same study population.

367 18 ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

Environmental Science & Technology

368

Supporting Information

369

Information on study populations (S1), abiotic parameters (S2), and weather variables (S3), full

370

details of the statistical analyses (S4), results for lifespan (S5), details of the path analyses (S6),

371

explanations for the simulation to estimate pesticide effect at the population level (S7).

372

Acknowledgments

373

We thank Selina Müller for her dedicated help during the experiment, and Geert Neyens for

374

setting up the insectary. Comments from four anonymous reviewers considerably improved our

375

manuscript. Financial support came from the Belspo project SPEEDY (IAP- project P7/04) and

376

research grants from the KU Leuven (PF/2010/07 and C16/17/002) and the FWO research

377

network EVENET.

378

References

379 380 381

(1)

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

382 383

(2)

Stehle, S.; Schulz, R. Agricultural insecticides threaten surface waters at the global scale. Proc. Natl. Acad. Sci. U. S. A. 2015, 112 (18), 5750–5755.

384 385 386

(3)

Wesner, J. S.; Kraus, J. M.; Schmidt, T. S.; Walters, D. M.; Clements, W. H. Metamorphosis enhances the effects of metal exposure on the mayfly, Centroptilum triangulifer. Environ. Sci. Technol. 2014, 48 (17), 10415–10422.

387 388 389

(4)

Debecker, S.; Dinh, K. V.; Stoks, R. Strong delayed interactive effects of metal exposure and warming: Latitude-dependent synergisms persist across metamorphosis. Environ. Sci. Technol. 2017, 51 (4), 2409–2417.

390 391

(5)

Forbes, V. E.; Calow, P. Population growth rate as a basis for ecological risk assessment of toxic chemicals. Philos. Trans. R. Soc. B Biol. Sci. 2002, 357 (1425), 1299–1306.

392 393

(6)

De Block, M.; Stoks, R. Fitness effects from egg to reproduction: Bridging the life history transition. Ecology 2005, 86 (1), 185–197.

394 395

(7)

Earl, J. E.; Whiteman, H. H. Are commonly used fitness predictors accurate? A metaanalysis of amphibian size and age at metamorphosis. Copeia 2015, 103 (2), 297–309.

396

(8)

Van Allen, B. G.; Briggs, V. S.; McCoy, M. W.; Vonesh, J. R. Carry-over effects of the 19 ACS Paragon Plus Environment

Environmental Science & Technology

larval environment on post-metamorphic performance in two hylid frogs. Oecologia 2010, 164 (4), 891–898.

397 398 399 400 401

(9)

Orizaola, G.; Dahl, E.; Nicieza, A. G.; Laurila, A. Larval life history and anti-predator strategies are affected by breeding phenology in an amphibian. Oecologia 2013, 171 (4), 873–881.

402 403 404

(10)

Stoks, R.; De Block, M.; Slos, S.; Van Doorslaer, W.; Rolff, J. Time constraints mediate predator-induced plasticity in immune function, condition, and life history. Ecology 2006, 87 (4), 809–815.

405 406 407

(11)

Orizaola, G.; Richter-Boix, A.; Laurila, A. Transgenerational effects and impact of compensatory responses to changes in breeding phenology on antipredator defenses. Ecology 2016, 97 (9), 2470–2478.

408 409 410

(12)

Boone, M. D.; Bishop, C. A.; Boswell, L. A.; Brodman, R. D.; Burger, J.; Davidson, C.; Gochfeld, M.; Hoverman, J. T.; Neuman-Lee, L. A.; Relyea, R. A.; et al. Pesticide regulation amid the influence of industry. Bioscience 2014, 64 (10), 917–922.

411 412

(13)

Tigreros, N. Linking nutrition and sexual selection across life stages in a model butterfly system. Funct. Ecol. 2013, 27 (1), 145–154.

413 414 415

(14)

Beketov, M. A.; Liess, M. Acute contamination with esfenvalerate and food limitation: Chronic effects on the mayfly, Cloeon dipterum. Environ. Toxicol. Chem. 2005, 24 (5), 1281–1286.

416 417 418 419

(15)

Spurlock, F.; Lee, M. Synthetic pyrethroid use patterns, properties, and environmental effects. In Synthetic Pyrethroids; Gan, J., Spurlock, F., Dendley, P., Weston, D. P., Eds.; ACS Symposium Series; American Chemical Society: Washington DC, 2008; Vol. 991, pp 3–25.

420 421 422

(16)

Dinh Van, K.; Janssens, L.; Therry, L.; Gyulavári, H. A.; Bervoets, L.; Stoks, R. Rapid evolution of increased vulnerability to an insecticide at the expansion front in a polewardmoving damselfly. Evol. Appl. 2016, 9 (3), 450–461.

423 424 425

(17)

Janssens, L.; Tüzün, N.; Stoks, R. Testing the time-scale dependence of delayed interactions: A heat wave during the egg stage shapes how a pesticide interacts with a successive heat wave in the larval stage. Environ. Pollut. 2017, 230, 351–359.

426 427

(18)

Dijkstra, K. B.; Lewington, R. Field guide to the dragonflies of Britain and Europe; British Wildlife Publishing: Gillingham, 2006.

428 429 430

(19)

Mikolajewski, D. J.; De Block, M.; Stoks, R. The interplay of adult and larval time constraints shapes species differences in larval life history. Ecology 2015, 96 (4), 1128– 1138.

431 432 433

(20)

Banks, M. J.; Thompson, D. J. Emergence, longevity and breeding area fidelity in Coenagrion puella (L.)(Zygoptera:Coenagrionidae). Odonatologica 1985, 14 (4), 279– 286.

434 435

(21)

Van Drooge, H. L.; Groeneveld, C. N.; Schipper, H. J. Data on application frequency of pesticide for risk assessment purposes. Ann. Occup. Hyg. 2001, 45 (1001), S95–S101. 20 ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

Environmental Science & Technology

436 437

(22)

Corbet, P. S. Dragonflies: behavior and ecology of Odonata; Cornell University Press: London, 1999.

438 439 440

(23)

Waringer, J. A.; Humpesch, U. H. Embryonic development, larval growth and life cycle of Coenagrion puella (Odonata: Zygoptera) from an Austrian pond. Freshw. Biol. 1984, 14 (4), 385–399.

441 442 443 444

(24)

Dabrowski, J. M.; Peall, S. K. C.; Reinecke, A. J.; Liess, M.; Schulz, R. Runoff-related pesticide input into the lourens river, South Africa: Basic data for exposure assessment and risk mitigation at the catchment scale. Water. Air. Soil Pollut. 2002, 135 (1–4), 265– 283.

445 446

(25)

EFSA (European Food Safety Authority). Conclusion on the peer review of the pesticide risk assessment of the active substance esfenvalerate. EFSA J. 2014, 12 (11), 3873.

447 448

(26)

Aarhus Amt. Pesticider i vandløb, kilder og søer i Aarhus Amt. Technical Report. Aarhus County, Denmark. 1999.

449 450 451

(27)

Cooper, C.; Smith, S.; Moore, M. Surface water, ground water and sediment quality in three oxbow lake watersheds in the Mississippi Delta Agricultural Region: Pesticides. Int. J. Ecol. Environ. Sci. 2003, 29, 171–184.

452 453 454 455 456

(28)

Werner, I.; Zalom, F. G.; Oliver, M. N.; Deanovic, L. A.; Kimball, T. S.; Henderson, J. D.; Wilson, B. W.; Krueger, W.; Wallender, W. W. Toxicity of storm-water runoff after dormant spray application in a French prune orchard, Glenn County, California, USA: Temporal patterns and the effect of ground covers. Environ. Toxicol. Chem. 2004, 23 (11), 2719–2126.

457 458 459 460

(29)

Bacey, J.; Spurlock, F.; Starner, K.; Feng, H.; Hsu, J.; White, J.; Tran, D. M. Residues and toxicity of esfenvalerate and permethrin in water and sediment, in tributaries of the Sacramento and San Joaquin rivers, California, USA. Bull. Environ. Contam. Toxicol. 2005, 74, 864–871.

461 462 463

(30)

Knäbel, A.; Stehle, S.; Schäfer, R. B.; Schulz, R. Regulatory FOCUS surface water models fail to predict insecticide concentrations in the field. Environ. Sci. Technol. 2012, 46 (15), 8397–8404.

464 465

(31)

Adelsbach, T. L.; Tjeerdema, R. S. Chemistry and fate of fenvalerate and esfenvalerate. Rev. Environ. Contam. Toxicol. 2003, 176, 137–154

466 467 468

(32)

Sternberg, K. Coenagrion puella. In Die Libellen Baden-Württembergs, Band 1: Allgemeiner Teil Kleinlibellen (Zygoptera); Sternberg, K.; Buchwald, R., Eds.; Verlag Eugen: Stuttgart, 1999; pp 278–287

469 470 471

(33) Thompson, D. J.; Hassall, C.; Lowe, C. D.; Watts, P. C. Field estimates of reproductive success in a model insect: Behavioural surrogates are poor predictors of fitness. Ecol. Lett. 2011, 14 (9), 905–913.

472 473

(34)

Banks, M. J.; Thompson, D. J. Lifetime mating success in the damselfly Coenagrion puella. Anim. Behav. 1985, 33 (4), 1175–1183.

474

(35)

Thompson, D. J. The effects of survival and weather on lifetime egg production in a 21 ACS Paragon Plus Environment

Environmental Science & Technology

model damselfly. Ecol. Entomol. 1990, 15 (4), 455–462.

475 476 477

(36)

Banks, M. J.; Thompson, D. J. Lifetime reproductive success of females of the damselfly Coenagrion puella. J. Anim. Ecol. 1987, 56 (3), 815–832.

478 479

(37)

Tammaru, T.; Esperk, T. Growth allometry of immature insects: Larvae do not grow exponentially. Funct. Ecol. 2007, 21 (6), 1099–1105.

480 481 482

(38)

Jackson, M. C.; Loewen, C. J. G.; Vinebrooke, R. D.; Chimimba, C. T. Net effects of multiple stressors in freshwater ecosystems: A meta-analysis. Glob. Chang. Biol. 2016, 22 (1), 180–189.

483 484

(39)

Grace, J. B. Structural equation modeling and natural systems; Cambridge University Press: Cambridge, 2006.

485 486

(40)

Shipley, B. Confirmatory path analysis in a generalized multilevel context. Ecology 2009, 90 (2), 363–368.

487 488

(41)

Lefcheck, J. S. piecewiseSEM: Piecewise structural equation modelling in r for ecology, evolution, and systematics. Methods Ecol. Evol. 2016, 7 (5), 573–579.

489 490

(42)

Shipley, B. The AIC model selection method applied to path analytic models compared using a d-separation test. Ecology 2013, 94 (3), 560–564.

491 492

(43)

Burnham, K. P.; Anderson, D. R. Model selection and Multimodel Inference: A Practical Information- Theoretic Approach; Springer: New York, 2002.

493 494

(44)

R Development Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing: Vienna, Austria 2013

495 496

(45)

Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67 (1), 1–48.

497 498

(46)

Pinheiro, J; Bates, D.; DebRoy, S.; Sarkar, D.; R Development Core Team. nlme: Linear and nonlinear mixed effects models. R package version 3.1-128, 2016.

499 500

(47)

Fox, J.; Weisberg, S. An {R} Companion to Applied Regression, 2nd ed.; Sage: Thousand Oaks, CA, 2011.

501

(48)

Barton, K. MuMIn: Multi-model inference. R package version 1.15.6, 2016.

502

(49)

Torchiano, M. effsize: Efficient effect size computation. R package version 0.7.1, 2017.

503 504 505

(50)

Rasmussen, J. J.; Reiber, L.; Holmstrup, M.; Liess, M. Realistic pesticide exposure through water and food amplifies long-term effects in a Limnephilid caddisfly. Sci. Total Environ. 2017, 580, 1439–1445.

506 507 508

(51)

Kunce, W.; Josefsson, S.; Örberg, J.; Johansson, F. Combination effects of pyrethroids and neonicotinoids on development and survival of Chironomus riparius. Ecotoxicol. Environ. Saf. 2015, 122, 426–431.

509 510 511

(52)

Rodrigues, A. C. M.; Gravato, C.; Quintaneiro, C.; Barata, C.; Soares, A. M. V. M.; Pestana, J. L. T. Sub-lethal toxicity of environmentally relevant concentrations of esfenvalerate to Chironomus riparius. Environ. Pollut. 2015, 207, 273–279. 22 ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

Environmental Science & Technology

512 513

(53)

Rowe, L.; Ludwig, D. Size and timing of metamorphosisin complex life cycles: Time constraints and variation. Ecology 1991, 72 (2), 413–427.

514 515 516

(54)

Abrams, P. A.; Leimar, O.; Nylin, S.; Wiklund, C. The effect of flexible growth rates on optimal sizes and development times in a seasonal environment. Am. Nat. 1996, 147 (3), 381–395.

517 518

(55)

Gotthard, K. Growth strategies of ectothermic animals in temperate environments. Environ. Anim. Dev. 2001, 287–304.

519 520 521

(56)

Laugen, A. T.; Laurila, A.; Merilä, J. Latitudinal and temperature-dependent variation in embryonic development and growth in Rana temporaria. Oecologia 2003, 135 (4), 548– 554.

522 523

(57)

De Block, M.; Stoks, R. Life-history variation in relation to time constraints in a damselfly. Oecologia 2004, 140 (1), 68–75.

524 525 526

(58)

Lowe, C. D.; Harvey, I. F.; Watts, P. C.; Thompson, D. J. Reproductive timing and patterns of development for the damselfly Coenagrion puella in the field. Ecology 2009, 90 (8), 2202–2212.

527 528 529

(59)

Dańko, M. J.; Dańko, A.; Golab, M. J.; Stoks, R.; Śniegula, S. Latitudinal and age-specific patterns of larval mortality in the damselfly Lestes sponsa: Senescence before maturity? Exp. Gerontol. 2017, 95, 107–115.

530 531 532

(60)

Lancaster, L. T.; Morrison, G.; Fitt, R. N. Life history trade-offs, the intensity of competition, and coexistence in novel and evolving communities under climate change. Philos. Trans. R. Soc. B Biol. Sci. 2017, 372 (1712), 20160046.

533 534

(61)

Harrison, X. A.; Blount, J. D.; Inger, R.; Norris, D. R.; Bearhop, S. Carry-over effects as drivers of fitness differences in animals. J. Anim. Ecol. 2011, 80 (1), 4–18.

535 536 537

(62)

O’Connor, C.; Norris, D. R.; Crossin, G. T.; Cooke, S. J. Biological carryover effects: linking common concepts and mechanisms in ecology and evolution. Ecosphere 2014, 5 (3), art28.

538 539 540

(63)

Forbes, V. E.; Cold, A. Effects of the pyrethroid esfenvalerate on life-cycle traits and population dynamics of Chironomus riparius—importance of exposure scenario. Environ. Toxicol. Chem. 2005, 24 (1), 78–86.

541 542 543

(64)

Liess, M.; Schulz, R. Chronic effects of short-term contamination with the pyrethroid insecticide fenvalerate on the caddisfly Limnephilus lunatus. Hydrobiologia 1996, 324 (2), 99–106.

544 545 546

(65)

Palmquist, K. R.; Jepson, P. C.; Jenkins, J. J. Impact of aquatic insect life stage and emergence strategy on sensitivity to esfenvalerate exposure. Environ. Toxicol. Chem. 2008, 27 (8), 1728.

547 548 549

(66)

Schulz, R.; Liess, M. Toxicity of aqueous-phase and suspended particle-associated fenvalerate: Chronic effects after pulse-dosed exposure of Limnephilus lunatus (Trichoptera). Environ. Toxicol. Chem. 2001, 20 (1), 185–190.

23 ACS Paragon Plus Environment

Environmental Science & Technology

550 551 552

(67)

Śniegula, S.; Janssens, L.; Stoks, R. Integrating multiple stressors across life stages and latitudes: Combined and delayed effects of an egg heat wave and larval pesticide exposure in a damselfly. Aquat. Toxicol. 2017, 186, 113–122.

553 554 555 556

(68)

Hued, A. C.; Oberhofer, S.; De Los Ángeles Bistoni, M. Exposure to a commercial glyphosate formulation (Roundup®) alters normal gill and liver histology and affects male sexual activity of Jenynsia multidentata (Anablepidae, cyprinodontiformes). Arch. Environ. Contam. Toxicol. 2012, 62 (1), 107–117.

557 558 559

(69)

Frederick, P.; Jayasena, N. Altered pairing behaviour and reproductive success in white ibises exposed to environmentally relevant concentrations of methylmercury. Proc. R. Soc. B 2011, 278 (1713), 1851–1857.

560 561 562 563

(70)

Crawley, S. E.; Gordon, J. R.; Kowles, K. A.; Potter, M. F.; Haynes, K. F. Impact of sublethal exposure to a pyrethroid-neonicotinoid insecticide on mating, fecundity and development in the bed bug Cimex lectularius L. (Hemiptera: Cimicidae). PLoS One 2017, 12 (5), 1–14.

564 565 566

(71)

Wei, H.; Huang, Y.; Du, J. Sex pheromones and reproductive behavior of Spodoptera litura (Fabricius) moths reared from larvae treated with four insecticides. J. Chem. Ecol. 2004, 30 (7), 1457–1466.

567 568

(72)

Cedergreen, N.; Ritz, C.; Streibig, J. C. Improved empirical models describing hormesis. Environ. Toxicol. Chem. 2005, 24 (12), 3166–3172.

569 570 571 572

(73)

Duncan, J.; Hinchcliffe, A.; Palmquist, K. Evidence of potential long term effects in (aquatic and terrestrial) invertebrates after short term pulsed exposure, Literature reviews on ecotoxicology of chemicals with special focus on plant protection products. 2009. Reference: FT/EFSA/PPR/2008/01. http://www.efsa.europa.eu/en/supporting/pub/17e.htm

573 574 575

(74)

Schmidt, T. S.; Kraus, J. M.; Walters, D. M.; Wanty, R. B. Emergence flux declines disproportionately to larval density along a stream metals gradient. Environ. Sci. Technol. 2013, 47 (15), 8784–8792.

576

(75)

Roff, D. A. Life history evolution; Sinauer Associates: Sunderland, MA, 2002.

577 578 579

(76)

Rohr, J. R.; Palmer, B. D. Aquatic herbicide exposure increases salamander desiccation risk eight months later in a terrestrial environment. Environ. Toxicol. Chem. 2005, 24 (5), 1253–1258.

580 581 582

(77)

Goedkoop, W.; Spann, N.; Åkerblom, N. Sublethal and sex-specific cypermethrin effects in toxicity tests with the midge Chironomus riparius Meigen. Ecotoxicology 2010, 19 (7), 1201–1208.

583 584 585

(78)

Janssens, L.; Stoks, R. Fitness effects of chlorpyrifos in the damselfly Enallagma cyathigerum strongly depend upon temperature and food level and can bridge metamorphosis. PLoS One 2013, 8 (6), e68107.

586 587 588

(79)

Debecker, S.; Sommaruga, R.; Maes, T.; Stoks, R. Larval UV exposure impairs adult immune function through a trade-off with larval investment in cuticular melanin. Funct. Ecol. 2015, 29 (10), 1292–1299. 24 ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

Environmental Science & Technology

589

25 ACS Paragon Plus Environment

Environmental Science & Technology

590

Figure legends

591

Figure 1. Candidate path models used to test whether the effect of a predictor on lifetime mating

592

success (LMS) of the damselfly Coenagrion puella is (a) direct, (b) indirect, or (c) both direct

593

and indirect. The predictor in the models was selected based on linear mixed effect analyses of

594

treatment means (see main text for details). Double-headed arrows are partial correlations

595

between traits.

596

Figure 2. Larval survival (a), growth rate (b), and post-pesticide exposure survival (c) of the

597

damselfly Coenagrion puella as a function of pesticide treatment and hatching period. Growth

598

rate is controlled for initial mass. Note that larval survival during exposure (a) was calculated

599

based on larvae that survived the winter (79%, see Results). Given are least square means ± 1

600

SE.

601

Figure 3. Metamorphic traits (age at emergence (a) and mass at emergence (b)) for the females

602

of the damselfly Coenagrion puella as a function of pesticide treatment and hatching period.

603

Both traits are controlled for larval density in the microcosms. Given are least square means ± 1

604

SE.

605

Figure 4. Lifetime mating success (LMS) as a function of pesticide treatment (a,b) and path

606

diagrams depicting the effects of the pesticide treatment on life-history traits and LMS (c,d) for

607

the females of the damselfly Coenagrion puella. Separately shown are plots and path diagrams

608

for the early hatching period (a,c) and late hatching period (b,d). Given are least square means ±

609

1 SE for the upper plots. Pesticide treatment was significant in the early period (*), but not in the

610

late period (NS). Dashed grey lines indicate non-significant paths, whereas double-headed

611

arrows are partial correlations between traits. Binary coding for pesticide treatment in the path

26 ACS Paragon Plus Environment

Page 26 of 31

Page 27 of 31

Environmental Science & Technology

612

model: control = 0, pesticide = 1. Standardized path coefficients are given next to the arrows (see

613

Table S6 for details).

27 ACS Paragon Plus Environment

Environmental Science & Technology

614 615

Figure 1.

28 ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

Environmental Science & Technology

616 617

Figure 2.

29 ACS Paragon Plus Environment

Environmental Science & Technology

618 619

Figure 3

30 ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31

Environmental Science & Technology

620

621 622

Figure 4

31 ACS Paragon Plus Environment