High-Throughput Effect-Directed Analysis Using ... - ACS Publications

Mar 16, 2018 - M. H. A.; Andersson, P. L.; Legler, J.; Brouwer, A. In vitro profiling of the endocrine-disrupting potency of brominated flame retardan...
0 downloads 2 Views 797KB Size
Subscriber access provided by Queen Mary, University of London

Ecotoxicology and Human Environmental Health

High-throughput Effect-Directed Analysis using downscaled in vitro reporter gene assays to identify endocrine disruptors in surface water Nick Zwart, Shan Li Nio, Corine J Houtman, Jacob De Boer, Jeroen Kool, Timo Hamers, and Marja H. Lamoree Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b06604 • Publication Date (Web): 16 Mar 2018 Downloaded from http://pubs.acs.org on March 20, 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 31

Environmental Science & Technology

1

1

High-throughput Effect-Directed Analysis using downscaled in vitro reporter

2

gene assays to identify endocrine disruptors in surface water

3 4

1

1

2

1

3

1

Nick Zwart , Shan Li Nio , Corine J. Houtman , Jacob de Boer , Jeroen Kool , Timo Hamers , Marja H. Lamoree

1

5 6

1

Department Environment & Health, VU University, Amsterdam, The Netherlands 2

7 8

The Water Laboratory, Haarlem, The Netherlands

3

Biomolecular Analysis group, VU University, The Netherlands

9 10

E-mail contact: [email protected]

11 12

Effect-Directed Analysis (EDA) is a commonly used approach for effect-based identification of endocrine

13

disruptive chemicals in complex (environmental) mixtures. For routine toxicity assessment of e.g. water

14

samples however, current EDA approaches are considered time consuming and laborious. We achieved

15

faster EDA and identification, by downscaling of sensitive cell-based hormone reporter gene assays and

16

increasing fractionation resolution to allow testing of smaller fractions with reduced complexity. The

17

high-resolution EDA approach is demonstrated by analysis of four environmental passive sampler

18

extracts. Downscaling of the assays to 384-well format allowed analysis of 64 fractions in triplicate (or

19

192 fractions without technical replicates) without affecting sensitivity compared to the standard 96-

20

well format. Through a parallel exposure method, agonistic and antagonistic androgen- and estrogen

21

receptor activity could be measured in a single experiment following a single fractionation. From sixteen

22

selected candidate compounds, identified through non-targeted analysis, thirteen could be confirmed

23

chemically and ten were found biologically active of which the most potent non-steroidal estrogens

24

were identified as oxybenzone and piperine. The increased fractionation resolution and the higher

ACS Paragon Plus Environment

Environmental Science & Technology

2 25

throughput that downscaling provides, allows for future application in routine high-resolution screening

26

of large numbers of samples in order to accelerate identification of (emerging) endocrine disruptors.

27 28

Introduction

29

Endocrine disruption is an important endpoint in toxicological and environmental screening as well as

30

the water quality control of the drinking water production process. Endocrine disruptive chemicals

31

(EDCs) like estrogens1 and androgens2 have been detected in the aquatic environment as pollutants. A

32

large portion of EDCs are emitted into the aquatic environment through urban (steroid hormones, flame

33

retardants, plasticizers) or industrial wastewater,3 agricultural runoff4 or deposition after combustion

34

(PAHs).5 Subsequent exposure can lead to disrupted signaling of endogenous hormones. Further liver

35

metabolism can enhance receptor binding potency of (inactive) pollutants and increase the risk of

36

endocrine disruption.6 While well characterized potent EDCs are actively monitored, unknown

37

compounds with endocrine disruptive potency, including their metabolites, transformation products and

38

degradation products, remain to be discovered. The ability to detect and identify relevant and yet

39

unknown EDCs is essential to efforts aimed at reducing their presence in the aquatic environment and

40

reducing human exposure.

41

Cell-based reporter gene assays have been used in Effect-Directed Analysis (EDA) for the identification of

42

(emerging) EDCs in environmental samples.7–9 Via EDA, compounds not analyzed by routine (chemical)

43

analysis are identified based on their biological activity in reporter bioassays. Activity measured by the

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Environmental Science & Technology

3 44

reporter gene assays in particular fractions collected during chromatographic separation can be

45

correlated to mass spectrometry data. A response in one or more fractions can direct efforts to identify

46

the compound responsible for the observed activity to a limited number of corresponding masses on the

47

mass chromatogram.

48

Reducing fraction complexity through high-resolution fractionation decreases the number of

49

compounds and masses to be identified per fraction, but increases the total number of fractions.

50

Luciferase reporter gene cell lines, while showing high sensitive towards their respective (ant)agonists,

51

are usually performed in a 96-well plate format which limits the number of samples or fractions that can

52

