Athabasca Oil Sands Petcoke Extract Elicits ... - ACS Publications

Apr 28, 2017 - Heterocyclic Aromatics in Petroleum Coke, Snow, Lake Sediments, and Air Samples from the Athabasca Oil Sands Region. Environmental ...
1 downloads 0 Views 900KB Size
Subscriber access provided by Eastern Michigan University | Bruce T. Halle Library

Article

Athabasca oil sands petcoke extract elicits biochemical and transcriptomic effects in avian hepatocytes Doug Crump, Kim Williams, Suzanne Chiu, Yifeng Zhang, and Jonathan W. Martin Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 28 Apr 2017 Downloaded from http://pubs.acs.org on May 3, 2017

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 35

Environmental Science & Technology

1

Athabasca oil sands petcoke extract elicits biochemical and transcriptomic effects in avian

2

hepatocytes

3 4

Doug Crumpa*, Kim L. Williamsa, Suzanne Chiua, Yifeng Zhangb, Jonathan W. Martinb

5 6

a

7

National Wildlife Research Centre, Carleton University, Ottawa, ON, Canada K1A 0H3

8

b

9

Canada, T6G 2G3

Ecotoxicology and Wildlife Health Division, Environment and Climate Change Canada,

Department of Laboratory Medicine and Pathology, University of Alberta, Edmonton, Alberta,

10 11

*To whom correspondence should be addressed:

12

Doug Crump, Environment and Climate Change Canada, National Wildlife Research Centre,

13

1125 Colonel By Drive, Raven Road, Ottawa, Ontario, Canada K1A 0H3 or K1S 5B6 (courier).

14

Fax: +1 (613) 998-0458. E-mail: [email protected]

15 16 17 18 19

Keywords: petroleum coke; polycyclic aromatic compounds; in vitro screening; EROD; avian;

20

PCR array; mRNA expression

1 ACS Paragon Plus Environment

Environmental Science & Technology

21 22

Page 2 of 35

Abstract Petroleum coke or “petcoke” is a granular carbonaceous material produced during the

23

upgrading of heavy crude oils, including bitumen. Petcoke dust was recently reported as an

24

environmental contaminant in the Athabasca oil sands region, but the ecotoxicological hazards

25

posed by this complex bitumen-derived material – including those to avian species – have not

26

been characterized. In this study, solvent extracts (x) of delayed and fluid petcoke (xDP and xFP)

27

were prepared and dissolved in dimethylsulfoxide. A water accommodated fraction of delayed

28

petcoke (waDP) was also prepared. Graded concentrations of xDP, xFP, and waDP were

29

administered to chicken and double-crested cormorant hepatocytes to determine effects on 7-

30

ethoxyresorufin-O-deethylase (EROD) activity, porphyrin accumulation, and mRNA expression.

31

Polycyclic aromatic compounds (PACs) were characterized and xDP, xFP, and waDP had total

32

PAC concentrations of 93000, 270, and 5.3 ng/ml. The rank order of biochemical and

33

transcriptomic responses was xDP>xFP>waDP (e.g. EROD EC50s were lower for xDP compared

34

to xFP and waDP). A total of 22, 18, and 4 genes were altered following exposure to the highest

35

concentrations of xDP, xFP, and waDP, respectively, using a chicken PCR array comprising 27

36

AhR-related genes. To provide more exhaustive coverage of potential toxicity pathways being

37

impacted, two avian ToxChip PCR arrays – chicken and double-crested cormorant – were

38

utilized and xDP altered the expression of more genes than xFP. Traditional PAC-related toxicity

39

pathways and novel mechanisms of action were identified in two avian species following

40

petcoke extract exposure. Extrapolation to real world exposure scenarios must consider the

41

bioavailability of the extracted PACs compared to those in exposed organisms.

42 43 44 2 ACS Paragon Plus Environment

Page 3 of 35

45

Environmental Science & Technology

Introduction

46

Petroleum coke or “petcoke” is a by-product of petroleum upgrading. The process

47

involves thermally cracking the complex hydrocarbons and separating the non-volatile and heavy

48

hydrocarbons (petcoke) from the lighter and valuable hydrocarbons used for synthetic crude oil.

49

In the Athabasca oil sands region, two such coking methods are used to produce marketable

50

synthetic crude oil – delayed coking and continuous fluid coking.1,2 In Alberta, approximately

51

15% of extracted bitumen is converted to petcoke, and this material cannot be used for fuel in

52

most markets due to its high sulfur content.3 This granular, carbonaceous, hydrophobic coal-like

53

industrial by-product can be stored in large piles near processing facilities where it is susceptible

