Subscriber access provided by HACETTEPE UNIVERSITESI KUTUPHANESI
Article
Metamorphosis affects metal concentrations and isotopic signatures in a mayfly (Baetis tricaudatus): Implications for the aquatic-terrestrial transfer of metals Jeff S. Wesner, David M Walters, Travis S Schmidt, Johanna M. Kraus, Craig A. Stricker, William H. Clements, and Ruth E. Wolf Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b05471 • Publication Date (Web): 12 Jan 2017 Downloaded from http://pubs.acs.org on January 19, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
Environmental Science & Technology
1
Metamorphosis affects metal concentrations and isotopic signatures in a mayfly (Baetis
2
tricaudatus): Implications for the aquatic-terrestrial transfer of metals
3 4
Jeff S. Wesner*,1, David M. Walters2, Travis S. Schmidt3, Johanna M. Kraus2, Craig A. Stricker2,
5
William H. Clements,4 Ruth E. Wolf5
6 7
1
Department of Biology, University of South Dakota, Vermillion, SD 57069, USA
8
2
U. S. Geological Survey, Fort Collins Science Center, Fort Collins, CO 80526, USA
9
3
U. S. Geological Survey, Colorado Water Science Center, Denver, CO 80225, USA
10
4
Department of Fish, Wildlife & Conservation Biology and Graduate Degree Program in
11
Ecology, Colorado State University, Fort Collins, CO 80523, USA
12
5
13
*Corresponding Author: 605-677-6178, Email:
[email protected] 14
word count: ~6875
15
Abstract
Perkin Elmer, Inc. San Jose, CA 95134
16
Insect metamorphosis often results in substantial chemical changes that can alter
17
contaminant concentrations and fractionate isotopes. We exposed larval mayflies (Baetis
18
tricaudatus) and their food (periphyton) to an aqueous zinc gradient (3-340 µg Zn/l) and
19
measured zinc concentrations at different stages of metamorphosis: larval, subimago, and imago.
20
We also measured changes in stable isotopes (δ15N and δ13C) in unexposed mayflies. Larval zinc
21
concentrations were positively related to aqueous zinc, increasing 9-fold across the exposure
22
gradient. Adult zinc concentrations were also positively related to aqueous zinc, but were 7-fold
23
lower than larvae. This relationship varied according to adult substage and sex. Tissue
1 ACS Paragon Plus Environment
Environmental Science & Technology
24
concentrations in female imagoes were not related to exposure concentrations, but the converse
25
was true for all other stage-by-sex combinations. Metamorphosis also increased δ15N by ~0.8‰,
26
but not δ13C. Thus, the main effects of metamorphosis on insect chemistry were large declines in
27
zinc concentrations coupled with increased δ15N signatures. For zinc, this change was largely
28
consistent across the aqueous exposure gradient. However, differences among sexes and stages
29
suggest that caution is warranted when using nitrogen isotopes or metal concentrations measured
30
in one insect stage (e.g. larvae) to assess risk to wildlife that feed on subsequent life stages (e.g.
31
adults).
32 33 34
Introduction Aquatic insects are nearly ubiquitous in freshwater ecosystems, where they form key
35
links between aquatic and terrestrial food webs when emerging adults enter terrestrial food webs
36
via consumer1-3 or detrital pathways4. In addition to transferring nutrients and energy to
37
terrestrial ecosystems3, 5, adult aquatic insects may also transfer aquatic-derived contaminants6-11.
38
For example, tetragnathid spiders that fed on emerging chironomids from a lake contaminated
39
with polychlorinated biphenyls (PCBs) had concentrations that exceeded wildlife values for
40
arachnivorous birds12. While there is clear potential for adult aquatic insects to transfer
41
contaminants from water to land, not all are transferred equally. For example, contaminants that
42
biomagnify in food chains (e.g. organochlorines and methyl mercury) tend to be retained during
43
metamorphosis, while those that do not biomagnify (e.g. metals) are typically lost during
44
metamorphosis13-15. Maternal transfer of metals to eggs represents another elimination
45
mechanism. Kim et al.13 found that gravid mayflies (Neocloeon triangulifer) had ~27% higher
46
zinc body burdens than postpartum females, indicating substantial transfer of zinc to eggs.
2 ACS Paragon Plus Environment
Page 2 of 33
Page 3 of 33
47
Environmental Science & Technology
There is evidence that the loss of contaminants is dependent on exposure concentrations.
48
In a meta-analysis, Kraus et al.14 found that organisms exposed to high concentrations of
49
contaminants lost proportionally more contaminants during metamorphosis than insects exposed
50
to lower concentrations. However, this result was based on categorical comparisons between
51
“high” versus “low” concentrations, rather than a formal test across an experimental contaminant
52
gradient. Such tests are needed to understand how metal loss in insects varies across
53
environmentally relevant exposure conditions. In addition to contaminant loss, stables isotopes
54
can fractionate during metamorphosis14, 16-18. In particular, δ15N typically increases in adults
55
relative to larvae due to protein metabolism during metamorphosis and the subsequent loss of
56
metabolic waste that is depleted in 15N14, 18. This finding is not universal, however, and δ15N can
57
also be reduced in adults relative to larval insects in some cases16. δ13C, which is used as a diet
58
tracer, can also fractionate during metamorphosis, potentially affecting the estimation of carbon
59
sources in food webs16.
60
Understanding how metamorphosis alters food web tracers like stable isotopes and the
61
relationship between contaminant loss and exposure concentration is critical for predicting the
62
flux of contaminants from water to land and for assessing risks to wildlife that feed on adult
63
aquatic insects14, 19, 20. For example, if the magnitude of contaminant loss is consistent across
64
exposure concentrations (assuming an element is not limiting), then it should be relatively
65
straightforward to predict adult concentrations and the risk to predators of adults if larval
66
concentrations are known, and vice versa. Likewise, isotope fractionation during metamorphosis
67
can alter source information, predictions of trophic position, and estimates of bioaccumulation
68
for organisms that eat larval (e.g. fish) versus adult (e.g. birds) insects14, 16, 18. For example, since
69
δ15N is typically increased in adults versus larvae, the trophic level of a hypothetical organism
3 ACS Paragon Plus Environment
Environmental Science & Technology
70
that eats only adult aquatic insects would be underestimated if it was determined from larval
71
values of δ15N14. Such underestimation could in turn alter estimates of biomagnification based on
72
trophic level14.
73
To experimentally test for concentration-dependent loss of zinc during metamorphosis,
74
we exposed the mayfly Baetis tricaudatus to a gradient of aqueous zinc in the lab. Mayflies are
75
hemimetabolous insects that do not have a pupal form, but rather transform directly from larvae
76
to adult (i.e. incomplete metamorphosis). The first adult life stage, the subimago, is immature,
77
and must complete an additional molt to enter the sexually mature, or imago, adult stage. We
78
measured zinc concentrations in insects at each of these life-stages (larva, subimago, imago) to
79
test the hypothesis that metal loss during metamorphosis changes as a linear function of exposure
80
concentrations. Specifically, we created an exposure gradient of dissolved aqueous zinc that
81
varied over two orders of magnitude at environmentally relevant concentrations. We then tested
82
1) whether tissue concentrations in larvae, subimagos, and imagos were positively related to
83
exposure concentrations, 2) whether the magnitude of zinc lost during metamorphosis varied
84
across the exposure gradient, and 3) whether these relationships differed for male and female
85
insects. We also measured stable nitrogen and carbon isotopes in insects to test the hypothesis
86
that fractionation occurs between larvae, subimagoes, and imagoes.
87 88
Methods
89
Experimental design
90
Three days prior to the start of the experiment, we collected ~800 late-instar B.
91
tricaudatus larvae and ~70 periphyton-covered rocks from Spring Creek in Fort Collins, CO,
92
USA. We identified late instars as mayflies that were relatively large and similarly sized, but
4 ACS Paragon Plus Environment
Page 4 of 33
Page 5 of 33
Environmental Science & Technology
93
without wing pads (mayflies in the final instar have visible wing pads). Aqueous zinc
94
concentration at Spring Creek at the time of collection was 6.2 µg/L and hardness was 147 mg/L.
95
Spring creek is a cold water urban stream with cobble substrates. Larvae and rocks were
96
transported to holding tanks at the Aquatic Experimental Laboratory (AXL) at the Fort Collins
97
Science Center (U.S. Geological Survey (USGS), Fort Collins, CO). Holding tanks consisted of
98
20 1L plastic cups filled with dechlorinated tap water and suspended in a water bath within a
99
Living Stream™ recirculating tank (model LS-900) equipped with a water chiller (model D1-33)
100
supplied by Frigid Units, Inc. (Toledo, OH). Water in the living stream was maintained at a
101
temperature of ~15ºC with a 16:8 light:dark cycle. Mayflies began emerging from the holding
102
tanks after 1-2 days, and we used the emerging individuals and larvae collected on the same days
103
from the holding tanks for stable isotope analyses. Sample sizes for δ15N were 10 larvae, 6
104
subimagoes, and 8 imagoes. Samples sizes for δ13C were 11 larvae, 8 subimagoes, and 9
105
imagoes. The remaining larvae were used in the zinc experiment described below.
106
After three days of acclimation, 40 larvae and 5 rocks were placed in each of 12 replicate
107
flow-through microcosms (Fig. S1). A microcosm consisted of a 15 L bucket submerged in a
108
water bath in two Living Streams (6 buckets per stream) to which we added 5 L of soft water
109
reconstituted from reverse-osmosis filtered water. Flow-through conditions were created by a
110
mesh-covered PCV pipe, which allowed surface water to flow out of the bucket. Clean water
111
(reconstituted as above) was continuously gravity-fed to each bucket from a common head tank
112
above each living stream, creating a flow-through rate of 5 L/day. Stream temperature was held
113
at a constant 16.2 +/- 0.2℃ (mean +/- sd) within the microcosms. Flow was maintained via
114
magnetic pumps attached to the bottom of the bucket with PVC. The pumps only circulated
5 ACS Paragon Plus Environment
Environmental Science & Technology
115
water within each bucket and did not introduce water from the water bath. We fit mesh netting
116
above each bucket to capture emerging insects.
117 118
Water quality
119
On the same day that larvae were added (day 0) we created a zinc gradient by adding
120
dissolved ZnSO4 at nominal concentrations ranging from 0 to 567 µg Zn/L via constant drip
121
from twelve 20L Mariotte bottles positioned 1.5 m above each microcosm (Fig. S1). This
122
gradient was intended to span the range of Zn concentrations typical of nearby streams, while
123
also spanning below and above the U.S. EPA aquatic life criterion for zinc (120 µg Zn/L).
124
Microcosms were randomly assigned a nominal zinc concentration within each living stream,
125
which were treated as a block. Because zinc concentrations were not replicated, we alternated the
126
assignment of each successive concentration among blocks to ensure that blocks contained a
127
similar range of zinc concentrations. We did not include block in the final analyses because
128
boxplots revealed no difference in any response variable between blocks.
129
We sampled aqueous zinc concentrations on days 1, 6, and 18. For each sample, we
130
filtered 50 mL of water through a 0.45 µm Acrodisc syringe filter with Supor membrane into a
131
50mL Falcon™ tube. Water samples were acidified with ultrapure HNO3 and analyzed as
132
described below. We used the mean of these three measurements (day 1, 6, and 18) as the basis
133
for the data analysis below. Temperature and conductivity were measured on days 2 and 6; pH
134
was measured on day 6 (HQ40, Hach Company, Loveland, CO).
135 136
Insect and algal sampling
6 ACS Paragon Plus Environment
Page 6 of 33
Page 7 of 33
137
Environmental Science & Technology
Each morning during the experiment, all nets were inspected for recently emerged adults,
138
noting their developmental stage (subimago or imago) and sex at the time of collection. Adults
139
collected on day one were excluded from analysis to ensure that all insects had at least 24 hours
140
of exposure (Fig. S2). We collected between 1 and 4 larvae from each microcosm on days 4 and
141
12. The experiment ended on day 19, following three days with no emergence. On this day, we
142
collected algal samples by scraping algae into a centrifuge tube from at least two rocks per
143
microcosm. All samples were placed in individually labeled, factory-clean centrifuge tubes (2.5
144
mL) and stored frozen at -20°C. All samples were later dried at 60°C for > 48 hours and weighed
145
to the nearest 0.01 mg.
146 147
Stable isotope and metals analysis
148
Samples for stable isotopes were oven dried, ground to a fine powder, and massed into 4
149
x 6 mm tin capsules. Insects were pooled when necessary to achieve a target mass of ~1 mg dry
150
mass, and were analyzed for δ15N and δ13C by continuous flow isotope ratio mass spectrometry
151
using a Carlo Erba NC2500 interfaced to a Micromass Optima mass spectrometer21. Isotopic data
152
were normalized to air and V-PDB with USGS 40 (δ15N = -4.52 ‰, δ13C = -26.24 ‰) and USGS
153
41 (δ15N = 47.57 ‰, δ13C = 37.76 ‰). Quality control and accuracy were assessed using internal
154
standards and primary standards respectively; precision was ±0.2 ‰ for both isotopes.
155
For metals analysis, insect samples were pooled across collection dates to achieve a target
156
mass of 1 mg dry mass per sample. Composites contained a mix of 3-6 insects collected between
157
days 4-12 (larvae) or days 2-10 (subimagoes and imagoes) (Fig. S2), with the exception of one
158
treatment (567 µg/l nominal zinc) in which larvae were only collected on day 4. Larvae from this
159
treatment were not included in the analysis to ensure that all larval samples had similar lengths of
7 ACS Paragon Plus Environment
Environmental Science & Technology
160
exposure. Prior to analysis we identified two clear outliers. One was a larval composite with
161
24,200 µg Zn/g. The other was a subimago composite with 4,270 µg Zn/g. These values were an
162
order of magnitude higher than any other larval or subimago samples, respectively (Fig. S3). We
163
assumed they were contaminated or mislabeled (e.g. a larval sample mislabeled as an adult
164
sample) and excluded them from the analysis. All other samples were used to estimate the
165
average concentration of zinc in insects or algae for each microcosm.
166
We measured zinc concentrations in larval (n = 11), subimago and imago (n = 12 each)
167
samples at the end of the experiment. These estimates represent means of 2-6 composite samples
168
per microcosm. Because most of our samples contained larvae and adults pooled across different
169
dates, we were unable to reliably track concentration changes over time. However, for emerging
170
mayflies, the mean exposure time varied from 2 to 8 days among samples (mean 3.7 days). That
171
is, a composite sample with two imagoes collected on days two and four would represent three
172
days of exposure on average, and a sample of two imagoes collected on days seven and nine
173
would represent a mean of eight days of exposure. However, we found no evidence of an
174
interaction between mean length of exposure and aqueous zinc concentrations in explaining adult
175
tissue concentrations (exposure length x aqueous Zn interaction, mean and (95% credible
176
intervals): 0.0003 (-0.00088, 0.00095, Table S1), indicating that pooling across time likely had
177
minimal effects.
178
Most of our imago composites (46/55 composites) contained either all females or all
179
males. We used these samples to test whether the relationship between zinc tissue concentrations
180
and aqueous zinc differed between sexes for both subimagoes and imagoes. Not all combinations
181
were present in each microcosm, so the analyses are based on the following sample sizes: male
8 ACS Paragon Plus Environment
Page 8 of 33
Page 9 of 33
Environmental Science & Technology
182
subimagoes (n = 8 microcosms), female subimagoes (n = 7), male imagoes (n = 11), and female
183
imagoes (n = 9). Composited samples were acid digested using a microwave digester (MarsXpress; CEM
184 185
Corporation, Matthews, NC, USA) and a mixture of 4mL HNO3, 2 mL H2O2, and 4 mL H2O.
186
The resultant solutions, along with water samples, were analyzed for total zinc by inductively
187
coupled plasma-mass spectrometry (ICP-MS) at the Crustal Geophysics and Geochemical
188
Science Center (USGS, Denver, CO). Measured zinc concentrations were normalized to the dry
189
mass of insects used in each digestion. All zinc tissue concentrations are reported in units of raw
190
or log-transformed µg Zn/g dry mass. Quality analysis and control samples included blanks and
191
Standard Reference Material (50 mg). Blanks for tissues consisted of nitric acid, DI water, and
192
peroxide run through the digestion procedure exclusive of tissue. In order to compare digestion
193
blank values to tissue concentrations directly the blank values were converted to µg/g using a
194
maximum sample weight of 4 mg. Tissue blanks (n = 8) contained 55 ± 43 µg Zn/g (mean ± sd,
195
range: 8-108 µg/g), representing < 3% of mean insect tissue concentrations (~1900 µg/g).
196
Recovery of zinc in standard reference material for metals in biological tissue (DOLT3 –
197
National Research Council Canada) was 125% of reported values. For water analyses, blanks (DI
198
only) were below detection (99% probability that the slope was >0 (Table S2, Fig. 1a; slope (95%
298
CrI): 0.62 (0.42, 0.8)). There was a positive relationship between predicted values generated
299
from the posterior distribution and observed values (r2 = 0.87, Fig. S5b), indicating good model
300
fit.
301
Larval zinc concentrations averaged 3027 µg Zn/g (95% CrI: 2185-4104 µg Zn/g), and
302
were positively related to aqueous zinc concentrations with a 98% probability that the slope was
303
>0 (Table S2, Fig. 1b). Zinc concentrations in adults (subimagoes + imagoes) declined by >7-
304
fold relative to larvae, averaging 417 (CrI: 332-510 µg Zn/g), but were still positively related to
305
aqueous zinc (Fig. 1b). However, the positive relationship between adult zinc and aqueous zinc
306
differed between subimagoes and imagoes. Subimago zinc concentrations were positively related
307
to aqueous zinc (Fig. 2b), with a slope that was nearly identical to that of larvae [subimago slope:
308
0.38 (0.06, 0.69), larva slope: 0.39 (0.06, 0.72); Table S2; Fig. 3)]. In contrast, imago zinc
309
concentrations were unrelated to aqueous zinc [slope: 0.08 (-0.42, 0.69)] with only a 62%
310
probability that the slope was > 0 (Table S2, Fig. 2c).
311
Changes in the relationship between zinc in insects and aqueous zinc during
312
metamorphosis are reflected in the fold-change of zinc between life-stages. Adults had 2-17 fold
313
lower zinc concentrations than larvae, but this change was not related to aqueous zinc
314
concentrations for either larvae vs adults (Fig. 2d) or larvae vs. subimagoes (Fig. 2e), with only a
315
68 or 65% probability, respectively, that the slope was negative. In contrast, the fold-change
316
between subimagoes and imagoes ranged from -2.3 to 2.9 and was negatively related to aqueous
14 ACS Paragon Plus Environment
Page 14 of 33
Page 15 of 33
Environmental Science & Technology
317
zinc with a probability of 99%, indicating that imagoes exposed to higher aqueous zinc
318
continued to lose zinc relative to subimagoes (Fig. 2f).
319
Male and female subimagoes (Fig. 3a and 3b) and male imagoes (Fig. 3c) were positively
320
related to aqueous zinc, and the model slopes (0.36, 0.45, and 0.42, respectively; Table S2) were
321
similar to those for larvae [0.39 (0.12, 0.70), Table S2]. In contrast, there was no relationship
322
between zinc in female imagoes and aqueous zinc [slope: -0.05 (-0.35, 0.25); Table S2, Fig. 3d],
323
with a probability of 99% that the slopes between male and female imagoes differed (Table S2).
324 325 326
Isotopes Both subimagoes and imagoes were higher in δ15N relative to larvae. The mean increase
327
was 0.8‰ (0.3-1.3) for subimagoes and 0.9‰ (0.4-1.3) for imagoes with a >99% probability that
328
these changes were >0 (Table S2; Fig. 4). In contrast, while there was a 92% probability that
329
δ13C differed between larva and subimago, the mean difference between these stages was small
330
and may not be ecologically important [0.5‰ (-0.2, 1.2); Fig. 4]. There was only a 66%
331
probability that δ13C differed between imagoes and subimagoes, with imagoes declining by
332
0.2‰ (0.9, -0.6) (Fig. 4), which is within analytical uncertainty (0.2 per mil).
333 334 335
Discussion The change in zinc concentrations during metamorphosis in B. tricaudatus, and the
336
subsequent relationship between exposure concentrations and tissue concentrations varied as a
337
function of adult stage (subimago vs imago) and sex (male vs female). Despite a 7-fold reduction
338
in zinc concentrations of adults compared to larvae, there was still a positive relationship to
15 ACS Paragon Plus Environment
Environmental Science & Technology
339
aqueous zinc concentrations. This relationship was retained in male subimagoes, male imagoes,
340
and female subimagoes, but not in female imagoes.
341
In field studies (and most laboratory studies), insects are typically characterized only as
342
larvae or adults, rather than by adult substage or sex. In that sense, our finding that adult mayflies
343
overall retained a positive relationship between zinc in tissues and zinc in the water are
344
consistent with findings from a broad field survey in which zinc concentrations in adult mayflies
345
(but not most other insects) were positively correlated to aqueous zinc concentrations9.
346
Moreover, zinc concentrations in riparian tetragnathid spiders were positively correlated with
347
concentrations in their mayfly prey, particularly Baetis spp., in a field study.9 However, Kraus et
348
al.9 did not measure zinc concentrations in larval insects and did not distinguish between
349
subimagoes and imagoes or by sex. While our study shows that distinguishing subimagoes and
350
imagoes in mayflies can reveal important differences in zinc concentrations, it is limited to
351
mayflies, which are the only aquatic insect that molt as adults, while other aquatic insect orders
352
have only a single adult stage28. Thus, the differences between males and females in this study
353
are likely more generalizable to other taxa than differences in adult substage (mayflies only), and
354
deserve further study. For example, our results, along with those of Kim et al.13, suggest that
355
using zinc concentrations of adult aquatic insects as proxies for measuring aquatic contamination
356
is possible, but only after accounting for zinc loss in females. This is supported by the striking
357
similarity in the mean slopes of aqueous concentrations and tissue concentrations for larvae,
358
male subimagoes, female subimagoes, and male imagoes (0.39, 0.36, 0.45, 0.42, respectively;
359
Table S2), compared to female imagoes (-0.05)
360 361
Regardless of whether adult insect concentrations reflect aqueous concentrations, it is clear that adult insects lose the majority of their zinc during metamorphosis from larva to
16 ACS Paragon Plus Environment
Page 16 of 33
Page 17 of 33
Environmental Science & Technology
362
subimago. However, the amount of zinc lost appears variable among taxa, which may simply be
363
attributable to variation in the initial amount contained in larvae. For example, in a meta-
364
analysis, Kraus et al.14 found that adult insects from multiple taxa (which includes results from B.
365
tricaudatus in this study) had ~2-fold lower average zinc concentrations than larvae,
366
considerably smaller than the 7-fold reduction in our study. However, final average zinc
367
concentrations in adults (~450 µg/g) in this study were similar to averages of the five mayfly
368
species in Kraus et al.14 (~380 µg/g). In contrast, larval B. tricaudatus averaged 3632 µg/g in this
369
study, which is more than 2-fold higher than the five other mayfly species analyzed in Kraus et
370
al. (2014b), only one of which contained larval concentrations >1000 µg/g.
371
The mechanism explaining the reduction in metal concentrations between larvae and
372
adults is unclear. One route of metal loss for insects is through the exuvium, which is shed
373
during molts. However, metal loss in the exuvium is often a small fraction (50% of metal loss. However, while the
378
meconium seems to be a likely pathway of metal loss in our study, we are unaware of studies
379
documenting the presence of meconium in aquatic insects. Nor have we observed obvious signs
380
of it during our collections. Resolving the mechanisms behind metal loss during metamorphosis
381
in aquatic insects remains a critical next step in understanding how metals are removed during
382
this critical life-history event.
383 384
In addition to changes in metal concentration during metamorphosis, we also found differences between male and female imagoes, in which males retained a positive relationship
17 ACS Paragon Plus Environment
Environmental Science & Technology
385
between aqueous zinc and zinc tissue concentrations, but females did not. A likely explanation
386
for this difference is that excess zinc in females was transferred to eggs and lost during
387
oviposition. Kim et al.13 found that eggs contained ~27% of the total zinc concentration in
388
imagoes of the parthenogenetic mayfly, Neocloeon triangulifer. We did not collect eggs nor
389
distinguish between pre-partum and post-partum females in our experiment, however, so this
390
mechanism is speculative.
391
Fractionation of stable nitrogen isotopes occurred entirely during the larva to subimago
392
molt. The ~0.8‰ increase in δ15N is consistent with other studies showing higher δ15N in adults
393
relative to larvae16, 18, 30. It is also ecologically relevant, representing ~1/4 of estimated
394
fractionation during trophic transfer (i.e. trophic enrichment factor = 3.4‰31). The lack of
395
additional fractionation between subimagoes and imagoes indicates that nitrogen stable isotope
396
signatures in B. tricaudatus adults, and perhaps other mayflies, can reliably be used regardless of
397
whether subimagoes or imagoes are consumed. When combined with large changes in metal
398
concentration, it is clear that estimates of trophic magnification and exposure risk for predators
399
of adult aquatic insects need to be estimated directly from adults or corrected for changes in
400
nitrogen isotopes and metal concentrations during metamorphosis from larva to adult14, 16.
401
Failure to do so can result in misleading estimates of trophic position and exposure risk for
402
riparian wildlife that rely on aquatic-terrestrial subsidies14.
403
Survival to emergence was low overall, but unrelated to aqueous zinc. This was
404
somewhat surprising, given that Baetis tricaudatus are often among the most sensitive aquatic
405
insects in mesocosms studies of metal exposure32-35, even though they persist or quickly recover
406
in response to metal contamination in field studies36. In a previous experiment with the baetid
407
mayfly Neocloeon triangulifer, Wesner et al.37 found that survival to emergence was reduced at
18 ACS Paragon Plus Environment
Page 18 of 33
Page 19 of 33
Environmental Science & Technology
408
aqueous exposures below the EPA chronic criterion (85 µg Zn/l – hardness adjusted) due to
409
mortality during metamorphosis. In the current experiment, the hardness-adjusted criterion was
410
56 µg Zn/l. Average zinc concentrations in 10 of our 12 treatments in this experiment exceeded
411
that level, so it is possible mortality during metamorphosis affected insects in most of our
412
treatments, though we did not measure this directly.
413
Metamorphosis is a critical period in insect life-cycles that alters survival and insect
414
chemistry14, 37, 38, with strong implications for linked freshwater-terrestrial food webs39. Recent
415
work indicates that metamorphosis can alter survival estimates37, contaminant concentrations14,
416
and isotopic signatures14, 16, 18 of organisms exposed to contaminants. Our study adds to this
417
knowledge by showing that a) the majority of metal loss and stable isotope change in mayflies
418
occurs during the molt from larva to adult, b) the amount of metal lost in this molt is not
419
concentration-dependent, and c) sex-specific differences in metal concentrations can lead to
420
different relationships between zinc in adult insects and exposure concentrations.
421 422
Acknowledgements
423
We thank John Simon for help in the construction and maintenance of the artificial streams,
424
Robert Zuellig for help in locating and identifying B. tricaudatus, and Lauren Hargis for
425
sampling help. Statistical analyses were improved by travel support for JSW to attend a Bayesian
426
workshop at Colorado State University (Training in Bayesian Modeling for Practicing
427
Ecologists, NSF Award #1145200). Primary funding for the experiments was from the U.S.
428
Geological Survey through a research grant to WHC. This research was subjected to USGS
429
review and approved for publication. Any use of trade, product, or firm names is for descriptive
430
purposes only and does not imply endorsement by the U.S. Government.
19 ACS Paragon Plus Environment
Environmental Science & Technology
431
Supporting Information
432
Description of Bayesian analysis, figures with experimental details and statistical results, and a
433
table with full model results. The data used in analyses are available at this link:
434
https://dx.doi.org/10.5066/F72V2D85. This information is available free of charge via the
435
Internet at http://pubs.acs.org.
436 437
References
438
1.
439
invertebrate prey link streams and riparian zones. Freshw. Biol. 2005, 50, (2), 201-220.
440
2.
441
habitats of a tallgrass prairie stream. Am. Midl. Nat 1993, 129, 288-300.
442
3.
443
terrestrial and aquatic food webs. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, (1), 166-170.
444
4.
445
subsidies: experimental addition of aquatic insects increases terrestrial arthropod densities.
446
Ecology 2011, 92, (11), 2063-2072.
447
5.
448
bird by altering aquatic‐insect subsidies. Ecology 2010, 91, (8), 2406-2415.
449
6.
450
emergent insect-mediated flux of methyl mercury across a gradient of contamination. Environ.
451
Sci. Technol. 2013, 47, (3), 1614-1619.
452
7.
453
from an aquatic insect to terrestrial insect predators. PloS ONE 2013, 8, (6), e67817.
Baxter, C. V.; Fausch, K. D.; Carl Saunders, W. Tangled webs: reciprocal flows of
Gray, L. J. Response of insectivorous birds to emerging aquatic insects in riparian
Nakano, S.; Murakami, M. Reciprocal subsidies: dynamic interdependence between
Hoekman, D.; Dreyer, J.; Jackson, R. D.; Townsend, P. A.; Gratton, C. Lake to land
Epanchin, P. N.; Knapp, R. A.; Lawler, S. P. Nonnative trout impact an alpine‐nesting
Tweedy, B. N.; Drenner, R. W.; Chumchal, M. M.; Kennedy, J. H. Effects of fish on
Mogren, C. L.; Walton, W. E.; Parker, D. R.; Trumble, J. T. Trophic transfer of arsenic
20 ACS Paragon Plus Environment
Page 20 of 33
Page 21 of 33
Environmental Science & Technology
454
8.
Cristol, D. A.; Brasso, R. L.; Condon, A. M.; Fovargue, R. E.; Friedman, S. L.; Hallinger,
455
K. K.; Monroe, A. P.; White, A. E. The movement of aquatic mercury through terrestrial food
456
webs. Science 2008, 320, (5874), 335-335.
457
9.
458
Cross‐ecosystem impacts of stream pollution reduce resource and contaminant flux to riparian
459
food webs. Ecol. Appl. 2014, 24, (2), 235-243.
460
10.
461
cross‐system subsidy: Chronic stream pollution controls riparian spider populations. Ecology
462
2011, 92, (9), 1711-1716.
463
11.
464
export organic contaminants to riparian predators. Ecol. Appl. 2008, 18, (8), 1835-1841.
465
12.
466
from Contaminated Sediments to Terrestrial Ecosystems and Potential Risks to Arachnivorous
467
Birds†. Environ. Sci. Technol. 2009, 44, (8), 2849-2856.
468
13.
469
dynamics in the mayfly Centroptilum triangulifer. Ecotoxicology 2012, 21, (8), 2288-2296.
470
14.
471
Metamorphosis alters contaminants and chemical tracers in insects: Implications for food webs.
472
Environ. Sci. Technol. 2014, 48, (18), 10957-10965.
473
15.
474
chironomids (Diptera, Chironomidae). Environ. Pollut. 1989, 62, (1), 73-85.
Kraus, J. M.; Schmidt, T. S.; Walters, D. M.; Wanty, R. B.; Zuellig, R. E.; Wolf, R. E.
Paetzold, A.; Smith, M.; Warren, P. H.; Maltby, L. Environmental impact propagated by
Walters, D. M.; Fritz, K. M.; Otter, R. R. The dark side of subsidies: adult stream insects
Walters, D. M.; Mills, M. A.; Fritz, K. M.; Raikow, D. F. Spider-Mediated Flux of PCBs
Kim, K.; Funk, D.; Buchwalter, D. Dietary (periphyton) and aqueous Zn bioaccumulation
Kraus, J. M.; Walters, D. M.; Wesner, J. S.; Stricker, C. A.; Schmidt, T. S.; Zuellig, R. E.
Timmermans, K. R.; Walker, P. A. The fate of trace metals during the metamorphosis of
21 ACS Paragon Plus Environment
Environmental Science & Technology
475
16.
Alp, M.; Peckarsky, B. L.; Bernasconi, S. M.; Robinson, C. T. Shifts in isotopic
476
signatures of animals with complex life-cycles can complicate conclusions on cross-boundary
477
trophic links. Aquat. Sci. 2013, 75, (4), 595-606.
478
17.
479
level and metamorphosis on discrimination of hydrogen isotopes in a plant-herbivore system.
480
PLoS ONE 2012, 7, (3), e32744.
481
18.
482
diet: an ontogenetic shift in δ15N during insect metamorphosis. Funct. Ecol. 2008, 22, (1), 109-
483
113.
484
19.
485
transport to terrestrial arthropod consumers in a multiuse river system. Environ. Pollut. 2016,
486
213, 53-62.
487
20.
488
tree swallows: implications for trace-element exposure after habitat remediation. Arch. Environ.
489
Contam. Toxicol. 2013, 65, (3), 575-587.
490
21.
491
coupled. delta. 13C and. delta. 15N measurements. Anal. Chem. 1992, 64, (3), 288-291.
492
22.
Team, R. C. R: A language and environment for statistical computing. 2013.
493
23.
McElreath, R. Statistical Rethinking: A Bayesian Course with Examples in R and Stan.
494
CRC Press: Boca Raton, U.S.A. 2016; Vol. 122.
495
24.
496
Princeton University Press: Princeton, U.S.A., 2015.
497
25.
Peters, J. M.; Wolf, N.; Stricker, C. A.; Collier, T. R.; del Rio, C. M. Effects of trophic
Tibbets, T. M.; Wheeless, L. A.; Del Rio, C. M. Isotopic enrichment without change in
Alberts, J. M.; Sullivan, S. M. P. Factors influencing aquatic-to-terrestrial contaminant
Beck, M. L.; Hopkins, W. A.; Jackson, B. P. Spatial and temporal variation in the diet of
Fry, B.; Brand, W.; Mersch, F.; Tholke, K.; Garritt, R. Automated analysis system for
Hobbs, N. T.; Hooten, M. B. Bayesian Models: a Statistical Primer for Ecologists.
Team, S. D. Stan: a C++ library for probability and sampling, version 1.3. 2013.
22 ACS Paragon Plus Environment
Page 22 of 33
Page 23 of 33
Environmental Science & Technology
498
26.
Buerkner, P. brms: Bayesian Regression Models using Stan. R package version 0.6. 0
499
2015.
500
27.
501
P. Bayesian Data Analysis in Ecology Using Linear Models with R, BUGS, and Stan. Academic
502
Press: Cambridge, U.S.A. 2015.
503
28.
Lancaster, J.; Downes, B. J. Aquatic Entomology. OUP Oxford: Oxford, U.K. 2013.
504
29.
Yasunobu, A.; Suzuki, K. T. Excretion of cadmium and change in the relative ratio of
505
iso-cadmium-binding proteins during metamorphosis of fleshfly (Sarcophaga peregrina). Comp.
506
Biochem. Physiol. Part C: Comparative Pharmacology 1984, 78, (2), 315-317.
507
30.
508
and gut analysis. Amer. Midl. Nat. 1994, 131, (1), 146-155.
509
31.
510
assumptions. Ecology 2002, 83, (3), 703-718.
511
32.
512
Ginés, M. a. J. s. The use of field and mesocosm experiments to quantify effects of physical and
513
chemical stressors in mining-contaminated streams. Environ. Sci. Technol. 2016, 50, (14), 7825-
514
7833.
515
33.
516
in stream microcosms: understanding differences between single species tests and field
517
responses. Environ. Sci. Technol. 2013, 47, (13), 7506-7513.
518
34.
519
assemblage from a Rocky Mountain stream in experimental microcosms. J. North Am.
520
Benthological Soc. 1994, 13, (4), 511-523.
Korner-Nievergelt, F.; Roth, T.; von Felten, S.; Guélat, J.; Almasi, B.; Korner-Nievergelt,
Mihuc, T.; Toetz, D. Determination of diets of alpine aquatic insects using stable isotopes
Post, D. M. Using stable isotopes to estimate trophic position: models, methods, and
Cadmus, P.; Clements, W. H.; Williamson, J. L.; Ranville, J. F.; Meyer, J. S.; Gutiérrez
Clements, W. H.; Cadmus, P.; Brinkman, S. F. Responses of aquatic insects to Cu and Zn
Kiffney, P. M.; Clements, W. H. Effects of heavy metals on a macroinvertebrate
23 ACS Paragon Plus Environment
Environmental Science & Technology
521
35.
Mebane, C. A.; Schmidt, T. S.; Balistrieri, L. S. Larval aquatic insect responses to
522
cadmium and zinc in experimental streams. Environ. Toxicol. Chem. 2016.
523
36.
524
stream ecosystem. Elementa: Science of the Anthropocene 2015, 3, (1), 000042.
525
37.
526
Metamorphosis enhances the effects of metal exposure on the mayfly, Centroptilum triangulifer.
527
Environ. Sci. Technol. 2014, 48, (17), 10415-10422.
528
38.
529
between larval stress, adult asymmetry and individual quality. Funct. Ecol. 2008, 22, (2), 271-
530
277.
531
39.
532
disproportionately to larval density along a stream metals gradient. Environ. Sci. Technol. 2013,
533
47, (15), 8784-8792.
Mebane, C. A.; Eakins, R. J.; Fraser, B. G.; Adams, W. J. Recovery of a mining-damaged
Wesner, J. S.; Kraus, J. M.; Schmidt, T. S.; Walters, D. M.; Clements, W. H.
Campero, M.; De Block, M.; Ollevier, F.; Stoks, R. Metamorphosis offsets the link
Schmidt, T. S.; Kraus, J. M.; Walters, D. M.; Wanty, R. B. Emergence flux declines
534 535
24 ACS Paragon Plus Environment
Page 24 of 33
Page 25 of 33
Environmental Science & Technology
536
Table 1. Water quality in each microcosm. Mean Zn, Ca, Mg, and Hardenss are averaged over three measurements (Day 1, Day 6, and Day 18), conductivity and temperature (°C) are averaged over two measurements, and pH is a single measurement on day 6. Nominal Mean Zn Zn (µg/l) (µg/l) 0 52 103 155 206 258 309 361 412 464 515 567
3 24 59 78 75 141 100 199 340 248 264 264
Day 1 (µg/l)
Day 6 (µg/l)
Day 18 (µg/l)
3 17 13 38 27 32 29 8 135 38 102 60
2 20 57 89 64 141 117 257 306 279 300 282
3 34 106 108 135 251 155 332 579 426 390 450
pH Cond Temp.
6.7 6.7 6.7 6.65 6.64 6.62 6.59 6.56 6.6 6.69 6.58 6.57
150 160 154 160 153 160 164 147 149 158 160 153
16 16.3 16.2 16.5 16 16.2 16.4 16 15.9 16.5 16.7 16.3
537
25 ACS Paragon Plus Environment
Ca Mg Hardness (mg/l) (mg/l) (mg/L) 7.6 8.1 7.5 7.8 8.3 6.1 7.6 8.8 8.7 8.8 8.9 8.1
5.5 5.4 5.5 5.4 5.3 5.4 5.6 5.6 5.9 5.5 5.5 5.7
41.5 45.9 42.7 44.8 42.6 44.7 44.6 37.4 41.3 42.2 43.5 41.7
Environmental Science & Technology
538
Figure Captions
539
Figure 1. Relationship between aqueous zinc exposure (x-axis) and zinc tissue concentrations (y-
540
axis) in a) algae and b) larval and adult insects. Circles are raw data. Lines show the regression
541
line (straight line) and 95% credible intervals (curved lines). Pr is the probability that the slope is
542
greater than 0. Values for the mean and 95% credible intervals of the slope are as follows: a)
543
0.61 (0.42, 0.8), b) Larvae: 0.39 (0.02, 0.77), Adults: 0.25 (-0.12, 0.59).
544 545
Figure 2. Relationship between aqueous zinc exposure (x-axis) and zinc tissue concentrations
546
(top row) for B. tricaudatus sampled as a) larvae, b) subimagoes, and c) imagoes. The bottom
547
row shows the fold change in concentrations between d) larvae and all adults (subimagoes +
548
imagoes), e) larvae and only subimagoes, and f) subimagoes and imagoes. Circles are raw data.
549
Lines show the regression line (straight line) and 95% credible intervals (curved lines). Pr is the
550
probability that the slope is greater than 0 (for zinc concentration) or less than 0 (for fold
551
change). Values for the mean and 95% credible intervals of the slope are as follows: a) 0.39
552
(0.06, 0.72), b) 0.38 (0.06, 0.69), c) 0.08 (-0.42, 0.69), d) -0.75 (-4.29, 2.8), e) 0.58 (-2.85, 4.22),
553
f) -1.16 (-2.13, -0.21).
554 555
Figure 3. Relationship between aqueous zinc exposure (x-axis) and zinc tissue concentrations (y-
556
axis) for B. tricaudatus males and females sampled as subimagoes or imagoes. Circles are raw
557
data. Lines show the regression line (straight line) and 95% credible intervals (curved lines). Pr
558
is the probability that the slope is greater than 0. Values for the mean and 95% credible intervals
559
of the slope are as follows: a) 0.36 (0, 0.71), b) 0.45 (-0.07, 0.97), c) 0.42 (0.12, 0.7), d) -0.05, (-
560
0.35, 0.25).
26 ACS Paragon Plus Environment
Page 26 of 33
Page 27 of 33
Environmental Science & Technology
561 562
Figure 4. Fractionation of δ15N and δ13C during metamorphosis in the mayfly, Baetis tricaudatus.
563
Open circles are the raw values of individual insects. Black circles and error bars are the
564
posterior mean and 95% credible intervals.
565 566 567
27 ACS Paragon Plus Environment
Environmental Science & Technology
Figure 1 a) Algae, Pr = 1
b)
Zinc concentration (µg Zn/g dry mass)
568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589
Aqueous zinc (µg/L)
28 ACS Paragon Plus Environment
Larvae , Pr = 0.98 Adults, Pr = 0.92
Page 28 of 33
Page 29 of 33
Environmental Science & Technology
Pr = 0.99
b)
Pr = 0.99
c)
Pr = 0.62
d)
Pr = 0.68
e)
Pr = 0.65
f)
Pr = 0.99
Fold-change in tissue zinc concentrations
Zinc concentration (µg Zn/g dry mass)
a)
Aqueous zinc (µg/L) 590 591
Figure 2
29 ACS Paragon Plus Environment
Environmental Science & Technology
Figure 3
a) male subimagoes
Zinc concentration (µg Zn/g dry mass)
592 593 594 595 596 597 598 599 600 601 602 603 604 605
b) female subimagoes
Pr = 0.97
Pr = 0.96
c) male imagoes
d) female imagoes
Pr = 0.36
Pr = 1
Aqueous zinc (µg/L)
30 ACS Paragon Plus Environment
Page 30 of 33
Page 31 of 33
Figure 4 12
-30 -31
11 δ13C
-32 δ15N
606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626
Environmental Science & Technology
10
-33 9
8
-34 -35
0.5 larvae 1 1.5subimagoes 2 2.5 imagoes 3 3.5
0.5 larvae 1 1.5subimagoes 2 2.5 imagoes 3 3.5
627
31 ACS Paragon Plus Environment
Environmental Science & Technology
166x141mm (150 x 150 DPI)
ACS Paragon Plus Environment
Page 32 of 33
Page 33 of 33
Environmental Science & Technology
c)
b)
a)
10 R² = 0.1306
9 8.5 8 7.5
e)
10 R² = 0.8425
fold-change larvae to adultsobserved
ln (µg Zn/g)observed
7
9 8 7 6 5
-12 -17 -11
-9 fold-change larvae to adultspredicted
R² = 0.4286 -2
0
2
larvae
subimagos
-7
13Cobserved
6 5
3
7 ln (µg Zn/g)predicted
R² = 0.0149
-7 -12 -17
larvae
subimagos
-7
imagos
11 10 9 R² = 0.4141
8 10
5 5.5 6 6.5 ln (µg Zn/g)predicted
-9 -8 fold-change larvae to subimagospredicted
12
R² = 0.4017
5
9
-2
i)
7
imagos
-31 -32 -33 -34
R² = 0.0332 -35 -32.2 -32 13Cpredicted
7
-10
female imagos male imagos female subimagos male subimagos
9
fold-change subimago to imagopredicted
-30
-7
7
8
f)
-2
h)
4 3 2 1 0 -1 -2 -3
9
10
3 R² = 0.0276
9
ln (µg Zn/g)observed
fold-change subimago to imagoobserved
g)
7 ln (µg Zn/g)predicted
8 9 ln (µg Zn/g)predicted
adults
R² = 0.8425
5
16
fold-change larvae to subimagosobserved
13 14 15 survivalpredicted
5
j)
R² = 0.868
7
5 12
ln (µg Zn/g)observed
15
9.5
10
15Nobserved
ln (µg Zn/g)observed
survivalobserved
20
d)
larvae
10
25
-31.8
ACS Paragon Plus Environment
10.5 15Npredicted
11