Larval Red Drum - ACS Publications - American Chemical Society

Aug 2, 2017 - Larval Red Drum (Sciaenops ocellatus) Sublethal Exposure to. Weathered Deepwater Horizon Crude Oil: Developmental and. Transcriptomic Co...
3 downloads 8 Views 2MB Size
Subscriber access provided by Binghamton University | Libraries

Article

Larval red drum (Sciaenops ocellatus) sublethal exposure to weathered Deepwater Horizon crude oil: Developmental and transcriptomic consequences Elvis Genbo Xu, Alexis Khursigara, Jason Magnuson, Edward Starr Hazard, Gary T Hardiman, Andrew Esbaugh, Aaron P. Roberts, and Daniel Schlenk Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b02037 • Publication Date (Web): 02 Aug 2017 Downloaded from http://pubs.acs.org on August 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 37

Environmental Science & Technology

1

Larval red drum (Sciaenops ocellatus) sublethal exposure to weathered Deepwater Horizon

2

crude oil: Developmental and transcriptomic consequences

3

Elvis Genbo Xu1*, Alex J. Khursigara2, Jason Magnuson3, E. Starr Hazard4,5, Gary Hardiman5,6,

4

Andrew J. Esbaugh2, Aaron. P. Roberts3, Daniel Schlenk1

5

6

1

Department of Environmental Sciences, University of California, Riverside, CA 92521

7

2

Marine Science Institute, University of Texas at Austin, TX 78373

8 9 10 11 12 13 14

3

Department of Biological Sciences & Advanced Environmental Research Institute, University of North Texas, TX 76203 4

Center for Genomics Medicine, Medical University of South Carolina, Charleston, SC 29403

5

Computational Biology Resource Center, Medical University of South Carolina, Charleston, SC 29403 6

Departments of Medicine & Public Health Sciences, Medical University of South Carolina, Charleston, SC 29403

15

16

*Corresponding author

17

Corresponding author: Elvis Genbo Xu

18

Corresponding email: [email protected]

19

Corresponding address: Department of Environment Sciences, University of California,

20

Riverside, CA 92521, USA

21

Corresponding Tel.: 1-951-313-7643

22

Fax No. 1-951-827-3993 1 ACS Paragon Plus Environment

Environmental Science & Technology

23

24

Abstract

25

The Deepwater Horizon (DWH) incident resulted in extensive oiling of the pelagic zone and

26

shoreline habitats of many commercially important fish species. Exposure to water

27

accommodated fraction (WAF) of oil from the spill causes developmental toxicity through

28

cardiac defects in pelagic fish species. However, few studies have evaluated the effects of the oil

29

on near-shore estuarine fish species such as red drum (Sciaenops ocellatus). Following exposure

30

to a certified weathered slick oil (4.74 µg/L ∑PAH50) from the DWH event, significant sub-

31

lethal impacts were observed ranging from impaired nervous system development (average 17%

32

and 22% reduction of brain and eye area at 48 hpf, respectively) to abnormal cardiac morphology

33

(100% incidence at 24, 48 and 72 hpf) in red drum larvae. Consistent with the phenotypic

34

responses, significantly differentially expressed transcripts, enriched gene ontology, and altered

35

functions and canonical pathways predicted adverse outcomes in nervous and cardiovascular

36

systems, with more pronounced changes at later larval stages. Our study demonstrated that the

37

WAF of weathered slick oil of DWH caused morphological abnormalities predicted by a suite of

38

advanced bioinformatic tools in early developing red drum, and provided also the basis for a

39

better understanding of molecular mechanisms of crude oil toxicity in fish.

40 41 42 43 44 45 46 47 48

2 ACS Paragon Plus Environment

Page 2 of 37

Page 3 of 37

Environmental Science & Technology

49 50 51

52 53 54 55 56 57 58 59

Table of Contents (TOC) Art

Introduction The Deepwater Horizon (DWH) incident in 2010 was the largest marine oil spill in U.S. history,

60

resulting in the release of approximately 700 million liters of crude oil, and extensively oiling the

61

pelagic zone and shoreline habitats.1 Previous studies have investigated the developmental

62

toxicity of crude oil to different fish and identified a variety of abnormalities in cardiac function,

63

formation of the craniofacial skeleton, nervous system as well as reduced swimming

64

performance.2-9 Most of the previous studies were performed on pelagic species, but little is

65

known about the effects on local fast-developing estuarine fish species such as red drum

66

(Sciaenops ocellatus). The red drum inhabits the Southern Atlantic and Gulf of Mexico coasts, 3 ACS Paragon Plus Environment

Environmental Science & Technology

67

and generally spend their initial first 3-4 years within estuarine and nearshore waters.

68

Importantly, this species is a highly valued recreational species in Gulf of the Mexico, which

69

represents significant economic importance for Gulf coast communities.10 Red drum develops

70

faster than many of the studied pelagic species and other estuarine species,6,8 which may enhance

71

susceptibility to the developmental impacts of oil. Unlike mahi-mahi (Coryphaena hippurus),

72

red drum embryos were less sensitive to naturally weathered oil than oil directly obtained from

73

source well.11 Consequently, further investigation on molecular events and pathways responsible

74

for developmental toxicity is essential to understand different mechanisms of crude oil toxicity to

75

the different fish species.

76

The morphological and functional impacts of crude oil on fish embryonic development have

77

been intensively investigated. However, molecular initiating events preceding these defects are

78

not well understood.2,5-7,9 Studies in scombrid fish species have indicated inhibition of plasma

79

membrane potassium ion channels and intracellular calcium transporters in cardiac myocytes.12

80

Reductions in expression of these genes and others controlling intracellular calcium have also

81

been noted in different fish species.3,13 However, it is unclear how oil or PAHs within oil impact

82

the regulation and function of these channels and transporters. In addition to cardiotoxicity, our

83

previous study by RNA sequencing and physiological assessment in mahi-mahi demonstrated

84

that weathered oil also resulted in significant perturbations in metabolism, steroid biosynthesis,

85

vision, and AhR pathway suggesting other targets in addition to the heart may be involved in the

86

developmental toxicity of DWH oil.3 A more recent study in Atlantic haddock (Melanogrammus

87

aeglefinus) explored the transcriptional basis for cardiac formation, craniofacial development,

88

ionoregulation and cholesterol homeostasis, confirming a key role of intercellular calcium

89

cycling and excitation-contraction coupling after crude oil exposure.14 While cardiotoxicity of

4 ACS Paragon Plus Environment

Page 4 of 37

Page 5 of 37

Environmental Science & Technology

90

crude oil and PAHs have long been studied, central nervous system development is also affected

91

by PAH mixtures, resulting in behavioral dysfunction, alterations in locomotion, impaired visual

92

acuity, and reduced foraging efficiency.15,16 However, the underlying molecular mechanisms for

93

these effects have yet to be identified for early life stages of fish. Given that crude oil is a

94

complex chemical mixture, there likely are multiple targets involving multiple interacting

95

molecular mechanisms. High throughput sequencing (HTS) allows relatively unbiased

96

quantification of expression levels of transcripts with a high sensitivity and broad genome

97

coverage at one time of assessment, compared to other methods such as microarray and targeted

98

PCR, and downstream bioinformatic assessment has the potential to predict phenotypic outcomes

99

during developmental processes after toxicant exposure.

100

Here, we exposed red drum larvae to weathered DWH slick oil and evaluated morphological

101

anchoring transcriptional effects at three different critical windows (24, 48 and 72 hour post

102

fertilization (hpf)) of larval development using a HTS approach coupled with advanced

103

bioinformatic tools.3,4 Based on genome-wide differentially expressed transcripts, the most

104

enriched gene ontology, impacted biological processes, and canonical pathways in nervous,

105

cardiovascular and ocular systems, as well as upstream regulators were identified. These

106

molecular responses were compared at different developmental stages with phenotypical

107

measurements, providing novel insights into the mechanisms of DWH oil-induced

108

developmental toxicities in this rapidly growing estuarine species. This study also provides

109

potential phenotype-specific biomarkers through linking DEGs to toxicity endpoints, which

