Effect-Directed Analysis of Aryl Hydrocarbon ... - ACS Publications

Sep 18, 2016 - ABSTRACT: The construction of the Three Gorges Dam. (TGD) in the Yangtze River raises great concern in ecotoxicological research since ...
0 downloads 0 Views 812KB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Effect-directed analysis of aryl hydrocarbon receptor agonists in sediments from the Three Gorges Reservoir, China Hongxia Xiao, Martin Krauss, Tilman Floehr, Yan Yan, Arnold Bahlmann, Kathrin Eichbaum, Markus Brinkmann, Xiaowei Zhang, Xingzhong Yuan, Werner Brack, and Henner Hollert Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b03231 • Publication Date (Web): 18 Sep 2016 Downloaded from http://pubs.acs.org on September 20, 2016

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 34

Environmental Science & Technology

1

Effect-directed analysis of aryl hydrocarbon receptor agonists in sediments

2

from the Three Gorges Reservoir, China

3

Hongxia Xiao1, Martin Krauss2, Tilman Floehr1, Yan Yan1, Arnold Bahlmann2, Kathrin

4

Eichbaum1, Markus Brinkmann1,3, Xiaowei Zhang4, Xingzhong Yuan5, Werner Brack1,2,

5

Henner Hollert1,4,5,6*

6

1

7

Biology and Biotechnology, RWTH Aachen University, Aachen 52074, Germany

8

2

9

UFZ, Leipzig 04318, Germany

Department of Ecosystem Analysis, Institute for Environmental Research, ABBt – Aachen

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

10

3

11

Canada

12

4

13

Nanjing University, Nanjing 210046, China

14

5

15

China

16

6

17

Shanghai 200092, China

School of Environment and Sustainability, University of Saskatchewan, Saskatoon S7N 5B3

State Key Laboratory of Pollution Control & Resource Reuse, School of the Environment,

College of Resources and Environmental Science, Chongqing University, Chongqing 400030,

Key Laboratory of Yangtze Water Environment, Ministry of Education, Tongji University,

18 19

*Corresponding author:

20

Henner Hollert

21

Worringerweg 1, 52074 Aachen, Germany

22

Phone: +49 (0)241 – 80 / 26669, Fax: +49 (0)241 – 80 / 22182

23

E-mail: [email protected] (H.Hollert);

24 25 26 1 ACS Paragon Plus Environment

Environmental Science & Technology

27

ABSTRACT

28

The construction of the Three Gorges Dam (TGD) in the Yangtze River raises great concern

29

in ecotoxicological research since large amounts of pollutants enter the Three Gorges

30

Reservoir (TGR) water bodies after TGD impoundment. In this work, effect-directed analysis

31

(EDA) combining effect assessment, fractionation procedure, target and non-target analyses,

32

was used to characterize aryl hydrocarbon receptor (AhR) agonists in sediments of the TGR.

33

Priority polycyclic aromatic hydrocarbons (PAHs) containing four to five aromatic rings were

34

found to contribute significantly to the overall observed effects in the area of Chongqing. The

35

relatively high potency fractions in the Kaixian area were characterized by PAHs and

36

methylated derivatives thereof and heterocyclic polycyclic aromatic compounds (PACs) such

37

as dinaphthofurans. Benzothiazole and derivatives were identified as possible AhR agonists in

38

the Kaixian area based on non-target liquid chromatography-high resolution mass

39

spectrometry (LC-HRMS). To our knowledge, this study is the first one applying the EDA

40

approach and identifying potential AhR agonists in TGR.

41

Key words: Effect-directed analysis (EDA); aryl hydrocarbon receptor agonists; Three

42

Gorges Reservoir (TGR); non-target analysis; QSAR

43

2 ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

44

Environmental Science & Technology

TOC/Abstract art:

45 46

Introduction

47

The Three Gorges Reservoir (TGR), created in consequence of the Yangtze River’s

48

impoundment by the Three Gorges Dam (TGD), spreads over a distance of 663 km between

49

the town of Sandouping, Hubei Province, and the Jiangjin district of Chongqing Municipality.

50

Besides the obvious benefits of the TGD, such as hydroelectricity, navigation, and flood

51

control, construction of the TGD also poses great challenges to the unique ecosystem,1

52

particularly

53

industrialization, and intensified shipping activities.2 Furthermore, construction of the dam

54

reduced the river’s flow velocity,3,

55

suspend particles and adhering contaminants.5

when

facing

numerous

4

anthropogenic

impacts,

e.g.,

overpopulation,

and consequently increased the sedimentation rate of

56

Sediments are considered as the final sink of persistent and lipophilic pollutants in the

57

environment.6 They can become a potential source of pollutants through resuspension of

58

particulate matter, e.g., during flood events.7-9 Frequently occurring floods in the water

59

fluctuation zone of TGR may increase the pollutants bioavailability through remobilization

60

and direct exposure to benthic organisms.9, 10 To ensure the environmental and public health,

61

suitable and effective monitoring strategies of sediments in TGR are urgently demanded.

62

Extensive chemical-analytical research has been performed on TGR, indicating that the

63

sediments in TGR are polluted with a mixture of persistent organic pollutants (POPs), such as

64

polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs) and polycyclic aromatic

65

hydrocarbons (PAHs).3,14,15 Most of these compounds are known to induce cytochrome P450

66

1A (CYP1A) by ligand-activation of the aryl hydrocarbon receptor (AhR).11, 12 It has been

67

reported that the binding of xenobiotics to the AhR can trigger a broad spectrum of adverse 3 ACS Paragon Plus Environment

Environmental Science & Technology

68

effects, such as on biochemistry, physiology and reproduction in many organisms.12-14

69

Furthermore, exposure of early life stage of fish to AhR agonists may cause increased

70

mortality in further developmental stages15-18 and adverse outcomes in wild fish

71

populations.19-23. This is of particular concern, as the Yangtze River plays an important role

72

for fishery production in China. 24, 25

73

As recently shown in a review,26 only limited research was done on bioassays to

74

determine AhR-mediated activity in this area. Moreover, chemical analysis of target

75

compounds provides only limited information on adverse biological effects of complex

76

mixture.2,

77

fractionation, chemical analysis and bioassays, has been demonstrated to be a suitable tool for

78

the identification of causative toxicants in complex environmental samples. This approach is

79

directed by the biotests (in vitro or in vivo assays), through the assignment of toxicity to

80

several groups of toxicants by separation steps including extraction, clean up, and

81

fractionation. The aim is to remove compounds without significant contribution to sample

82

toxicity and to identify the predominant toxicants using chemical analytical tools. EDA can be

83

particularly useful for addressing the effects on ecological health and environmental

84

management of contamination source, thus further supporting the prioritization and regulation

85

of environmental contaminants.29-31 Although this approach has been successfully used for

86

toxicant identification in various environmental matrices,28, 32-37 the application of EDA in

87

China is still quite limited.38-40, and few research has been done in the TGR.

27

The concept of effect-directed analysis (EDA)28, featuring a combination of

88

In our previous studies, sampling sites along the TGR were screened according to the triad

89

approach 41 to achieve a comprehensive perspective on ecotoxicological status of this area.2, 42

90

Two sites – close to the cities of Chongqing and Kaixian – were identified as regional “hot-

91

spots” with respect to dioxin-like activity and mutagenicity. The present study aimed to

4 ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34

Environmental Science & Technology

92

characterize and identify individual AhR agonist in the sediments extracts in support of the

93

prioritization and regulation of environmental contaminants present in sediments of the TGR.

94

95

Material and methods

96

Sample Collection and Preparation

97

Sediment samples were collected using a Van-Veen sampler in September 2011. Three

98

samples were collected at the mainstream close to Chongqing – upstream (CNG-U),

99

downstream (CNG-D) and directly at the tributary’s inlet (CNG-T), as well as one sample

100

from Hanfeng Lake (HF-L) in Kaixian. For detailed information see Support Information (SI,

101

Figure S1).

102

All sediments were freeze-dried, sieved (≤ 0.2 cm), and thoroughly homogenized by

103

using pestle and mortar. Thereafter, the applied EDA strategy followed the flowchart as

104

shown in Figure 1. Sediments (20 g each) were extracted per run, with acetone: hexane (1:1;

105

v:v) in a pressurized liquid extractor (PLE) (Speed Extractor E-916, Büchi Labortechnik AG)

106

at 100°C and 120 bar in two cycles (heat up 1 min; hold 10 min; discharge 2 min; flush with

107

solvent 1 min; flush with gas 4 min), as described in detail in a previous study.42 Parent

108

extracts were rotary-evaporated close to dryness and re-dissolved in acetone: hexane (1:1; v:v )

109

to a final concentration of 20 g sediment equivalents (SEQ) per mL solvent. The parent

110

extracts were kept in amber glass at 4°C for fractionation.

111 112

Fractionation of Samples

113

To reduce the interference of non-toxic substances such as minerals, salts and large

114

biogenic organic molecules, parent extracts were solvent exchanged to dichloromethane 5 ACS Paragon Plus Environment

Environmental Science & Technology

115

(DCM), and purified by gel-permeation chromatography (GPC) (Biobeads SX3, Bio Rad) as

