Oil Sands Derived Naphthenic Acids Are Oxidative Uncouplers and

Aug 13, 2018 - Canadian Rivers Institute, Department of Biomedical Science, Atlantic Veterinary College, University of Prince Edward Island , Charlott...
0 downloads 0 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

Ecotoxicology and Human Environmental Health

Oil Sands-Derived Naphthenic Acids are Oxidative Uncouplers and Impair Electron Transport in Isolated Mitochondria Kate Rundle, Mahmoud Sharaf, Don Stevens, Collins Kamunde, and Michael R. Van Den Heuvel Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b02638 • Publication Date (Web): 13 Aug 2018 Downloaded from http://pubs.acs.org on August 14, 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 32

Environmental Science & Technology

TOC art only 84x47mm (300 x 300 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Oil Sands-Derived Naphthenic Acids are Oxidative Uncouplers and Impair Electron Transport in Isolated Mitochondria

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28

Kate I. Rundle1, Mahmoud S. Sharaf2, Don Stevens2, Collins Kamunde2, and Michael R. van den Heuvel1,3* 1

Canadian Rivers Institute, Department of Biomedical Science, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Canada, C1A 4P3 2

Department of Biomedical Science, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Canada, C1A 4P3 3

Canadian Rivers Institute, Department of Biology, University of Prince Edward Island, Charlottetown, Canada, C1A 4P3

* Author to whom correspondence should be addressed Department of Biology University of Prince Edward Island 550 University Avenue Charlottetown, Canada C1A 4P3 ph: 1-902-388-0895 fax: 1-902-566-0740 [email protected]

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

Environmental Science & Technology

Rundle et al. 2

29 30 31

ABSTRACT

32

cause numerous acute and chronic effects in fishes. However, the mechanism of toxicity underlying these

33

effects has not been fully elucidated. Due to their carboxylic acid moiety and the reported disruption of

34

cellular bioenergetics by similar structures, we hypothesized that NAs would uncouple mitochondrial

35

respiration with the resultant production of reactive oxygen species (ROS). Naphthenic acids were

36

extracted and purified from 17-year-old oil sands tailings waters yielding an extract of 99% carboxylic

37

acids with 90% fitting the classical O2-NA definition. Mitochondria were isolated from rainbow trout

38

liver and exposed to different concentrations of NAs. Mitochondrial respiration, membrane potential, and

39

ROS emission were measured using the Oroboros fluorespirometry system. Additionally, mitochondrial

40

ROS emission and membrane potential were evaluated with real-time flow cytometry. Results showed

41

NAs uncoupled oxidative phosphorylation, inhibited respiration, and increased ROS emission. The

42

effective concentration (EC50) and inhibition concentration (IC50) values for the endpoints measured

43

ranged from 21.0 to 157.8 mg/L, concentrations similar to tailings waters. For the same endpoints,

44

EC10/IC10 values ranged from 11.8 to 66.7 mg/L, approaching concentrations found in the environment.

45

These data unveil mechanisms underlying effects of NAs that may contribute to adverse effects on

46

organisms in the environment.

Naphthenic acids (NAs) are predominant compounds in oil sands-influenced waters. These acids

47 48

Keywords: oil sands, naphthenic acid, mitochondria, ROS, membrane potential

49 50 51 52 53 54 55

ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 32

Rundle et al. 3

56 57

INTRODUCTION The Athabasca oil sands in Alberta, Canada are among the largest oil deposits in the world. In

58

2016, approximately 2.5 million barrels of crude oil were being produced every day and an estimated 165

59

billion barrels remained in oil reserves.1 Most of Alberta's crude oil is produced using the Clarke hot

60

water extraction method, where large volumes of hot water and a caustic soda (NaOH) are added to

61

promote the separation of bitumen that is further processed into synthetic crude oil.2,3 Large volumes of

62

waste water, commonly referred to as oil sands process-affected water (OSPW) are produced during this

63

process and are stored in tailings ponds. At the end of 2013, these tailings covered an area of

64

approximately 220 km2 and housed over 975 million liters of OSPW.1,2 OSPW contains residual bitumen,

65

dissolved salts, minerals, and a complex mixture of inorganic and organic compounds.2,4

66

Naphthenic acids (NAs) make up a large proportion of the organic compounds present in OSPW

67

and have been implicated for much of the observed OSPW toxicity of various species.2,4 These acids are a

68

mixture of alkyl-substituted acyclic and cycloaliphatic carboxylic acids that are natural components of

69

petroleum.4,5 Classical NAs refer to compounds that follow the general formula CnH2n+ZO24, where n is the

70

carbon number and z is the hydrogen deficiency, though this definition has been broadened to include

71

carboxylic acids with additional oxygens and compounds containing N or S, the definition we used

72

herein. Organisms exposed to NAs extracted from OSPW exhibit inhibition of growth;6,7 disruption of the

73

endocrine system and reproductive physiology;8 increased developmental abnormalities;9,10 and exposure

74

is correlated with opportunistic disease in fishes.11 NAs are considered bioavailable and have been

75

suggested to act via narcosis or physical disruption of the cell membrane.12,13 However, because of the

76

complexity of NAs, a model of multiple mechanisms of toxic action also has been suggested.7,14

77

While no specific mechanistic studies have been conducted on NAs, their carboxylic acid moiety

78

suggests that NAs could have the potential to disrupt mitochondrial membrane permeability and

79

consequently uncouple oxidative phosphorylation (OXPHOS). Supporting this theme, another group of

80

carboxylic acids prevalent in pulp and paper effluents, resin acids, were found to alter mitochondrial

81

membrane integrity, uncouple OXPHOS, and decrease ATP production.15,16 Uncoupling of OXPHOS is

ACS Paragon Plus Environment

Page 5 of 32

Environmental Science & Technology

Rundle et al. 4

82

the disconnection of two reactions. The first reaction is the movement of electrons through a series of

83

protein complexes in the inner mitochondrial membrane, referred to as the electron transport system

84

(ETS), causing the creation of a hydrogen ion gradient known as the proton motive force. The second is a

85

phosphorylation event by ATP synthase powered by the proton motive force that is responsible for

86

converting ADP to ATP.17 When OXPHOS is uncoupled, the proton motive force is dissipated without

87

being used to produce ATP.17 While mild uncoupling has been associated with a decrease in production of

88

reactive oxygen species (ROS), uncoupling that occurs in mitochondria where the ETS has been inhibited

89

increases production of ROS.18 Reactive oxygen species play a role in cell signaling,19 but in excess, can

90

impose oxidative stress, leading to the release of cytochrome c from mitochondria resulting in apoptosis.20

91

We hypothesized that NAs can uncouple OXPHOS and alter the overall mitochondrial function.

92

This hypothesis was tested by examining the effects of NAs on the respiration of mitochondria isolated

93

from rainbow trout liver. Based on this, we predicted that NAs would uncouple OXPHOS and inhibit the

94

ETS, and that these effects would lead to oxidative stress. These predictions were tested using NAs

95

extracted from aged OSPW using acid precipitation and by measuring oxygen consumption, ROS

96

emission, and membrane potential of mitochondria isolated from rainbow trout liver using respirometry

97

and real-time flow cytometry.

98 99 100 101

MATERIALS AND METHODS Chemicals and Reagents Sucrose, potassium phosphate, EGTA, bovine serum albumin (BSA), aprotinin, malate,

102

glutamate, adenosine diphosphate (ADP), and hydrogen peroxide were supplied by Sigma Aldrich

103

(Oakville, ON). Tris-hydrogen chloride was from EMD Chemicals Inc. (Gibbstown, NJ, USA). The

104

fluorescent dyes, Amplex red and H2DCFDA, as well as horseradish peroxidase (HRP) were from

105

Invitrogen (Oregon, USA). The membrane potential indicator lipophilic cation tetraphenylphosphonium

106

chloride (TPP+), was from Alfa Aesar (Heysham, England) and the fluorescent dye JC-1 was from

107

Biotium Inc. (California, USA). Sulfuric acid and diethylaminoethyl (DEAE) cellulose used in the NA

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 32

Rundle et al. 5

108

extraction were from Sigma Aldrich (Oakville, ON) and hydrochloric acid, dichloromethane (DCM),

109

ammonium hydroxide, and acetonitrile were supplied by Caledon Laboratories Chemicals (Georgetown,

110

ON). Oligomycin-A, carbonyl cayanide-4-(trifluoromethoxy)phenylhydrazone (FCCP), and antimycin-A

111

were from Sigma Aldrich (Oakville, ON).

112

Naphthenic Acid Extraction

113

Naphthenic acids were extracted from 17-year-old tailings water from Syncrude Canada Ltd.'s

114

experimental pond 10 based on previously described methods.21,22 In the present study, the water used

115

was from the same source sample used in previous studies22,23,33 Briefly, 4,000 L of tailings water were

116

acidified using sulfuric acid and the precipitate was removed, concentrated by centrifugation (17,000 × g)

117

and dissolved in 0.1 M NaOH. The sample was then centrifuged again at 17,000 × g to remove any clay

118

or insoluble material and the supernatant was vacuum filtered through cleaned DEAE cellulose to remove

119

any humic material. The filtered supernatant was extracted using DCM in glass liquid-liquid extractors

120

for 72 h to remove neutral compounds. The sample was again acidified with sulfuric acid and centrifuged

121

to obtain NAs. The NA extract was washed several times with distilled water to remove residual sulfuric

122

acid and was then freeze dried to obtain a solid material. All concentrations in this study are based on

123

weighing out the dried solid material. As there are no analytical standards for NAs, this represents the

124

most accurate concentration available.

125

Naphthenic Acid Characterization

126

Mass spectrometry to quantify the mixture was performed by direct injection of a solution of the

127

extracted naphthenic acids into a Thermo Velos Orbitrap high-resolution mass spectrometer equipped

128

with an electrospray ionization interface in negative mode. Naphthenic acids dissolved in a 10 mM 1:1

129

NH4OH:acetonitrile solution were introduced into the mass spectrometer by direct infusion at 3 µL/min.

130

Ion intensity data at 100,000 mass resolution were evaluated against the predicted molecular weights of

131

naphthenic acids according to the formula CnH2n+ZO2. Further chemical characterization using NMR and

132

IR spectrometry was previously conducted on an earlier sample of NAs extracted by identical methods

ACS Paragon Plus Environment

Page 7 of 32

Environmental Science & Technology

Rundle et al. 6

133

and is consistent with a pure mixture of carboxylic acids.23 The characteristics reflected in those analyses

134

would be expected to be similar to the sample used.

135

A sample of naphthenic acid extract was sent to Galbraith Laboratories (Knoxville, TN) for

136

analysis of percent carbon, hydrogen, nitrogen, sulfur, and oxygen. Briefly, carbon, nitrogen, hydrogen,

137

and sulfur were determined by combustion of the sample in pure oxygen at high temperature (920 – 1350

138

°C) and the combustion products CO2, H2O, N2, and SO2 were automatically analyzed in a self-

139

integrating, steady state thermal conductivity analyzer (PerkinElmer 2400 Series II CHNS/O

140

Analyzer).24,25 Oxygen percentage was determined by pyrolysis of the sample in a helium atmosphere.

141

After pyrolysis, N, H, and CO are produced and passed through an absorption filter where halogenated

142

compounds were retained and N, H, and CO were separated via a chromatographic column. The amount

143

of CO was automatically analyzed the PerkinElmer 2400 Series II CHNS/O Analyzer.26 Detection limits

144

for C, H, N, and O were 0.5% and 0.1% for sulfur.

145

Mitochondria Isolation

146

Fish were obtained from Ocean Farms (Brookvale, PE) and were approximately 1 year old. They

147

were maintained in a 250 L tank in the Atlantic Veterinary College Aquatics Facility and fed 1% of their

148

weight daily. Mitochondria were isolated from rainbow trout (Oncorhynchus mykiss; mean fish weight

149

135±11 g) livers using differential centrifugation.27 Mitochondria were re-suspended in mitochondrial

150

respiration buffer (MRB: 10 mM Tris HCl, 25 mM KH2PO4, 100 mM KCl, 1 mg/mL BSA, 2 mg/mL

151

aprotinin, pH 7.3) and immediately used for experiments. Protein concentrations of the mitochondrial

152

suspensions were measured by the method of Bradford protein assay28 using BSA as a standard.

153

Measurement of Mitochondrial Respiration, H2O2 Emission, and Mitochondrial Membrane

154

Potential (∆Ψmt)

155

Oxygen consumption, H2O2 emission (representing a portion of ROS emission), and

156

mitochondrial membrane potential (∆Ψmt) were measured simultaneously using the Oxygraph-2k

157

fluorespirometer (Oroboros Instruments, Innsbruck, Austria) as described previously.29 Briefly, oxygen

158

concentration and changes in concentration per mg of mitochondrial protein were recorded in real-time

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 32

Rundle et al. 7

159

(continuous measurement). The ∆Ψmt was obtained by measuring TPP+ uptake by mitochondria using the

160

TPP+ ion selective electrode with respect to a reference electrode. Calibration with a known TPP+

161

concentration was done at the beginning of each experiment. Lastly, the H2O2 concentration and rate of

162

emission were measured using the Amplex red - HRP assay. Amplex red reacts with H2O2 in a 1:1

163

stoichiometric reaction giving rise to the fluorescent product resorufin, which is measured by the

164

Oroboros. A calibration curve was obtained using known concentrations (0-0.45 µM) of H2O2 at the

165

beginning of each experiment.

166

Isolated mitochondria (0.75 mg protein/mL) were first added to the Oroboros chambers, followed

167

by the substrates for ETS complex I, malate and glutamate. State 3 respiration, where ATP is actively

168

being produced, was initiated by the addition of ADP (6.25 µM). After all ADP was depleted,

169

mitochondria entered state 4 respiration, when no ATP production takes place. To investigate the effects

170

of NAs on different mitochondrial energetic states, dose response curves consisting of concentrations

171

ranging from 20 to 320 mg/L NA were obtained during state 3 and then again during state 4 respiration.

172

Here, NAs were added sequentially in 20 or 40 mg/L doses during each run to achieve predetermined

173

cumulative target concentrations. Sequential addition was required in order to complete a full dose

174

response will within the 4 h viability of mitochondria and to limit the number of live fish required. At

175

least three biological replicates for each state were performed.

176

Real-time Flow Cytometry Measurement of ∆Ψmt and ROS Emission

177

Real time ∆Ψmt and mitochondrial oxidation state (representing a portion of ROS emission)

178

measurements were done using fluorescent dyes JC-1 and H2DCFDA, respectively, as previously

179

described.27 Isolated mitochondria were added to MRB containing 5 mM of malate and 5 mM of

180

glutamate (state 2 respiration, high membrane potential, low respiration most comparable to state 4).

181

During ∆Ψmt measurement, mitochondria were incubated with JC-1 for 10 min prior to beginning of the

182

experiment. During ROS emission measurements, H2DCFDA was added approximately 200 s after

183

beginning the experiment. A dose-response assessment was conducted using the same NA concentrations

184

as in the Oroboros experiments. Doses of NA (20 or 40 mg/L) were added approximately 160 s apart and

ACS Paragon Plus Environment

Page 9 of 32

Environmental Science & Technology

Rundle et al. 8

185

a sample size of 3 was obtained for each parameter. Data were gated and presented using Flowjo 7.6.5

186

(TreeStar Inc., Oregon, USA). Three biological replicates for each parameter were performed.

187

Statistical Analysis

188

In order to obtain more consistent data for curve fitting, numbers were converted to a percentage

189

of the control, which was represented by state 3, state 4, or state 2 depending on the measurement. The

190

control measurement was used as 100% and subsequent measurements were converted to a percent by

191

dividing the measured value by the control and multiplying by 100. Dose-response curves were prepared

192

using Statistica v.13 (Dell, Oklahoma, USA) using a four-parameter logistic equation (three parameters

193

where the effect variable started at zero) using a Quasi-Newton estimation method in order to determine

194

the IC50 (inhibition concentration) or EC50 (effective concentration) for all endpoints. Curve fitting

195

employed a standard loss function, (observed-predicted)2. Estimates of IC10 or EC10 were calculated using

196

the equation of the curve. Curves were estimated for every replicate and the mean and 95% confidence

197

intervals calculated from the three replicates of each endpoint.

198 199

RESULTS AND DISCUSSION

200

Characterization of Naphthenic Acid Extract

201

Examination of NAs by negative ion mode direct injection HRMS revealed that our sample was

202

99% carboxylic acids with 90.4% of the ion intensity conforming to the NA structural formula,

203

CnH2n+ZO2. The most abundant ions had either 15 (30%) or 16 carbons (28%) and Z-values of -4 (40%) or

204

-6 (34%) (Figure 1). Other studies23,30 have observed similar carbon numbers and the Z values as found in

205

the sample examined herein. In the present study and previous studies, over 90% of compounds were in

206

the O2 category and fit the classical NA formula.23,30 This was expected as the extraction method used in

207

this study is specific to extracting amphipathic NAs that fit the classical NA formula.21 It is worth noting

208

that when comparing purified NAs with those directly extracted from the same tailings in a previous

209

study, the authors found that low and high molecular weight NAs tend to be depleted in the present

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 32

Rundle et al. 9

210

extraction process due to the multiple acid precipitation steps and the dichloromethane solvent extraction

211

steps, respectively.21,22

212

While the majority of negative ions in the extract fit within the classical NA formula, NA

213

compounds outside of the classical formula were also found. Of the ions identified, 2.0 % were O3

214

structures that are potentially carboxylic acids with an alcohol (or ether) functional group. A total of 7.0

215

% of the mixture was accounted for by O4 compounds that are likely dominated by compounds containing

216

two carboxylic acid groups. Based on this, the mixture is estimated to be 99.4% carboxylic acids at a

217

minimum, but may be higher due to traces of O5 compounds also identified in the mixture. Oxidized O3-,

218

O4-, and even O5-NAs have been reported in various OSPW mixtures.31,32 However, these compounds are

219

more water soluble than O2-NAs and therefore are less likely to precipitate out during extraction.22,23 A

220

study examining the same tailings water that was directly extracted using C18 methods found that O3 and

221

O4 compound accounted for 40% of the mixture.33

222

Elemental analysis of our extract showed that 70.4%, 18.4%, and 10.6% of the sample by mass

223

was C, O, and H, respectively. Nitrogen was not detectable in the mixture (< 0.5%), and sulfur was

224

present at 0.25% indicating a relatively pure hydrocarbon mixture. Accounting for the O2, O3, and O4

225

compounds, the elemental analysis suggested that the average carbon number was 11.14, which is lower

226

than the weighted average carbon number from the mass spectrometry analysis of 15.5. This may have

227

occurred if the sample was not completely dry at the time of analysis. An elemental analysis of seven NA

228

extracts from a previous study showed a range of 59-71% carbon, 15-20% oxygen, 6-9% hydrogen, 0-

229

0.7% nitrogen, and 1-6% sulfur.34 This analysis is considered consistent with our sample because previous

230

studies have documented that OSPW and NA extract samples from different areas have varying chemical

231

compositions and individual acid proportions.14,35 Overall, the results from all of the chemical analyses

232

suggest that the compound extracted in the present study is a NA mixture with > 99% purity this mixture

233

could be reliably used to estimate effects of OSPW-derived NAs on energy metabolism.

234

NAs Uncouple Oxidative Phosphorylation

ACS Paragon Plus Environment

Page 11 of 32

Environmental Science & Technology

Rundle et al. 10

235

An example of the effects of cumulative additions of NAs on mitochondrial oxygen consumption,

236

∆Ψmt, and ROS emission are shown in Figure 2, where tracings of each of the three parameters during a

237

representative experiment of state 4 respiration are overlaid. There is a period of instability for

238

approximately 1 min, particularly for the ROS endpoint as seen by large spikes, and all readings were

239

taken after the system returned to stability (2 min). The time of addition of mitochondria, necessary

240

substrates for respiration, and sequential doses of NAs are indicated. After the addition of substrates, the

241

mitochondria enter state 3 respiration, where ADP is being converted to ATP. The resultant movement of

242

protons by the ETS to drive ATP synthase results in an increase of mitochondrial respiration (red line) as

243

complex IV of the electron transport chain consumes oxygen to create proton motive force (blue line).

244

When all added ADP is converted to ATP, the mitochondria enters into state 4, or a resting state, which is

245

indicated by the reduction in respiration, and increase in ∆Ψm to a steady state. In state 4, the sequential

246

doses of NA reduce ∆Ψm and this reduction in ∆Ψm triggers an increase in respiration to compensate.

247

However, this was not consistently observed in the present study at higher doses of NA, there is an

248

ultimate decrease in respiration as the electron transport chain presumably ceases to function at a level

249

sufficient to sustain an increase respiration. Similarly, there is an increase in ROS emission (green lines,

250

large spikes are only a transient effect after NA addition), but at higher doses, ROS emissions also

251

decrease.

252

To demonstrate the functionality of the system to predict mitochondrial dysfunction and provide

253

examples of how effect at different levels of the ETS can be detected, a number of reference toxicants

254

were used, including oligomycin, an inhibitor of ATP synthase, carbonyl cayanide-4-

255

(trifluoromethoxy)phenylhydrazone (FCCP), a model uncoupler of OXPHOS, and antimycin-A, which

256

binds to cytochrome C and inhibits electron transport (Figure S1). Oligomycin completely eliminates

257

respiration (as ATP synthase is the only step that uses oxygen) with an increase in ∆Ψmt as ATP synthase

258

is not utilizing the proton motive force, and causes no change in the ROS emission. FCCP decreases ∆Ψmt

259

as it uncouples the proton gradient with an increase in respiration to compensate, and a slight increase in

260

ROS. Antimycin-A almost completely stops respiration as it impairs the ETS with only a slight drop in

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 32

Rundle et al. 11

261

∆Ψmt as the upstream complexes are still active. The subsequent back up of electrons causes substantive

262

ROS generation.

263

Dose-effect curves for all experiments combined show that successive additions of NAs resulted

264

in consistent decreases in respiration during state 3 with almost complete cessation of respiration at the

265

highest dose (ATP being produced; Figure 3a). When respiration is already low during state 4 there was a

266

trend towards a 25% increase in respiration, however the effect was not consistent in all replicates as

267

indicated by the high variability in the effect, the IC50 values were not calculated (ATP not being

268

produced; Figure 3b). NAs showed a trend to dissipate ∆Ψmt by less than 10% during state 3 when

269

membrane potential is already relatively low as compared to state 4 (ATP being produced; Figure 4a).

270

However, this effect was not consistent and IC50 values could not be calculated for all replicates. During

271

state 4, there was a consistent 20% drop in ∆Ψmt (ATP not being produced; Figure 4b). Real-time flow

272

cytometry measurement was used to validate these results and a stronger (50%) dose-dependent

273

dissipation of ∆Ψmt was observed (Figure 4c, d), though the effective concentrations using the Oroboros

274

and flow cytometry methods were similar (Table 1). The uncoupling of mitochondrial respiration and

275

dissipation of ∆Ψmt that were observed using two independent methods show that NAs can act as

276

mitochondrial uncouplers. The implications of sequential additions of NAs were evaluated by comparison

277

to the sequential addition of solvent controls (Figures S2 and S3) demonstrating the over time the system

278

remained relatively stable (< 5% drift) with sequential solvent additions. Furthermore, the effects of only

279

a single dose, equivalent to two sequential additions, of NA was investigated, and were found to have

280

similar effects as sequential doses (Figure S1).

281

The respective result for ∆Ψmt and respiration in state 3 and state 4 would be expected in the case

282

of a mitochondrial uncoupler that effects both the proton gradient, and the ETS. In state 3 respiration,

283

oxygen is being consumed at its maximum rate to maintain membrane potential to drive ATP synthesis,

284

the dissipation of ∆Ψmt alone should not result in a decrease in respiration. Thus, these results indicate that

285

the ETS is being inhibited independently of the reduction in ∆Ψmt. As there was almost complete

286

reduction of respiration, and only partial reduction of ∆Ψmt, this also suggests that the inhibition of the

ACS Paragon Plus Environment

Page 13 of 32

Environmental Science & Technology

Rundle et al. 12

287

ETS may be a more important mechanism. In contrast, in state 4, membrane potential is high, and oxygen

288

consumption is low, thus dissipation of the proton gradient by NAs causes an increased consumption of

289

oxygen that is indicative of the ETS attempting to restore the gradient. While disruption of the ETS may

290

be occurring in state 4, it does not become apparent until higher doses of NAs where respiration comes to

291

a plateau, and then appears to drop off (Figure 3b).

292

The observed uncoupling effect of NAs are likely due to the ability of NAs to act as

293

protonophores in which they are protonated in the intermembrane space and transported across the inner

294

mitochondrial membrane into the mitochondrial matrix where they release the proton in the alkaline

295

environment, thus dissipating the proton gradient (i.e., uncoupling OXPHOS).36,37 While and impact on

296

the ETS is evident from the decrease in state 3 respiration, the precise cause of this parallel effect is not

297

clear. Since there are multiple steps in the ETS, which are dependent on previous reactions, NAs could be

298

having inhibitory effects at any of these potential steps that drive the ETS, including direct interference

299

with ETS complexes.

300

NAs Stimulate ROS Emission

301

During both state 3 and 4 respiration dose responses using the Oroboros system, H2O2

302

concentration increased steadily for the first half of the dose range, followed by a general plateau or

303

decrease (Figure 5a, b). Real-time flow cytometry measurements were again used to validate the results

304

and a similar trend was observed (Figure 5c, d). The calculated EC50 and EC10 values were similar

305

between both energetic states and measurement techniques, ranging from 11.8 to 36.7 mg/L (Table 1).

306

The finding of stimulation of ROS emission by NAs measured by both Oroboros (H2O2 emissions) and

307

flow cytometry (mitochondrial oxidation state) is consistent with earlier reports that OSPW induced the

308

expression of genes related to oxidative stress in multiple species.38,39 The increased ROS emission

309

measured lends mechanistic support to an earlier report that correlated oxidative stress due to OSPW

310

exposure with developmental malformations observed in fathead minnow (Pimephales promelas)

311

embryos.40 Thus, mitochondria, a major source of cellular ROS, may be an important target for NAs.7,39

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 32

Rundle et al. 13

312

The underlying cause of increased ROS production by NAs is unknown. Generally, toxicants

313

have been shown to increase ROS production by interacting with the ETS via two different mechanisms.41

314

The first is by blocking electron transport thus increasing the reduction level of the ETS complexes

315

upstream of the inhibition.41 This results in increased leakage of electrons from the ETS complexes with

316

partial reduction of oxygen to form ROS.41,42 The second is that a toxicant may accept an electron from

317

the ETS complexes and transfer it to molecular oxygen. This mechanism stimulates ROS production

318

without inhibiting the ETS.41 We propose the first explanation is the most plausible mechanism of

319

increased ROS emission by NAs due to the observed inhibition of state 3 respiration (Figure 3a, 4a). The

320

inhibition of state 3 is indicated by the impairment of the ETS43, which is associated with an immediate

321

increase in ROS and dissipation of ∆Ψmt after addition of NAs. ROS emissions were the most sensitive

322

endpoints with EC50/IC50 values 3 to 6-fold lower than the other endpoints. In addition 95% confidence

323

intervals for the median effective concentration for state 3 respiration and ROS generation using the

324

Oroboros did not overlap, suggesting that ROS production can be detected well before effects on

325

OXPHOS are apparent. The ROS emission increased up to a NA concentration of ~ 160 mg/L on the

326

Oroboros system and ~ 120 mg/L during real-time flow cytometry, after which there was a drop in the

327

rate of ROS emission. The drop of ROS emission is likely due to impaired mitochondria function at

328

higher doses of NAs shown by almost complete inhibition of state 3 respiration by the 240 mg NA/L dose

329

(Figure 3a, 5a).

330

Environmental Relevance

331

Environmental concentrations of NAs are typically much lower than the EC50/IC50 values we

332

calculated, but EC10/IC10 estimates, particularly for ROS emission approach the high end of the range of

333

concentrations found in the receiving environment. Any comparison to environmental concentrations is

334

confounded by the many methods to determine such concentrations, including what is used as a standard,

335

as well as what the investigators consider to be NAs. In the present study, we were able to weigh

336

quantities of the NA mixture, which represented accurate concentrations given the mixture's high purity

337

as carboxylic acids. While our mixture consists of mostly O2-NAs, we believe toxicity of the source

ACS Paragon Plus Environment

Page 15 of 32

Environmental Science & Technology

Rundle et al. 14

338

material would be relatively unchanged based on another study from our laboratory that exposed isolated

339

mitochondria to that source material (unpublished data). Recent estimates of surface water NA

340

concentrations using Fourier-transform infrared (FTIR) spectroscopy methods in local lakes and

341

experimental ponds, including Beaver Creek, Demonstration Pond, Bison Pond, Mildred Lake, and

342

Horizon Lake, showed NA levels ranging from < 1 to 14 mg/L.35 The FTIR method overestimates

343

measured total carboxylic acids as compared to HMRS by a factor of two.33 A study conducted on

344

different sites, Muskeg River and Jackpine Mines, found surface water concentrations > 10 mg/L for O2

345

NAs and groundwater levels in the order of 5 mg/L using ultra performance liquid chromatography-time

346

of flight-mass spectrometry (UPLC-TOF-MS) methods. Although given different standard mixtures, it is

347

difficult to compare to HRMS methods compared by our laboratory.44 Levels of NAs in the Athabasca

348

River using liquid chromatography-high resolution mass spectrometry (LC-HRMS) techniques showed

349

concentration of total NA (O2-species only, thus an under estimate of all NAs) ranging from 0.7 to 52

350

µg/L.45 In general, concentrations in most affected waters are well within the EC10/IC10 values determined

351

in this study suggesting they may cause at least 10% changes in mitochondrial function.

352

The effective concentration of NAs experienced by aquatic organism can also be expected to be

353

dictated by the potential of NAs to bioconcentrate. Partitioning of oil sands organic compounds in

354

polydimethylsiloxane showed that O2-species have apparent KOW (pH 8.4) that ranging from negligible to

355

as high as 100 46; with KOW increasing with carbon number and decreasing with double bond equivalents.

356

However, as most NAs are present in their naphthenate form at environmental pH (> 7.5),

357

bioaccumulation is likely to be low. In a study that exposed rainbow trout to 3 mg/L at pH 8.2,

358

bioaccumulation factor (BAF) was reported as ~2.47 However, that study was done with commercially

359

available NAs, known to be rich in free fatty acids, unlike the aged samples utilized herein that may have

360

higher BAFs. However, yellow perch (Perca flavescens) exposed to 13 mg/L NA (by FTIR) aged tailings

361

water in an experimental pond were found to contain non-detectable levels of NAs (< 1 mg/kg) in muscle

362

tissue.48 Those fish accumulated and excreted NAs as indicated by biliary NA levels as high as 200

363

mg/L.48 There is support for bioaccumulation of NAs in fish from the Athabasca River (where levels of

ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 32

Rundle et al. 15

364

NAs are much lower) that had tissue concentrations ranging from non-detectable to as high as 2.8 mg/kg.

365

NAs tend to concentrate at the highest levels in gills and liver49, but data on compound-specific BAF in

366

fish are sparse and inconsistent. There has not been any work regarding the subcellular partitioning of

367

NAs and further work is needed to estimate and characterize the concentration of NAs that could be

368

present in the mitochondria of living organisms.

369

The question of how much reduction in ATP production would result in cellular or organismal

370

impacts is also critical to the environmental relevance of the present results. Cells must produce ATP at

371

the same rate they consume it. Therefore, a cell that is producing 10% less ATP than the cell demands can

372

deplete its ATP within two minutes.50 Typically, full or partial uncoupling of OXPHOS and the

373

dissipation of ∆Ψmt decrease or prevent the production of ATP. Cytotoxicity studies using mouse and

374

human cells have reported rapid (< 48 h) partial cell death occurs with 30% reduction in ATP levels51 and

375

greater than 50% cell death occurs when ATP levels decline by half.51,52 At the organismal level, this may

376

limit available energy for the organisms to carry out normal physiological processes, possibly resulting in,

377

or at least contributing to, the toxic effects reported in various studies. The NA mixture in the present

378

study was observed to cause complete mortality in threespine stickleback (Gasterosteus aculeatus) at 15

379

mg/L, comparable to EC10/IC10 values reported here (authors unpublished data), suggesting that energetic

380

mechanisms may play a factor in mortality. While caution must be taken when interpreting in vitro results

381

in vivo, the EC10/IC10 values determined herein may represent concentrations that are relevant to sub-

382

lethal physiological effects of NAs in trout and perhaps in aquatic organisms in general. As such, we

383

conclude that oxidative uncoupling is a likely mechanism of sub-lethal toxicity that occurs at

384

concentrations of NAs present in the oil sands region of Canada.

385 386

ACKNOWLEDGEMENTS

387

This research was funded through NSERC Discovery grants to MRV and CK. The work was made

388

possible through NSCER RTI and CFI grants to MRV and CK. The authors wish to thank Gillian

389

MacDonald, Dr. Rabin Bissesseur, and Dr. Hebelin Correa for their assistance in completing this work.

ACS Paragon Plus Environment

Page 17 of 32

Environmental Science & Technology

Rundle et al. 16

390 391

SUPPORTING INFORMATION

392

Raw Oroboros data traces for respiration, ROS, and membrane potential for positive controls oligomycin,

393

FCCP, and antimycin-A, solvent controls in state 3 respiration, solvent controls in state 4 respiration (S1-

394

S3).

395 396

REFERENCES

397 398

(1) About Oil Sands; http://www.energy.alberta.ca/OilSands/585.asp

399 400

(2) Allen, E.W. Process water treatment in Canada's oil sands industry: I. Target pollutants and treatment

401

objectives. J. Environ. Eng. Sci. 2008, 7, 123-138.

402 403

(3) Headley, J.V.; McMartin, D.W. A review of the occurence and fate of naphthenic acids in aquatic

404

environments. J. Environ. Sci. Health, Part A. 2004, 39 (8), 1989-2010.

405 406

(4) Clemente, J.S.; Fedorak, P.M. A review of the occurrence, analysis, toxicity, and biodegradation of

407

naphthenic acids. Chemosphere. 2005, 60, 585-600.

408 409

(5) Madill, R.E.A.; Orzechowski, M.T.; Chen, G.; Brownlee, B.G.; Bunce, N.J. Preliminary risk

410

assessment of the wet landscape option for reclamation of oil sands mine tailings: bioassays with mature

411

fine tailings pore water. Environ. Toxicol. 2001, 16 (3), 197-208

412 413

(6) Melvin, S.D.; Lanctôt, C.M.; Craig, P.M.; Moon, T.W.; Peru, K.M.; Headley, J.V.; Trudeau, V.L.

414

2013. Effects of naphthenic acid exposure on development and liver metabolic processes in anuran

415

tadpoles. Environ. Pollut. 2013, 177, 22-27.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 32

Rundle et al. 17

416 417

(7) Wiseman, S.B.; Anderson, J.C.; Liber, K.; Giesy, J.P. Endocrine disruption and oxidative stress in

418

larvae of Chironomus dilutus following short-term exposure to fresh and aged oil sands process-affected

419

water. Aquat. Toxicol. 2013a, 142-143, 414-421.

420 421

(8) Kavanagh, R.J.; Frank, R.A.; Solomon, K.R.; Van Der Kraak, G. Reproductive and health assessment

422

of fathead minnows (Pimphales promelas) inhabiting a pond containing oil sands process-affected water.

423

Aquat. Toxicol. 2013, 130-131, 201-209.

424 425

(9) Marentette, J.R.; Frank, R.A.; Bartlett, A.J.; Gillis, P.L.; Hewitt, L.M.; Peru, K.M.; Headley, J.V.;

426

Brunswick, P.; Shang, D.; Parrorr, J.L. Toxicity of naphthenic acid fraction components extracted from

427

fresh and aged oil sands process-affected waters, and commercial naphthenic acid mixtures, to fathead

428

minnow (Pimephales promelas) embryos. Aquat. Toxicol. 2015, 164, 108-117.

429 430

(10) Mohseni, P.; Hahn, N.A.; Frank, R.A.; Hewitt, L.M.; Hajibabaei, M.; Van Der Kraak, G. Naphthenic

431

acid mixtures from oil sands process-affected water enhance differentiation of mouse embryonic stem

432

cells and affect development of the heart. Environ. Sci. Technol. 2015, 49, 10165-10172.

433

(11) Hogan, N.S., K.L. Thorpe, and M.R. van den Heuvel. Opportunistic disease in yellow perch in

434

response to decadal changes in the chemistry of oil sands-affected waters. Environ. Poll. 2018, 234, 769-

435

778.

436

(12) Frank, R.A.; Kavanagh, R.; Burnison, B.K.; Arsenault, G.; Headley, J.V.; Peru, K.M.; Van Der

437

Kraak, G.; Solomon, K.R. Toxicity assessment of collected fractions from an extracted naphthenic acid

438

mixture. Chemosphere. 2008, 72, 1309-1314.

439

ACS Paragon Plus Environment

Page 19 of 32

Environmental Science & Technology

Rundle et al. 18

440

(13) He, Y.; Wiseman, S.B.; Zhang, X.; Hecker, M.; Jones, P.D.; Gamal El-Din, M.; Martin, J.W.; Giesy,

441

J.P. Ozonation attenuates the steroidogenic disruptive effects of sediment free oil sands process water in

442

the H295R cell line. Chemosphere. 2010, 80, 578-584

443 444

(14) Reinardy, H.C.; Scarlett, A.G.; Henry, T.B.; West, C.E.; Hewitt, L.M.; Frank, R.A.; Rowland, S.J.

445

Aromatic naphthenic acids in oil sands process-affected water, resolved by GCxGC-MS, only weakly

446

induce the gene for vitellogenin production in zebrafish (Danio rerio) larvae. Environ. Sci. Technol. 2013,

447

47, 6614-6620.

448 449

(15) Rabergh, C.M.I.; Isomaa, B.; Eriksson, J.E. The resin acids dehydroabietic acid and isopimaric acid

450

inhibit bile acid uptake and perturb potassium transport in isolated hepatocytes from rainbow trout

451

(Oncorhynchus mykiss). Aquat. Toxicol. 1992, 23, 169-180.

452 453

(16) Rissanen, E.; Krumschnabel, G.; Nikinmaa, M. Dehydroabietic acid, a major component of wood

454

industry effluents, interferes with cellular energetics in rainbow trout hepatocytes. Aquat. Toxicol. 2003,

455

62, 45-53.

456 457

(17) Terada, H. Uncouplers of oxidative phosphorylation. Environ. Health Perspect. 1990, 87, 213-218.

458 459

(18) Boveris, A., Chance, B. The mitochondrial generation of hydrogen peroxide: general properties and

460

effect of hyperbaric oxygen. Biochem. J. 1973, 134, 707-716.

461 462

(19) Collins, Y.; Chouchani, E.T.; James, A.M.; Menger, K.E.; Comchemé, H.M.; Murphy, M.P.

463

Mitochondrial redox signalling at a glance. J. Cell Sci. 2012, 125, 801-806.

464

ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 32

Rundle et al. 19

465

(20) Kowaltowski, A.J.; Vercesi, A.E. Mitochondrial damage induced by conditions of oxidative stress.

466

Free Radical Biol. Med. 1999, 26 (3/4), 463-471.

467 468

(21) Frank, R.A.; Kavanagh, R.; Burnison, B.K.; Headley, J.V.; Peru, K.M.; Ven Der Kraak, G.;

469

Solomon, K.R. Diethylaminoethyl-cellulose clean-up of a large volume naphthenic acid extract.

470

Chemosphere. 2006, 64, 1346-1352.

471 472

(22) Macdonald, G.Z., Hogan, N.S., Köllner, Thorpe, K.L., Phalen, L.J., Wagner, B.D., van den Heuvel,

473

M.R. 2013. Immunotoxic effects of oil sands-derived naphthenic acids to rainbow trout. Aquat. Toxicol.

474

2013, 126, 95-103.

475 476

(23) Leclair, L.A.; Pohler, L.; Wiseman, S.B.; He, Y.; Arens, C.J.; Giesy, J.P.; Scully, S.; Wagner, B.D.;

477

van den Heuvel, M.R., Hogan, N.S. In vitro assessment of endocrine disrupting potential of naphthenic

478

acid fractions derived from oil sands-influenced waters. Environ. Sci. Technol. 2015, 49, 5743-5752.

479 480

(24) ASTM D4239-83, Sulfur in the Analysis Sample of Coal and Coke Using High Temperature Tube

481

Furnace Combustion Methods, Method C. In: Annual Book of ASTM Methods, Vol. 05.05; ASTM 1992.

482 483

(25) Bance, S. Handbook of Practical Organic Microanalysis; Ellis Horwood, Ltd. 1980.

484 485

(26) Thermo Fischer Scientific, FlashEA™ 1112 Elemental Analyzer Operating Manual.

486

2016

487 488

(27) Sharaf, M.S.; van den Heuvel, M.R.; Stevens, D.; Kamunde, C. Zinc and calcium modulate

489

mitochondrial redox state and morphofunctional integrity. Free Radical Biol. Med. 2015, 84, 142-153.

490

ACS Paragon Plus Environment

Page 21 of 32

Environmental Science & Technology

Rundle et al. 20

491

(28) Bradford, M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein

492

utilizing the principle of protein–dye binding. Anal. Biochem. 1976, 72, 248–254.

493 494

(29) Sharaf, M.S., Stevens, D., Kamunde, C. Mitochondrial transition ROS spike (mTRS) results from

495

coordinated activities of complex I and nicotinamide nucleotide trans- hydrogenase. BBA - Bioenergetics.

496

2017, 1858, 955-965.

497

(30) Headley, J.V.; Peru, K.M.; Fahlman, B.; Colodey, A.; McMartin, D.W. Selective solvent extraction

498

and characterization of the acid extractable fraction of Athabasca oils sands process waters by Orbitrap

499

mass spectrometry. Int. J. Mass Spectrom. 2013, 345-347, 104-108.

500 501

(31) Huang, R.; McPhedran, K.N.; Gamal El-Din, M. Ultra performance liquid chromatography ion

502

mobility time-of-flight mass spectrometry characterization of naphthenic acids species from oil sands

503

process-affected water. Environ. Sci. Technol. 2015, 49, 11737-11745.

504 505

(32) Wang, C.; Huang, R.; Klamerth, N.; Chelme-Ayala, P.; Gamal El-Din, M. Positive and negative

506

electrospray ionization analyses of the organic fractions in raw and oxidized oil sands process-affected

507

water. Chemosphere. 2016, 165, 239-247.

508 509

(33) Leclair, L.A.; MacDonald, G.Z.; Phalen, L.J.; Köllner, B.; Hogan, N.S.; van den Heuvel, M.R. The

510

immunological effects of oil sands surface waters and naphthenic acids on rainbow trout (Oncorhynchus

511

mykiss). Aquat. Toxicol. 2013, 142-143, 185-194.

512 513

(34) Grewer, D.M.; Young, R.F.; Whittal, R.M.; Fedorak, P.M. Naphthenic acids and other acid-

514

extractables in water samples from Alberta: What is being measured? Sci. Total Environ. 2010, 408,

515

5997-6010.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 32

Rundle et al. 21

516 517

(35) Arens, C.J.; Hogan, N.S.; Kavanagh, R.J.; Mercer, A.G.; Van Der Kraak, G.J.; van den Heuvel, M.R.

518

Sublethal effects of aged oil sands-affect water on white sucker (Catostomus commersonii). Environ.

519

Toxicol. Chem. 2015, 34(3), 589-599.

520 521

(36) Kadenbach, B. Intrinsic and extrinsic uncoupling of oxidative phosphorylation. Biochim. Biophys.

522

Acta. 2003, 1604, 77-94.

523 524

(37) Rottenberg, H.; Hashimoto, K. Fatty acid uncoupling of oxidative phosphorylation in rat liver

525

mitochondria. Biochemistry. 1986, 25, 1747-1755.

526 527

(38) Gagné, F.; Douville, M.; André, C.; Debenest, T.; Talbot, A.; Sherry, J.; Hewitt, L.M.; Frank, R.A.;

528

McMaster, M.E.; Parrott, J.; Bickerton, G. Differential changes in gene expression in rainbow trout

529

hepatocytes exposed to extracts of oil sands process-affected water and the Athabasca River. Comp.

530

Biochem. Physiol., Part C. 2012, 155, 551-559.

531 532

(39) Wiseman, S.B.; He, Y.; Gamal-El Din, M.; Martin, J.W.; Jones, P.D.; Hecker, M.; Giesy, J.P.

533

Transcriptional responses of male fathead minnows exposed to oil sands process-affected water. Comp.

534

Biochem. Physiol. C: Pharmacol. Toxicol. 2013b, 157, 227-235.

535 536

(40) He, Y.; Patterson, S.; Wang, N.; Hecker, M.; Martin, J.W.; Gamal El-Din, M.; Giesy, J.P.; Wiseman,

537

S.B. Toxicity of untreated and ozone-treated oil sands process-affected water (OSPW) to early life stages

538

of the fathead minnow (Pimephales promelas). Water Res. 2012, 46, 6359-6368.

539 540

(41) Turrens, J.F. Mitochondrial formation of reactive oxygen species. J. Physiol. 2003, 552 (2), 335-344.

541

ACS Paragon Plus Environment

Page 23 of 32

Environmental Science & Technology

Rundle et al. 22

542

(42) Barrientos, A.; Moraes, C.T. Titrating the effects of mitochondrial complex I impairment in the cell

543

physiology. J. Biol. Chem. 1999, 23 (4), 16188-16197.

544 545

(43) Kalbácová, M.; Vrbacky, M.; Drahota, Z.; Melkova, Z. Comparison of the effect of mitochondrial

546

inhibitors on mitochondrial membrane potential in two different cell lines using flow cytometry and

547

spectrofluorometry. Cytometry Part A. 2003, 52A, 110-116.

548

(44) Huang, R.; Chen, Y.; Meshref, M.N.A.; Chelme-Ayala, P.; Dong, S.; Ibrahim, M.D.; Wang, C.;

549

Klamerth, N.; Hughes, S.A.; Headley, J.V.; Peru, K.M.; Brown, C.; Mahaffey, A.; Gamal El-Din, M.

550

Characterization and determination of naphthenic acids species in oil sands process-affected water and

551

groundwater from oil sands development area of Alberta, Canada. Water Res. 2018, 128, 29-137.

552 553

(45) Arens, C.J.; Arens, J.C.; Hogan, N.A.; Kavanagh, R.J.; Berrue, F.; Van Der Kraak, G.J.; van den

554

Heuvel, M.R. Population impacts in white sucker (Catostomus commersonii) exposed to oil sands-derived

555

contaminants in the Athabasca River. Environ. Tox. Chem. 2017, 36 (8), 2058-2067.

556 557

(46) Zhang, K.; Pereira, A.S.; Martin, J.W. Partitioning for thousands of dissolved organic species in oil

558

sands process-affected water. Environ. Sci. Technol. 2015, 49, 8907-8913.

559 560

(47) Young, R.F; Wismer, W.V.; Fedorak, P.M. Estimating naphthenic acids concentrations in

561

laboratory-exposed fish and in fish from the wild. Chemosphere. 2008, 73, 498–505.

562 563

(48) van den Heuvel, M.R.; Hogan, N.S.; MacDonald, G.Z.; Berrue, F.; Young, R.F.; Arens, C.J.; Kerr,

564

R.G.; Fedorak, P.M. Assessing accumulation and biliary excretion of naphthenic acids in yellow perch

565

exposed to oil sands-affected waters. Chemosphere. 2014, 95, 691-627.

566

ACS Paragon Plus Environment

Environmental Science & Technology

Page 24 of 32

Rundle et al. 23

567

(49) Young, R.F.; Martínez Michael, L.; Fedorak, P.M. Distribution of naphthenic acids in tissues of

568

laboratory-exposed fish and in wild fishes from near the Athabasca oil sands in Alberta, Canada.

569

Ecotoxicol. Environ. Saf. 2011, 74, 889–896.

570 571

(50) Moyes, C.D.; Schulte, P.M. Principles of Animal Physiology, 2nd ed.; Pearson Education Inc.: San

572

Fransico, U.S., 2008

573 574

(51) Lieberthal, W.; Menza, S.A.; Levine, J.S. Graded ATP depletion can cause necrosis or apoptosis of

575

cultured mouse proximal tubular cells. Am. J. Renal Physiol. 1998, 274 (2), F315-F327.

576

(52) Rogers, C.; Davis, B.; Neufer, P.D.; Murphy, M.P.; Anderson, E.J.; Robidoux, J. A transient increase

577

in lipid peroxidation primes preadipocytes for delayed mitochondrial inner membrane permeabilization

578

and ATP depletion during prolonged exposure to fatty acids. Free Radic. Biol. Med. 2014, 67, 330-341.

579 580 581 582 583 584 585

ACS Paragon Plus Environment

Page 25 of 32

Environmental Science & Technology

Rundle et al. 24

586

Table 1. IC50 or EC50 and IC10 or EC10 values mg/L (95% Confidence Interval) calculated for each

587

parameter based on Oroboros and flow cytometry measurements. Curve fits were conducted on individual

588

replicates and parameters averaged. The mean variance explained by the curve fits for each parameter are

589

indicated. Parameters for which dose-effect curves could not be fitted for all replicates are indicated as no

590

consistent effect.

591 Parameter

Mean EC50/IC50

Mean EC10/IC10

Mean % Variance

(mg/L)

(mg/L)

Explained

Oroboros State 3

104.7 (34.2-175.1)

35.1 (13.9-56.2)

99.0

Oroboros State 4

no consistent effect

no consistent effect

Oroboros State 3

no consistent effect

no consistent effect

Oroboros State 4

157.8 (58.4-257.2)

66.7 (25.5-107.9)

99.2

Flow Cytometry State 2

88.7 (61.5-115.8)

29.4 (0-85.5)

96.7

Oroboros State 3

21.0 (15.0-27.0

15.0 (1.1-29.0)

88.5

Oroboros State 4

36.7 (1-109.8)

24.1 (1.2-47.1)

87.0

25.6 (12.4-38.7)

11.8 (0.4-23.2)

66.8

Oxygen Consumption

Mitochondrial Membrane Potential

Reactive Oxygen Species

Flow Cytometry State 2

592

ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 32

Rundle et al. 25

593

Figure 1. Relative abundance of molecular ions present in naphthenic acid extract obtained using negative

594

ion direct injection HRMS. The Z-value represents the hydrogen deficiency of the compound. The

595

absolute value of Z divided by 2 gives the number of double bonds or rings in the compounds. Each Z-

596

value reported in the extract is indicated with a different colour.

597 598

Figure 2. Representative Oroboros tracing of a state 4 respiration dose response. This tracing includes all

599

three parameters, oxygen consumption (red line), mitochondrial membrane potential (∆Ψmt; blue line),

600

and ROS represented by H2O2 production (green line) overlaid. The units for H2O2 production are

601

pmol/mg protein/min and are represented in the right axis, while the ∆Ψmt (mV) are arbitrary for the

602

purposes of this figure. The transient spikes of H2O2 production (green line) are an artifact of addition of

603

NAs or substrates to the experiment chamber. At the beginning of the experiment, isolated mitochondria

604

and the necessary substrates malate and glutamate (MG) and ADP were added as shown with the first

605

three arrows. During the appropriate respiration state, NAs were added sequentially 3 minutes apart to

606

achieve a cumulative concentration of NAs, which is indicated by arrows (eg., NA20 indicates 20 mg/L

607

NAs). All data for analysis were collected once steady states were reached.

608 609

Figure 3. (a) Dose-effect curve for oxygen consumption during state 3 respiration using the Oroboros

610

system with the calculated IC50. (b) Dose-effect curve for oxygen consumption during state 4 respiration

611

using the Oroboros system (EC50 could not be consistently calculated for state 4). Error bars are standard

612

error of the mean. N=3 for each point.

613 614

Figure 4. (a) Dose-effect curve for mitochondrial membrane potential (∆Ψmt) during state 3 respiration

615

measured by the Oroboros system (IC50 could not be consistently calculated for state 3). (b) Dose-effect

616

curve for ∆Ψmt during state 4 respiration measured by the Oroboros system during state 4 respiration with

617

calculated IC50. (c) Dose response curve for ∆Ψmt during state 2 respiration measured by real-time flow

618

cytometry with calculated IC50. (d) Representative tracing for the effects of different doses of NA on ∆Ψmt

ACS Paragon Plus Environment

Page 27 of 32

Environmental Science & Technology

Rundle et al. 26

619

during state 2 respiration measured by real-time flow cytometry where the x-axis is time and the y-axis is

620

the amount of fluorescence. The colour represents the density of the distribution of fluorescence. The

621

additions of the NA doses are indicated on the figure with arrows and are labeled to each corresponding

622

dose. Error bars are standard error of the mean. N=3 for each point.

623 624

Figure 5. (a) Overall dose-effect curve for H2O2 production during state 3 respiration measured by the

625

Oroboros system with calculated EC50. (b) Overall dose-effect curve for H2O2 production during state 4

626

respiration measured by the Oroboros system with calculated EC50. (c) Dose-effect curve for

627

mitochondria ROS during state 2 respiration measured real-time flow cytometry with calculated EC50.

628

EC50 valcues were calculated from individual replicates and then averaged. (d) Representative tracing for

629

the effects of different doses of NA on ROS emission during state 2 respiration measured by real-time

630

flow cytometry where the x-axis is time and the y-axis is the amount of fluorescence. The colour

631

represents the density of the distribution of fluorescence. The additions of the NA doses are indicated on

632

the figure with arrows and are labeled to each corresponding dose. Error bars are standard error of the

633

mean. N=3 for each point.

634

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 1 692x563mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

Environmental Science & Technology

Figure 2 1058x595mm (96 x 96 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 3 1044x404mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

Environmental Science & Technology

Figure 4 789x578mm (96 x 96 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 5 901x581mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 32 of 32