Molecular Effects, Speciation, and Competition of Inorganic and

Sci. Technol. , 2018, 52 (15), pp 8876–8884. DOI: 10.1021/acs.est.8b02124. Publication Date (Web): July 9, 2018. Copyright © 2018 American Chemical...
0 downloads 0 Views 576KB Size
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