110

could be useful in ecological risk assessment of oil pollution.

111 112 113

Materials and Methods

5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 37

114

Animals and DWH oil exposure

115

Embryonic red drum were collected from brood stock tanks at the Texas Parks and Wildlife –

116

CCA Marine Development Center in Corpus Christi, Texas and transported under constant

117

aeration to the University of Texas Marine Science Institute. Embryos were subsequently treated

118

with formalin during aeration. Embryos were then rinsed with sterilized seawater and checked

119

for buoyancy and coloration using a Nikon SMZ2800N microscope. Spawns with low

120

fertilization rates or poor egg quality were not used. The oil was a weathered oil collected from

121

a slick in the Gulf of Mexico on June 29th, 2010 from the hold of barge number CT02404. Oil

122

exposures were generated according to standard protocols for high energy water accommodated

123

fractions (HEWAF), as previously described.7 Oil loading rate was 1 g per 1 liter of seawater (35

124

ppt). Red drum spawned at night and eggs were collected in the next morning. Treatments

125

started at 12 hpf, and red drum eggs typically hatch at 24 hpf. Time course exposure to WAF of

126

weathered slick oil (2.5%; 4.74 µg/L PAHs) was at 24, 48 and 72 hpf with four replicates. The

127

test was performed in an environmental control chamber set at 25 °C and 30 ppt salinity with a

128

14:10 h light: dark photoperiod. Survival was assessed daily. A minimum of a 70% hatching

129

success was required followed by a 80% survival for all tests. All experiments were approved by

130

the institutional animal care and use committee (IACUC) at the University of Texas at Austin

131

(AUP-2014-00375).

132

133

Water chemistry analysis

134

A sub-sample of the initial diluted HEWAF was collected in a 250 mL amber bottle and stored at

135

4°C

until

analyzed

by

ALS

Environmental

(Kelso,

6 ACS Paragon Plus Environment

WA)

for

PAHs

using

gas

Page 7 of 37

Environmental Science & Technology

136

chromatography/mass spectrometry-selective ion monitoring (GC/MS-SIM; based on EPA

137

method 8270D). Individual measurements (µgL-1) for 76 PAHs were quantified from initial

138

samples collected from HEWAF dilutions. Reported ΣPAH values represent the sum of 50 select

139

PAH analytes. Temperature, pH, dissolved oxygen, and salinity were measured at test start (0

140

hour), and subsequent days (24, 48, and 72 hpf). A summary of all measured water quality

141

parameters and ΣPAH concentrations are provided in Supporting Information Tables S1, S2.

142

143

Morphological Characteristics

144

Initial range finding morphological measurements were conducted on 24, 48 and 72 hpf larvae

145

treated with the WAF of slick oil (15.8 µg/L PAHs) as well as controls with one replicate (n = 10)

146

to characterize developmental features. Given more pronounced morphological changes were

147

observed at 48 hpf, a second treatment was carried out on larvae treated with the slick oil (4.74

148

µg/L PAHs) and controls collected at 48 hpf from four replicate beakers.

149

anesthetized using 250 mgL-1 of MS222 (buffered with 500 mgL-1 NaHCO3). Individuals were

150

then mounted in left lateral view onto 3 % methylcellulose in a Petri dish for image collection.

151

All images were collected using a Nikon SMZ800N microscope and Nikon Digital Sight DS U-3

152

and associated software. Still frames were imported into Image J to quantify individual

153

morphological endpoints, including brain area, eye area, iris area, pericardial area spine length

154

and total body length.

155

156

RNA sequencing and de novo assembly

7 ACS Paragon Plus Environment

Larvae were

Environmental Science & Technology

157

A detailed description of RNA isolation, cDNA library construction and sequencing is presented

158

in Xu et al.3 Larvae collected in triplicate at 24, 48 and 72hpf were used for sequencing. Briefly,

159

the 25 red drum larvae from each replicate were pooled and total RNA was isolated and purified

160

with RNeasy Mini Kit (Qiagen, Valencia, California). The 200 ng of total RNA were used to

161

prepare RNA-Seq libraries using the TruSeq RNA Sample Prep kit following the Illumina

162

protocol (Illumina, San Diego, CA). Single read sequencing (1X50; 50 million reads per sample)

163

was performed on an Illumina HiSeq 2500. Data were subjected to Illumina quality control (QC)

164

procedures (>80% of the data yielded a Phred score of 30). The read data were deposited in the

165

NCBI database (Accession Number: GSE90113). Adapter sequences were trimmed off from the

166

raw sequences and filtered using Trimmomatic (version 0.33).17 Trinity (version 2.2.0) was used

167

for de novo assembly.18 Trans-decoder (version 3.0.0) for coding sequence prediction, and

168

Trinotate (2.0.2) for functional annotation. A detailed assembly and annotation protocol can be

169

found at Xu et al.4

170

171

Time-course transcriptomic analysis

172

A time-course (24, 48 and 72 hpf) transcriptomic analysis was carried out on an OnRamp

173

Bioinformatics Genomics Research Platform (OnRamp Bioinformatics, San Diego, CA).

174

Basically, OnRamp mapped reads to the Fugu (Takifugu rubripes) transcriptome (FUGU4) using

175

BLASTX, generated gene-level count data, and performed differential expression analysis with

176

DEseq2. The protein FASTA sequences from Ensembl for Fugu were compared using Ensembl's

177

homology to create an Entrez gene list that mapped via Fugu to red drum. A detailed description

178

of these methods is presented in Xu et al.3 Statistical bioinformatics analysis was carried out

8 ACS Paragon Plus Environment

Page 8 of 37

Page 9 of 37

Environmental Science & Technology

179

using the sorted gene list by DAVID Bioinformatics Resources, Advaita Pathway Guide and

180

Ingenuity Pathway Analysis (IPA, Qiagen, Valencia, CA). The gene lists were first examined on

181

GO terms against Fugu reference using DAVID Bioinformatics Resources. The analysis was

182

performed

183

(GOTERM_CC_DIRECT)

184

KEGG_PATHWAY using the functional annotation tool with a corrected Fisher exact p-value

185

(EASE score) < 0.1. The rationale behind using Advaita and IPA is that this approach has been

186

demonstrated to improve functional analysis of fish genes with a more sensitive systems level

187

interrogation, by providing access to the best-annotated databases for human/mouse/rat

188

models,19,20 while limitations of the mapping due to the extra genome duplication events in

189

teleost fish and species differences in gene functions still exist. Biological process, microRNA

190

and disease were enriched and predicted by Advaita. IPA was used to examine the toxicity

191

pathways and reveal biological pathways/mechanisms underlying toxicity-specific phenotypes,

192

as well as to provide insights into phenotype-specific biomarkers through linking DEGs to their

193

known role in specific toxicity endpoints.4 Canonical pathways were further overlapped and

194

clustered using IPA to understand how pathways may impact function, and whether there are any

195

linkages among significantly disturbed pathways.

on

biological

process

(GOTERM_BP_DIRECT),

molecular

function

cellular

component

(GOTERM_MF_DIRECT)

and

196

197

Results

198

Chemical composition of weathered slick oil

199

The weathered slick oil predominantly consisted of 3-ring (72%) and 4-ring (22%) PAHs (Fig.

200

S1). The most abundant compounds were Phenanthrenes and Anthracenes, following by 9 ACS Paragon Plus Environment

Environmental Science & Technology

201

Dibenzothiophene. The profile is similar to those obtained from HEWAF preparations using

202

other sources of slick oil from the DWH spill.3,5

Page 10 of 37

203 204

Morphological abnormalities from slick oil exposure

205

Developmental stages of red drum are similar to those of many warm water marine fishes with

206

early cell divisions essentially complete by 1.5 hours and an early embryo formed by 12 h. At a

207

rearing temperature of 25°C, red drum embryos began hatching at 21 hpf, and most hatched at 24

208

hpf. At 48 hpf, the liver, stomach and intestine were well-differentiated, and the eyes were

209

pigmented. The mouth was functional and the yolk-sac was lost between 72 and 96 hpf. The

210

