Toxicogenomic Assessment of 6-OH-BDE47-Induced Developmental

Oct 17, 2016 - Using the averaged fold change values of 437 whole DEGs at each measured concentration, a gamma model was fitted to draw dose–respons...
2 downloads 6 Views 1MB Size
Subscriber access provided by BOSTON UNIV

Article

Toxicogenomic Assessment of 6-OH-BDE47 Induced Developmental Toxicity in Chicken Embryo Ying Peng, Pu Xia, Junjiang Zhang, Daniel L. Villeneuve, Jiamin Zhang, Zhihao Wang, Si Wei, Hongxia Yu, and Xiaowei Zhang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b04467 • Publication Date (Web): 17 Oct 2016 Downloaded from http://pubs.acs.org on October 19, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 45

Environmental Science & Technology

1

Toxicogenomic Assessment of 6-OH-BDE47 Induced

2

Developmental Toxicity in Chicken Embryo

3 4

Ying Peng†#, Pu Xia†#, Junjiang Zhang†, Daniel L. Villeneuve‡, Jiamin Zhang†,

5

Zhihao Wang†, Si Wei†, Hongxia Yu†, Xiaowei Zhang†*

6

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

7

Environment, Nanjing University, Nanjing, 210023, PR China

8



United States Environmental Protection Agency, Mid-Continent Ecology Division,

9

Duluth, MN, 55804, USA.

10

# These authors contributed equally to this paper

11

*Correspondence:

12

Xiaowei Zhang, PhD, Prof

13

School of the Environment

14

Nanjing University

15

Nanjing, 210089, China;

16

Tel.: 86-25-89680623

17

Fax: 86-25-89680623

18

E-mail: [email protected]

19

[email protected]

20 21

ACS Paragon Plus Environment

Environmental Science & Technology

22

Abstract:

23

Hydroxylated polybrominated diphenyl ethers (OH-PBDEs) are analogs of PBDEs

24

with hundreds of possible structures and are frequently detected in the environment.

25

However, the in vivo evidence on the toxicity of OH-PBDEs is still very limited. Here

26

the developmental toxicity of 6-OH-BDE47, a predominant congener of OH-PBDEs

27

detected in the environment, in chicken embryos was assessed using a toxicogenomic

28

approach. Fertilized chicken eggs were dosed via in ovo administration of 0.006 to

29

0.474 nmol 6-OH-BDE47/g egg followed by 18-days incubation. Significant embryo

30

lethality (LD50 = 1.940 nmol/g egg) and increased hepatic somatic index (HSI) were

31

caused by 6-OH-BDE47 exposure. The functional enrichment of differentially

32

expressed genes (DEGs) was associated with oxidative phosphorylation, generation of

33

precursor metabolites and energy, and electron transport chain, which suggest that

34

6-OH-BDE47 exposure may disrupt the embryo development by altering the function

35

of energy production in mitochondria. Moreover, aryl hydrocarbon receptor (AhR)

36

mediated responses including up-regulation of CYP1A4 was observed in the livers of

37

embryos exposed to 6-OH-BDE47. Overall, this study confirmed the embryo lethality

38

by 6-OH-BDE47 and further improved the mechanistic understanding on OH-PBDEs

39

caused toxicity. Ecological risk assessment via application of both no observed effect

40

level (NOEL) and the sensitive NOTEL (transcriptional NOEL) suggested that

41

OH-PBDEs might cause ecological risk to wild birds.

42

ACS Paragon Plus Environment

Page 2 of 45

Page 3 of 45

Environmental Science & Technology

43

Introduction

44

Hydroxylated

and

methoxylated

polybrominated

diphenyl

ethers

45

(OH-/MeO-PBDEs) are analogs of PBDEs which have been used as flame retardants

46

and are widely distributed in the environment.1 OH-/MeO-PBDEs are chemicals with

47

hundreds of possible structures and many of them have been detected in abiotic

48

matrices or biota including surface water, fish, birds, polar bear, marine mammals and

49

human blood.2-6 Many in vitro studies have shown that some OH-PBDEs are more

50

potent than their postulated precursor PBDEs and corresponding MeO-PBDEs for

51

several toxicologically relevant molecular endpoints, including binding to steroid

52

hormone receptors,7-10 toxicity to neural cells,11 cytotoxic effects,12 and alteration of

53

DNA repair pathways.13 However, the in vivo evidence concerning the toxicity of

54

OH-/MeO-PBDEs remains limited.

55

6-OH-BDE47 was one of the predominant congeners of OH-PBDEs detected in

56

biota.14-16 Our previous study using an in vitro AhR- luciferase reporter gene (LRG)

57

test suggested that OH-/MeO-PBDEs can activate aryl hydrocarbon receptor (AhR) in

58

both mammalian and avian species, and 6-OH-BDE47 was one of the most potent

59

congeners.17,

60