116

described elsewhere.43 The purified extracts were evaporated close to dryness, and re-

117

dissolved in hexane: DCM (9:1; v:v) for fractionation.

118

The fractionation step was performed on three coupled normal-phase high

119

performance liquid chromatography (NP-HPLC) columns, including cyanopropyl (CN),

120

nitrophenyl (NO) and porous graphitized carbon (PGC), which has been described in detail in

121

Lübcke-von Varel et. al.44 Polar polycyclic aromatic compounds (PACs) were trapped on and

122

eluted from CN column, while nonpolar compound groups were retained on NO and PGC

123

columns. The fractionation windows and eluting compounds in respective fractions according

124

to Lübcke-von Varel et al.44 are given in SI (Table S1). Extract amounts of 40 — 60 g SEQ

125

were fractioned per run. Each fraction was rotary-evaporated close to dryness and re-

126

dissolved in acetone: hexane (1:1; v:v) to a final concentration of 20 g SEQ per mL solvent.

127

Subsequently, 1 mL sub-fractions were solvent exchanged with dimethyl sulfoxide (DMSO)

128

and tested for their AhR-mediated activities as described below. The remaining sub-fractions

129

were stored in amber glass vials at –20°C for chemical analysis. It should be noted here that

130

the whole sample preparation and fractionation procedure is dedicated towards hydrophobic

131

and polyaromatic compound. Thus, this study might potentially exclude other, more polar

132

AhR agonists present in sediments.

133 134

Ethoxyresorufin-O-deethylase (EROD) induction assay

135

In order to avoid cytotoxic effects in the EROD assay, the neutral red retention assay

136

(NR) was carried out with RTL-W1 cells according to the method described by Babich and

137

Borenfreund.45 The AhR-mediated activity was determined by the EROD induction assay

138

using rainbow trout (Oncorhynchus mykiss) liver cells (RTL-W1) according to previously

6 ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

Environmental Science & Technology

139

described methods by Gustavsson et al.46 with slight modifications.2 Samples were analyzed

140

in three replicates, each with three internal replicates serially diluted with medium in seven

141

1:2 steps. The well characterized substance 2,3,7,8-tetrachlorodibenzo-p-dioxin (2,3,7,8-

142

TCDD) was used in a test concentration range from 3.13 to 100 pM on each plate as a

143

reference standard. EROD activity was determined fluorometrically via a multiwell plate

144

reader (Infinite M200, Tecan Austria GmbH, Grödig). The artificial substrate 7-

145

ethoxyresorufin (7-EXT) was deethylated to resorufin by cytochrome P450 (CYP) enzymes.

146

Fluorescence of resorufin was measured with excitation/emission wavelength of 544/590 nm,

147

whereas the amount of protein was detected by fluorescence with excitation/emission wave-

148

length of 360/460 nm. EROD activity was determined based on the quantity of produced

149

resorufin per total amount protein and reaction time. Dose–response curves for EROD

150

induction were computed by log concentration (agonist) vs. response with variable slope

151

using Prism 6.0 (GraphPad Software Inc., San Diego). In order to provide comparability with

152

other research, bioassay-derived TCDD equivalents (BEQs) were calculated on the basis of a

153

fixed effect level of EC25 of the maximum response caused by the samples to those of the

154

standard compound 2,3,7,8-TCDD, using the formula BEQsample (pgTCDD/gSEQ) =

155

EC25TCDD(pgTCDD/mL) / EC25sample(gSEQ/mL). Moreover, the measured BEQs in

156

fractions were compared to the TCDD equivalents (TEQs), which were obtained by

157

multiplying the concentration of each compound in each sample by its relative potency (REP),

158

and summing up these values. This comparison enables an estimation of the contribution of

159

the analyzed chemical contaminants to the measured AhR-mediated activity.11

160 161

Targeted Analysis of PACs by GC-MS

162

A list of 32 PACs including priority PAHs, alkylated PAHs and heterocycle

163

polyromatic chemicals (oxygen, and sulfur-containing PAHs) was selected for measurement 7 ACS Paragon Plus Environment

Environmental Science & Technology

164

based on their occurrence in sediments (SI, Table S5). The identification and quantification of

165

PACs were performed with a GC (Agilent 6890) equipped with an auto sampler (Agilent

166

7683) and coupled with a mass spectrometer (Agilent 5973) with electron ionization source.

167

Aliquots of 1 µl of sample were injected in splitless mode at an injector temperature of 250°C.

168

After sample injection, the analytes were separated on capillary column (HP5MS,

169

35 m × 0.25 mm, i.d., 0.25 µm, Agilent) using helium as carrier gas at a constant flow rate of

170

1.3 mL/min. The oven temperature started at 60°C and ramped with 30°C/min until 150°C,

171

then with 6°C/min up to 186°C, followed by 4°C/min to 280°C held for 21.5 min. For priority

172

16 EPA-PAHs, quantification was performed in selected ion monitoring (SIM) mode using

173

external standards and corrected by means of the injection standards, i.e., 13C benzo[a]pyrene.

174

The concentration of substituted PAHs as well as heterocyclic aromatic compounds were

175

calculated based on methods described elsewhere.47

176 177

Non-target Analysis in polar fractions by LC-HRMS

178

The liquid chromatography separation was performed with an Agilent 1200 system

179

equipped with a Kinetex Core-Shell C18 column (100 mm × 3.0 mm; 2.6 µm; Phenomenex).

180

A linear gradient elution with water and methanol both containing 0.1 % formic acid at a flow

181

rate of 0.2 mL/min was used. The LC system was connected to an ion trap-Orbitrap hybrid

182

instrument (LTQ Orbitrap XL, Thermo Scientific). Analytes were ionized by electrospray

183

ionization (ESI) and atmospheric pressure chemical ionization (APCI) in separate runs, both

184

in positive and negative ion mode, respectively. Details on instrument settings are given in the

185

supporting information (SI, Table S1). Full scan spectra were recorded at nominal resolving

186

power of 100,000 (referenced to m/z 400) at a mass range of m/z 100 – 1000. High resolution

187

product ion spectra (HRMS/MS) were acquired data-dependent for the two most intense

188

precursor ions of the full scan at a nominal resolving power of 15,000, using an isolation 8 ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

Environmental Science & Technology

189

width of 1.3 m/z, a minimum precursor ion intensity of 50,000 and a dynamic exclusion time

190

of 20 s. Both collision induced dissociation (CID) at 35% and higher-energy collisional

191

dissociation (HCD) at 100% were used for fragmentation.

192 193

Processing of HRMS data for compound identification

194

The MZmine 2.10 software48 was used for peak detection and peak lists (consisting of m/z

195

values, retention times, and signal intensities) obtained from full scan chromatograms of the

196

samples, solvent and processing blanks. The processing steps and settings of MZmine were

197

given in SI (SI, Table S2). The peak lists were further processed using a R script to remove

198

background contamination peaks occurring in solvent and processing blanks (intensity ratio of

199

sample to blank < 10), and those originating from background signals not resembling

200

Lorentzian or Gaussian peak shapes (details are given in Hug et al.49 ).

201

To further narrow down the peak lists to containing only compounds potentially being

202

hydrophobic or polyaromatic AhR agonists, which were targeted by the sample preparation

203

procedure, the following filtering steps were applied:

204

(i)

205

all peaks with retention times > 15 minutes were kept, as the active sediment fractions of interest should comprise compounds exerting certain hydrophobicity;

206

(ii)

207

all peaks were kept, which showed an at least five times higher intensity in the active fractions F13 to F15 than in any of the non-active fractions F16 to F18;

208

(iii)

all peaks were kept, which had a mass defect below an upper boundary defined by

209

the molecular formula CnH (2n-6-0.5n) (n is an even number), as most AhR

210

receptor agonists should be of polyaromatic nature or contain at least a certain

211

number of double bonds. This upper boundary corresponds to 6.5 double-bond

212

equivalents (DBEs) for n = 10 and 7.0 for n = 12, etc. 9 ACS Paragon Plus Environment

Environmental Science & Technology

213

The filtered peak lists of potential AhR agonists of the active fractions were processed

214

using the R “nontarget” package.40, 41 This allowed searching for bounds of isotope peaks (13C,

215

15

216

M+Na+, M+K+, M+NH4+ in positive mode; M-H− [M+formate]− in negative mode)

217

considering single and double charged ions. Peaks were finally grouped into components,

218

encompassing the monoisotopic peak and its associated isotope or adduct peaks representing

219

an individual chemical compound. The package’s settings are given in the SI (Table S3).

N,

34

S,

37

Cl, and

81