54

to wind erosion; a potential vector for widespread environmental and human exposure to petcoke

55

dust. Zhang et al.1 recently reported that delayed petcoke was a major source of polycyclic

56

aromatic hydrocarbon (PAH) deposition to the environment around oil sands development. PAH

57

concentrations of 25-28 µg/g dry weight were reported in authentic delayed petcoke from the

58

Athabasca oil sands region.1 Moreover, metals, including nickel and vanadium, have been

59

measured in petcoke at concentrations greater than 100 ppm.1,4 PAHs are known mutagens and

60

carcinogens that are included on the Toxic Substances List (Schedule 1) of the Canadian

61

Environmental Protection Act and can elicit biochemical and molecular-level effects via aryl

62

hydrocarbon receptor (AhR) activation.5,6

63

Although PAHs are classified as toxic substances, the screening-level hazard

64

characterization of petcoke (both green [CASRN 64741-79-3] and calcined [CASRN 64743-05-

65

1]) conducted as part of the US Environmental Protection Agency’s (EPA) High Production

66

Volume Challenge Program did not classify petcoke as a hazardous waste material based on the

67

analysis of potential human health and environmental impacts.7 In fact, most toxicity analyses

68

reported a low potential of petcoke to cause adverse effects on aquatic or terrestrial environments 3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 35

69

as well as a low health hazard potential in humans, with no observed carcinogenic, reproductive,

70

or developmental effects. In rats and monkeys exposed via inhalation to dust aerosol

71

concentrations of delayed process petcoke (10.2 and 30.7 mg/m3) for 5 days/week over 2 years,

72

Klonne et al. reported no significant adverse effects with the exception of discolouration of lungs

73

(due to deposition of petcoke) and increased lung weight; no mortality was observed up to the

74

highest concentration tested.8 Klonne et al. also conducted a nose-only inhalation study with

75

micronized petcoke for 52 days to assess reproductive/developmental toxicity and reported no

76

effects up to the highest concentration tested; ~0.3 mg/L.8 No mortality occurred in fathead

77

minnows (Pimephales promelas) exposed to water accommodated petcoke fractions under semi-

78

static conditions for 96 h up to a maximum loading rate of 1000 mg/L; however, petcoke has

79

negligible water solubility and attempts to measure unalkylated PAHs, metals and sulphur

80

revealed that concentrations were below detection limits.7 However, it is not clear if these former

81

studies are relevant to petcoke produced in the Canadian oil sands region. Furthermore, to the

82

best of our knowledge, no studies have assessed the potential toxic effects of petcoke in birds,

83

which could be exposed, for example, via the diet or by inhalation of airborne petcoke dust

84

particles. In general, more comprehensive toxicity tests such as in vivo and in vitro assays using

85

birds, fish, invertebrates, and mammals are warranted to assess the toxicity of petcoke in

86

complex organisms, cell lines, and/or primary cell cultures.

87

The objectives of the present study were to: 1) prepare solvent extracts from delayed and

88

fluid petcoke (xDP and xFP) and a water accommodated fraction of delayed petcoke (waDP); 2)

89

characterize the concentration of a suite of polycyclic aromatic compounds (PACs) in xDP, xFP

90

and waDP; and 3) utilize a well-established, high-throughput avian in vitro assay (i.e. primary

91

embryonic hepatocyte assay), two well-characterized biochemical assays, and custom-designed

4 ACS Paragon Plus Environment

Page 5 of 35

Environmental Science & Technology

92

avian PCR arrays to elucidate biochemical and transcriptomic effects of petcoke extract exposure

93

in both a domestic and wild avian species (i.e. chicken and double-crested cormorant). The

94

inclusion of two avian species permits the comparison of end points in a lab model species that is

95

typically considered highly sensitive to a wide range of environmental pollutants, with a wild

96

species that could realistically be exposed to petcoke in the environment. Ultimately, the goal

97

was to determine if subtle biochemical and molecular biological end points could be used to rank

98

the petcoke extracts in a manner that corresponded to the actual PAC concentrations and

99

determine predicted (i.e. AhR-mediated end points) and novel mechanisms of petcoke toxicity in

100

avian cells.

101 102

Experimental section

103

Sample preparation for contaminant analysis and administration to hepatocytes

104

For the DP and FP extract, 10 g of authentic DP and FP solids from the Athabasca oil

105

sands region in Canada were extracted twice with 250 mL Optima-grade dichloromethane

106

(DCM) and then with 250 mL Optima-grade methanol (MeOH). The combined solvent was

107

filtered through a GB-140 glass fiber membrane filter (90 mm diameter, pore size = 0.4 µm).