growth in total length was minimal from 24 to 72 hpf, but significant growth in area of eye and

211

brain during 24-72 hpf indicated the larvae were undergoing accelerated head development (see

212

controls in Table S3). At 48 hpf, the slick oil treatment significantly increased pericardial area

213

(45%), and decreased brain area (17%) and eye area (22%) in larvae compared to controls. The

214

total body length was also significantly reduced by 12% (Table 1). The slick oil treatment

215

increased pericardial area and decreased tectum area in larvae at 24, 48 and 72 hpf compared to

216

controls, respectively (one replicate; Table S3). The slick oil treatment decreased the length of

217

body and notochord at 24, 48 and 72 hpf, respectively, compared to time-matched controls.

218

Notably, eyes were smaller in oil-treated larvae indicated by reduced lens diameter at 24, 48 and

219

72 hpf (Table S3). The pupil diameter was also reduced at 48 and 72 hpf, and the effect was

220

more pronounced at 48 hpf (Table S3).

221

222

De novo assembly of red drum transcriptome and annotation

10 ACS Paragon Plus Environment

Page 11 of 37

Environmental Science & Technology

223

A total of 2,168,925,036 Illumina HiSeq reads from red drum larvae were generated. After

224

trimming the adapters, 181,576,332 bases were assembled with Trinity resulting in 257,434

225

transcript contigs with an average length of 705 bp and an N50 of 1,071 bases. 30,950

226

HMMER/PFAM protein domains (Pfam), 9,219 predicted transmembrane regions (TmHMM),

227

34,080 non-supervised orthologous groups of genes (eggNOG), and 38,647 GO_blast were

228

determined using Trinotate pipeline (Table S4). The final transcriptome assembly provided a

229

valuable resource for wider and deeper genetic research on red drum.

230

231

Transcriptional responses at 24 hpf

232

In 24 hpf animals exposed to slick oil, 104 genes were significantly differentially expressed at an

233

false discovery rate (FDR) < 0.4, with 78% up-regulated (Fig. S2a). The significant Gene

234

Ontology (GO) terms (biological process, cellular component, molecular function) and Kyoto

235

Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by DAVID against Fugu

236

background (Fig. S3). Biological processes and pathways changed by oil exposure included

237

transcripts involved in fatty acid metabolism, lipid homeostasis, cell redox homeostasis, ECM-

238

receptor interaction, focal adhesion, metabolic pathways (tyrosine, arachidonic acid, glutathione,

239

carbon), and ubiquinone biosynthesis. For the ontology of molecular function, Flavin adenine

240

dinucleotide binding and oxidoreductase activity were the predominant transcripts altered by oil

241

exposure. Consistent with DAVID assessments, Advaita analysis revealed metabolic pathways,

242

protein digestion, metabolism of xenobiotics by cytochrome P450, extracellular matrix (ECM)-

243

receptor interaction, steroid biosynthesis, and focal adhesion as the most highly ranked bio-

244

pathways (Table S5). IPA predicted significant activation of aryl hydrocarbon receptor (AhR)

11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 37

245

signaling, NRF2-mediated oxidative stress responses, fatty acid metabolism, hepatic fibrosis,

246

xenobiotic metabolism, cholesterol biosynthesis, and transmembrane potential of mitochondria

247

by slick oil at 24 hpf (Fig. 1).

248

249

Transcriptional responses at 48 hpf

250

Compared to 24 hpf, the expression of more genes (769 genes) was significantly altered at 48 hpf

251

after slick oil exposure. The enriched GO terms were generally different between 24 hpf and 48

252

hpf after slick oil exposure, although some common terms were enriched at both 24 and 48 hpf,

253

such as lipid homeostasis, proteinaceous extracellular matrix, collagen trimmer, mitochondrion,

254

oxidoreductase activity, ECM-receptor interaction, and carbon metabolism. The most

255

representative GO terms by slick oil exposure at 48 hpf included proteolysis and proteasome-

256

related terms. Notably, some neuron-associated terms started to be highly enriched at this stage,

257

such as axonogenesis, neural crest cell migration, neuron projection, and synapse, suggesting the

258

disruption of nervous system development. Similar to DAVID, Advaita revealed proteasome,

259

protein digestion, and metabolic pathways as the highest ranked bio-pathways. Several genes

260

(e.g. rho, grm6, gnat1,gnat2, cnga3, cngb3, cacna1f) that were significantly differentially

261

downregulated after slick oil exposure were linked to eye diseases, e.g. retinitis pigmentosa,

262

congenital stationary night blindness, cataracts, and cone-rod dystrophy

263

nervous system development and function (174 genes) was the highest ranked physiological

264

phenotype at 48 hpf predicted by IPA. The highest ranked toxicity predictions included NRF-

265

mediated oxidative stress response (ranked No.1; 26 genes; p-value = 1.65E-08), cardiac

12 ACS Paragon Plus Environment

(Table S5).

The

Page 13 of 37

Environmental Science & Technology

266

hypertrophy, mitochondrial dysfunction, hepatic fibrosis, xenobiotic metabolism, AhR signaling,

267

increased heart failure, liver proliferation, and PPARα/RXRα activation (Fig. 1).

268

269

Transcriptional responses at 72 hpf

270

The number of significantly differentially expressed genes (FDR < 0.4) greatly increased from

271

769 at 48 hpf to 2,564 genes at 72 hpf after slick oil exposure. Compared to 24 and 48 hpf, the

272

significantly enriched GO terms were vastly different indicated by DAVID in 72 hpf after slick

273

oil exposure. Protein translation, cytosolic ribosomal subunits, structural constituents of the

274

ribosome, and ribosome were the most significantly enriched biological processes, cellular

275

components, molecular functions and KEGG pathways in 72 hpf after slick oil exposure,

276

respectively (Fig. S3). Ribosome and metabolic pathways were the highest ranked bio-pathways

277

identified by Advaita in 72 hpf animals after slick oil exposure (Table S5). Comparing DAVID

278

and Advaita Pathway Guide analyses, both approaches were consistent with regard to the

279

enriched GO pathways identified, with ribosomal and metabolic pathways being the highest

280

ranked pathways (Fig. 1). Advaita also predicted that other biological processes and molecular

281

functions were altered, including phototransduction, steroid biosynthesis, retinol metabolism as

282

well as the disease pathways, diamond-blackfan anemia (faulty ribosome biogenesis), retinitis

283

pigmentosa, and cytochrome c oxidase deficiency. Notably, at 72 hpf after slick oil exposure,

284

the genes involved in ribosomal pathway were all significantly downregulated. IPA analysis

285

further predicted the activation of cardiac hypertrophy (ranked No.1; 95 genes; p-value = 1.49E-

286

12), mitochondrial dysfunction, and xenobiotic metabolism, renal necrosis, TR/RXR activation,

287

cardiac necrosis, FXR/RXR activation, AhR signaling, hypoxia-inducible factor signaling, liver

13 ACS Paragon Plus Environment

Environmental Science & Technology

288

proliferation, increases heart failure, cholesterol biosynthesis, and decreases in transmembrane

289

potential of mitochondria by slick oil exposure (Fig. 1).

Page 14 of 37

290

291

Upstream analysis

292

The upstream analysis by IPA was used to identify the activation state of upstream regulators

293

that could explain the observed gene expression profile alterations in slick oil-treated red drum

294

larvae at 24, 48 and 72 hpf. The predicted top activated or down-regulated regulators in red drum

295

larvae treated with slick oil are shown in Fig. 2. More similar regulators were observed between

296

48 and 72 hpf groups with fewer at 24 hpf. Nuclear factor nfe2l2 was predicted as the most

297

relevant activated upstream regulator across all developmental stages. nfe2l2 regulates genes

298

which contain antioxidant response elements (ARE) in their promoters, and many of these genes

299

are known to encode proteins involved in cellular responses to oxidative stress, immune system

300

processes, metabolism of xenobiotics, oxidoreductase activity and apoptosis. Examples include

301

sqstm1, gsto1, prdx1, hspa9, creg1 and osgin1, which were all up-regulated in the present study.

302