Br) with a rule-based algorithm and peaks of relevant adducts (M+H+,

220

Starting from these annotated component lists, a determination of molecular formulas was

221

done based on the raw data file using the QuanBrowser of the Xcalibur software (Thermo

222

Scientific). For obtained molecular formulas, the Chemspider compound database50 was

223

searched to retrieve candidate structures. Plausible candidate compounds were identified

224

based on the number of references in Chemspider as an indicator of human use and

225

commercial importance. For these compounds reference standards were obtained if possible

226

and used for confirmation based on retention times and MS/MS spectra. For compounds

227

without reference compound, a tentative identification was done based on an interpretation of

228

MS/MS spectra in comparison to those of known compounds.

229 230

QSAR modeling using VirtualToxLab

231

The AhR binding affinities of candidate compounds, which were identified from both

232

target and non-target analysis were simulated by VirtualToxLab. The VirtualToxLab is an in

233

silico tool for predicting the toxic potential of chemicals by simulating and quantifying their

234

interactions towards a series of proteins, which are known to trigger biological effects using

235

automated, multi-dimensional QSAR.51

236 10 ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

237

Environmental Science & Technology

Statistical analysis

238

Statistical analyses were performed using SigmaPlot 12.0 (Systat Software Inc). Non-

239

parametric Shapiro-Wilk ANOVA on ranks followed by Holm-Sidak’s post hoc test (p ≤ 0.05)

240

was used to determine significant differences in the EROD induction of fractions compared to

241

the process control.

242

Results and discussion

243

Bioassay-Derived induction equivalent quantities

244

AhR-mediated activity of all sub-fractions, parent- (par) and reconstituted extracts (rec), and

245

arithmetic sum (sum) of all fractions was examined in the EROD assay with RTL-W1 cells

246

(Figure 2). BEQs showed significant fraction-specific differences, but followed the same

247

relative distribution pattern in all four samples. The observation was in accordance with other

248

studies which applied the same automated on-line fractionation method as the present study.47,

249

52, 53

250

naphthalenes (PCNs), coplanar PCBs, polychlorinated dibenzo-p-dioxins (PCDDs) and

251

dibenzofurans (PCDFs), caused no or very low AhR-mediated activity, suggesting low levels

252

of these compounds in TGR sediments. This result was in accordance with previous findings

253

that only low levels of PCBs and PCDD/Fs were detected in water and sediments of the

254

TGR.4,

255

fractions F7 – F10, characterized by PAHs with four to six aromatic rings, as well as fractions

256

F13 – F15, mostly characterized by intermediately polar to polar PACs (SI, Table S1). The

257

highest BEQs were detected in HF-L, ranging from 42 to 177 pg BEQ/SEQ g and 69 to

258

109 pg BEQ/SEQ g in fractions F7-F10 and F13-F15, respectively.

The fractions F3-F5, co-eluting with typical AhR agonists such as polychlorinated

42, 54

Significantly greater EROD inducing potency (p ≤ 0.05) was detected in

11 ACS Paragon Plus Environment

Environmental Science & Technology

259

In order to determine eventual losses of activity during sample processing, all fractions

260

were combined to form reconstituted extracts and were tested with the same procedure.

261

Recovery was estimated on the basis of reconstituted fractions of the parent extracts, and was

262

67% for HF-L, 57% for CNG-T, and 93% for CNG-D. Due to limited sample amount, BEQs

263

of the parent extract CNG-U were not available, and thus no recovery could be determined.

264

Our results indicated good chemical recoveries, in which typical chemically analyzed

265

standard compound concentration recoveries of 60 – 80% were reported.55

266

The total response to exposure with these fractions was site-specific. The sample from

267

the mainstream (CNG-U) exhibited higher levels of contaminants than the samples originating

268

from tributaries (CNG-T) and downstream regions (CNG-D), but with similar fractionation

269

patterns. Thus, higher levels of contamination in the upstream of Yangtze river at Chongqing

270

area may be diluted by sediments of the less contaminated tributary, the Jialing River. This

271

results in lower concentrations of contaminants downstream of Chongqing, which is in

272

agreement with previous findings.2, 56 The artificial Hanfeng Lake at Kaixian exhibited the

273

highest AhR-mediated activity, which could be attributed to a lower dilution of discharged

274

contamination compared to the Chongqing area, and a potential accumulation of

275

contamination in the lake.2

276 277

Target Chemical Analysis

278

The PAH fractions (F7 – F10), which exhibited significant EROD induction activity

279

(Figure 2), were selected for priority PAH analysis. Moreover, the highest potency fractions

280

HF-L-F9 – F10 were scanned for methylated PAHs as well as heterocyclic aromatic

281

compounds (SI, Table S4). The fractions in the Chongqing area exhibited similar PAHs

282

profiles, which were characterized by four– and five– ringed PAHs. In addition, the

12 ACS Paragon Plus Environment

Page 12 of 34

Page 13 of 34

Environmental Science & Technology

283

concentration of CNG-U showed a higher contamination than the tributaries (CNG-T) and

284

downstream (CNG-D), which is in accordance with the bio-analytical results of the present

285

study. The fractions from the Kaixian area showed a contamination with pyrogenic five– and

286

six– ringed PAHs, together with some methylated PAHs, and their oxygen (O-) and sulfur (S-)

287

heterocyclic aromatic compounds.

288

The total priority EPA-PAHs ranged from 9 — 101 ng/g sediment (SI, Table S4).

289

Among the parent PAHs, fluoranthene was measured at the highest concentrations

290

(56.6 ng/g sediment in CNG-U-F7). Levels of PAHs with EROD induction capacity,57

291

including

292

indeno[1,2,3-c,d]pyrene, benzo[b]fluoranthene, chrysene, and benzo[a]anthracene were

293

detected up to 40.2 ng/g sediment. The concentrations of heterocyclic aromatic compounds

294

were detected up to 5 ng/g sediment. The analytical results revealed that relatively low levels

295

of priority PAHs in TGR, which are one to two orders of magnitude lower in comparison to

296

levels detected in the Mulde River58, Danube River52, and Oslo harbour. 59

benzo[k]fluoranthene,

dibenzo[a,h]anthracene,

benzo[a]pyrene,

297 298

Contribution of PACs to the EROD induction potency of fractions

299

To estimate the fraction of activity explained by priority PAHs, we compared TEQs,

300

calculated using the REP values given by Bols, et al.57, to BEQs from EROD assay (SI,

301

Figure S2). The results revealed that the contribution of priority PAHs contributed up to 43%

302

of the AhR-mediated activity of responding fractions. Low contributions of chemically

303

derived TEQs were found in F7 and F8 of CNG-U, where priority PAHs such as fluoanthene

304

and pyrene were detected, which have been reported as non-inducers and weak AhR-agonists

305

in previous studies.57,

306

chrysene,

60

Four– to six– ringed parent PAHs such as benzo[a]anthracene,

benzo[b+k]fluoranthene,

benzo[a]pyrene

were

13 ACS Paragon Plus Environment

commonly

found

in

F9,

Environmental Science & Technology

Page 14 of 34

307

indeno[1,2,3-cd]pyrene, dibenz[a,h]anthracene were found in F10 of CNG-U and HF-L.

308

Especially the four– to six– ringed PAHs, which have been confirmed as AhR agonists,35, 48, 49

309

accounted for 43% of the BEQs in CNG-U-F9, thus likely being the main contributors.

310

Similar results were observed in CNG-D-F9 (28%) and CNG-T-F9 (40%). In the significant

311

high potency of fractions HF-L-F9 and HF-L-F10, the priority PAHs accounted for only about

312

9% and 12% of BEQs. The remainder of the activity was investigated and the potentially

313

active chemicals presented below.

314

In the faction of HF-L-F10, we hypothesize that me-PAHs and heterocycles detected

315

may account for a part of the observed effects in the fraction. Dinaphtho[1,2-b;1’,2’-d]furan

316

and dinaphtho[1,2-b;2’,3’-d]furan were identified as potent inducers of EROD activity in the

317

RTL-W1 cell line58 and rat hepatoma H4IIEGud.Luc 1.1 cell line.61 Based on the REP values

318

from

319

8.7 pg TEQ / g SEQ, corresponding to 7% of the BEQ in HF-L-F10. Methylated

320

benzo[a]pyrenes were detected in the fraction. Due to too many isomers present in the sample,

321

concentrations of each compound were hard to quantify. Methylated PAHs were shown to

322

significantly induce AhR-mediated activity in reporter gene assays,62 and have been reported

323

to be responsible for large portions of AhR activity in sediments.52,

324

benzo[a]anthracenes were reported to be significantly more potent AhR inducers than their

325

parent compounds.63 Brack and Schirmer58 identified methylation of chrysene in the 1-

326

position and of benzo[a]anthracene in the 9-position to enhance the EROD induction potency

327

by one to two orders of magnitude. Thus, we hypothesize that methylated benzo[a]pyrenes

328

play a role of the AhR-mediated activity in the bioactive fraction. Due to a lack of respective

329

REP values, the AhR binding affinities of heterocyclic PACs were simulated using QSAR (SI,

330

Table S5). Most of the PACs were predicted to have a binding affinity to the AhR, thus taken

Brack

and

Schirmer,58

the

TEQ

of

dinaphthofurans

14 ACS Paragon Plus Environment

was

calculated

58

as

Methylated

Page 15 of 34

Environmental Science & Technology

331

as tentative compounds. To confirm the contribution of PACs, the relative potency of

332

heterocycles compounds should be an aim of further research.

333 334

Non-target analysis of polar fractions

335

In addition to the high potency fractions F7 to F10, the polar fractions, mainly eluting

336

with mono-nitro- PAHs (F13) as well as (hydroxyl-) quinones, keto-, dinitro-, hydroxyl-

337

PAHs and N-heterocycles (F14 – F15) showed relatively high AhR-mediated effects. Thus,

338

we also focused on characterizing the AhR agonists in polar fractions (F13 – F15) using LC-

339

HRMS. Because of the relatively low effects in samples CNG-D and CNG-T, we focused on

340

characterizing the AhR agonists in the higher effects of fraction F13 – F15 in the samples of

341

Chongqing upstream (CNG-U) and the Kaixian area (HF-L).

342

After peak detection using MZmine 2.10 and the removal of background peaks,

343

typically several hundred peaks remained in the peak lists for ESI+ and APCI+ mode, except

344

for fraction HF-L-F14, in which about 60 and 40 peaks could be detected (Table S5). In

345

APCI- and particularly ESI- mode, a significantly lower amount (< 90) of peaks remained.

346

The numbers of peaks were further reduced by an additional filtering step targeting likely

347

AhR agonists based on retention times, mass defects and absence in the non-active fraction

348

F16 – F18 (SI, Table S6). Nevertheless, in individual fractions more than 100 candidate peaks

349

remained. Thus, priority for identification was given to the most intense peaks in the active

350

fraction HF-L-F13.

351 352

Identification and selection of candidate AhR agonists

353

Based on accurate mass and isotope patterns, a plausible molecular formula could be

354

determined in ESI+ mode of all intense peaks (>106 a.u. intensity) in fraction HF-L-F13 (SI, 15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 34

355

Table S7). No additional compounds could be detected in APCI+ mode. While benzothiazole,

356

2-mercaptobenzothiazole and 2-(methylthio)benzothiazole were confirmed by reference

357

standards, MS/MS spectra of several other compounds suggests the occurrence of further

358

benzothiazole derivatives as indicated by fragments characteristic for a benzothiazole or a

359

mercaptobenzothiazole moiety (Table S7, Figure S3 – S6).

360

Benzothiazole derivatives are widely used as vulcanization accelerators and

361

antioxidants in rubber64,

65

and have been detected in water and sediments.49,

66, 67

362

compound at m/z 300.9923 and RT 27.4 was tentatively assigned as 2,2’-sulfanediylbis-(1,3-

363

benzothiazole), which could be a by-product of mercaptobenzothiazole use. Two isobaric

364

peaks at m/z 239.0669, RT 16.2 and 16.7, respectively, showed almost identical MS/MS

365

spectra with an intensive fragment at m/z 72.0803, which corresponds to a C4H10N group (SI,

366

Figure S5). This suggested N-t-butyl-2-benzothiazolesulphenamide, which is a vulcanization

367

accelerator as well, and a closely related isomer as further plausible candidates. He et al.

368

identified benzothiazole derivatives as AhR-active compounds in tire extracts by a toxicant-

369

identification-evaluation (TIE) approach.68 The study presented the ability of 2-

370

mercaptobenzothiazole to induce AhR-dependent gene expression, as well as benzothiazole to

371

be a weak AhR agonist in mammalian cell bioassays. In a previous study, both 2-

372

mercaptobenzothiazole and benzothiazole were shown as a relative potent AhR agonist of the

373

human AhR expressed in yeast.69

The

374

Another group of compounds in fraction HF-L-F13 were likely benzylamines, among

375

them the confirmed compound tribenzylamine, and several other compounds showing similar

376

MS/MS spectra with characteristic fragments C14H16N+ and C7H7+ (Table S7, Figure S7) were

377

found. The complete absence of the C6H7N+ ion (m/z 93.0573) suggests that no aniline

378

structure was present in these molecules.

16 ACS Paragon Plus Environment

Page 17 of 34

Environmental Science & Technology

379

Among the tentative identified toxicants, 2-mercaptobenzothiazole, benzothiazole

380

2-(methylthio)benzothiazole and tribenzylamine were predicted to have a binding affinity to

381

the AhR by QSAR (Table 1). Thus, we assume benzothiazole and its derivates to be important

382

pollutants in the Kaixian area.

383 384

Source of contamination and environmental significance

385

Through target chemical analysis, high molecular weight PAHs (four- to six- ringed PAHs)

386

(Figure. 3) contributed significant to the responding fractions, thus are taken as significant

387

AhR-active compounds inducing AhR-mediated activity in the Chongqing section. As

388

Chongqing is the most important industrial center in the TGR region, the origins of identified

389

pollutants are possibly from urban traffic emissions and runoff, coal combustion, as well as

390

intensified shipping activities since the impoundment of the reservoir3 Further research is

391

required to identify the causative toxicants for unexplained effects in this area. In the Kaixian

392

area, only a minor part of the EROD inducing potency of tested sediment extracts could be

393

explained by the priority PAHs. The bioactive fractions of PAHs were characterized by a

394

broad variety of heterocycles and methylated PAHs. Dinaphthofurans were identified as

395

pollutants to contribute part of the EROD induction potency in the bioactive fraction. The

396

methylated PAHs were hypothesized as contributors to the toxicity in the Kaixian area. PACs

397

originate from incomplete combustion and industrial processes or fossil fuels as well as from

398

natural sources like volcanic eruptions.70 Methylated benzo[a]pyrenes have been identified as

399

components from cigarette smoke, urban air particulates, gasoline engine, diesel exhaust, and

400

forest fire smoke, as well as in a variety of coal-derived liquids and tars.71-74 Biomass burning

401

such as wood fuel used for cooking or heating, or the burning of straw residues in the fields

402

has caused serious pollution in the southeast of China,75, 76 and should be considered since the

403

Kaixian area constitutes a rather rural area. Benzothiazole is an aromatic heterocyclic 17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 34

77

404

compound, which is used as parent compounds for the synthesis of larger structures

. The

405

compound 2-mercaptobenzothiazole (MBT) and its derivatives such as N-t-butyl-2-

406

benzothiazolesulfenamide (Figure. 3) are used in industry as accelerators for the vulcanization

407

of rubber.65, 78 It is documented that benzothiazoles are not be tightly bound in the rubber

408

matrix, and thus may leach into waste effluents, which enter sediments via runoff and reach

409

surface water via drainage systems.65 Evidence for large amount of MBT and other

410

benzothiazoles are released from vehicle tire rubber into the environment as a result of the

411

weathering and leaching of tire rubber,65, 79 the biological activity of benzothiazoles are still

412

limited. MBT was reported to show acute toxic to fish,80 micrograms,67 and cells.79 Industrial

413

wastewater and street runoff from impervious urban surface are taken as significant source of

414

benzothiazoles in the environment. As the location of sampling site located in a rather rural

415

area, we assume that benzothiazoles in our study most likely originated from an identified

416

rubber factory nearby. Tribenzylamine are well-known extractants in industrial wastewater.81-

417

83

418

Dibenzylamine has been reported to be detected in environment samples, such as surface

419

water84 and wastewater.85 It was reported as by-product during the rubber vulcanization

420

process.86

Thus, the compound was hypothesized to be associated with wastewater treatment.

421

The present investigation integrated a biological-chemical approach for the

422

characterization of AhR agonists in the sediment of TGR. To the best of our knowledge, this

423

is the first time that AhR agonists were analyzed in TGR in more detail using EDA. Our study

424

strongly supports that focus on prioritized pollutants may result in inadequate assessment of

425

complex

426

concentrations were detected in the present study, the vulnerable TGR ecosystem might still

427

be of concern for its absolute pollution mass.26 Long-term monitoring programs including the

428

causative compounds should be employed parallel to the proceeding economic and

environmental

mixtures.36,51,87

Although

relatively

18 ACS Paragon Plus Environment

low

target

chemical

Page 19 of 34

Environmental Science & Technology

429

demographic development due to the rapid industrialization and increased urbanization in

430

TGR. EDA is a suitable methodology to be included in monitoring programs to avoid

431

monitoring irrelevant compounds.

432

EDA has been proven to be useful for toxicant identification in the present study; however,

433

the approach is still not widely applied in environmental routine monitoring programs. With

434

respect to a tedious evaporation and solvent exchange steps, as well as large number of

435

fractions for bioassays, the relative laborious work largely limits a wider applicability of EDA.

436

Thereafter, workflows that can lead to a rapid assessment of the key toxicants are demanded30,

437

88

438

assessment and chemical identification are currently performed, with the aim to be applied as

439

a routine monitoring program, in order to support the prioritization of environmental

440

contaminants and the regulatory decisions in TGR.

. Correspondent studies in high throughput EDA integrating micro-fractionation, effect

441

442

Acknowledgements

443

This study has been carried out as part of the Yangtze-Hydro project (No. FKZ 02WT)

444

supported by the Germany Federal Ministry of Education and Research (BMBF),

445

SOLUTIONS project supported by the European Union Seventh Framework Programme

446

(FP7-ENV-2013-two-stage Collaborative project) under grant agreement No. 603437, and the

447

EDA-EMERGE project supported by the European Union Seventh Framework Programme

448

(FP7-PEOPLE-2011-ITN) under the grant agreement No. 290100. We also want to express

449

our gratitude to Dr. Niels Bols and Dr. Lucy Lee from the University of Waterloo, Canada,

450

who kindly provided the CYP1A expressing fibroblast-like permanent cell line RTL-W1 from

451

primary hepatocytes of rainbow trout (Oncorhynchus mykiss). We thank Marion Heinrich and

452

Lena Schinkel for conducting the analysis of polycyclic aromatic compounds by GC-MS. In 19 ACS Paragon Plus Environment

Environmental Science & Technology

453

addition, Hongxia Xiao received a personal grant supported by the scholarship program

454

Chinese Scholarship Council. Markus Brinkmann was supported by the German National

455

Academic Foundation (‘Studienstiftung des deutschen Volkes’) and is a Banting Fellow of the

456

Natural Sciences and Engineering Resarch Council of Canada (NSERC).

457 458

Supporting Information Available

459

Detailed information on the sampling information, LC-HRMS analysis, target chemical

460

analysis data, as well as the identification of non-target compounds is available in the

461

Supplementary Information. This information is available free of charge via the Internet at

462

http://pubs.acs.org.

20 ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34

Environmental Science & Technology

463

References

464

1.

Shen, G.; Xie, Z., Three Gorges Project: chance and challenge. Science 2004, 304, (5671), 681-681.

465 466 467 468

2.

Floehr, T.; Scholz-Starke, B.; Xiao, H.; Hercht, H.; Wu, L.; Hou, J.; Schmidt-Posthaus, H.; Segner, H.;

469 470 471

3.

472 473 474

4.

475 476 477

5.

478 479 480

6.

481 482 483

7.

484 485 486 487

8.

488 489 490

9.

491 492 493 494

10.

495 496

11.

Kammann, U.; Yuan, X.; Roß-Nickoll, M.; Schäffer, A.; Hollert, H., Linking Ah receptor mediated effects of sediments and impacts on fish to key pollutants in the Yangtze Three Gorges Reservoir, China — A comprehensive perspective. Science of The Total Environment 2015, 538, 191-211. Chen, J.; Chen, Z.; Xu, K.; Wei, T.; Li, M.; Wang, Z.; Watanabe, M., ADP-flow velocity profile to

interpret hydromorphological features of China’s Yangtze Three-Gorges valley. Chinese Science Bulletin 2005, 50, (7), 679-684. Wang, J.; Bi, Y.; Pfister, G.; Henkelmann, B.; Zhu, K.; Schramm, K.-W., Determination of PAH, PCB,

and OCP in water from the Three Gorges Reservoir accumulated by semipermeable membrane devices (SPMD). Chemosphere 2009, 75, (8), 1119-1127. Hu, B.; Yang, Z.; Wang, H.; Sun, X.; Bi, N.; Li, G., Sedimentation in the Three Gorges Dam and the

future trend of Changjiang (Yangtze River) sediment flux to the sea. Hydrology and Earth System Sciences 2009, 13, (11), 2253-2264. Hollert, H.; Keiter, S.; König, N.; Rudolf, M.; Ulrich, M.; Braunbeck, T., A new sediment contact assay

to assess particle-bound pollutants using zebrafish (Danio rerio) embryos. Journal of Soils and Sediments 2003, 3, (3), 197-207. Wölz, J.; Cofalla, C.; Hudjetz, S.; Roger, S.; Brinkmann, M.; Schmidt, B.; Schäffer, A.; Kammann, U.;

Lennartz, G.; Hecker, M.; Schüttrumpf, H.; Hollert, H., In search for the ecological and toxicological relevance of sediment re-mobilisation and transport during flood events. Journal of Soils and Sediments 2009, 9, (1), 1-5. Gerbersdorf, S.; Jancke, T.; Westrich, B., Sediment Properties for Assessing the Erosion Risk of

Contaminated Riverine Sites. An approach to evaluate sediment properties and their covariance patterns over depth in relation to erosion resistance. First investigations in natural sediments (11 pp). Journal of Soils and Sediments 2007, 7, (1), 25-35. Woelz, J.; Fleig, M.; Schulze, T.; Maletz, S.; Luebcke-von Varel, U.; Reifferscheid, G.; Kuehlers, D.;

Braunbeck, T.; Brack, W.; Hollert, H., Impact of contaminants bound to suspended particulate matter in the context of flood events. Journal of Soils and Sediments 2010, 10, (6), 1174-1185. Brinkmann, M.; Hudjetz, S.; Cofalla, C.; Roger, S.; Kammann, U.; Giesy, J.; Hecker, M.; Wiseman, S.;

Zhang, X.; Wölz, J.; Schüttrumpf, H.; Hollert, H., A combined hydraulic and toxicological approach to assess resuspended sediments during simulated flood events. Part I–multiple biomarkers in rainbow trout. Journal of Soils and Sediments 2010, 10, (7), 1347-1361. Van den Berg, M.; Birnbaum, L. S.; Denison, M.; De Vito, M.; Farland, W.; Feeley, M.; Fiedler, H.;

Hakansson, H.; Hanberg, A.; Haws, L., The 2005 World Health Organization reevaluation of human and

21 ACS Paragon Plus Environment

Environmental Science & Technology

497 498

mammalian toxic equivalency factors for dioxins and dioxin-like compounds. Toxicological sciences 2006, 93,

499

12.

Sorg, O., AhR signalling and dioxin toxicity. Toxicology Letters 2014, 230, (2), 225-233.

500 501 502

13.

Eichbaum, K.; Brinkmann, M.; Buchinger, S.; Reifferscheid, G.; Hecker, M.; Giesy, J. P.; Engwall, M.;

503 504 505 506

14.

507 508 509 510

15.

511

16.

512 513

(lsalvelinus fontinalis) during early development. Environmental toxicology and chemistry 1994, 13, (5), 817-

514 515 516

17.

517 518 519 520

18.

521 522

19.

523 524

20.

525 526

21.

527 528

22.

529 530 531 532

23.

(2), 223-241.

van Bavel, B.; Hollert, H., In vitro bioassays for detecting dioxin-like activity - Application potentials and limits of detection, a review. Science of the Total Environment 2014, 487, 37-48. Schiwy, A.; Brinkmann, M.; Thiem, I.; Guder, G.; Winkens, K.; Eichbaum, K.; Nüßer, L.; Thalmann,

B.; Buchinger, S.; Reifferscheid, G., Determination of the CYP1A-inducing potential of single substances, mixtures and extracts of samples in the micro-EROD assay with H4IIE cells. Nature protocols 2015, 10, (11), 1728-1741. Andreasen, E. A.; Spitsbergen, J. M.; Tanguay, R. L.; Stegeman, J. J.; Heideman, W.; Peterson, R. E.,

Tissue-specific expression of AHR2, ARNT2, and CYP1A in zebrafish embryos and larvae: effects of developmental stage and 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin exposure. Toxicological Sciences 2002, 68, (2), 403-419. Walker, M. K.; Peterson, R. E., Toxicity of 2, 3, 7, 8‐tetrachlorodibenzo‐P‐dioxin to brook trout

820. Di Paolo, C.; Groh, K. J.; Zennegg, M.; Vermeirssen, E. L. M.; Murk, A. J.; Eggen, R. I. L.; Hollert, H.;

Werner, I.; Schirmer, K., Early life exposure to PCB126 results in delayed mortality and growth impairment in the zebrafish larvae. Aquatic Toxicology 2015, 169, 168-178. Brinkmann, M.; Hudjetz, S.; Kammann, U.; Hennig, M.; Kuckelkorn, J.; Chinoraks, M.; Cofalla, C.;

Wiseman, S.; Giesy, J. P.; Schaffer, A.; Hecker, M.; Wolz, J.; Schuttrumpf, H.; Hollert, H., How flood events affect rainbow trout: evidence of a biomarker cascade in rainbow trout after exposure to PAH contaminated sediment suspensions. Aquatic Toxicology 2013, 128-129, 13-24. Gilbertson, M., PCB and dioxin research and implications for fisheries research and resource

management. Fisheries 1992, 17, (4), 26-27. Niimi, A., Biological and Toxicologtcal Effects of Environmental Contaminants in Fish and Their Eggs.

Canadian Journal of Fisheries and Aquatic Sciences 1983, 40, (3), 306-312. Van der Oost, R.; Beyer, J.; Vermeulen, N. P., Fish bioaccumulation and biomarkers in environmental

risk assessment: a review. Environmental toxicology and pharmacology 2003, 13, (2), 57-149. Whyte, J. J.; Jung, R.; Schmitt, C.; Tillitt, D., Ethoxyresorufin-O-deethylase (EROD) activity in fish as

a biomarker of chemical exposure. Critical reviews in toxicology 2000, 30, (4), 347-570. Ankley, G. T.; Bennett, R. S.; Erickson, R. J.; Hoff, D. J.; Hornung, M. W.; Johnson, R. D.; Mount, D.

R.; Nichols, J. W.; Russom, C. L.; Schmieder, P. K., Adverse outcome pathways: a conceptual framework to support ecotoxicology research and risk assessment. Environmental Toxicology and Chemistry 2010, 29, (3), 730-741.

22 ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34

Environmental Science & Technology

533 534

24.

535 536

25.

537 538 539

26.

540 541 542

27.

543 544

28.

545 546 547

29.

548 549 550

30.

551 552 553 554 555 556 557 558

31.

559 560 561 562

32.

563 564 565

33.

566 567 568

34.

Liu, J. K.; Cao, W. X., Fish resources of the Yangtze River basin and the tactics for their conservation.

Resources and Environment in the Yangtze Valley 1992, 1, (1), 17-22 (in Chinese). Liu, S. P.; Duan, X. B.; Chen, D. Q.; Liao, F. C.; Chen, W., Studies on status of fishery resources in the

Middle Reaches of the Yangtze River. Acta Hydrobiologica Sinica 2005, 29, (6), 708-711 (in Chinese). Floehr, T.; Xiao, H.; Scholz-Starke, B.; Wu, L.; Hou, J.; Yin, D.; Zhang, X.; Ji, R.; Yuan, X.;

Ottermanns, R.; Roß-Nickoll, M.; Schäffer, A.; Hollert, H., Solution by dilution?—A review on the pollution status of the Yangtze River. Environmental Science and Pollution Research 2013, 20, (10), 6934-6971. Wang, J.; Bovee, T. F.; Bi, Y.; Bernhöft, S.; Schramm, K.-W., Aryl hydrocarbon receptor (AhR)

inducers and estrogen receptor (ER) activities in surface sediments of Three Gorges Reservoir, China evaluated with in vitro cell bioassays. Environmental Science and Pollution Research 2014, 21, (4), 3145-3155. Brack, W., Effect-directed analysis: a promising tool for the identification of organic toxicants in

complex mixtures? Analytical and Bioanalytical Chemistry 2003, 377, (3), 397-407. Burgess, R. M.; Ho, K. T.; Biales, A. D.; Brack, W., Recent Developments in Whole Sediment Toxicity

Identification Evaluations: Innovations in Manipulations and Endpoints. In Effect-Directed Analysis of Complex Environmental Contamination, Brack, W., Ed. 2011; Vol. 15, pp 19-40. Brack, W.; Ait-Aissa, S.; Burgess, R. M.; Busch, W.; Creusot, N.; Di Paolo, C.; Escher, B. I.; Hewitt, L.

M.; Hilscherova, K.; Hollender, J., Effect-directed analysis supporting monitoring of aquatic environments—An in-depth overview. Science of The Total Environment 2016, 544, 1073-1118. Wernersson, A.-S.; Carere, M.; Maggi, C.; Tusil, P.; Soldan, P.; James, A.; Sanchez, W.; Dulio, V.;

Broeg, K.; Reifferscheid, G.; Buchinger, S.; Maas, H.; Van Der Grinten, E.; O’Toole, S.; Ausili, A.; Manfra, L.; Marziali, L.; Polesello, S.; Lacchetti, I.; Mancini, L.; Lilja, K.; Linderoth, M.; Lundeberg, T.; Fjällborg, B.; Porsbring, T.; Larsson, D.; Bengtsson-Palme, J.; Förlin, L.; Kienle, C.; Kunz, P.; Vermeirssen, E.; Werner, I.; Robinson, C.; Lyons, B.; Katsiadaki, I.; Whalley, C.; den Haan, K.; Messiaen, M.; Clayton, H.; Lettieri, T.; Carvalho, R.; Gawlik, B.; Hollert, H.; Di Paolo, C.; Brack, W.; Kammann, U.; Kase, R., The European technical report on aquatic effect-based monitoring tools under the water framework directive. Environmental Sciences Europe 2015, 27, (1), 1-11. Simon, E.; van Velzen, M.; Brandsma, S. H.; Lie, E.; Loken, K.; de Boer, J.; Bytingsvik, J.; Jenssen, B.

M.; Aars, J.; Hamers, T.; Lamoree, M. H., Effect-Directed Analysis To Explore the Polar Bear Exposome: Identification of Thyroid Hormone Disrupting Compounds in Plasma. Environmental Science & Technology 2013, 47, (15), 8902-8912. Luebcke-von Varel, U.; Bataineh, M.; Lohrmann, S.; Loeffler, I.; Schulze, T.; Flueckiger-Isler, S.; Neca,

J.; Machala, M.; Brack, W., Identification and quantitative confirmation of dinitropyrenes and 3nitrobenzanthrone as major mutagens in contaminated sediments. Environment International 2012, 44, 31-39. Radović, J. R.; Thomas, K. V.; Parastar, H.; Díez, S.; Tauler, R.; Bayona, J. M., Chemometrics-assisted

effect-directed Analysis of crude and refined oil using comprehensive two-dimensional gas chromatography– time-of-flight mass spectrometry. Environmental science & technology 2014, 48, (5), 3074-3083.

23 ACS Paragon Plus Environment

Environmental Science & Technology

569 570

35.

571 572

36.

573 574 575 576 577 578 579 580 581

37.

582 583 584

38.

585 586 587

39.

588

40.

589

hormonal activities in surface water and sediment from the Pearl River system, South China, using effect‐

590

directed in vitro bioassays. Environmental Toxicology and Chemistry 2011, 30, (10), 2208-2215.

591 592

41.

593 594 595 596

42.

597 598

43.

599 600 601

44.

602 603

45.

Hecker, M.; Hollert, H., Effect-directed analysis (EDA) in aquatic ecotoxicology: state of the art and

future challenges. Environmental Science and Pollution Research 2009, 16, (6), 607-613. Di Paolo, C.; Seiler, T.-B.; Keiter, S.; Hu, M.; Muz, M.; Brack, W.; Hollert, H., The value of zebrafish

as an integrative model in effect-directed analysis-a review. Environmental Sciences Europe 2015, 27, (1), 8. Brack, W.; Altenburger, R.; Schüürmann, G.; Krauss, M.; López Herráez, D.; van Gils, J.; Slobodnik, J.;

Munthe, J.; Gawlik, B. M.; van Wezel, A.; Schriks, M.; Hollender, J.; Tollefsen, K. E.; Mekenyan, O.; Dimitrov, S.; Bunke, D.; Cousins, I.; Posthuma, L.; van den Brink, P. J.; López de Alda, M.; Barceló, D.; Faust, M.; Kortenkamp, A.; Scrimshaw, M.; Ignatova, S.; Engelen, G.; Massmann, G.; Lemkine, G.; Teodorovic, I.; Walz, K.-H.; Dulio, V.; Jonker, M. T. O.; Jäger, F.; Chipman, K.; Falciani, F.; Liska, I.; Rooke, D.; Zhang, X.; Hollert, H.; Vrana, B.; Hilscherova, K.; Kramer, K.; Neumann, S.; Hammerbacher, R.; Backhaus, T.; Mack, J.; Segner, H.; Escher, B.; de Aragão Umbuzeiro, G., The SOLUTIONS project: Challenges and responses for present and future emerging pollutants in land and water resources management. Science of The Total Environment 2015, 503–504, 22-31. Qiao, M.; Chen, Y.; Zhang, Q.; Huang, S.; Ma, M.; Wang, C.; Wang, Z., Identification of ah receptor

agonists in sediment of meiliang bay, taihu lake, China. Environmental science & technology 2006, 40, (5), 1415-1419. Qu, G.; Shi, J.; Wang, T.; Fu, J.; Li, Z.; Wang, P.; Ruan, T.; Jiang, G., Identification of

tetrabromobisphenol A diallyl ether as an emerging neurotoxicant in environmental samples by bioassaydirected fractionation and HPLC-APCI-MS/MS. Environmental science & technology 2011, 45, (11), 5009-5016. Zhao, J. L.; Ying, G. G.; Yang, B.; Liu, S.; Zhou, L. J.; Chen, Z. F.; Lai, H. J., Screening of multiple

Chapman, P., Presentation and interpretation of Sediment Quality Triad data. Ecotoxicology 1996, 5, (5),

327-339. Floehr, T.; Scholz-Starke, B.; Xiao, H.; Koch, J.; Wu, L.; Hou, J.; Wolf, A.; Bergmann, A.; Bluhm, K.;

Yuan, X.; Roß-Nickoll, M.; Schäffer, A.; Hollert, H., Yangtze Three Gorges Reservoir, China: A holistic assessment of organic pollution, mutagenic effects of sediments and genotoxic impacts on fish. Journal of Environmental Sciences 2015, 38, 63-82. Streck, H.; Schulze, T.; Brack, W., Accelerated membrane-assisted clean-up as a tool for the clean-up of

extracts from biological tissues. Journal of Chromatography A 2008, 1196, 33-40. Luebcke-von Varel, U.; Streck, G.; Brack, W., Automated fractionation procedure for polycyclic

aromatic compounds in sediment extracts on three coupled normal-phase high-performance liquid chromatography columns. Journal of Chromatography A 2008, 1185, (1), 31-42. Babich, H.; Borenfreund, E., Neutral red assay for toxicology in vitro. In In Vitro Methods of

Toxicology, Watson, R., Ed. CRC Press, Boca Raton, Florida, USA: 1992; pp 237–251.

24 ACS Paragon Plus Environment

Page 24 of 34

Page 25 of 34

Environmental Science & Technology

604 605 606 607

46.

608 609 610 611

47.

612 613

48.

614 615

49.

616

50.

Royal Society of Chemistry http://www.chemspider.com

617 618

51.

Vedani, A.; Dobler, M.; Hu, Z.; Smieško, M., OpenVirtualToxLab—A platform for generating and

619 620 621

52.

622 623 624

53.

625 626 627

54.

628 629 630

55.

631 632 633

56.

634 635 636

57.

637 638 639

58.

Gustavsson, L.; Hollert, H.; Joensson, S.; Van Bavel, B.; Engwall, M., Reed beds receiving industrial

sludge containing nitroaromatic compounds - Effects of outgoing water and bed material extracts in the umu-C genotoxicity assay, DR-CALUX assay and on early life stage development in zebrafish (Danio rerio). Environmental Science and Pollution Research 2007, 14, (3), 202-211. Luebcke-von Varel, U.; Machala, M.; Ciganek, M.; Neca, J.; Pencikova, K.; Palkova, L.; Vondracek, J.;

Loeffler, I.; Streck, G.; Reifferscheid, G.; Flueckiger-Isler, S.; Weiss, J. M.; Lamoree, M.; Brack, W., Polar Compounds Dominate in Vitro Effects of Sediment Extracts. Environmental Science & Technology 2011, 45, (6), 2384-2390. Pluskal, T.; Castillo, S.; Villar-Briones, A.; Orešič, M., MZmine 2: modular framework for processing,

visualizing, and analyzing mass spectrometry-based molecular profile data. BMC bioinformatics 2010, 11, (1), 1. Hug, C.; Ulrich, N.; Schulze, T.; Brack, W.; Krauss, M., Identification of novel micropollutants in

wastewater by a combination of suspect and nontarget screening. Environmental Pollution 2014, 184, 25-32.

exchanging in silico toxicity data. Toxicology letters 2015, 232, (2), 519-532. Kaisarevic, S.; Luebcke-von Varel, U.; Orcic, D.; Streck, G.; Schulze, T.; Pogrmic, K.; Teodorovic, I.;

Brack, W.; Kovacevic, R., Effect-directed analysis of contaminated sediment from the wastewater canal in Pancevo industrial area, Serbia. Chemosphere 2009, 77, (7), 907-913. Wölz, J.; Schulze, T.; Lübcke-von Varel, U.; Fleig, M.; Reifferscheid, G.; Brack, W.; Kühlers, D.;

Braunbeck, T.; Hollert, H., Investigation on soil contamination at recently inundated and non-inundated sites. Journal of Soils and Sediments 2011, 11, (1), 82-92. Wolf, A.; Bergmann, A.; Wilken, R.-D.; Gao, X.; Bi, Y.; Chen, H.; Schüth, C., Occurrence and

distribution of organic trace substances in waters from the Three Gorges Reservoir, China. Environmental Science and Pollution Research 2013, 20, (10), 7124-7139. Brack, W.; Schirmer, K.; Kind, T.; Schrader, S.; Schuurmann, G., Effect-directed fractionation and

identification of cytochrome P4501A-inducing halogenated aromatic hydrocarbons in a contaminated sediment. Environmental Toxicology and Chemistry 2002, 21, (12), 2654-2662. Wang, J.; Bernhöft, S.; Pfister, G.; Schramm, K.-W., Water exposure assessment of aryl hydrocarbon

receptor agonists in Three Gorges Reservoir, China using SPMD-based virtual organisms. Science of The Total Environment 2014, 496, 26-34. Bols, N.; Schirmer, K.; Joyce, E.; Dixon, D.; Greenberg, B.; Whyte, J., Ability of polycyclic aromatic

hydrocarbons to induce 7-ethoxyresorufin-o-deethylase activity in a trout liver cell line. Ecotoxicology and Environmental Safety 1999, 44, (1), 118-128. Brack, W.; Schirmer, K., Effect-directed identification of oxygen and sulfur heterocycles as major

polycyclic aromatic cytochrome P4501A-Inducers in a contaminated sediment. Environmental Science & Technology 2003, 37, (14), 3062-3070.

25 ACS Paragon Plus Environment

Environmental Science & Technology

640 641 642

59.

643 644 645

60.

646 647 648

61.

649 650 651

62.

652 653 654

63.

655 656

64.

657 658

65.

659 660 661

66.

662 663

67.

664 665 666

68.

667 668 669

69.

670 671

70.

672 673

71.

674 675

72.

Grung, M.; Naes, K.; Fogelberg, O.; Nilsen, A. J.; Brack, W.; Luebcke-von Varel, U.; Thomas, K. V.,

Effects-Directed Analysis of Sediments From Polluted Marine Sites in Norway. Journal of Toxicology and Environmental Health-Part a-Current Issues 2011, 74, (7-9), 439-454. Machala, M.; Vondráček, J.; Bláha, L.; Ciganek, M.; Neča, J., Aryl hydrocarbon receptor-mediated

activity of mutagenic polycyclic aromatic hydrocarbons determined using in vitro reporter gene assay. Mutation Research/Genetic Toxicology and Environmental Mutagenesis 2001, 497, (1), 49-62. Vondráček, J.; Chramostová, K.; Plíšková, M.; Bláha, L.; Brack, W.; Kozubík, A.; Machala, M.,

Induction of aryl hydrocarbon receptor–mediated and estrogen receptor–mediated activities, and modulation of cell proliferation by dinaphthofurans. Environmental toxicology and chemistry 2004, 23, (9), 2214-2220. Machala, M.; Švihálková-Šindlerová, L.; Pěnčíková, K.; Krčmář, P.; Topinka, J.; Milcová, A.;

Nováková, Z.; Kozubík, A.; Vondráček, J., Effects of methylated chrysenes on AhR-dependent and -independent toxic events in rat liver epithelial cells. Toxicology 2008, 247, (2–3), 93-101. Marvanová, S. a.; Vondráček, J.; Pěnččíková, K. i.; Trilecová, L.; Krčmářř, P.; Topinka, J.; Nováková,

Z.; Milcová, A.; Machala, M., Toxic effects of methylated benz [a] anthracenes in liver cells. Chemical research in toxicology 2008, 21, (2), 503-512. Wik, A., Toxic components leaching from tire rubber. Bulletin of Environmental Contamination and

Toxicology 2007, 79, (1), 114-119. De Wever, H.; Verachtert, H., Biodegradation and toxicity of benzothiazoles. Water Research 1997, 31,

(11), 2673-2684. Kumata, H.; Sanada, Y.; Takada, H.; Ueno, T., Historical trends of N-cyclohexyl-2-benzothiazolamine,

2-(4-morpholinyl) benzothiazole, and other anthropogenic contaminants in the urban reservoir sediment core. Environmental science & technology 2000, 34, (2), 246-253. Kloepfer, A.; Jekel, M.; Reemtsma, T., Occurrence, sources, and fate of benzothiazoles in municipal

wastewater treatment plants. Environmental science & technology 2005, 39, (10), 3792-3798. He, H.; Hu, G.; Sun, C.; Chen, S.; Yang, M.; Li, J.; Zhao, Y.; Wang, H., Trace analysis of persistent

toxic substances in the main stream of Jiangsu section of the Yangtze River, China. Environmental Science and Pollution Research 2011, 18, (4), 638-648. Noguerol, T.-N.; Boronat, S.; Casado, M.; Raldúa, D.; Barceló, D.; Piña, B., Evaluating the interactions

of vertebrate receptors with persistent pollutants and antifouling pesticides using recombinant yeast assays. Analytical and bioanalytical chemistry 2006, 385, (6), 1012-1019. Laflamme, R.; Hites, R. A., The global distribution of polycyclic aromatic hydrocarbons in recent

sediments. Geochimica et cosmochimica Acta 1978, 42, (3), 289-303. Grimmer, G.; Böhnke, H., Profile analysis of polycyclic aromatic hydrocarbons and metal content in

sediment layers of a lake. Cancer letters 1975, 1, 75-83. Jensen, T. E.; Hites, R. A., Aromatic diesel emissions as a function of engine conditions. Analytical

Chemistry 1983, 55, (4), 594-599.

26 ACS Paragon Plus Environment

Page 26 of 34

Page 27 of 34

Environmental Science & Technology

676 677 678

73.

679 680

74.

681 682

75.

683 684

76.

685 686

77.

687 688 689

78.

690 691 692

79.

693 694 695

80.

696 697

81.

698 699

82.

700 701

83.

702 703 704

84.

705 706

85.

707 708 709

86.

Rice, J. E.; Rivenson, A.; Braley, J.; LaVoie, E. J., Methylated derivatives of pyrene and fluorene:

evaluation of genotoxicity in the hepatocyte/DNA repair test and tumorigenic activity in newborn mice. Journal of Toxicology and Environmental Health, Part A Current Issues 1987, 21, (4), 525-532. McMahon, C. K.; Tsoukalas, S. N., Polytiuclear Aromatic Hydrocarbons in Forest Fire Smoke. 1978, 3,

61-73. Pan, G.; Crowley, D.; Lehmann, J., Burn to air or burial in soil: the fate of China’s straw residues. IBI

web post 2011. Oregon State University Pollution controls used during China Olympics could save lives if continued;

Available at: http://oregonstate.edu/ua/ncs/node/13842. Ali, R.; Siddiqui, N., Biological aspects of emerging benzothiazoles: a short review. Journal of

Chemistry 2013, 2013. Engels, H. W.; Weidenhaupt, H. J.; Pieroth, M.; Hofmann, W.; Menting, K. H.; Mergenhagen, T.;

Schmoll, R.; Uhrlandt, S., Rubber, 4. Chemicals and additives. Ullmann's Encyclopedia of Industrial Chemistry 2004. He, G.; Zhao, B.; Denison, M. S., Identification of benzothiazole derivatives and polycyclic aromatic

hydrocarbons as aryl hydrocarbon receptor agonists present in tire extracts. Environmental Toxicology and Chemistry 2011, 30, (8), 1915-1925. Milanova, E.; Ellis, S.; Sithole, B., Aquatic toxicity and solution stability of two organic corrosion

inhibitors: 2-mercaptobenzothiazole and 1, 2, 3-benzotriazole. Nordic Pulp & Paper Research Journal 2001, 16, (3), 215-218. Rajesh, N.; Subramanian, M., A study of the extraction behavior of thallium with tribenzylamine as the

extractant. Journal of hazardous materials 2006, 135, (1), 74-77. Malvankar, P. L.; Shinde, V. M., Extraction of cobalt and nickel salicylates: application to steel and

industrial wastewater samples. Analyst 1990, 115, (11), 1473-1476. Kalidhasan, S.; Rajesh, N., Simple and selective extraction process for chromium (VI) in industrial

wastewater. Journal of Hazardous Materials 2009, 170, (2–3), 1079-1085. Botalova, O.; Schwarzbauer, J., Geochemical characterization of organic pollutants in effluents

discharged from various industrial sources to riverine systems. Water, Air, & Soil Pollution 2011, 221, (1-4), 7798. Botalova, O.; Schwarzbauer, J.; Sandouk, N. a., Identification and chemical characterization of specific

organic indicators in the effluents from chemical production sites. Water Research 2011, 45, (12), 3653-3664. Helmick, J. S.; Fiddler, W., Thermal decomposition of the rubber vulcanization agent, zinc

dibenzyldithiocarbamate, and its potential role in nitrosamine formation in hams processed in elastic nettings. Journal of Agricultural and Food Chemistry 1994, 42, (11), 2541-2544.

27 ACS Paragon Plus Environment

Environmental Science & Technology

710 711 712

87.

713 714

88.

715 716 717

89.

718 719

90.

Brack, W.; Schirmer, K.; Erdinger, L.; Hollert, H., Effect-directed analysis of mutagens and

ethoxyresorufin-O-deethylase inducers in aquatic sediments. Environmental Toxicology and Chemistry 2005, 24, (10), 2445-2458. Simon, E.; Lamoree, M. H.; Hamers, T.; de Boer, J., Challenges in effect-directed analysis with a focus

on biological samples. TrAC Trends in Analytical Chemistry 2015, 67, (0), 179-191. Schymanski, E. L.; Jeon, J.; Gulde, R.; Fenner, K.; Ruff, M.; Singer, H. P.; Hollender, J., Identifying

small molecules via high resolution mass spectrometry: communicating confidence. Environmental science & technology 2014, 48, (4), 2097-2098. Stephensen, E.; Adolfsson-Erici, M.; Hulander, M.; Parkkonen, J.; Förlin, L., Rubber additives induce

oxidative stress in rainbow trout. Aquatic toxicology 2005, 75, (2), 136-143.

720

28 ACS Paragon Plus Environment

Page 28 of 34

Page 29 of 34

Environmental Science & Technology

721

Table 1. Nontarget compounds identified in the sediments of Kaixian Compound

Identification level

89

Molecular formula

Chemspider

CAS

50

ID

Exact

Use/origin

AhR binding

Mass

Experimental

Predicted a

- 90

+

[amu] 2-mercaptobenzothiazole

1

C7H5NS2

608157

149-30-4

167.251

Rubber additive

(MBT)

+

79

benzothiazole

1

C7H5NS

6952

95-16-9

135.186

Rubber additive

+79 69

+

N-t-butyl-2-

2

C11H14N2S2

6960

95-31-8

238.372

Rubber additive

n.e

-

1

C8H7NS

11494

615-22-5

181.278

Rubber additive

+

2

C14H8N

746290

n.a

300.422

Byproduct

benzothiazolesulfenamide (TBBS) 2-(methylthio)benzothiazole

79

+

(MTBT) 2,2’-sulfanediylbis(1,3benzothiazole)

of

n.e

+

mercaptobenzothiazole

tribenzylamine

1

C21H21

22739

620-40-6

287.398

Electroplating extractant

n.e

+

dibenzylamine

2

C14H15

7373

103-49-1

197.276

Byproduct

n.e

-

of

rubber

vulcanization

722 723 724 725 726 727

Identification level based on Schymanski et al., 1= confirmed, 2=probable structure. n.a: no available. n.e. indicates that no literature was found related to the AhR-mediated activity of the compound. a Estimated by QSAR -: Not AhR agonist reported in literature or predicted to have no binding affinity to the AhR by QSAR modeling; +: AhR agonist reported in literature or predicted to have binding affinity to the AhR by QSAR modeling; A compound was considered as potentially AhR agonist, either reported in literature or predicted as positive by QSAR modeling.

29

ACS Paragon Plus Environment

Environmental Science & Technology

728

Figure captions

729

Figure 1: Flowchart of the applied EDA strategy. NP-HPLC: Normal Phase-High Performance

730

Liquid Chromatography; F: Fraction; GC/MS: Gas Chromatography-mass spectrometry; LC-HRMS:

731

Liquid Chromatography-High Resolution Mass Spectrometry; QSAR: Quantitative Structure-Activity

732

Relationships.

733

Figure 2: Aryl hydrocarbon receptor-mediated activities of parent (par), reconstituted (rec)

734

extracts, arithmetic sum of fraction BEQs (sum) and fraction activities (1-18) of four sediment

735

extracts. Asterisks denote significant differences between fractions and process control (Non-

736

parametric Shapiro-Wilk ANOVA on ranks with Holm-Sidak’s post hoc test, p≤0.05). dw: dry weight.

737

N.A.-No data available.

738

Figure 3: Structures of compounds identified and tentatively identified by target and nontarget

739

screening. 1# Benzo[a]anthracene; 2# Chrysene; 3# Benzo[b]fluoranthene: 4# Benzo[k]fluoranthene:

740

5# Indeno[1,2,3-cd]pyrene; 6# dibenz[a,h]anthracene: 7# benzo[g,h,i]perylene; 8# Dinaphtho[1,2-

741

b;1’,2’-d]furan; 9# Dinaphtho[1,2-b;2’,3’-d]furan; 10#2-mercaptobenzothiazole ; 11# benzothiazole;