108

Rotary evaporation was used to concentrate the extract, which was subsequently transferred to a

109

glass vial and blown to dryness under nitrogen gas. The final dry weights of pure solvent-

110

extracted material from 10 g DP and FP were 79.8 mg and 1.8 mg, respectively.

111

Given that the mass of total material extracted was so different between DP and FP,

112

concentrations of test materials for toxicity testing are reported here in two ways. First, to best

113

compare the toxic potency of the pure solvent-extracted material, concentrations are reported as

114

mg extract/volume DMSO. Stock solutions were prepared for both xDP and xFP at nominal

5 ACS Paragon Plus Environment

Environmental Science & Technology

115

concentrations of 0.01 mg extract/µL dimethyl sulfoxide (DMSO). Specifically, for delayed

116

petcoke, the 0.01 mg/µL stock was prepared by dissolving 3.9 mg of pure DP organic extract in

117

378 µL DMSO yielding a final concentration of 0.0103 mg/µL (Table 1). For fluid petcoke, the

118

entire extract (1.8 mg) was dissolved in 250 µL of DMSO, for a final concentration of 0.007

119

mg/µL (Table 1). The second method of reporting the concentrations is based on the equivalent

120

mass of original petcoke from which the solubilized test material originated. Thus, based on the

121

above stock solutions, the highest corresponding xDP equivalent concentration was 0.00126 g

122

DP/µL DMSO, and the highest corresponding xFP equivalent concentration was 0.04 g FP/µL

123

DMSO (Table 1).

Page 6 of 35

124

The water accommodated fraction (waDP) was prepared from solid DP (1 gram) that was

125

stirred in 10 mL Optima LC/MS water for 10 days at 23°C and then filtered as above. To remove

126

suspended particulate, this mixture was centrifuged at 15800 x g, transferred into a glass vial,

127

and serial dilutions were made from this stock for hepatocyte administration in Optima-grade

128

LC/MS water. Concentrations for toxicity testing are based on the total loading of DP in water,

129

thus, the highest equivalent concentration for waDP was 0.0001 g DP/µL water (1 g/10 000 µL;

130

Table 1).

131 132 133

Determination of PAC concentrations The polycyclic aromatic compounds that were analyzed are listed in the Supplemental

134

Information along with the GC-MS parameters. Briefly, 16 parent PAH standards (unlabelled)

135

and 16 internal standards (deuterium-labelled) were purchased from Wellington Laboratories

136

(Guelph, Canada) and the 4 alkyl-PAH and dibenzothiopene (DBT) standards were purchased

137

from Chiron (Norway). Accelerated solvent extraction was used for targeted analysis of xDP

6 ACS Paragon Plus Environment

Page 7 of 35

Environmental Science & Technology

138

and xFP in DMSO, while a liquid-liquid extraction method was used for analysis of waDP,

139

described in greater detail in the Supplemental Information. Solid phase extraction (SPE) was

140

used to clean up all three sample extracts prior to concentrating them in 200 µL of hexane for

141

quantitative GC-MS analysis. GC-MS analysis involved splitless injection of 1 µL from all

142

extracts. Chromatographic separation was achieved on a DB-5MS column (Agilent; 20 m x 0.18

143

mm x 0.18 µm). Additional details regarding the method and method detection limits (MDLs) for

144

all individual analytes are available in the Supplemental Information.

145 146 147

Preparation and dosing of chicken and double-crested cormorant hepatocytes Twenty-five fertilized, unincubated white leghorn chicken (Gallus domesticus) eggs were

148

obtained from the Canadian Food Inspection Agency and 15 double-crested cormorant (DCCO)

149

eggs were collected from nests containing one egg on May 20, 2016, at Doucet Rock, Lake

150