Another upstream regulator that was highly ranked was AhR. Its target genes, cyp1a1, ahrr,

303

cyp3a4 were among the most upregulated genes in our dataset; a number of collagens that play

304

an important role in maintaining the integrity of various tissues were significantly downregulated

305

also via AhR activation. One mature microRNA was also predicted as a top regulator, miR-124,

306

that targets ahrr mRNA in humans.21 A list of microRNAs was predicted based on differential

307

expression of their target mRNA by Advaita (Table S5), and deserve further investigation.

308

309

Discussion 14 ACS Paragon Plus Environment

Page 15 of 37

Environmental Science & Technology

310

The global changes in gene expression in embryonic exposures using transcriptomic and

311

bioinformatic tools can help target conventional morphological assessment in detecting

312

subsequent phenotypic effects in organisms exposed to various stressors. Assessments with

313

these tools may improve predictability, sensitivity and efficiency of toxicity assessments in

314

embryonic fish and other species.

315

opportunity to systematically identify most of the genes that are responsive in stressed embryos

316

or larvae. Embryonic gene expression is a highly dynamic process with unique changes at each

317

developmental stage. The patterns of gene expression across major developmental stages have

318

been known to have lasting impacts on the development of principal physiological systems in

319

vertebrate species.22-24 In particular, the embryo-to-larval stage is a crucial period in the life of

320

fish, and is more sensitive to environmental stress than adult stages.25 Our previous study in

321

mahi-mahi embryos and larvae using toxicogenomic tools demonstrated the toxicity of crude oil

322

is time- and oil type-dependent.3 To our knowledge, the present study is the first to de novo

323

assemble the transcriptome of red drum exposed to DWH crude oil. A list of phenotype-

324

associated DEGs for red drum was identified at different developmental stages (Table S6). This

325

DEG analysis illustrated similarities and differences in responses to DWH crude oil at different

326

developmental times at individual gene levels.

327

Consistent with studies for other fish species,2,3,13 our results also support evidence for the most

328

dominant malformations observed being abnormal development of nervous and cardiovascular

329

systems, which anchor changes in gene expression to predicted phenotypes by a suite of

330

bioinformatic tools (i.e., DAVID, Advaita and IPA). This agrees well with the studies in normal

331

developing fish species indicating that the primary differentially expressed functional groups

332

during the embryo-to-larval transition are nervous, muscular, and cardiovascular system

The advent of HTS technology provides an unbiased

15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 37

333

development.25 Similar to mahi-mahi,3 the transcriptional responses in red drum induced by oil

334

exposure is also time- dependent, and this could due to transcriptional changes associated with

335

normal embryonic developmental processes.26

336

PAHs have long been studied for their carcinogenic properties and cardiotoxicity, and more

337

recent studies have highlighted the neurotoxicity of individual PAHs and mixtures. Prenatal

338

exposure to benzo(a)pyrene (BaP) has been shown to induce neurological abnormalities like

339

cognitive impairment, learning difficulties, and loss of short-term memory in human.27 BaP

340

decreased

341

neurodegeneration in zebrafish.28 Exposure to PAH mixtures impaired neurodifferentiation and

342

caused corresponding neurobehavioral impairment in fish.16 A more recent study shows a

343

complex PAH mixture (mostly naphthalene, phenanthrene, fluoranthene and acenapthene,

344

similar to the composition of the slick oil in the present study) significantly impaired early-stage

345

neurodevelopment, resulting in a profound drop (up to 30%) in the glia-to-neuron ratio in

346

embryonic neural stem cells, which can affect brain homeostasis and function.29 In the present

347

study, nervous system development and function was predicted as the top-ranked physiological

348

phenotype in 48 hpf red drum exposed to slick oil, and included significant inhibition of synapse

349

development (p=1.43E-6, z score = -2.907), quantity of nervous tissue (p = 8.05E-05, z score = -

350

2.759), quantity of neurons (p = 1.56E-04; z score = -2.624), long term depression of synapse (p

351

= 5.72E-05; z score = -2.2), and impairment of synaptic transmission of nervous tissue (p =

352

1.44E-05, z sore = -2.138).

353

signaling pathways in the nervous system as top-ranked pathways in red drum larvae after slick

354

oil exposure at both 48 and 72 hpf (Fig. 3). These results are consistent with other studies in

355

embryonic fish species that demonstrated nervous system development and function were the

brain

mass,

locomotor

activity,

dopaminergic

neurons

and

resulted

in

Predictions by IPA also included impairment of a number of

16 ACS Paragon Plus Environment

Page 17 of 37

Environmental Science & Technology

356

most significantly enriched biological functions during normal embryo-to-larva transition.24,30 In

357

our study, the AhR pathway remained the top-ranked toxicity pathway in all three stages of red

358

drum larvae treated by slick oil, with the least significance (largest adjusted p-value) at 48 hpf

359

(Fig. 1), but a direct linkage from AhR pathways to neural signaling pathway and neurotoxicity

360

of slick was not identified in the present study. Activation of the AhR contributes to adverse

361

effects of PAHs on proliferation and migration of mouse but not human neural stem cells.31 The

362

molecular connections between AhR activation and neurotoxicity of oil is still unclear,

363

particularly considering PAHs with strong affinity to AhR were not effective in altering mouse

364

neural stem cells.32 Thus, it is likely that other mechanisms may contribute to the adverse effects

365

rather than just those mediated by the AhR. For example, the cholinergic system influences

366

cognition, anxiety, locomotion, and behavior by acting upon nicotinic acetylcholine receptors

367

(nAChRs), and diminishment of nAChRs can impair passive avoidance in mice.33 It has also

368

been established that the loss of nAChRs in the cerebral cortex and certain subcortical regions is

369

associated with brain dysfunction and degenerative diseases such as Alzheimer’s, Parkinson’s

370

and Lewy body diseases.34 We found three genes encoding nAChR subunits (chrna6, chrnb2 and

371

chrnd) that were all significantly downregulated in 48 hpf red drum larvae after slick oil

372

exposure (Fig. 4).

373

Cannabinoid receptor 1 (cnr1) is a G protein-coupled receptor located primarily in the central

374

and peripheral nervous systems. In medium-size spiny neurons from rat striatum, activation of

375

cnr1 is necessary for long term depression of excitatory synapses.35 In astroglial cells from

376

mouse brain, mutant mouse cnr1 knockout as well as Ca2+/ calmodulin dependent kinase camk2a

377

gene knockout decreases long term depression of synapses,36,37 and causes hyperactivity

378

behavior in mouse.38,39

In embryos of rockfish (Sebastiscus marmoratus), the decrease of

17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 37

379

camk2a mRNA expression might contribute to altered synaptic plasticity and neuronal survival

380

following exposure to the five-ring PAH, Benzo(a)pyrene.40 Neuronal signaling is also

381

diminished by glutamate ionotropic receptor gria1,41,42 and mutant gria1 gene knockouts

382

increased hyperactivity in mouse.43 The cnr1, camk2a, gria1and cholinergic receptor genes were

383

all significantly downregulated in 48 hpf red drum larvae after slick oil exposure (Fig. 4). This is

384

consistent with our previous study in mahi-mahi after slick oil exposure,3 suggesting

385

neurotoxicity of slick oil through depression of neurotrophin as well as changes in glutamate,

386

cholinergic, cannabinoid receptors and also Ca2+ homeostasis (Fig. 4). Since Ca2+ appears to also

387

play a significant role in cardiac toxicity, impairment of multiple targets may be responsible for a

388

number of phenotypic responses besides that of cardiotoxicity that may adversely affect

389

development at later life stages in red drum.

390

Although other targets may be possible, the present findings are consistent with previous studies

391

indicating that pericardial edema is a profound response in fish species after crude oil

392

exposure.3,7,11,13,44 The slick oil treatment significantly increased pericardial area in red drum

393

larvae at 24, 48 and 72 hpf compared to controls. Chemical blockade of IKr repolarizing

394

potassium currents and disruption of intra-cellular Ca2+ was suggested as the major initiating

395

events for cardiac defects in fish exposed to crude oil.12 A more recent study demonstrated a