be analyzed simultaneously.

53

This study focused on improvement of the current EDA approach by increasing throughput and

54

resolution to allow for faster identification as a step forward to a more routine application of EDA in

55

future (surface) water quality assessment.

56

Firstly, for high-resolution EDA of endocrine disruptive chemicals, an androgen receptor (AR) (AR-

57

EcoScreen),

58

EcoScreen GR-KO),11 aryl hydrocarbon receptor (AhR) receptor (DR-Luc)12 and estrogen receptor (ER)

59

(VM7Luc4E2, formerly known as BG1Luc4E2; further referred to as ER-Luc)13 reporter gene assay were

60

downscaled from 96- to 384-well plate format. Throughput is further improved by introducing a method

61

for parallel exposure of multiple endpoints with samples or fractions from a single source plate. In

62

addition, a metabolic system was incorporated in the downscaled assays, to allow for formation and

63

detection of active metabolites from inactive or less active pollutants.

64

Secondly, downscaled assays, using the parallel exposure method, were applied in an EDA approach to

65

analyze four passive sampler extracts which were separated by UPLC and collected as either 64 or 192

66

fractions. The mass spectra (recorded in parallel) were analyzed at retention times that correlated with

10

a recently developed AR-EcoScreen glucocorticoid receptor (GR) knockout mutant (AR-

ACS Paragon Plus Environment

Environmental Science & Technology

4 67

active fractions to select masses for identification. A qualitative non-targeted screening was performed

68

and selected candidates were confirmed chemically and biologically.

69

Materials and methods

70

Materials

71

Dulbecco's Modified Eagle Medium/Nutrient Mixture F-12 (DMEM/F12) medium with glutamax, phenol-

72

free DMEM/F12 medium with L-glutamine, low glucose phenol-free DMEM medium, and fetal bovine

73

serum were obtained from Gibco (Eggenstein, Germany), penicillin/streptomycin, G418, hygromycin,

74

zeocin, ATP, co-enzyme A, formic acid, acetonitrile (HPLC grade) and methanol (Chromasolv) from Sigma

75

(Zwijndrecht, The Netherlands), luciferin from Promega (Fitchburg, WI, USA), DTT (dithiothreitol) from

76

Duchefa (Haarlem, The Netherlands) and Aroclor 1254 induced rat liver S9 fraction was obtained from

77

MP Biomedicals (Santa Ana, USA). Water was purified on a Milli-Q Reference A+ purification system

78

(Millipore, Bedford, MA, USA). Reference compounds used for validation of the downscaled test

79

methods and candidate compounds for confirmation of hits were obtained from various suppliers (see

80

Supporting Information Table S1) and were dissolved in DMSO (Acros, Geel, Belgium).

81

Cell culture conditions

82

AR-EcoScreen (CHO-K1), exhibiting residual GR sensitivity, and AR-EcoScreen GR-KO (CHO-K1) cells, with

83

exclusive AR sensitivity, were maintained as described by Satoh et al.10 ER-Luc (MCF7 human breast

84

carcinoma) cells were maintained as described by Rogers and Denison.13 Briefly, cells were cultured at

85

37 °C with 5% CO2 in DMEM/F12 medium (with 10% fetal bovine serum and 1% penicillin/streptomycin),

86

further referred to as culture medium, and sub-cultured twice weekly. DR-Luc cells were maintained and

87

exposed as described in Supporting Information, Section S1.1.

88

Reporter gene assay protocols

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Environmental Science & Technology

5 89

Reporter gene assays were performed in transparent polystyrene CellStar 96-well plates (655160)

90

(Greiner Bio-One, Alphen aan den Rijn, The Netherlands) (as the standard format to compare against the

91

downscaled format) or in white µclear polystyrene CellStar 384-well plates (781098) (Greiner Bio-One,

92

Alphen aan den Rijn, The Netherlands) (downscaled format) using the same cell lines. To maintain the

93

cell-to-surface-area ratio during seeding and compensate for smaller well size, reaction volumes used in

94

the downscaled format were approximately 3-fold lower compared to the standard format (see

95

Supporting Information Table S2). Potential effects of medium evaporation on the measurement were

96

reduced by filling the outer ring with 100 µL (96-well plates) or two outer rings with 34 µL (384-well

97

plates) sterile milliQ water which resulted in 60 or 240 wells available for measurements, respectively.

98

Prior to exposure, trypsinated cells were resuspended in phenol-free DMEM/F12 L-glutamine (AR-

99

EcoScreen and AR-EcoScreen GR-KO), low glucose phenol-free DMEM (ER-Luc) medium (with 5%

100

charcoal stripped fetal bovine serum), further referred to as AR-EcoScreen or ER-Luc assay medium,

101

respectively, and seeded at 200.000 cells/mL in 100 µL (96-well plates) or 34 µL aliquots (384-well

102

plates). Plates were incubated for 24 hours at 37 °C and 5% CO2. In each experiment (n=1), seeded cells

103

were exposed at t=0 in triplicate to compounds dissolved in 100 µL (96-well plates) or 34 µL assay

104

medium (384-well plates), at a final DMSO concentration of 1 or 5 µL/mL for the AR-EcoScreen (and GR-

105

KO) or ER-Luc cells, respectively, using a single (96-well plates) or 8-channel (384-well plates) manual

106

pipette. Reference compound dilution series were prepared in DMSO and, prior to exposure, diluted to

107

exposure concentrations (see Supporting Information Table S1) in the respective assay medium for each

108

cell type at a final DMSO concentration of 1 or 5 µL/mL for the AR-EcoScreen (and GR-KO) or ER-Luc

109

cells, respectively. In antagonism experiments, cells were additionally exposed to compounds or

110

fractions in the presence of 200 pM DHT in AR-EcoScreen (and GR-KO) or 4 pM 17β-estradiol (E2) in ER-

111

Luc by spiking assay medium with 5 µM DHT and 0.1 µM E2 in DMSO, respectively. Cytotoxicity was

112

measured in the exposed cells by adding resazurin dissolved in PBS at t=22 h leading to a final

ACS Paragon Plus Environment

Environmental Science & Technology

6 113

concentration of 21 nM and measuring fluorescence (λex=530 nm; λem=590 nm) immediately after

114

addition and at t=24 h. Conversion rate of resazurin into resorufin in cells exposed to test substances

115

was compared to the conversion rate in cells exposed to vehicle (DMSO) only. At t=24 h, medium was

116

aspirated and cells were lysed for 10 minutes on an orbital shaker at 700 rpm with 50 µL (96-well plates)

117

or 17 µL lysis mix (25 mM TRIS (pH 7.8), 2 mM DDT, 2 mM 1,2-cyclohexylenedinitrilotetraacetic acid

118

(CDTA), 10% Glycerol and 1% Triton-X100) (384-well plates). Luminescence was measured on a

119

Varioskan luminometer (Thermo Fisher Scientific, Waltham, MA, USA) for one second after injection of

120

100 µL or 34 µL of glow-mix (2 mM Trycin, 1.07 mM C4H2Mg5O14, 2.67 mM MgSO4, 0.1 mM EDTA, 33.3

121

mM DTT, 270 mM Co-enzyme A, 470 mM Luciferin, and 530 mM ATP) followed by quenching of the

122

reaction with 100 µL (96-well plates) or 34 µL 0.1 M NaOH (384-well plates).

123

Biotransformation

124

Prior to exposure of cells seeded on a 384-well plate in the downscaled format, compounds or fractions

125

were incubated for 90 minutes at 37 °C in 50 µL DMEM phenol-free low glucose medium (with a final

126

volume of 3.6% DMSO (v:v) to ensure solubility of the concentrated compounds during metabolism) in

127

96-well plates, with or, as a control, without addition of 1.7 µL S9-mix (300 µL rat liver S9 fraction per

128

mL, 33 mM KCl, 8mM MgCl2∙6H2O, 6.5 mM glucose-6-phosphate, 4 mM NADP, 100 mM sodium

129

phosphate buffer pH 7.4), to the 50 µL reaction volume at a final concentration of 33 µL S9-mix per mL

130

reaction volume (0.2 mg protein per mL reaction mixture), for generation of metabolites. Incubations

131

were performed on single compounds at concentrations from 1.81 to 109 µM (BPA), 0.181 to 181 µM

132

(flutamide) or 1.81 to 181 µM (tamoxifen). After incubation, 216 µL DMEM phenol-free low glucose

133

medium was added to the reactions, reducing DMSO concentrations to 0.7% and S9-mix concentrations

134

to 6.4 µL/mL. Cells seeded on 384-well plates were prepared for exposure by adding 24 µL of their

135

respective assay medium to the aspirated cells. In each experiment (n=1) 10 µL of diluted

136

biotransformation reaction mixtures were added to the cells in triplicate to reach the final volume of 34

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Environmental Science & Technology

7 137

µL with a maximum DMSO concentration of 0.2% and an S9-mix concentration of 1.9 µL per mL

138

exposure medium.

139

Sample preparation

140

Adsorption-based Speedisk (SD) and partitioning-based silicone rubber (SR) passive samplers were

141

deployed for a period of six weeks between August and October 2014 in the river Meuse at Eijsden and

142

in the effluent stream of an ultra-low loaded activated sludge municipal wastewater treatment plant

143

(WWTP) with 25% industrial contribution serving a population equivalent of approximately 300.000, in

144

the Netherlands. Briefly, sample preparation consisted of dichloromethane extraction of SD and hexane

145

extraction of SR samplers. SD and SR extracts were solvent-exchanged from dichloromethane and

146

hexane, respectively, by evaporation under a flow of nitrogen at room temperature and redissolving the

147

samples in 1:2:3 (v:v:v) MeOH:ACN:H2O to 8.4 SD/mL and 60 g SR/mL, respectively.

148

LC-MS analysis and fractionation

149

Fractionation of SD and SR sample extracts was performed on a Kinetex C18 (100 x 2.1 mm, 1.7 µm

150

particle size) column using an Agilent Infinity 1290 UPLC pump and autosampler. Extracts were injected

151

(10 µL) at a flowrate of 250 µL/min in 80% mobile phase A (100% H2O + 0.1% formic acid) and 20%

152

mobile phase B (100% ACN + 0.1% formic acid). The solvent gradient increased to 80% mobile phase B

153

over 15 minutes and was subsequently kept as such for 13 minutes. Post-column, the flow was split in a

154

9:1 ratio on a Quicksplit adjustable flow splitter (ASI, Richmond, CA, USA) with 9 parts being diverted to

155

a nanofraction collector and 1 part to a microTOF II time-of-flight mass spectrometer (Bruker Daltonics,

156

Billerica, MA, USA). The MS was equipped with an ESI source set to positive mode and scanned masses

157

from 50 m/z to 3000 m/z at 10 Hz. Corona and capillary voltages were set to 500 and 4500 V

158

respectively. Nebulizer pressure was kept at 2 bar and nitrogen drying gas flow was kept at 6 L/min.

159

Exposure of fractionated extracts

ACS Paragon Plus Environment

Environmental Science & Technology

8 160

In each experiment (n=1), sample fractions were collected in either 96 (64 fractions collected at 26

161

second intervals) or 384-well plates (192 fractions collected at 9 second intervals) filled with 10 µL or 4

162

µL 10% DMSO in milliQ water, respectively, as keeper to increase recoveries.14 Collected fractions were

163

dried in a Centrivap concentrator (Labconco Corp., Kansas City, MO, USA) for 5 hours at 25°C and

164

redissolved in 50 µL or 12 µL ER-Luc assay medium supplemented with 3.6% DMSO added to improve

165

solubility of compounds, in 96- or 384-well plates respectively, for 10 minutes at 700 rpm. Redissolved

166

fractions were diluted with 216 µL or 50 µL ER-Luc assay medium, lowering the DMSO concentration to

167

0.7% (considering potential residual DMSO used as keeper negligible), in 96- or 384-well plates

168

respectively. Cells seeded on 384 well plates in the downscaled format were prepared for exposure to

169

fractions by aspirating the medium and adding 24 µL AR-EcoScreen assay medium to AR-EcoScreen (and

170

GR-KO) cells or 24 µL ER-Luc assay medium with 0.43% DMSO to ER-Luc cells. Cells were exposed to a

171

single concentration by adding 10 µL of the redissolved fraction to prepared cells, in triplicate for

172

measuring 64 fractions or in a single well for measuring 192 fractions, using a digital stepper- or manual

173

8-channel pipette. This resulted in a final DMSO concentration of 0.2% (2 µL/mL) in AR-EcoScreen (and

174

GR-KO) and 0.5% (5 µL/mL) in ER-Luc cells during exposure compared to 0.1% and 0.5% DMSO during

175

exposure, respectively, described in in the reporter gene assay protocol. As a result from the addition of

176

10 µL dissolved fraction in ER-Luc assay medium to 24 µL AR-EcoScreen assay medium, AR-EcoScreen

177

(and GR-KO) cells were exposed in 29.4% ER-Luc and 70.6% AR-EcoScreen assay medium. From each

178

fraction plate, multiple seeded 384-well plates were exposed to measure different endpoints in parallel.

179

Reference compound dilution series prepared in DMSO were diluted in the same assay medium

180

composition and at the same DMSO concentration at which cells were exposed to fractions. Cells were

181

exposed, in triplicate, to reference compounds by adding 34 µL of the diluted compounds to aspirated

182

cells.

183

Data analysis

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

Environmental Science & Technology

9 184

Bioassay results were analyzed in Prism 5.04 (Graphpad Software Inc, San Diego, CA). For each serial

185

dilution data set, D'Agostino-Pearson test was used to test data normality, and Levene’s test to test

186

homogeneity of variance (significance P < 0.05). Dose response curves of reference compounds and

187

candidate compounds were fitted with a four parametric logistic function [Y=A+(B-A)/(1+(x/C)^D)],

188

where A and B denote minimal and maximal response respectively, C is the EC50 or IC50, D is the

189

hillslope and x represents the tested concentration. Significant differences (P