showed that this chemical can also modulate other molecular pathways in zebrafish19,

61

20

62

molecular-level pathways were impacted following activation of AhR and potentially

63

other targets, by 6-OH-BDE47 in the embryo development remains unclear.

18

Although 6-OH-BDE47 is a known AhR activator,17,

18

evidence

and demonstrated thyroid hormone activity in in vitro Yeast hybrid assays.21 How

ACS Paragon Plus Environment

Environmental Science & Technology

64

Toxicogenomic approaches which provide global view of molecular pathways that

65

might be disturbed by the chemical exposure have been widely used to help elucidate

66

molecular mechanisms of chemical toxicity. In particular, RNA-seq-based

67

transcriptomics has increasingly been used to study gene transcript expression due to

68

its much broader genome coverage and higher sample throughput compared with the

69

previous microarray technologies.

70

The primary goal of this study was to test the hypothesis that 6-OH-BDE47 could

71

cause embryo toxicity in avian species. The developmental toxicity of 6-OH-BDE47

72

in chicken (Gallus gallus) was tested by in ovo administration. The hepatic

73

distribution of 6-OH-BDE47 was confirmed by instrumental analysis and the hepatic

74

transcriptomic profile was evaluated to examine the underlying molecular response in

75

6-OH-BDE47 induced embryo toxicity. Furthermore, we also assessed the potential

76

ecological risk of 6-OH-BDE47 on wild birds.

77

78

79

2. Materials & method

2.1 Chemicals and solutions

80

6-MeO-BDE4, 6-OH-BDE47 and other 10 MeO-PBDEs, 10 OH-PBDEs were

81

synthesized at City University of Hong Kong (purities > 98%).22 6-OH-BDE82,

82

6-MeO-BDE-82, 6-MeO-BDE85, 6-MeO-BDE87 were purchased from AccuStandard

83

(New Haven, CT, USA). 6-OH-BDE47 was dissolved in dimethyl sulfoxide (DMSO,

ACS Paragon Plus Environment

Page 4 of 45

Page 5 of 45

Environmental Science & Technology

84

Generay Biotech, Shanghai, China) to prepare stock solutions and then diluted in

85

DMSO to the desired test concentrations when injected to the chicken embryo.

86

2.2 Procedures of egg preparation, injection and incubation

87

Fertilized chicken eggs were purchased from local farm with no potential dioxin

88

exposure history in Nanjing, Jiangsu Province, China. In general, preparation of

89

injection solutions and egg injection procedures follow methodology described in

90

Powell et al. (1996)23 with minor modifications. Eggs were stored in a cooler for no

91

longer than one week at 13-15 °C until 24 h before injection. Eggs were weighed to

92

the nearest 0.1 g and then held to a bright light to detect subtle damage to the shell.

93

Undamaged eggs with mean weights of 42.08 ± 3.62 g (± SD) had the center of their

94

air cells marked with pencil to outline the injection site. A pipette with sterile

95

aciculiform tip (Eppendoff, Germany) was inserted horizontally through the air cell

96

after filing to thin the egg shell. Five doses (0.474, 0.158, 0.053, 0.018, 0.006 nmol

97

6-OH-BDE47 /g egg) with DMSO solution (5 µL) were injected into chamber of

98

chicken eggs (n=10 eggs per treatment). The air cell was chosen as the site of

99

injection because of ease and speed of delivery of the chemical into the egg.24 The site

100

of injection was then sealed using liquid paraffin wax. Controls included non-injected

101

and solvent vehicle control eggs. Concentration of DMSO in final test solutions did

102

not exceed 0.01%.

103

The ten eggs from each treatment were incubated in an automatic digital display

104

incubator (Wansheng, Nanjing, JS, China) until hatch or 18 d. After 18 d incubation,

ACS Paragon Plus Environment

Environmental Science & Technology

105

any eggs which had not hatched were opened and dead embryos were recorded.

106

Pipped embryos were removed and euthanized by decapitation. For all embryos that

107

pipped, the following measurements were recorded: embryo mass, liver mass, tarsus

108

length, head plus bill length (head + bill: back of head to tip of bill). Liver samples

109

were collected and quickly frozen with liquid nitrogen and four liver samples from

110

each treatment group were picked at random for further analysis. Each liver sample

111

was split with half of sample allocated for RNA extraction, subsequent RNA-seq and

112

qRT-PCR analysis (n=4). The other half of each liver was allocated for chemical

113

analysis and was stored at -80°C until extracted.

114

2.3 Chemical Analysis

115

Liver samples (approximately 0.3 g, ww) were mixed with anhydrous sodium

116

sulfate (5 times of liver sample weight). A surrogate recovery standard was added

117

(6-OH-BDE82 and 6-MeO-BDE85) and then samples were extracted with 50%

118

dichlormethane (DCM)/hexane (HEX) using an accelerated solvent extraction system

119