Huron (46° 8'22.09"N, 82°51'7.35"W). Eggs of both species were artificially incubated until ~2

151

days pre-hatch8,10,11 at which point, embryos were euthanized by decapitation and hepatocytes

152

were prepared by collagenase digestion and filtration and distributed into 48-well plates at ~40

153

µg/well, as described previously.9,10

154

Following incubation for 24 hours at 37.5°C and 5% CO2, hepatocytes were treated with

155

either the DMSO or water vehicle control (0.5% v/v) and serial dilutions of xDP, xFP or waDP

156

(n=3 wells/treatment for 7-ethoxyresorufin-O-deethylase [EROD], and n=6 wells/treatment for

157

PCR array). At the highest concentrations, final administered amounts (in g petcoke equivalents)

158

were 0.00315 g xDP, 0.1 g xFP and 0.00025 g waDP (Table 1). Serial dilutions covering 7 orders

159

of magnitude (1x10-3- 10-9 g petcoke equivalent/µL for xDP, 4x10-2-10-8 g petcoke equivalent/µL

160

for xFP, and 1x10-4-10-10 g petcoke equivalent/µL for waDP) and a 2,3,7,8-tetrachloro-dibenzo-

7 ACS Paragon Plus Environment

Environmental Science & Technology

161

p-dioxin (TCDD) positive control (100 nM)10 were administered to chicken embryonic

162

hepatocytes (CEH) for EROD activity determination. For PCR array analysis, serial dilutions

163

covering 6 orders of magnitude were used for CEH, and two orders of magnitude for double-

164

crested cormorant embryonic hepatocytes (DCEH). The number of concentrations and petcoke

165

extracts tested in DCEH (i.e. waDP was not screened in DCEH due to the minimal response

166

observed for EROD activity and gene expression in CEH) was lower due to the limited pool of

167

cormorant hepatocytes available. The limited pool of DCEH also meant that the EROD and

168

porphyrin assays were only conducted in chicken cells. Hepatocytes of both species were

169

incubated for 24 h following extract administration until processing for subsequent RNA

170

isolation, EROD/porphyrin determination or cell viability.10,11

Page 8 of 35

171 172

Cell viability determination

173

Viability of CEH following exposure to concentrations of xDP, xFP, and waDP covering

174

4 orders of magnitude from the most concentrated solution was estimated with the ViaLight Plus

175

kit (Lonza), according to the manufacturer’s instructions. DCEH viability was not determined

176

due to the limited pool of cells. Untreated cells were used as a negative control and cells dosed

177

with a final, in-well concentration of 300 µM tris(1,3-dichloro-2-propyl) phosphate were used as

178

a positive control.11 Luminescence data were analyzed and interpreted as described previously.10

179 180 181 182

EROD and porphyrin assays EROD and total porphyrin assays were conducted as described previously9,12 and reagents were obtained from Sigma-Aldrich (St. Louis, MO). EROD activity, total protein

8 ACS Paragon Plus Environment

Page 9 of 35

Environmental Science & Technology

183

concentration and total porphyrin concentration were measured using a fluorescence plate-reader

184

(Beckman Coulter DTX880 multimode detector, Brea, CA, USA), as previously described.9,12

185 186 187

PCR array Total RNA was extracted from CEH and DCEH (n=3 technical replicates/treatment

188

group) using the Qiagen RNeasy 96-kit (Qiagen, Mississauga, ON), reverse transcribed (~200

189

ng) to cDNA and then added directly to the RT2 SYBR Green Mastermix (Qiagen).10,13

190

Mastermix (25 µL) was aliquoted to each well of the PCR array containing primers at

191

preoptimized concentrations.

192

Two custom chicken RT2 Profiler PCR Arrays were built by Qiagen/SABiosciences

193

(Frederick, MD) according to our specifications: 1) a PCR array containing 27 genes14

194

(Catalog#: CAPG12083) to determine the effects of petcoke extract exposure on the AhR

195

pathway; and 2) a chicken ToxChip PCR array comprising 43 target genes, which permits a more

196