396

transcriptional cascade that is tightly linked to defects in cardiomyocyte intracellular calcium

397

cycling and heart chamber growth through bmp10 in Atlantic haddock.14 However, mRNA of

398

bmp 10 was not changed in the current study in red drum, nor in previous studies with mahi-

399

mahi treated with slick oil that caused abnormal cardiac phenotypes. Some 4- and 6-ring PAHs

400

and oxygen-substituted PAHs can inappropriately activate the AhR in developing

401

cardiomyocytes, leading to primary defects in cardiac morphogenesis in fish species,45-47 and this

18 ACS Paragon Plus Environment

Page 19 of 37

Environmental Science & Technology

402

form of cardiac toxicity is dependent on the AhR and is prevented by AhR gene knockdown.47,48

403

In contrast, the cardiotoxicity of weathered crude oil and single non-alkylated tricyclic PAHs

404

occurs without activation of the AhR in fish cardiomyocytes,49,50 and is not prevented by AhR

405

gene knockdown.51 Since the slick DHW oil is a complex mixture, multiple toxic mechanisms

406

were possibly involved in the cardiotoxicity, including both AhR-dependent and AhR-

407

independent pathways. Our data demonstrated several AhR-independent canonical pathways

408

associated with heart development and function were modified at the transcriptional level. For

409