(Dionex ASE 200). The column extraction eluent was concentrated to 10 mL and 1

120

mL (10% of the sample) was removed for gravimetric lipid determination, the

121

remaining volume of extract was dissolved with HEX and then partitioned with

122

potassium hydroxide (0.5M in 50% ethanol). The neutral fraction used for BDEs and

123

MeO-BDEs analysis was treated with 2 mL concentrated sulfuric acid to remove

124

lipids, then cleaned on a multilayer silica gel column (inner diameter 0.8 cm, packed

125

with 8 cm neutral silica (3% deactivated, w/w) and 8 cm sulfuric acid silica (2:1

ACS Paragon Plus Environment

Page 6 of 45

Page 7 of 45

Environmental Science & Technology

126

w/w)) and eluted with 35 mL DCM: HEX (1:1, v/v). Finally, the extract was blown to

127

near-dryness under a gentle stream of nitrogen, and 10 µL BDE-77 and

128

6-MeO-BDE87 were added respectively as internal standards and made up to 100 µL

129

prior to GC/MS analysis. The alkaline phase used for OH-BDE analysis was acidified

130

to pH < 2 with hydrochloric acid , then extracted two times with HEX: MTBE (1:1,

131

v/v, 6 mL and 3 mL, respectively) and dried by anhydrous sodium sulfate. Samples

132

were derivatized with diazomethane overnight to methylate the hydroxyl group. After

133

methylation, the samples were purified on a silica gel column (inner diameter 0.8 cm)

134

packed with 8 cm neutral silica (3% deactivated, w/w) and 8 cm sulfuric acid silica

135

(2:1 w/w) and eluted with 35 mL DCM: HEX (1:1, v/v). Each extract was blown to

136

near dryness and reconstituted in 50 mL iso-octane. 6-MeO-BDE87 was added as an

137

internal standard for GC/MS analysis. The analysis of MeO-PBDEs and BDEs was

138

performed using a TSQ Quantum GC/MS (Thermo Scientific, USA) coupled with

139

electron capture negative ionization (ECNI) source in the selected ion monitoring

140

mode (SIM). Concentration and potential biotransformation of 10 OH-PBDEs, 10

141

MeO-PBDEs were also determined. Detailed protocols for instrumental conditions

142

and quality assurance and quality control (QA/QC) are provided in the Supporting

143

Information.

144

2.4 RNA Sequencing

145

To analyze transcriptomic response in the liver of embryos after exposure to

146

6-OH-BDE47, samples from the 0.01% DMSO vehicle treatment group, 0.018 nmol/g

ACS Paragon Plus Environment

Environmental Science & Technology

147

egg, 0.052 nmol/g egg, 0.158 nmol/g egg and 0.474 nmol/g egg treatment groups

148

were used for RNA-seq. Total RNA of chicken livers was isolated by RNeasy mini

149

kit, followed by digestion of genomic DNA (QIANGEN, GmbH, Hilden). Total RNA

150

samples were stored at -80°C. The concentrations of total RNA were measured by

151

Synergy H4 Hybrid Take3 reader (BioTek Instruments, Winooski, VT), followed by

152

determination of the RNA quality using Agilent 2100 bioanalyzer (Agilent

153

technologies, Santa Clara, CA, US). 16 total RNA samples, including 4 control

154

samples and 12 treatment samples (N=4 for 0.018 nmol/g and 0.158 nmol/g,

155

respectively; N=2 for 0.052 nmol/g and 0.474 nmol/g, respectively), were submit to

156

the libraries construction using Dynabeads mRNA DIRECT Micro Kit (Life

157

technologies, AS, Oslo, Norway) and Ion Total RNA-Seq Kit v2 (Life technologies,

158

Austin). The prepared libraries were then submitted to high throughput sequencing on

159

Ion Torrent Proton (Life technologies, Carlsbad, CA, USA) following manufacturer’s

160

protocol.

161

2.5 RNA-seq Data Analysis

162

Raw sequences were mapped to Gallus gallus 4.75 (Ensembl Genes 8325) using the

163

Torrent Mapping Alignment Program (TMAP). The aligned files were submitted to R

164

package, Genomic Alignments (Bioconductor version 1.6.3)26 for gene expression

165

analysis. Normalization of gene expression was performed to account for library

166

preparation effects. The normalized expressions of each gene across control to

ACS Paragon Plus Environment

Page 8 of 45

Page 9 of 45

Environmental Science & Technology

167

treatment groups (measured concentrations) were submitted to analysis of

168

differentially expressed genes (DEGs).

169

Two approaches were used to identify DEGs according to the measured hepatic

170

chemical concentrations. First, linear regression analysis on expression of genes was

171

conducted using basic function ‘lm’ in R.27 The p-values of linear regression of each

172

gene were submitted to multiple hypothesis test (Benjamini & Hochberg method) to

173

get adjusted-p value. Significant linear regression genes (adjusted-p value