exhaustive coverage in terms of potential toxicity pathways being impacted (see Table S2 for a

197

list of genes and pathways; Catalog#: CAPG13553). The double-crested cormorant ToxChip

198

PCR array was custom-designed and built in our lab and is described in detail elsewhere.10 Each

199

96-well PCR array contained either three identical sets of 27 target genes or two identical sets of

200

43 target genes plus 5 control genes. The 5 control genes included 2 internal control genes, a

201

positive PCR control (tests the efficiency of the PCR among/within plates using a pre-dispensed

202

artificial DNA sequence and the primer set that detects it), a reverse transcription control, and a

203

genomic DNA contamination control. Arrays were run using the Stratagene MX3005P PCR

204

system (Agilent Technologies, Santa Clara, USA) as described previously.10,13,14

205

9 ACS Paragon Plus Environment

Environmental Science & Technology

206 207

Page 10 of 35

Data analysis A one-way analysis of variance (ANOVA) with Tukey’s honestly significant differences

208

(HSD) post hoc test was used to determine significant differences in CEH viability between

209

various concentrations of the 3 petcoke extract preparations and the vehicle control in GraphPad

210

Prism v5.02 (San Diego, CA, USA).

211

EROD activity data were fit to a modified Gaussian curve as described previously9,10 and

212

EC50 values were calculated if possible (i.e. if a maximal EROD induction was reached).

213

ECthreshold (ECthr) values were determined if a maximal EROD induction was not reached and for

214

porphyrin concentration data, as described previously.10

215

PCR array data were analyzed using MxPro v4.10 software (Agilent Technologies) and

216

an example of raw cycle threshold data is provided in Table S3. Cycle threshold (Ct) data were

217

normalized to 2 internal control genes (elongation factor 1-alpha [Eef1a1] and ribosomal protein

218

L4 [Rpl4]) and the fold change of target gene mRNA abundance relative to the vehicle control

219

was calculated using the 2-∆∆Ct method.15 The PCR Arrays demonstrate strong correlations across

220

technical replicates with average correlation coefficients > 0.99 ensuring reliable detection of

221

differences in expression between biological samples (Qiagen/SABiosciences). Significant fold

222

change differences (pBaP>Chr~BaA~benz[ghi]perylene (BghiP).6 The most potent PAH, DahA, comprised

315

3.5% of the ∑14PAHs in the delayed petcoke extract, while Chr and BaP (4th and 5th most potent

316

in Head et al.6) comprised 40 and 11% of the ∑14PAHs, respectively. Therefore, it can be

317

hypothesized that specific PAHs are disproportionally responsible for the observed EROD and

318

porphyrin induction, but the potential for agonism, synergism, and/or antagonism among the

319

mixture of PACs cannot be discounted.

14 ACS Paragon Plus Environment

Page 15 of 35

Environmental Science & Technology

320 321

PCR array

322

i)

323

A custom-designed chicken PCR array was utilized to confirm the AhR pathway effects

324

described above with those at the transcriptome level. The expression levels of the two internal

325

control genes, Eef1a1 and Rpl4, were not affected by treatment with any of the petcoke extracts.

326

There was no amplification observed in the genomic DNA contamination control and the reverse

327

transcription control and positive PCR control met the appropriate quality control/assurance

328

criteria. Unprocessed mRNA expression fold change data for the various petcoke extract

329

concentrations are available in Table S4.

330

AhR pathway PCR array

The hierarchical clustering of xDP, xFP, and waDP extracts based on their transcriptional

331

alterations in CEH resulted in three main branches: 1) extracts that significantly altered ≥ 18/27

332

genes; 2) extracts that significantly altered 12-13/27 genes; and 3) extracts that significantly

333

altered ≤ 4 genes (Figure 2). To compare gene profiles of the petcoke extracts to a known AhR

334

agonist (TCDD), data from a previous study from our laboratory14 were included in the

335

hierarchical cluster analysis. The 3 most concentrated xDP extracts (1x10-3-10-5 g petcoke

336

equivalent/µL) elicited the most gene expression changes (22-25/27 genes) and clustered with

337

the most concentrated xFP extract (4x10-2 g petcoke equivalent/µL), which altered 18/27 genes,

338

and TCDD (1 nM), which altered 19/27 genes (Figure 2). Exposure of CEH to xDP (1x10-7) and

339

xFP (4x10-3) led to the modulation of 13 and 12 genes, respectively, and these extracts formed a

340

distinct branch. Finally, 4 extract preparations led to the alteration of ≤ 4 genes and included an

341

extremely diluted xDP extract (1x10-8, which only altered one gene), xFP (4x10-4), and waDP

342

(1x10-4 and 10-6 g petcoke equivalent/µL) (Figure 2). These transcriptomic responses for AhR

15 ACS Paragon Plus Environment

Environmental Science & Technology

343

pathway-specific gene targets are in agreement with the EROD and porphyrin assays

344

demonstrating concordance at multiple levels of biological organization.

345

Page 16 of 35

The AhR pathway PCR array has been used previously to determine gene expression

346

profiles of TCDD, natural AhR agonists and photodegraded flame retardants.14,22 We observed

347

82% directional concordance between CEH exposed to the petcoke extracts in the highly-

348

responsive cluster and TCDD (Figure 2, Table S4). Cyp1a4 was the most responsive transcript

349

on the PCR array and was up-regulated hundreds of fold in the highest 3 xDP treatment groups.

350

In CEH exposed to 1 nM TCDD, Cyp1a4 was induced 1976-fold.14 In vitro measures of avian

351

AhR activation (e.g. Cyp1a4 induction) have been shown to be good predictors of overt DLC

352

toxicity in avian embryos.23 Due to the large Cyp1a4 induction observed following exposure of

353

CEH to xDP and to a lesser extent xFP, future whole animal egg injection studies are warranted

354

to explore the toxicological implications of these transcriptomic results.

355

Regucalcin (Rgn) was significantly up-regulated (≥1.5-fold, p