742

12#2-(methylthio)benzothiazole; 13#2,2’-sulfanediylbis(1,3-benzothiazole); 14# tribenzylamine.

743

30 ACS Paragon Plus Environment

Page 30 of 34

Page 31 of 34

Environmental Science & Technology

744 745 746 747

Figure 1. Flowchart of the applied EDA strategy. NP-HPLC: Normal Phase-High Performance Liquid Chromatography; F: Fraction; GC/MS: Gas Chromatography-mass spectrometry; LC-HRMS: Liquid Chromatography-High Resolution Mass Spectrometry; QSAR: Quantitative Structure-Activity Relationships.

748

31 ACS Paragon Plus Environment

Environmental Science & Technology

749 750 751 752 753

Figure 2. Aryl hydrocarbon receptor-mediated activities of parent (par), reconstituted (rec) extracts, arithmetic sum of fraction BEQs (sum) and fraction activities (1-18) of four sediment extracts. Asterisks denote significant differences between fractions and process control (Non-parametric Shapiro-Wilk ANOVA on ranks with Holm-Sidak’s post hoc test, p≤0.05). dw: dry weight. N.A.-No data available.

754

32 ACS Paragon Plus Environment

Page 32 of 34

Page 33 of 34

Environmental Science & Technology

755 756

757 758 759 760 761 762 763

Figure 3: Structures of compounds identified and tentatively identified by target and nontarget screening. 1# Benzo[a]anthracene; 2# Chrysene; 3# Benzo[b]fluoranthene: 4# Benzo[k]fluoranthene: 5# Indeno[1,2,3cd]pyrene; 6# dibenz[a,h]anthracene: 7# benzo[g,h,i]perylene; 8# Dinaphtho[1,2-b;1’,2’-d]furan; 9# Dinaphtho[1,2-b;2’,3’-d]furan; 10#2-mercaptobenzothiazole ; 11# benzothiazole; 12#2(methylthio)benzothiazole; 13#2,2’-sulfanediylbis(1,3-benzothiazole); 14# tribenzylamine.

33 ACS Paragon Plus Environment

Environmental Science & Technology

TOC 254x190mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 34 of 34