Subscriber access provided by UNIV OF DURHAM
Ecotoxicology and Human Environmental Health
Molecular effects, speciation and competition of inorganic and methyl mercury in the aquatic plant Elodea nuttallii Rebecca Beauvais-Flück, Vera I Slaveykova, Ulf Skyllberg, and Claudia Cosio Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b02124 • Publication Date (Web): 09 Jul 2018 Downloaded from http://pubs.acs.org on July 12, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
Environmental Science & Technology
1
Molecular effects, speciation and competition of inorganic and methyl mercury in the
2
aquatic plant Elodea nuttallii
3
Rébecca Beauvais-Flück1, Vera I. Slaveykova1, Ulf Skyllberg2, Claudia Cosio*,1,†
4 5
1
6
Sciences, Faculty of Sciences, University of Geneva, Switzerland.
7
2
8
Sciences, Umeå, Sweden.
Department F.-A. Forel for environmental and aquatic sciences, Earth and Environmental
Department of Forest Ecology and Management, Swedish University of Agricultural
9 10
†
11
UMR-I 02 (SEBIO), Université de Reims Champagne Ardenne, F-51687 Reims, France.
present address: Unité Stress Environnementaux et BIOSurveillance des Milieux Aquatiques
12 13
*Corresponding author:
[email protected],
[email protected] 1 ACS Paragon Plus Environment
Environmental Science & Technology
14
Abstract
15
Mercury (Hg) remains hazardous in aquatic environments, because of its biomagnification in
16
food webs. Nonetheless, Hg uptake and impact in primary producers is still poorly
17
understood. Here, we compared the cellular toxicity of inorganic and methyl Hg (IHg; MeHg)
18
in the aquatic plant Elodea nuttallii. IHg and MeHg regulated contigs involved in similar
19
categories (e.g. energy metabolism, development, transport, secondary metabolism), but
20
MeHg regulated more contigs, supporting a higher molecular impact than IHg. At the
21
organism level, MeHg induced antioxidants, while IHg decreased chlorophyll content.
22
The uptake of Hg and expression of a subset of contigs was subsequently studied in complex
23
media. Measured uptake pointed to a contrasted impact of cell walls and copper (Cu) on IHg
24
and MeHg. Using speciation modeling, differences in uptake were attributed to the differences
25
in affinities of IHg and MeHg to organic matter in relation to Cu speciation. We also
26
identified a distinct gene expression signature for IHg, MeHg and Cu, further supporting
27
different molecular toxicity of these trace elements.
28
Our data provided fundamental knowledge on IHg and MeHg uptake in a key aquatic primary
29
producer and confirmed the potential of transcriptomic to assess Hg exposure in
30
environmentally realistic systems.
31 32
Keywords: copper, dissolved organic matter (DOM), pigment, redox, RNA-sequencing,
33
speciation.
2 ACS Paragon Plus Environment
Page 2 of 33
Page 3 of 33
Environmental Science & Technology
34
Introduction
35
Mercury (Hg) is one of the most hazardous metals, mostly because of the biomagnification of
36
its methylated form (MeHg) in aquatic food webs [1]. Freshwater ecosystems remain
37
currently affected by Hg pollution worldwide, despite past and ongoing efforts to reduce
38
anthropogenic emissions of Hg [2]. Aquatic plants have a key ecological role in shallow water
39
ecosystems and are instrumental for Hg transfer to top consumers in food webs [3-6]. For
40
example, the macrophytes Elodea nuttallii and Typha domingensis showed high Hg
41
accumulation from both sediments and water column [7-9]. Hg accumulated in plants can
42
affect their fitness and life cycle. In macrophytes, frequently reported toxic effects of Hg
43
include the generation of oxidative stress, such as increased lipid peroxidation and altered
44
activity of enzymes involved in the regulation of reactive oxygen species (ROS) [10-13]. At
45
the organism level, inorganic Hg (IHg) inhibited root growth in rice seedlings, [14]. Similarly,
46
exposure 24h to 70 ng·L-1 IHg reduced E. nuttallii root growth, while 23 ng·L-1 MeHg
47
increased root growth [4]. Histology and proteomics further revealed an increased
48
lignification of cell walls in response to IHg exposure [4]. Few studies analyzed the impact of
49
Hg at the transcriptome level in macrophytes. In rice seedlings, IHg affected genes involved
50
in amino acid metabolism, cell-wall formation, chemical detoxification, secondary
51
metabolism, signal transduction, abiotic stress response and water transport [14]. In E.
52
nuttallii, a strong impact of IHg on genes involved in energy metabolism and metal
53
homeostasis was reported [15]. Similarly, transcriptomics of E. nuttallii exposed to MeHg
54
revealed a significant impact on the expression level of genes involved in energy metabolism
55
and water transport [16]. However, most studies were conducted in simplified media, lacking
56
organic matter, which limits their environmental relevance. Moreover, to our knowledge no
57
previous study aimed at comparing uptake and effects of IHg and MeHg with comparable
58
experimental conditions on the transcriptome of aquatic plants was performed. Detailed
3 ACS Paragon Plus Environment
Environmental Science & Technology
59
characterization of molecular and cellular mechanisms of toxicity and uptake as well as
60
tolerance mechanisms in aquatic plant in realistic conditions is a prerequisite to understand the
61
fate and the impact of IHg and MeHg in the aquatic environment.
62
This study aimed first to compare the molecular toxicity pathways triggered by 2h-long
63
exposures to a large range of IHg and MeHg concentrations (0.01 to 15 µg·L-1) in the
64
macrophyte Elodea nuttallii using whole transcriptome analysis (RNA-Sequencing; RNA-
65
Seq) and physiological endpoints including pigments content (chlorophyll and anthocyanin)
66
and the activity of anti-oxidative enzymes. Subsequently, we investigated the Hg uptake and
67
expression level of a subset of genes in E. nuttallii exposed to mixtures, containing IHg,
68
MeHg, dissolved organic matter (DOM) and copper (Cu) to gain fundamental knowledge on
69
Hg uptake and tolerance mechanisms in a key aquatic primary producer in more realistic
70
exposure scenarios.
71 72
Material and Methods
73
Plant exposure conditions
74
All experiments were done in triplicates. Elodea nuttallii were grown in the laboratory [7].
75
For all treatments, shoots (~10 cm long without roots) were exposed 2h in 500 mL of an
76
artificial medium (91.0 mg·L-1 CaCl2.2H2O, 43.2 mg·L-1 MgSO4.7H2O, 23.5 mg·L-1
77
NaHCO3, 13.6 mg·L-1 KH2PO4 and 0.4 mg·L-1 NH4NO3, pH 6.9 ± 0.1) in controlled
78
conditions (16:8h light:dark, 1000 lux, 20°C). This short time of exposure was chosen to
79
target early-response genes to IHg and MeHg, before more general stress response appears.
4 ACS Paragon Plus Environment
Page 4 of 33
Page 5 of 33
Environmental Science & Technology
80
For the first set of experiments (whole transcriptome analysis): treatments included 0.015,
81
0.15 and 15 µg·L-1 IHg (Hg(NO3)2) and 0.01, 0.10, 1 or 10 µg·L-1 MeHg (MeHgCl). For
82
comparison, in Europe, the current environmental quality standard (EQS) for THg in surface
83
water is 0.07 µg·L-1.
84
For the second set of experiments (targeted transcriptome analysis): To improve
85
environmental realism, treatments included 0.15 (low) and 15 µg·L-1 IHg (high), 0.25 µg·L-1
86
MeHg, 4 µg·L-1 Cu (CuSO4), low IHg + MeHg, high IHg + MeHg, low IHg + 4 µg·L-1 Cu,
87
70 µg·L-1 IHg (high) + 4 µg·L-1 Cu and MeHg + 17 µg·L-1 Cu. All these conditions were in
88
addition tested in presence of 1 or 10 mg·L-1 Suwanee River humic acid (SRHA). Thereafter
89
treatments without IHg, MeHg, Cu, SRHA are named Control, while treatments with SRHA
90
alone (in the absence of metal) are named 1 mg·L-1 SRHA and 10 mg·L-1 SRHA according to
91
their respective concentration.
92
Total RNA extraction
93
After 2h of exposure, 3 shoots per treatment were snap-frozen in liquid nitrogen, ground and
94
total RNA was
95
(http://www.thermofisher.com/order/catalog/product/AM9738). RNA quality was determined
96
by Qubit and bioanalyzer.
97
Whole transcriptome analysis (RNA-Seq)
98
See details in Supplementary text. Briefly, libraries of cDNA were prepared and sequenced on
99
Illumina HiSeq 2000. Differential gene expression (DGE) was computed with the EdgeR
100
Bioconductor package. We defined differently expressed contigs as significant at a threshold
101
of false discovery rate corrected p-values (FDR) 99 % of IHg and MeHg to be complexed by SRHA RS
335
functional groups, thus an extremely low abundance of IHg and MeHg complexed by the
336
inorganic ligands Cl- and OH- (Table 3 and S7). At 1 mg·L-1 of SRHA, its RS functional
337
groups were abundant enough to complex >99% of IHg and MeHg at their lowest addition
338
levels, but when IHg was increased, SRHA RS groups were more than saturated by Hg2+ and
339
inorganic complexes involving Cl- and OH- dominated the speciation of both IHg and MeHg
340
(Table 3). Notably, Cu(II) could not compete with MeHg for DOM associated RS groups.
341
However, if Cu(II) was significantly transformed to Cu(I) in our experiments, Cu(I)
342
outcompeted the bonding of MeHg to RS groups, with MeHgOH and MeHgCl dominating the
343
speciation of MeHg at 1 mg·L-1 SRHA (Table 3, SI text).
344
We made a calculation assuming all Cu added was reduced to Cu(I) and that the bonding to
345
SRHA-RS groups was similar in strength to RS encountered in aquatic environments [55]. A
346
transformation of all Cu(II) to Cu(I) is obviously an exaggeration, but formation of Cu(I) in
347
our experiment is not unlikely, given the experimental conditions and uptake mechanisms of
348
Cu as Cu(I), which is conserved through Eukaryotes including terrestrial and aquatic plants
349
[7, 56, 57]. In principal our calculation shows that Cu(I) would compete efficiently with
350
MeHg in presence of 1 mg·L-1 SRHA (Table 3) and in essence occupy all RS groups. In
15 ACS Paragon Plus Environment
Environmental Science & Technology
351
presence of 10 mg·L-1 SRHA, the RS groups in the SRHA should be abundant enough to bind
352
in principal all Cu(I) and MeHg, but because of differences in the relative binding strengths
353
between Cu(I) and MeHg, the concentration of Cu complexes would be substantially smaller
354
than the concentration of inorganic MeHg complexes. This suggests that the competition for
355
plant uptake gets in favor of MeHg in presence of 10 mg·L-1 SRHA. The fact that the uptake
356
of MeHg indeed follows the pattern predicted by the chemical speciation model may be taken
357
as support for Cu2+ being reduced to Cu+ in presence of SRHA under our experimental
358
conditions. Because IHg forms relatively stronger complexes than MeHg with DOM-RS
359
groups [52], Cu(I) have less of a competitive effect towards Hg in the IHg + Cu treatments.
360
Concomitantly, in absence of SRHA, Cu + IHg/MeHg significantly decreased Cu uptake as
361
compared to Cu addition alone. Further, the Cu uptake significantly increased when 1 mg·L-1
362
SRHA was added in mixtures, but not in Cu treatment only (t-test, p-value = 0.06; Figure S2).
363
With 10 mg·L-1 SRHA, [Cu]intra was not significantly different vs without SRHA, except for
364
Cu + MeHg, for which Cu uptake significantly decreased (1.4×), confirming the interaction
365
between Cu and IHg/MeHg, and previously observed competition of Cu on [THg]intra uptake
366
as well as increased [THg]intra in MeHg + Cu mixtures.
367
These observations of different effects of Cu on IHg and MeHg uptake, together with
368
differences of cell wall binding mentioned above are congruent with the 2-4× higher
369
bioaccumulation of MeHg than IHg of field studies [3] and highlight that binding to DOM
370
and competition with other metals are central factors for IHg and MeHg bioavailability in
371
natural waters. To better understand and predict Hg uptake and trophic transfer in aquatic
372
systems, it is desirable to in-depth study the effect of parameters controlling Hg uptake in
373
primary producers at the base of food webs. To reach this aim, new tools to measure
374
bioavailability are highly desirable.
16 ACS Paragon Plus Environment
Page 16 of 33
Page 17 of 33
Environmental Science & Technology
375
The expression level of selected contigs was analyzed by hierarchical clustering
376
encompassing the expression level of all contigs (Figure 3) to assess its potential as
377
ecotoxicogenomic tool of exposure in complex media. Indeed such tools to measure
378
bioavailability of IHg and MeHg in the field would be of great interest. Transcriptomic is seen
379
as sensitive and specific to assess the impact of low contaminant concentrations in non-model
380
organisms [15, 30, 58-61]. Contigs were selected based on their functions in the main
381
impacted biological pathways and dose-dependent response with increasing IHg or MeHg
382
exposure. After quality control and removal of contigs showing a log2FC >0.5 in 1 and 10
383
mg·L-1 SRHA, 16 contigs only passed quality control (Table S2) and included a catalase, the
384
PSI subunit, two pectinesterases, 4-coumarate-ligase and 2 S-adenosylmethionine (SAM)
385
enzymes, involved in the biosynthesis of secondary metabolites.
386
Without SRHA, contigs signature discriminated well the samples exposed to Cu from other
387
treatments, and IHg vs MeHg (Figure 3). In response to IHg exposure, 9 genes were up-
388
regulated (e.g. 4-coumarate ligase, DNA polymerase) at low and high IHg, and 4 were up-
389
regulated only at higher IHg concentration (catalase, SAM carrier 1, pectinesterases). This
390
result confirms the RNA-Seq results (e.g. cell wall biosynthesis). The IHg-MeHg mixtures
391
were clustered closer to MeHg than IHg treatments, suggesting a higher sensitivity of the
392
selected contigs to MeHg than IHg. On the opposite, MeHg + Cu was clustered next to Cu,
393
certainly because of the down-regulation of Locus_8381 and the receptor kinase which
394
showed specificity to Cu exposure.
395
The presence of 1 mg·L-1 SRHA did not influence the clustering of samples (Figure 3),
396
confirming its non-significant effect on the uptake. At 10 mg·L-1 SRHA, the signature of high
397
IHg treatments was significantly distinct from the other treatments (Figure 3) forming 2
398
clusters: IHg and IHg + MeHg samples clustered together, suggesting that the MeHg-
399
sensitivity was mitigated compared without DOM.
17 ACS Paragon Plus Environment
Environmental Science & Technology
400
Without SRHA, Cu samples clustered according to measured Cu uptake, ranging from 2.20 ±
401
0.02 µg·g-1dw in high IHg + Cu to 6.37 ± 0.13 µg·g-1dw in Cu treatment (Figure 3). At 1 mg·L-1
402
SRHA, Cu uptake increased (vs in the absence of SRHA), but the gene signature of high IHg
403
+ Cu was closer to IHg/MeHg treatments and 1 mg·L-1 SRHA (Figure 3). Similarly, at 10
404
mg·L-1 SRHA, although Cu uptake decreased in Cu + MeHg vs Cu, non-significant change in
405
the gene expression signature was observed and Cu treatments were identically clustered than
406
for other SRHA concentrations (Figures 3 and S2).
407
Data suggested that the subset of selected genes responded more specifically to IHg and
408
MeHg than Cu. However, gene signature successfully discriminated IHg, MeHg and Cu
409
treatments in presence of DOM. Similarly, a previous study in E. nuttallii was able to
410
correlate gene signature with gradients of Hg concentrations in natural water both in situ and
411
in the laboratory, nevertheless it was unable to differentiate IHg from Cu, Cd or IHg + Cd [15,
412
30]. Previous studies showed transcriptomic efficient for analysis of short-term exposure, to
413
correlate with gradients of contaminants in natural waters, to be more sensitive than classical
414
bioassays (e.g. bioaccumulation or physiological effects), as well as to be able to identify
415
toxicant-specific signatures [15, 30, 62-64]. These properties seem particularly interesting for
416
in-situ analysis characterized by a cocktail of different metals with DOM [21, 65-69], but a
417
subset of genes targeting more diverse pollutants need to be identified. Here we identified 16
418
contigs that can complete the previous set of contigs targeting metals available in E. nuttallii
419
[7]. Our data confirmed the hypothesis that transcriptomic is a very sensitive approach, and is
420
promising for deriving biomarkers of exposure in environmentally relevant scenarios, for
421
example, in mixtures at sublethal concentrations and should be further explored for its
422
suitability for ecotoxicological testing.
423 424
Acknowledgments
18 ACS Paragon Plus Environment
Page 18 of 33
Page 19 of 33
Environmental Science & Technology
425
The Swiss National Science Foundation (contracts 205321_138254, 200020_157173) for
426
financial support, Drs. Mylène Docquier and Didier Chollet for nCounter. RBF acknowledges
427
the support of funds Constantin Topali.
428 429
Supporting Information
430
Supplementary text of Material and Methods for RNAseq, cellular endpoints, mercury
431
analysis, chemical speciation modelling. Supplementary text of Results and Discussion for
432
chemical speciation modelling. Supplementary tables of differential gene expression in
433
simplified and complex media, effective concentrations in media, stability constants for
434
chemical speciation modelling, uptake and number of reads and contigs, Mapman ontology,
435
speciation modelling. Supplementary figures of Biological pathways, Cu uptake, Relationship
436
between Cu and humic acids.
437 438
References
439
1.
440
of gold mining and damming on mercury contamination in fish. Environ Sci Technol 2005,
441
39, (8), 2448-2454.
442
2.
443
Environmental Transport. In UNEP, Geneva, Switzerland: 2013.
444
3.
445
E.; Ungureanu, V.-G.; Dominik, J., Extremely elevated methyl mercury levels in water,
446
sediment and organisms in a Romanian reservoir affected by release of mercury from a chlor-
447
alkali plant. Water Res. 2014, 49, 391-405
448
4.
449
proteomic changes suggest an important role of cell walls in the high tolerance to metals of
450
Elodea nuttallii. J Hazard Mater 2013, 263, (2), 575-583.
451
5.
452
V. I.; Cosio, C., Role of cellular compartmentalization in the trophic transfer of mercury
453
species in a freshwater plant-crustacean food chain. J. Hazard. Mater. 2016, 320, 401-407.
454
6.
Boudou, A.; Maury-Brachet, R.; Coquery, M.; Durrieu, G.; Cossa, D., Synergic effect
UNEP, Global Mercury Assessment 2013: Sources, Emissions, Releases and
Bravo, A. G.; Cosio, C.; Amouroux, D.; Zopfi, J.; Chevalley, P.-A.; Spangenberg, J.
Larras, F.; Regier, N.; Planchon, S.; Pote, J.; Renaut, J.; Cosio, C., Physiological and
Beauvais-Fluck, R.; Chaumot, A.; Gimbert, F.; Queau, H.; Geffard, O.; Slaveykova,
Beauvais-Fluck, R.; Gimbert, F.; Mehault, O.; Cosio, C., Trophic fate of inorganic and 19 ACS Paragon Plus Environment
Environmental Science & Technology
455
methyl-mercury in a macrophyte-chironomid food chain. J Hazard Mater 2017, 338, 140-147.
456
7.
457
Mercury bioaccumulation in the aquatic plant Elodea nuttallii in the field and in microcosm:
458
accumulation in shoots from the water might involve copper transporters. Chemosphere 2013,
459
90, (2), 595-602.
460
8.
461
Cosio, C., Elodea nuttallii roots effect on bacterial communities and MMHg proportion in a
462
Hg polluted sediment. Plos one 2012, 7, (9), e45565, 1-8.
463
9.
464
Doimeadios, R. C. R.; Millan, R., Mercury species accumulation and distribution in Typha
465
domingensis under real field conditions (Almaden, Spain). Environ Sci Poll Res 2018, 1-7.
466
10.
467
M.; Singh, S. P., Mercury bioaccumulation induces oxidative stress and toxicity to submerged
468
macrophyte Potamogeton crispus L. Bull Environ Contam Toxicol 2000, 65, (5), 573-582.
469
11.
470
antioxidant mechanisms in the marine phanerogam Posidonia oceanica. Dis Aquat Org 2002,
471
50, (2), 157-160.
472
12.
473
L. C.; Huang, H. J., Mercury-induced biochemical and proteomic changes in rice roots. Plant
474
Physiol Biochem 2012, 55, 23-32.
475
13.
476
of mercury in aquatic systems. Environ Toxicol Chem 2014, 33, (6), 1225-1237.
477
14.
478
Huang, T. L.; Lin, C. Y.; Huang, H. J., Transcriptome profiling and physiological studies
479
reveal a major role for aromatic amino acids in mercury stress tolerance in rice seedlings. Plos
480
One 2014, 9, (5), e95163, 1-11.
481
15.
482
Elodea nuttallii transcriptome in response to mercury and cadmium pollution: development of
483
sensitive tools for rapid ecotoxicological testing. Environ Sci Technol 2013, 47, (15), 8825-
484
8834.
485
16.
486
environmental and high concentrations of methylmercury on the transcriptome of the
487
macrophyte Elodea nuttallii. Aquat Toxicol 2018, 194, 103-111.
Regier, N.; Larras, F.; Bravo, A. G.; Ungureanu, V. G.; Amouroux, D.; Cosio, C.,
Regier, N.; Frey, B.; Converse, B.; Roden, E.; Grosse-Honebrick, A.; Bravo, A. G.;
Lominchar, M. A.; Sierra, M. J.; Jimenez-Moreno, M.; Guirado, M.; Martin-
Ali, M. B.; Vajpayee, P.; Tripathi, R. D.; Rai, U. N.; Kumar, A.; Singh, N.; Behl, H.
Ferrat, L.; Romeo, M.; Gnassia-Barelli, M.; Pergent-Martini, C., Effects of mercury on
Chen, Y. A.; Chi, W. C.; Huang, T. L.; Lin, C. Y.; Thi, T. Q. N.; Hsiung, Y. C.; Chia,
Cosio, C.; Flück, R.; Regier, N.; Slaveykova, V. I., Effects of macrophytes on the fate
Chen, Y. A.; Chi, W. C.; Trinh, N. N.; Huang, L. Y.; Chen, Y. C.; Cheng, K. T.;
Regier, N.; Baerlocher, L.; Munsterkotter, M.; Farinelli, L.; Cosio, C., Analysis of the
Beauvais-Fluck, R.; Slaveykova, V. I.; Cosio, C., Effects of two-hour exposure to
20 ACS Paragon Plus Environment
Page 20 of 33
Page 21 of 33
Environmental Science & Technology
488
17.
Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton,
489
S.; Cooper, A.; Markowitz, S.; Duran, C.; Thierer, T.; Ashton, B.; Meintjes, P.; Drummond,
490
A., Geneious Basic: an integrated and extendable desktop software platform for the
491
organization and analysis of sequence data. Bioinformatics 2012, 28, (12), 1647-1649.
492
18.
493
Speleman, F., Accurate normalization of real-time quantitative RT-PCR data by geometric
494
averaging of multiple internal control genes. Genome biol 2002, 3, (7), Research0034, 1-12.
495
19.
496
Fernie, A. R.; Stitt, M.; Usadel, B., Mercator: a fast and simple web server for genome scale
497
functional annotation of plant sequence data. Plant Cell Environ 2014, 37, (5), 1250-8125.
498
20.
499
development stages indicates a set of class III peroxidase genes potentially involved in pod
500
shattering in Arabidopsis thaliana. Bmc Genomics 2010, 11, 16, 1-16.
501
21.
502
to mercury exposure under enhanced ultraviolet radiation: Effects on bioaccumulation,
503
transcriptome, pigment content and oxidative stress. Aquat Toxicol 2016, 180, 218-226.
504
22.
505
Protection Agency 2007.
506
23.
507
US Environmental Protection Agency 2002.
508
24.
509
and chemical speciation of mercury and methylmercury driven by organic thiols and
510
micromolar sulfide concentrations in boreal wetland soils. Environ Sci Technol 2017, 51, (7),
511
3678-3686.
512
25.
513
Muller, L. A.; Rhee, S. Y.; Stitt, M., MAPMAN: a user-driven tool to display genomics data
514
sets onto diagrams of metabolic pathways and other biological processes. Plant j cell molec
515
biol 2004, 37, (6), 914-939.
516
26.
517
to using MapMan to visualize and compare Omics data in plants: a case study in the crop
518
species, Maize. Plant Cell Environ 2009, 32, (9), 1211-1229.
519
27.
520
data. Bioinformatics 2002, 18, 207-208.
Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.;
Lohse, M.; Nagel, A.; Herter, T.; May, P.; Schroda, M.; Zrenner, R.; Tohge, T.;
Cosio, C.; Dunand, C., Transcriptome analysis of various flower and silique
Regier, N.; Beauvais-Flück, R.; Slaveykova, V. I.; Cosio, C., Elodea nuttallii exposure
USEPA, Aquatic life ambient freshwater quality criteria- copper. US Environmental
USEPA, Method 1631: Mercury in Water by Oxidation, Purge and Trap, and CVAFS.
Liem-Nguyen, V.; Skyllberg, U.; Bjorn, E., Thermodynamic modeling of the solubility
Thimm, O.; Blasing, O.; Gibon, Y.; Nagel, A.; Meyer, S.; Kruger, P.; Selbig, J.;
Usadel, B.; Poree, F.; Nagel, A.; Lohse, M.; Czedik-Eysenberg, A.; Stitt, M., A guide
Sturn, A.; Quackenbush, J.; Trajanoski, Z., Genesis: cluster analysis of microarray
21 ACS Paragon Plus Environment
Environmental Science & Technology
521
28.
Bravo, A. G.; Le Faucheur, S.; Monperrus, M.; Amouroux, D.; Slaveykova, V. I.,
522
Species-specific isotope tracers to study the accumulation and biotransformation of mixtures
523
of inorganic and methyl mercury by the microalga Chlamydomonas reinhardtii. Environ Poll
524
2014, 192, (0), 212-215.
525
29.
526
distribution in willow. J environ qual 2004, 33, (5), 1779-1785.
527
30.
528
Worms, I. A. M.; Petit, B.; Civic, N.; Docquier, M.; Slaveykova, V. I., Transcriptomic
529
approach for assessment of the impact on microalga and macrophyte of in-situ exposure in
530
river sites contaminated by chlor-alkali plant effluents. Water Res. 2017, 121, 86-94.
531
31.
532
spermidine on polyamine metabolism in water hyacinth leaves under mercury stress. Plant
533
Growth Regul 2010, 60, (1), 61-67.
534
32.
535
the antioxidant behavior of spermine in cadmium- and copper-treated wheat leaves. Biometals
536
2007, 20, (2), 185-195.
537
33.
538
molecules in plant responses and adaptation to heavy metal stress. J Exp Bot 2006, 57, (4),
539
711-726.
540
34.
541
M. T.; Frey, B.; Pauly, M.; Bacic, A.; Ringli, C., Arabidopsis leucine-rich repeat extensin
542
(LRX) proteins modify cell wall composition and influence plant growth. BMC Plant Biol
543
2015, 15, 1-11.
544
35.
545
Martinoia, E.; Lee, Y., The phytochelatin transporters AtABCC1 and AtABCC2 mediate
546
tolerance to cadmium and mercury. Plant J 2012, 69, (2), 278-288.
547
36.
548
Metal toxicity in Chlamydomonas reinhardtii. Effect on sulfate and nitrate assimilation.
549
Biomol Engin 2003, 20, (4-6), 199-203.
550
37.
551
during and after mercury exposure. J Hazard Mater 2012, 217, 271-278.
552
38.
553
Anthocyanins, thiols, and antioxidant scavenging enzymes are involved in Lemna gibba
Wang, Y.; Greger, M., Clonal differences in mercury tolerance, accumulation, and
Dranguet, P.; Cosio, C.; Le Faucheur, S.; Beauvais-Fluck, R.; Freiburghaus, A.;
Ding, C. X.; Shi, G. X.; Xu, X. Y.; Yang, H. Y.; Xu, Y., Effect of exogenous
Groppa, M. D.; Tomaro, M. L.; Benavides, M. P., Polyamines and heavy metal stress:
Sharma, S. S.; Dietz, K.-J., The significance of amino acids and amino acid-derived
Draeger, C.; Fabrice, T. N.; Gineau, E.; Mouille, G.; Kuhn, B. M.; Moller, I.; Abdou,
Park, J.; Song, W. Y.; Ko, D.; Eom, Y.; Hansen, T. H.; Schiller, M.; Lee, T. G.;
Mosulen, S.; Dominguez, M. J.; Vigara, J.; Vilchez, C.; Guiraum, A.; Vega, J. M.,
Wu, Y.; Wang, W. X., Thiol compounds induction kinetics in marine phytoplankton
Leao, G. A.; de Oliveira, J. A.; Felipe, R. T. A.; Farnese, F. S.; Gusman, G. S.,
22 ACS Paragon Plus Environment
Page 22 of 33
Page 23 of 33
Environmental Science & Technology
554
tolerance to arsenic. J Plant Inter 2014, 9, (1), 143-151.
555
39.
556
on amino acid and anthocyanin content of potato tubers. BMC Plant Biol 2008, 8, 1-10.
557
40.
558
of heavy metals on plant growths and pigment contents in Arabidopsis thaliana. Plant Pathol
559
J 2012, 28, (4), 446-452.
560
41.
561
kinetically controls mercury bioavailability to Bacteria. Environ Sci Technol 2014, 48, (6),
562
3153-3161.
563
42.
564
matter: the role of the mercury-to-DOM concentration ratio. Environ Sci Technol 2002, 36,
565
(16), 3564-3570.
566
43.
567
review. Chemosphere 2004, 55, (3), 319-331.
568
44.
569
Meech, J. A., Induced plant uptake and transport of mercury in the presence of sulphur-
570
containing ligands and humic acid. New Phytol 2005, 166, (2), 445-454.
571
45.
572
Plant tolerance to mercury in a contaminated soil is enhanced by the combined effects of
573
humic matter addition and inoculation with arbuscular mycorrhizal fungi. Environ sci poll res
574
intern 2016, 23, (11), 11312-11322.
575
46.
576
bioavailability of copper for lettuce sprouts. Environ Intern 2005, 31, (4), 603-608.
577
47.
578
matter: importance of the copper-to-dissolved organic matter ratio and implications for the
579
biotic ligand model. Environ Sci Technol 2012, 46, (18), 9948-9955.
580
48.
581
reduced organic sulfur in humic acid as affected by S/Hg ratio. Environ Sci Technol 2001, 35,
582
(13), 2741-2745.
583
49.
584
mercury(II) in soil organic matter: EXAFS evidence for linear two-coordination with reduced
585
sulfur groups. Environ Sci Technol 2006, 40, (13), 4174-4180.
586
50.
Dancs, G.; Kondrak, M.; Banfalvi, Z., The effects of enhanced methionine synthesis
Baek, S. A.; Han, T.; Ahn, S. K.; Kang, H.; Cho, M. R.; Lee, S. C.; Im, K. H., Effects
Chiasson-Gould, S. A.; Blais, J. M.; Poulain, A. J., Dissolved organic matter
Haitzer, M.; Aiken, G. R.; Ryan, J. N., Binding of mercury(II) to dissolved organic
Ravichandran, M., Interactions between mercury and dissolved organic matter - a
Moreno, F. N.; Anderson, C. W.; Stewart, R. B.; Robinson, B. H.; Ghomshei, M.;
Cozzolino, V.; De Martino, A.; Nebbioso, A.; Di Meo, V.; Salluzzo, A.; Piccolo, A.,
Inaba, S.; Takenaka, C., Effects of dissolved organic matter on toxicity and
Craven, A. M.; Aiken, G. R.; Ryan, J. N., Copper(II) binding by dissolved organic
Hesterberg, D.; Chou, J. W.; Hutchison, K. J.; Sayers, D. E., Bonding of Hg(II) to
Skyllberg, U.; Bloom, P. R.; Qian, J.; Lin, C. M.; Bleam, W. F., Complexation of
Drexel, R. T.; Haitzer, M.; Ryan, J. N.; Aiken, G. R.; Nagy, K. L., Mercury(II) 23 ACS Paragon Plus Environment
Environmental Science & Technology
587
sorption to two Florida Everglades peats: evidence for strong and weak binding and
588
competition by dissolved organic matter released from the peat. Environ Sci Technol 2002,
589
36, (19), 4058-4064.
590
51.
591
methyl mercury to reduced sulfur groups in soil and stream organic matter as determined by
592
x-ray absorption spectroscopy and binding affinity studies. Geochim Cosmochim Acta 2002,
593
66, (22), 3873-3885.
594
52.
595
Hg and MeHg in wetland soils and sediments under suboxic conditions: Illumination of
596
controversies and implications for MeHg net production. J Geophys Res: Biogeosci 2008,
597
113, G00C03 1-14.
598
53.
599
Modeling of the structure-specific kinetics of abiotic, dark reduction of Hg(II) complexed by
600
O/N and S functional groups in humic acids while accounting for time-dependent structural
601
rearrangement. Geochim Cosmochim Acta 2015, 154, 151-167.
602
54.
603
complexation between mercury and dissolved organic matter in a contaminated environment.
604
Environ Sci Technol 2009, 43, (22), 8548-8553.
605
55.
606
copper, complexing ligands and thiol compounds in shallow lagoon waters. Mar Environ Res
607
2009, 67, (1), 17-24.
608
56.
609
tightrope. Biochim Biophys Acta Molec Cell Res 2006, 1763, (7), 737-746.
610
57.
611
transporter family in Arabidopsis thaliana. Plant Mol Biol 2003, 51, (4), 577-587.
612
58.
613
P.; Budzinski, H.; Durrieu, G.; Sarraco, J.; Elie, P.; Couture, P.; Baudrimont, M.; Bernatchez,
614
L., Transcriptome profile analysis reveals specific signatures of pollutants in Atlantic eels.
615
Ecotoxicology 2015, 24, (1), 71-84.
616
59.
617
expression profiling in Daphnia magna, part II: Validation of a copper specific gene
618
expression signature with effluent from two copper mines in California. Environ Sci Technol
619
2008, 42, (16), 6257-6263.
Qian, J.; Skyllberg, U.; Frech, W.; Bleam, W. F.; Bloom, P. R.; Petit, P. E., Bonding of
Skyllberg, U., Competition among thiols and inorganic sulfides and polysulfides for
Jiang, T.; Skyllberg, U.; Wei, S.; Wang, D.; Lu, S.; Jiang, Z.; Flanagan, D. C.,
Miller, C. L.; Southworth, G.; Brooks, S.; Liang, L.; Gu, B., Kinetic controls on the
Chapman, C. S.; Capodaglio, G.; Turetta, C.; van den Berg, C. M., Benthic fluxes of
Balamurugan, K.; Schaffner, W., Copper homeostasis in eukaryotes: Teetering on a
Sancenón, V.; Puig, S.; Mira, H.; Thiele, D.; Peñarrubia, L., Identification of a copper
Baillon, L.; Pierron, F.; Coudret, R.; Normendeau, E.; Caron, A.; Peluhet, L.; Labadie,
Poynton, H. C.; Zuzow, R.; Loguinov, A. V.; Perkins, E. J.; Vulpe, C. D., Gene
24 ACS Paragon Plus Environment
Page 24 of 33
Page 25 of 33
Environmental Science & Technology
620
60.
Simon, D. F.; Davis, T. A.; Tercier-Waeber, M. L.; England, R.; Wilkinson, K. J., In
621
situ evaluation of cadmium biomarkers in green algae. Environ. Pollut. 2011, 159, (10), 2630-
622
2636.
623
61.
624
profiling of Chlamydomonas reinhardtii exposed to trace levels of free cadmium. Environ
625
Toxicol Chem 2008, 27, (8), 1668-1675.
626
62.
627
Pankratz, M.; Jakel, J.; Strahle, U., Transcriptional profiling reveals barcode-like
628
toxicogenomic responses in the zebrafish embryo. Genome Biol 2007, 8, (10), R227, 1-17.
629
63.
630
B.; Varshavsky, J.; Loguinov, A. V.; Vulpe, C. D.; Perkins, E. J., Biomarker discovery and
631
transcriptomic responses in Daphnia magna exposed to munitions constituents. Environ Sci
632
Technol 2009, 43, (11), 4188-4193.
633
64.
634
transcriptomic response of Arabidopsis thaliana to polymetallic contamination: implications
635
for the identification of potential biomarkers of metal exposure. Metallomics 2016, 8, (5),
636
518-531.
637
65.
638
Design and calibration of microarrays as universal transcriptomic environmental biosensors.
639
Comp Funct Genomics 2005, 6, (3), 132-137.
640
66.
641
of a pesticide/heavy metal mixture in marine bivalves: a transcriptomic assessment. BMC
642
Genomics 2011, 12, 195, 1-17.
643
67.
644
signatures in Chlamydomonas reinhardtii as Cd biomarkers in metal mixtures. Aquat Toxicol
645
2010, 100, (1), 120-127.
646
68.
647
Patarnello, T.; Bargelloni, L., Transcriptomic resources for environmental risk assessment: a
648
case study in the Venice lagoon. Environ Pollut 2015, 197, 90-98.
649
69.
650
Makynen, E. A.; Durhan, E. J.; Kahl, M. D.; Thomas, L. M.; Perkins, E. J.; Ankley, G. T.,
651
Ecotoxicogenomics to support ecological risk assessment: a case study with bisphenol A in
652
fish. Environ Sci Technol 2012, 46, (1), 51-59.
Simon, D. F.; Descombes, P.; Zerges, W.; Wilkinson, K. J., Global expression
Yang, L.; Kemadjou, J. R.; Zinsmeister, C.; Bauer, M.; Legradi, J.; Muller, F.;
Garcia-Reyero, N.; Poynton, H. C.; Kennedy, A. J.; Guan, X.; Escalon, B. L.; Chang,
Gomez-Sagasti, M. T.; Barrutia, O.; Ribas, G.; Garbisu, C.; Becerril, J. M., Early
Almeida, J. S.; McKillen, D. J.; Chen, Y. A.; Gross, P. S.; Chapman, R. W.; Warr, G.,
Dondero, F.; Banni, M.; Negri, A.; Boatti, L.; Dagnino, A.; Viarengo, A., Interactions
Hutchins, C. M.; Simon, D. F.; Zerges, W.; Wilkinson, K. J., Transcriptomic
Milan, M.; Pauletto, M.; Boffo, L.; Carrer, C.; Sorrentino, F.; Ferrari, G.; Pavan, L.;
Villeneuve, D. L.; Garcia-Reyero, N.; Escalon, B. L.; Jensen, K. M.; Cavallin, J. E.;
25 ACS Paragon Plus Environment
Environmental Science & Technology
653 654
Table 1. Proportion (%) of regulated contigs in main MapMan categories (FDR < 0.1 %)
655
after 2h exposure of E. nuttallii to a range of IHg and MeHg concentrations. IHg (µg·L-1) Category 0.015 0.15 15 amino acid metabolism 0.5 2.6 3.1 cell structure 6.5 8.5 10.3 energy metabolism 5.7 10.3 11.6 gene regulation 21.2 25.6 22.9 other 3.5 8.6 8.4 redox 0.9 1.2 1.4 secondary metabolism 1.0 4.2 2.7 stress 2.0 3.7 3.0 transport 2.3 8.7 8.2 unknown 56.3 26.4 28.1
MeHg (µg·L-1) 0.01 0.1 1 10 2.1 1.9 1.3 1.4 8.4 8.6 7.1 7.4 7.1 6.8 5.0 5.4 28.0 30.0 29.0 28.4 6.7 6.9 6.6 7.7 1.4 1.2 1.1 1.2 2.3 2.3 1.4 1.7 3.0 3. 5 2.5 2.1 4.6 5.4 3.5 3.8 35.9 33.1 42.4 40.9
656 657
26 ACS Paragon Plus Environment
Page 26 of 33
Page 27 of 33
Environmental Science & Technology
658 659
Table 2. Chlorophyll content, anthocyanin content and activities of anti-oxidative stress
660
enzymes (POD, class III peroxidases; SOD, superoxide dismutases) in E. nuttallii exposed 2h
661
to IHg or MeHg (results are presented as % of control; bold characters indicate significant
662
difference with control, mean ± SD, n = 3; t-test p-value < 0.05). Treatment (µg·L-1) IHg
MeHg
control 0.015 0.15 15 control 0.01 0.1 1 10
chlorophyll (%) 100.0 ± 15.2 77.3 ± 16.2 70.5 ± 9.6 66.7 ± 9.9 100.0 ± 18.7 82.4 ± 35.3 91.8 ± 21.3 96.0 ± 5.7 113.0 ± 11.7
anthocyanin (%)
100.0 ± 10.3 764.2 ± 103.2 502.5 ± 227.9 315.2 ± 89.0 483.7 ± 108.1
POD activity (%)
SOD activity (%)
100 ± 67.9 156.2 ± 54.4 264.0 ± 8.5 20.8 ± 10.7 100.0 ± 20.62 280.0 ± 132.5 328.0 ± 69.5 240.9 ± 56.4 151.6 ± 33.1
100 ± 8.5 117.1 ± 24.8 167.9 ± 27.4 84.26 ± 12.6 100 ± 8.9 86.3 ± 9.0 87.0 ± 5.9 66.8 ± 22.4 94.0 ± 5.9
663 664
27 ACS Paragon Plus Environment
Environmental Science & Technology
665 666 667
668
Table 3. Fraction (M·M-1) of inorganic chemical species of IHg, MeHg and Cu in different treatments.
Treatment Low IHg Low IHg + MeHg Low IHg + Cu
no SRHA IHg: 1.0 IHg: 1.0 MeHg: 1.0 IHg: 1.0 Cu(II): 1.0
1 mg·L-1 SRHA IHg: 4 × 10-15 IHg: 4 × 10-15 MeHg: 4 × 10-6 IHg: 4 × 10-15 Cu(II): 0.73
10 mg·L-1 SRHA IHg: 4 × 10-17 IHg: 4 × 10-17 MeHg: 4 × 10-7 IHg: 4 × 10-17 Cu(II): 0.03
High IHg High IHg + MeHg High IHg + Cu
IHg: 1.0 IHg: 1.0 MeHg: 1.0 IHg: 1.0 Cu(II): 1.0
IHg: 0.83 IHg: 0.83 MeHg: 0.99 IHg: 0.83 Cu(II): 0.93
IHg: 2.3 × 10-16 IHg: 2 × 10-16 MeHg: 1 × 10-6 IHg: 2 × 10-15 Cu(II): 0.46
MeHg: 5 × 10-6 MeHg: 9 × 10-6 Cu(II): 0.89 MeHg: 1.0 Cu(I): 0.77
MeHg: 5 × 10-7 MeHg: 4 × 10-7 Cu(II): 0.27 MeHg: 2 × 10-8 Cu(I): 9 × 10-12
MeHg MeHg: 1.0 MeHg MeHg: 1.0 + Cu Cu(II): 1.0 *MeHg MeHg: 1.0 + Cu(I) Cu(I): 1.0 * Calculated for Cutot = Cu(I)
669 670
28 ACS Paragon Plus Environment
Page 28 of 33
Page 29 of 33
671
Environmental Science & Technology
Figure legends
672 673 674 675 676
Figure 1. Uptake of IHg (white) and MeHg (black) in washed E. nuttallii measured as total Hg (THg = IHg + MeHg, [THg]intra) vs effective concentration of IHg and MeHg in media ([THg]med) (n=3) (A); Number of contigs significantly up- (triangle) or down-regulated (down-pointing triangle; FDR < 0.1 %) vs bioaccumulated ([THg]intra; dw, dry weight) (B).
677 678 679 680 681 682 683
Figure 2. Effect of Cu or MeHg on IHg uptake and of Cu on MeHg uptake in absence (A) and presence of SRHA (B). Uptake was measured as THg (= IHg + MeHg) concentration in washed plants ([THg]intra) and divided by effective THg concentration in medium ([THg]med) at the beginning of the exposure. Asterisks indicate a significant difference with the respective treatment without competitor (A) or with the respective treatment without SRHA normalized to 100% (dashed line) (mean ± SD, n = 3, t-test, p-value