example, IPA analysis predicted the significant activation of NFAT (nuclear factor of activated

410

T-cells) in cardiac hypertrophy (Rank #25 at 48 hpf and Rank #6 at 72 hpf; Fig. 3). Previous

411

studies have recognized the importance of Ca2+-sensitive signaling molecules, including

412

calcineurin, CalmK, and MAPK in hypertrophic pathways, in which NFAT plays a critical role,

413

being the best-characterized target for the development of cardiac hypertrophy.52 Prolonged

414

cardiac hypertrophy is associated with arrhythmia, sudden death, decompensation, and dilated

415

cardiomyopathy.53 The hypertrophic response is orchestrated by growth factors and cytokines

416

acting through several interdependent signaling cascades whose molecules include G-proteins

417

such GNAI1, GTPases such as SOS Ras, and kinases such as MAPK and CAMK, as well as

418

transcription factors HAND1, HDAC, which are essential for cardiac development and altered by

419

slick oil exposure (Fig. S4).

420

In addition to pericardial edema, impaired eye formation was the most noted indicator at 48 and

421

72 hpf. Forty-four significantly differentially expressed genes were identified in the pathway of

422

eye formation and retinal degeneration uniquely associated with the 48 hpf larvae after slick oil

423

exposure (Fig. S5).

424

oxidation producing numerous oxygenated PAHs.54 Eye malformation was induced in zebrafish

DWH slick oil has been weathered and undergone significant photo-

19 ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 37

425

embryos treated with oxygenated PAHs (such as 1,9-benz-10-anthrone ), and mRNAs associated

426

with visual perception were also altered.46,55 IPA predicted inhibition of phototransduction and

427

retinal degeneration by slick oil at 48 hpf (Fig. S5), and retinitis pigmentosa (RP) was the top

428

ranked disease by Advaita at both 48 and 72 hpf (Table S5). Photoreceptors respond to light by

429

the closure of a cyclic nucleotide-gated (CNG) channel (activated by binding of cGMP/cAMP),

430

causing hyperpolarization of the plasma membrane and decrease of the synaptic glutamate

431

release. CNG channels provide the only source for Ca2+ influx into rod and cone outer segments

432

for calcium, the lack of this channel would result in a decrease of cellular calcium

433

concentration.56,57 At low calcium concentrations, the membrane guanylate cyclase may be

434

stimulated permanently, leading to elevated levels of cGMP, which may trigger the degeneration

435

of photoreceptors.58 The absence of functional CNG channels is equivalent to the permanent

436

closure of channels which occur under continuous bright light conditions. Continuous exposure

437

of experimental animals to light has been shown to result in photoreceptor degeneration.59,60

438

Consistent with previous study in mahi-mahi, cyclic nucleotide-gated channel alpha 3 (cnga3),

439

retinoid isomerohydrolase (rpe65), guanylate cyclase 2D (gucy2d), transducin (gnat2), and their

440

upstream regulator orthodenticle homeobox 1 (otx1) were significantly downregulated, which

441

could lead to the perturbation of rhodopsin regeneration and phototransduction by slick oil (Fig.

442

S5). Similarly, both BaP and crude oil caused ocular toxicity and adverse affected visual-

443

associated behavior in fish species.61,62

444

exposure by disturbing 14 different pathways that interact with phototransduction pathways (Fig.

445

S6), including relaxin signaling, a-adrenergic signaling, CXCR4 signaling, CREB signaling in

446

neurons and androgen signaling through molecular connections such as G protein transducins

447

and guanine nucleotide binding protein (Fig. 4). The observed eye morphology would cause

Such ocular toxicity may be induced by slick oil

20 ACS Paragon Plus Environment

Page 21 of 37

Environmental Science & Technology

448

severe disadvantages during fish life cycles by affecting feeding, competing, escaping, migration

449

responses and survival.

450

In conclusion, this is the first study to de novo assemble the transcriptome and investigate the

451

time-course transcriptomic responses in estuarine fish larvae exposed to DWH oil. Comparisons

452

of phenotypic responses with transcriptomic profiles were consistent with significantly

453

differentially expressed genes, enriched gene ontology, functions and canonical pathways mostly

454

involved in nervous and cardiovascular systems, both of which were validated at the phonotypic

455

level using morphometric analysis, and with more pronounced changes at later larval stages.

456

Overall, the present study has provided valuable molecular resources for expanding our

457

understanding of the developmental toxicity of DWH oil in fish. The developed de novo and

458

comparative toxicity pathway methods will also be useful for toxicogenomic study on other non-

459

model fish species. This study calls for follow-up studies to confirm the initial transcriptomic

460

results with additional methods (e.g. in situ expression in specific organ and tissue, dose-

461

response assays) at an expanded range of oil concentrations, and to assess the later life

462

consequences of larval exposure at the population level.

463 464

Supporting Information

465

Further information is available that provides PAH measurements (Fig. S1; Table S2), water

466

quality measurements (Table S1), plots showing relative expression of genes (Fig. S2),

467

significantly enriched gene ontology (GO) terms and KEGG pathways, diseases and miRNAs

468

(Fig. S3; Table S5), morphological measurements (Table S3), statistics of transcriptome

469

assembly and annotation (Table S4), plot showing the role of NFAT in cardiac hypertrophy (Fig.

470

S4), inhibition of eye formation and induced retinal degeneration ribosome biosynthesis (Fig. S5), 21 ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 37

471

overlapping canonical pathways (Fig. S6), the top predicted canonical pathways with

472

corresponding genes (Table S6).

473 474

Acknowledgments

475

This research was made possible by a grant from The Gulf of Mexico Research Initiative. Grant

476

No: SA-1520; Name: Relationship of Effects of Cardiac Outcomes in fish for Validation of

477

Ecological Risk (RECOVER). The authors are grateful to the Gulf of Mexico Research Initiative

478

Information and Data Cooperative (GRIIDC) for supporting data management system to store

479

the data generated ( doi:10.7266/N7S180KJ). G. Hardiman acknowledges Medical University of

480

South Carolina College of Medicine start-up funds. This research was supported in part by the

481

Genomics Shared Resource, Hollings Cancer Center.

482

483

References

484

(1)

485

fact and conclusions of law: Phase Two trial. In re: Oil spill by the oil rig “Deepwater Horizon”

486

in the Gulf of Mexico, on April 20, 2010, No. MDL 2179, 2015 WL 225421 (LA. E.D. Jan. 15,

487

2015). (Doc. 14021). U.S. District Court for the Eastern District of Louisiana.

488

(2)

489

Larsen, M.; Stekoll, M.S.; Rice, S.D.; Collier, T.K; Scholz, N.L. Sublethal exposure to crude oil

490

during embryonic development alters cardiac morphology and reduces aerobic capacity in adult

491

fish. Proc. Natl. Acad. Sc. 2011, 108 (17): 7086−7090.

United States of America v. BP Exploration & Production, Inc., et al. 2015. Findings of

Hicken, C.E.; Linbo, T.L.; Baldwin, D.H.; Willis, M.L.; Myers, M.S.; Holland, L.;

22 ACS Paragon Plus Environment

Page 23 of 37

Environmental Science & Technology

492

(3)

Xu, E.G.; Mager, E.M.; Grosell, M.; Pasparakis, C.; Schlenker, L.S.; Stieglitz, J.D.;

493

Benetti, D.; Hazard, E.S.; Courtney, S.M.; Diamante, G; Freitas, J. Time-and oil-dependent

494

transcriptomic and physiological responses to Deepwater Horizon oil in mahi-mahi (Coryphaena

495

hippurus) embryos and larvae. Environ. Sci. Technol. 2016, 50(14):7842-7851.

496

(4)

497

transcriptome assembly and comparative toxicity pathway analysis in mahi-mahi (Coryphaena

498

hippurus) embryos and larvae exposed to Deepwater Horizon oil. Sci. Rep., 2017, 7: 44546.

499

(5)

500

Scholz, N.L.; Benetti, D.D.; Grosell, M. Acute embryonic or juvenile exposure to Deepwater

501

Horizon crude oil impairs the swimming performance of mahi-mahi (Coryphaena

502

hippurus). Environ. Sci. Technol. 2014, 48(12): 7053-7061.

503

(6)

504

Stieglitz, J. D.; French, B. L.; Labenia, J. S.; Laetz, C. a.; et al. Deepwater Horizon crude oil

505

impacts the developing hearts of large predatory pelagic fish. Proc. Natl. Acad. Sci. 2014, 111

506

(15), E1510–E1518.

507

(7)

508

Linbo, T.L.; Lay, C.; Forth, H.; Scholz, N.L.; Incardona, J.P. The effects of weathering and

509

chemical dispersion on Deepwater Horizon crude oil toxicity to mahi-mahi (Coryphaena

510

hippurus) early life stages. Sci. Total Environ. 2016, 543: 644-651.

511

(8)

512

Genomic, and Developmental Effects of the Deepwater Horizon Oil Spill on Resident Gulf Killi

513

fish (Fundulus grandis). Environ. Sci. Technol. 2013, 47, 5074–5082.

Xu, E.G.; Mager, E.M.; Grosell, M.; Hazard, E.S.; Hardiman, G; Schlenk, D. Novel

Mager, E.M.; Esbaugh, A.J.; Stieglitz, J.D.; Hoenig, R.; Bodinier, C.; Incardona, J.P.;

Incardona, J. P.; Gardner, L. D.; Linbo, T. L.; Brown, T. L.; Esbaugh, a. J.; Mager, E. M.;

Esbaugh, A.J.; Mager, E.M.; Stieglitz, J.D.; Hoenig, R.; Brown, T.L.; French, B.L.;

Dubansky, B.; Whitehead, A.; Miller, J. T.; Rice, C. D.; Galvez, F. Multitissue Molecular,

23 ACS Paragon Plus Environment

Environmental Science & Technology

Page 24 of 37

514

(9)

Johansen, J.L.; Esbaugh, A.J. Sustained impairment of respiratory function and swim

515

performance following acute oil exposure in a coastal marine fish. Aquat. Toxicol. 2017, 187,

516

82-89.

517

(10)

518

U.S. Dept. of Commerce, NOAA Tech. Memo. NMFS-F/SPO-163, 237p

519

(11)

520

function and survival are affected by crude oil in larval red drum, Sciaenops ocellatus. Sci. Total

521

Environ. 2017, 579, 797-804.

522

(12)

523

impairs cardiac excitation-contraction coupling in fish. Science, 2014, 343(6172): 772-776.

524

(13)

525

Meeren, T.; Thorsen, A.; Thorbjørnsen, M.; Jentoft, S.; Edvardsen, R.B. Crude oil exposures

526

reveal roles for intracellular calcium cycling in haddock craniofacial and cardiac development.

527

Sci. Rep. 2016, 6: 31058.

528

(14)

529

Edvardsen, R.B.; Jentoft, S. Novel adverse outcome pathways revealed by chemical genetics in a

530

developing marine fish. eLife, 2017, 6, e20707.

531

(15)

532

of naphthalene and phenanthrene on visual and prey capture endpoints during early stages of the

533

dourado Salminus Brasiliensis. Mar. Environ. Res. 2008, 66(1), 205-207.

National Marine Fisheries Service. 2016. Fisheries Economics of the United States, 2014.

Khursigara, A.J.; Perrichon, P.; Bautista, N.M.; Burggren, W.W.; Esbaugh, A.J. Cardiac

Brette, F.; Machado, B.; Cros, C.; Incardona, J.P.; Scholz, N.L.; Block, B.A. Crude oil

Sørhus, E.; Incardona, J.P.; Karlsen, Ø.; Linbo, T.; Sørensen, L.; Nordtug, T.; van der

Sørhus, E.; Incardona, J.P.; Furmanek, T.; Goetz, G.W.; Scholz, N.L.; Meier, S.;

Carvalho, P.S.; Daniel da Costa; B.K.; Novelli, G.A.; Bainy, A.C.; Fraga, A.P.M. Effects

24 ACS Paragon Plus Environment

Page 25 of 37

Environmental Science & Technology

534

(16)

Brown, D.R.; Bailey, J.M.; Oliveri, A.N.; Levin, E.D.; Di Giulio, R.T. Developmental

535

exposure to a complex PAH mixture causes persistent behavioral effects in naive Fundulus

536

heteroclitus (killifish) but not in a population of PAH-adapted killifish. Neurotoxicol. Teratol.

537

2016, 53, 55–63.

538

(17)

539

sequence data. Bioinformatics, 2014, 30 (15): 2114-2120.

540

(18)

541

Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; Chen, Z. Full-length transcriptome assembly

542

from RNA-Seq data without a reference genome. Nature Biotechnol. 2011, 29(7), 644-652.

543

(19)

544

Global transcriptome analysis of Atlantic cod (Gadus morhua) liver after in vivo methylmercury

545

exposure suggests effects on energy metabolism pathways. Aquat. Toxicol. 2013, 126:314–25.

546

(20)

547

Martin, S.A. Differential responses of the gut transcriptome to plant protein diets in farmed

548

Atlantic salmon. BMC genomics, 2016, 17(1), 1.

549

(21)

550

D.P.; Linsley, P.S.; Johnson, J.M. Microarray analysis shows that some microRNAs

551

downregulate large numbers of target mRNAs. Nature. 2005, 433(7027):769-773.

552

(22)

553

transition in mammals. Molecular reproduction and development. 2009, 76(9):805-818.

Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: a flexible trimmer for Illumina

Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.;

Yadetie, F.; Karlsen, O.A.; Lanzén, A.; Berg, K.; Olsvik, P.; Hogstrand, C.; Goksøyr, A.

Król, E.; Douglas, A.; Tocher, D.R.; Crampton, V.O.; Speakman, J.R.; Secombes, C.J.;

Lim, L.P.; Lau, N.C.; Garrett-Engele, P.; Grimson, A.; Schelter, J.M.; Castle, J.; Bartel,

Evsikov, A.V.; Marín de Evsikova, C. Gene expression during the oocyte‐to‐embryo

25 ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 37

554

(23)

Vesterlund, L.; Jiao, H.; Unneberg, P.; Hovatta, O.; Kere, J. The zebrafish transcriptome

555

during early development. BMC Dev. Biol. 2011, 11(1):1.

556

(24)

557

Exploring the larval transcriptome of the common sole (Solea solea L.). BMC genomics. 2013,

558

14(1):1.

559

(25)

560

understanding marine fish larval development 1 1 This review is part of a virtual symposium on

561

current topics in aquaculture of marine fish and shellfish. Can. J. Zool. 2011, 89(7):599-611.

562

(26)

563

Schlenk, D. Developmental transcriptomic analyses for mechanistic insights into critical

564

pathways involved in embryogenesis of pelagic mahi-mahi (Coryphaena hippurus). PlosOne.

565

2017, 12(7): e0180454.

566

(27)

567

Ramesh, A.; Aschner, M.; Campbell, D.; Levitt, P.; Hood, D.B. Prenatal polycyclic aromatic

568

hydrocarbon exposure leads to behavioral deficits and downregulation of receptor tyrosine

569

kinase, MET. Toxicol. Sci. 2010, 118(2): 625–634.

570

(28)

571

pyrene level causes neurodegenerative disease-like syndromes in zebrafish (Danio rerio). Aquat.

572

Toxicol. 2015, 167: 200-208.

573

(29)

574

differences in the developmental neurotoxicity of an environmental polycylic aromatic

Ferraresso, S.; Bonaldo, A.; Parma, L.; Cinotti, S.; Massi, P.; Bargelloni, L.; Gatta, PP.

Mazurais, D.; Darias, M.; Zambonino-Infante, J.L.; Cahu, C.L. Transcriptomics for

Xu, E.G.B.; Mager, E.M.; Grosell, M.; Stieglitz, J.D.; Hazard, E.S.; Hardiman, G.;

Sheng, L.; Ding, X.X.; Ferguson, M.; McCallister, M.; Rhoades, R.; Maguire, M.;

Gao, D.; Wu, M.; Wang, C.; Wang, Y.; Zuo, Z. Chronic exposure to low benzo [a]

Slotkin, T.A.; Skavicus, S.; Card, J.; Di Giulio, R.T.; Seidler, F.J. In vitro models reveal

26 ACS Paragon Plus Environment

Page 27 of 37

Environmental Science & Technology

575

hydrocarbon mixture compared to benzo [a] pyrene: Neuronotypic PC12 Cells and embryonic

576

neural stem cells. Toxicology, 2017, 377, 49-56.

577

(30)

578

Developmental transcriptomics in Atlantic haddock: Illuminating pattern formation and

579

organogenesis in non-model vertebrates. Dev. Biol. 2016, 411(2):301-313.

580

(31)

581

N.; Rockel, T.D.; Schreiber, T.; Fritsche, E. Species-specific differential AhR expression

582

protects human neural progenitor cells against developmental neurotoxicity of PAHs. Environ.

583

Health Perspect. 2011, 118, 1571–1577.

584

(32)

585

binding affinities and AHH induction potencies of 29 polynuclear aromatic hydrocarbons.

586

Toxicol. Lett. 1986, 34, 67–74.

587

(33)

588

neuronal nicotinic acetylcholine receptors. Neuroscience, 2004, 129(3), 575-582.

589

(34)

590

V.; Buxbaum, J.D.; Nãsland, J.; Davis, K.; Gotti, C. Nicotinic receptor subtypes in human brain

591

ageing, Alzheimer and Lewy body diseases. Eur. J. Pharmacol. 2000, 393(1), 215-222.

592

(35)

593

striatal long-term depression. J. Neurosci. 2007, 27(19):5260-5264.

594

(36)

595

forms of synaptic plasticity. Curr. Biol. 1994, 4(8):687-693.

Sørhus, E.; Incardona, JP.; Furmanek, T.; Jentoft, S.; Meier, S.; Edvardsen, R.B.

Gassmann, K.; Abel, J.; Bothe, H.; Haarmann-Stemmann, T.; Merk, H.F.; Quasthoff, K.

Piskorska-Pliszczynska, J.; Keys, B.; Safe, S.; Newman, M.S. The cytosolic receptor

Marubio, L.M.; Paylor, R. Impaired passive avoidance learning in mice lacking central

Perry, E.; Martin-Ruiz, C.; Lee, M.; Griffiths, M.; Johnson, M.; Piggott, M.; Haroutunian,

Singla, S.; Kreitzer, A.C.; Malenka, R.C. Mechanisms for synapse specificity during

Stevens, C.F.; Tonegawa, S.; Wang, Y. The role of calcium-calmodulin kinase II in three

27 ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 37

596

(37)

Han, J.; Kesner, P.; Metna-Laurent, M.; Duan, T.; Xu, L.; Georges, F.; Koehl, M.;

597

Abrous, D.N.; Mendizabal-Zubiaga, J.; Grandes, P.; Liu, Q.; Bai, G.; Wang, W.; Xiong, L.; Ren,

598

W.; Marsicano, G.; Zhang, X. Acute cannabinoids impair working memory through astroglial

599

CB1 receptor modulation of hippocampal LTD. Cell. 2012, 148(5):1039-50.

600

(38)

601

Imperato, A.; Pedrazzini, T.; Roques, B.P.; Vassart, G.; Fratta, W.; Parmentier, M.

602

Unresponsiveness to cannabinoids and reduced addictive effects of opiates in CB1 receptor

603

knockout mice. Science. 1999, 283(5400):401-404.

604

(39)

605

children's development (4): Brain-behavior phenotypying of genetically-engineered mice using a

606

comprehensive behavioral test battery on research of neuropsychiatric disorders. J. Toxicol. Sci.

607

2009;34 Suppl 2:SP293-305.

608

(40)

609

influences

610

marmoratus). Neurotoxicology 2012, 33(4): 758-762.

611

(41)

612

Garner, C.C.; Shenolikar, S.; Ehlers, M.D. Neurabin/protein phosphatase-1 complex regulates

613

dendritic spine morphogenesis and maturation. Mol. Biol. Cell. 2005, 16(5):2349-2362.

614

(42)

615

converge on a common protein synthesis-dependent long-term depression that is affected in

616

fragile X syndrome mental retardation. J. Neurosci. 2007, 27(43):11624-11634.

Ledent, C.; Valverde, O.; Cossu, G.; Petitet, F.; Aubert, J.F.; Beslot, F.; Böhme, G.A.;

Takao, K.; Miyakawa, T. 2009. Intrauterine environment-genome interaction and

He, C.; Wang, C.; Zhou, Y.; Li, J.; Zuo, Z. Embryonic exposure to benzo (a) pyrene neural

development

and

function

in

rockfish

(Sebastiscus

Terry-Lorenzo, R.T.; Roadcap, D.W.; Otsuka, T.; Blanpied, T.A.; Zamorano, P.L.;

Volk, L.J.; Pfeiffer, B.E.; Gibson, J.R.; Huber, K.M. Multiple Gq-coupled receptors

28 ACS Paragon Plus Environment

Page 29 of 37

Environmental Science & Technology

617

(43)

Fitzgerald, P.J.; Barkus, C.; Feyder, M.; Wiedholz, L.M.; Chen, Y.C.; Karlsson, R.M.;

618

Machado-Vieira, R.; Graybeal, C.; Sharp, T.; Zarate, C.; Harvey-White, J.; Du, J.; Sprengel, R.;

619

Gass, P.; Bannerman, D.; Holmes, A. Does gene deletion of AMPA GluA1 phenocopy features

620

of schizoaffective disorder? Neurobiol. Dis. 2010, 40(3):608-21.

621

(44)

622

Esbaugh, A.J.; Mager, E.M.; Stieglitz, J.; Hoenig, R.; Benetti D. Corresponding morphological

623

and molecular indicators of crude oil toxicity to the developing hearts of mahi mahi. Sci. Rep.

624

2015, 5, 17326.

625

(45)

626

teratogenesis of polycyclic aromatic hydrocarbons and PCB-126 in Atlantic killifish (Fundulus

627

heteroclitus). Aquat. Toxicol. 2010, 99, 232–240

628

(46)

629

ring polycyclic aromatic hydrocarbons in zebrafish is differentially dependent on AH receptor

630

isoforms and hepatic cytochrome P4501A metabolism. Toxicol. Appl. Pharmacol. 2006, 217,

631

308–321.

632

(47)

633

Tanguay, R.L. Ligand-specific transcriptional mechanisms underlie aryl hydrocarbon receptor-

634

mediated developmental toxicity of oxygenated PAHs. Toxicol. Sci. 2015, 147(2): 397–411.

635

(48)

636

toxicity and XRE- and ARE-associated gene induction in zebrafish (Danio rerio). Toxicol. Appl.

637

Pharmacol. 2011, 254, 280–287

Edmunds, R.C.; Gill, J.A.; Baldwin, D.H.; Linbo, T.L.; French, B.L.; Brown, T.L.;

Clark, B. W.; Matson, C. W.; Jung, D.; Di Giulio, R. T. AHR2 mediates cardiac

Incardona, J. P.; Day, H. L.; Collier, T. K.; Scholz, N. L. Developmental toxicity of 4-

Goodale, B.C.; La Du, J.; Tilton, S.C.; Sullivan, C.M.; Bisson, W.H.; Waters, K.M.;

Van Tiem, L. A.; Di Giulio, R. T. AHR2 knockdown prevents PAH-mediated cardiac

29 ACS Paragon Plus Environment

Environmental Science & Technology

Page 30 of 37

638

(49)

Incardona, J.P.; Carls, M.G.; Teraoka, H.; Sloan, C.A.; Collier, T.K.; Scholz, N.L. Aryl

639

hydrocarbon receptor-independent toxicity of weathered crude oil during fish development.

640

Environ. Health Perspect. 2005, 113, 1755–1762.

641

(50)

642

J.P. Geologically distinct crude oils cause a common cardiotoxicity syndrome in developing

643

zebrafish. Chemosphere 2013, 91,1146–1155.

644

(51)

645

malformation is an early response to TCDD in embryonic zebrafish. Toxicol. Sci. 2005, 84, 368–

646

377.

647

(52)

648

NFATc1 mediates vascular endothelial growth factor-induced proliferation of human pulmonary

649

valve endothelial cells. J. Biol. Chem. 2003, 278(3), 1686-1692.

650

(53)

651

in coordination with the MAPKs. Cardiovasc. Res. 2004, 63(3), 467-475.

652

(54)

653

M.C.; Valentine, D.L.; Reddy, C.M. Oil weathering after the Deepwater Horizon disaster led to

654

the formation of oxygenated residues. Environ. Sci. Technol. 2012, 46(16), 8799-8807.

655

(55)

656

Anderson, K.A.; Waters, K.M.; Tanguay, R.L. Comparative developmental toxicity of

657

environmentally relevant oxygenated PAHs. Toxicol. Appl. Pharmacol. 2013, 271(2), 266-275.

Jung, J.H.; Hicken, C.E.; Boyd, D.; Anulacion, B.F.; Carls, M.G.; Shim, W.J.; Incardona,

Antkiewicz, D. S.; Burns, C. G.; Carney, S. A.; Peterson, R. E.; Heideman, W. Heart

Johnson, E.N.; Lee, Y.M.; Sander, T.L.; Rabkin, E.; Schoen, F.J.; Kaushal, S.; Bischoff, J.

Molkentin, J.D. Calcineurin–NFAT signaling regulates the cardiac hypertrophic response

Aeppli, C.; Carmichael, C.A.; Nelson, R.K.; Lemkau, K.L.; Graham, W.M.; Redmond,

Knecht, A.L.; Goodale, B.C.; Truong, L.; Simonich, M.T.; Swanson, A.J.; Matzke, M.M.;

30 ACS Paragon Plus Environment

Page 31 of 37

Environmental Science & Technology

658

(56)

Hodgkin, A.L.; McNaughton, P.A.; Nunn, B.J. The ionic selectivity and calcium

659

dependence of the light-sensitive pathway in toad rods. J. Physiol. 1985, 358, 447.

660

(57)

661

extended family with diverse functions. Annu. Rev. Physiol. 1996, 58(1), 395-426.

662

(58)

663

coupling of a Ca2+/calmodulin-dependent nitric oxide synthase and a soluble guanylyl cyclase in

664

vertebrate photoreceptor cells. EMBO J. 1994, 13(14), 3312.

665

(59)

666

equivalent-light hypothesis from arrestin knockout mice. Invest. Ophthalmol. Vis. Sci. 1999,

667

40(12), 2770-2772.

668

(60)

669

induced photoreceptor apoptosis and neuroprotection for retinal degeneration. Prog. Retin. Eye

670

Res. 2005, 24(2), 275-306.

671

(61)

672

Kitamura, S.I.; Takata, H.; Handoh, I.C.; Nakayama, K.; Murakami, Y. Nervous system

673

disruption and concomitant behavioral abnormality in early hatched pufferfish larvae exposed to

674

heavy oil. Environ. Sci. Pollut. Res. 2012, 19(7): 2488-2497.

675

(62)

676

developmental toxicity in zebrafish embryos and the possible mechanisms involved. J. Hazard

677

Mater. 2013, 261: 172-180.

Finn, J.T.; Grunwald, M.E.; Yau, K.W. Cyclic nucleotide-gated ion channels: an

Koch, K.W.; Lambrecht, H.G.; Haberecht, M.; Redburn, D.; Schmidt, H.H. Functional

Fain, G.L.; Lisman, J.E. Light, Ca2+, and photoreceptor death: new evidence for the

Wenzel, A.; Grimm, C.; Samardzija, M.; Remé, C.E. Molecular mechanisms of light-

Kawaguchi, M.; Sugahara, Y.; Watanabe, T.; Irie, K.; Ishida, M.; Kurokawa, D.;

Huang, L.; Wang, C.; Zhang, Y.; Wu, M.; Zuo, Z. Phenanthrene causes ocular

678

31 ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 37

679

Figure legends

680

Fig. 1 Toxicity predictions determined from differentially expressed genes between controls

681

versus slick oil-treated red drum larvae for 24, 48 and 72 hpf. The left Y-axis indicates the –log

682

(p-value; red dots in upper line) determined by Fisher's exact test, and the right Y-axis indicates

683

the ratio (lower line), which is calculated as the number of genes found in each given pathway

684

divided by the total number of genes in that pathway.

685

Fig. 2 Activated and inhibited upstream regulators in slick oil treated red drum larvae predicted

686

by Upstream Analysis in IPA. The prediction of activation is based on the global direction of

687

changes of the modulated genes. The activation Z-score indicates whether the observed gene

688

responses to upstream regulators agree with expected changes derived from the literature that

689

accrued in IPA database, was used to predict the activation state. Z-scores ≥ 2 or ≤ − 2 indicates

690

that the upstream regulator was predicted to be activated or inhibited, respectively. A Fisher’s

691

Exact Test was used to determine the significance of the overlap between the regulator and the

692

responsive genes.

693

Fig. 3 Canonical Pathways for 48 hpf (a) and 72 hpf (b). The y-axis displays the -log of p-value

694

which is calculated by Fisher's exact test right-tailed. The numbers on the right indicate the

695

number of genes affiliated with that pathway. Red bar, number of up-regulated genes; green bar:

696

down-regulated; gray bar, unchanged genes.

697

Fig. 4 Predicted mechanisms through Ingenuity Pathway Analysis showing how slick oil may

698

lead to the suppression of long term depression of synaptic transmission, decreased quantity of

699

nervous tissue, and hyperactive behavior.

700

32 ACS Paragon Plus Environment

Page 33 of 37

Environmental Science & Technology

701 702

Fig. 1 33 ACS Paragon Plus Environment

Environmental Science & Technology

703 704

Fig. 2

705

34 ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37

Environmental Science & Technology

706

707

Fig. 3

708

709

710

35 ACS Paragon Plus Environment

Environmental Science & Technology

711 712

Fig. 4

713

36 ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37

Environmental Science & Technology

714

Table 1 Morphological measurements: brain area, eye area, pericardial area, Iris area, spine and

715

total body length in control and oil treated red drum larvae at 48 hpf. Data were analyzed by t-

716

test and a two tailed p- value indicates the significant effect of slick oil exposure at p < 0.05 (N =

717

4). Morphological

Control

Slick oil

p-value

Brain Area (µm2)

15057 ± 429

12556 ± 828

0.029

Eye Area (µm2)

13138 ± 282

10296 ± 788

0.001

Iris Area (µm2)

1703 ± 63

1370 ± 100

0.090

Pericardial Area (µm2)

6444 ± 464

9343 ± 1200