Whole-Life-Stage Characterization in the Basic Biology of Daphnia

Nov 8, 2017 - Animal Maintenance and Growth, Reproduction, and Survival Curves ... Additionally, the number of dead D. magna in each beaker was ...
0 downloads 0 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Whole-Life-Stage Characterization in the Basic Biology of Daphnia magna and Effects of TDCIPP on Growth, Reproduction, Survival and Transcription of Genes Han Li, Siliang Yuan, Guanyong Su, Meng Li, Qiangwei Wang, Guonian Zhu, Robert J. Letcher, Yufei Li, Zhihua Han, and Chunsheng Liu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04569 • Publication Date (Web): 08 Nov 2017 Downloaded from http://pubs.acs.org on November 9, 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 34

Environmental Science & Technology

1

Whole-Life-Stage Characterization in the Basic Biology of Daphnia

2

magna and Effects of TDCIPP on Growth, Reproduction, Survival

3

and Transcription of Genes

4 5

Han Li¶, Siliang Yuan¶,ф, Guanyong Suξ, Meng Li§, Qiangwei Wang§, Guonian

6

Zhu§, Robert J. Letcher‡, Yufei LiЖ, Zhihua Han$, Chunsheng Liu¶,†,*

7 8



9



College of Fisheries, Huazhong Agricultural University, Wuhan 430070, China Collaborative Innovation Centre for Efficient and Health Production of Fisheries in

10

Hunan Province, Changde 415000, China

11

ф

12

Control and Prevention, Wuhan 430070, China

13

ξ

14

of Environmental and Biological Engineering, Nanjing University of Science and

15

Technology, Nanjing 210094, China

16

§

17

310058, China

18



19

5B6, Canada

20

Ж

21

of PR China, Beijing 100045, China

22

$

Hubei Engineering Technology Research Center for Aquatic Animal Diseases

Jiangsu Key Laboratory of Chemical Pollution Control and Resources Reuse, School

Institute of Pesticide and Environmental Toxicology, Zhejiang University, Hangzhou

Departments of Chemistry and Biology, Carleton University, Ottawa, Ontario K1S

China Rural Technology Development Centre, Ministry of Science and Technology

Nanjing Institute of Environmental Science, MEP, Nanjing 210042, Jiangsu, China

23 24

*Author for correspondence:

25

Tel: 86 27 87282113

26

Fax: 86 27 87282114.

27

Email: [email protected] or [email protected] (C. Liu) 1

ACS Paragon Plus Environment

Environmental Science & Technology

28

Page 2 of 34

ABSTRACT

29

Toxicity tests of chemicals have mainly focused on the partial life-cycle

30

evaluation of model animals. Limited information is available for the evaluation of

31

effects of chemicals from a whole-life-stage exposure perspective. The objective of

32

this study was to perform a whole-life-stage characterization in the basic biology of

33

Daphnia magna (D. magna) and evaluate the effects of a known organophosphate

34

ester (OPE) contaminant, tris (1,3-dichloro-2-propyl) phosphate (TDCIPP), on growth,

35

reproduction,

36

characterization in growth, reproduction and survival of D. magna was conducted,

37

and representative sampling time points for the three developmental stages were

38

identified (day 6, day 32 and day 62). Transcriptomic profiles for these three stages

39

were compared and stage-specific PCR arrays of D. magna were developed. The

40

whole-life-stage exposure to environmentally relevant or greater concentrations of

41

TDCIPP significantly inhibited growth and reproduction of D. magna and decreased

42

survival at the later stage of the exposure experiment (≥ 32 days). Such adverse

43

effects were not observed in the early stage of the exposure (< 32 days), suggesting

44

that short-term toxicity tests, such as the standard 21-day test, might underestimate the

45

environmental risk of TDCIPP. Furthermore, expressions of genes selected at day 6,

46

day 32 and day 62 were significantly changed after TDCIPP exposure, and the

47

changes in the expressions of partial genes were correlated to the inhibitory effects on

48

growth, reproduction and survival.

survival

and

transcription

of

genes.

49 50 51

2

ACS Paragon Plus Environment

The

whole-life-stage

Page 3 of 34

Environmental Science & Technology

52

INTRODUCTION

53

Previous studies that have assessed the toxicity of exogenous chemicals, were

54

mainly conducted based on partial life-cycle evaluations of selected model animals.1-3

55

For example, experiments in 48-h acute and 21-day chronic exposure studies of

56

flowback and produced wastewater (FPW) experiments were conducted using

57

Daphnia magna (D. magna) to characterize toxicity of FPW.4 In another study,

58

medaka fish (Oryzias latipes) were selected and a 28-day chronic exposure was

59

performed to study cellular and molecular effects of Microcystin-LR on the liver.5

60

Zebrafish were exposed to the organophosphate ester (OPE), triphenyl phosphate, for

61

11-14 days to study the bioconcentration and tissue distribution of the chemical, thus

62

evaluating its fate and potential toxicity.6 Collectively, these studies centralized the

63

toxic effects to a certain life stage of a test organism although they usually had two

64

typical shortcomings. Firstly, in these tests, the lowest observed effect concentrations

65

of chemicals were relatively high and mostly not environmentally relevant due to the

66

short exposure period. Therefore, it was difficult to accurately predict risks of

67

chemicals at environmentally relevant concentrations. Secondly, in these standard

68

toxicity tests, the evaluation of toxic effects of chemicals was performed in a partial

69

life-cycle exposure period, and thus underestimated the environmental impacts of the

70

chemical since in natural environments some wildlife might be exposed to some

71

chemical pollutants throughout their life span and especially for those organisms with

72

short life spans such as D. magna.

73

To address this knowledge gap, in the present study a whole-life-stage toxicity 3

ACS Paragon Plus Environment

Environmental Science & Technology

74

test protocol with D. magna was developed and a exposure experiment was performed

75

using the known OPE environmental contaminant, tris(1,3-dichloro-2-propyl)

76

phosphate (TDCIPP) . D. magna is a model organism in ecotoxicology studies

77

because it is an imperative animal link in the aquatic food chain.7,8 Furthermore, it has

78

a short life span, and the whole-life-stage exposure of chemical pollutants in natural

79

aquatic environment might be a frequent occurrence. Therefore, D. magna may be a

80

suitable test animal for evaluating the toxic effects of chemicals after a

81

whole-life-stage exposure. TDCIPP is a current-use flame retardant which has been

82

added to a diverse array of commercial products.9,10 Being an additive flame retardant,

83

TDCIPP can easily be released into the surrounding environment. Environmental

84

monitoring has demonstrated that TDCIPP is extensively distributed in various

85

environmental media,11-15 and the highest concentration reported in natural waters was

86

377 ng/L.16 Previous studies reported that exposure to high or relatively low

87

concentrations of TDCIPP could cause neurological, developmental and reproductive

88

damage,17-27 however, to the best of our knowledge, all these studies were performed

89

based on partial life stage tests with organisms. It remains unknown whether

90

whole-life-stage exposure could cause greater environmental hazards compared with

91

those of partial life stage tests.

92

The objective of this study was to perform a whole-life-stage characterization in

93

the basic biology of D. magna and evaluate the effects of a known organophosphate

94

ester (OPE) contaminant, TDCIPP, on growth, reproduction, survival and

95

transcription of genes. Briefly, in this study, the whole-life-stage growth, reproduction 4

ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34

Environmental Science & Technology

96

and survival curves of D. magna were determined. Appropriate sampling time points

97

for toxicity testing of the chemical TDCIPP were also determined. RNA-Seq

98

technology was then used to obtain stage-specific genes of D. magna for PCR array

99

development. Finally, TDCIPP was selected to perform a whole-life-stage exposure

100

experiment and its effects at environmentally relevant or greater concentrations on

101

growth, reproduction, survival and transcription of genes were evaluated.

102

MATERIALS AND METHODS

103

Chemicals and Reagents

104

Tris (1,3-dichloro-2-propyl) phosphate (TDCIPP) was purchased from Sigma (St.

105

Louis, MO, USA), and it was diluted with dimethyl sulfoxide (DMSO) using a serial

106

dilution approach. TRIzol reagent, reverse transcription and SYBR Green kits were

107

purchased from Takara (Dalian, Liaoning, China). d15-TDCIPP was purchased from

108

Toronto Research Chemical (Toronto, Canada). All the other reagents used in this

109

study were of analytical grade.

110

Normal Culture of D. magna

111

Animal maintenance and growth, reproduction and survival curves. D. magna

112

was cultured as has been reported elsewhere.9 Culture medium was aerated tap-water,

113

which was pretreated by an Angel water cleaner (J1205-ROB12, Shenzhen, China).

114

The culture medium was changed completely every two days. To get a comprehensive

115

insight into its whole-life-stage characteristics, the growth, reproduction and survival

116

curves of D. magna were determined. Here, two experiments were run in parallel. In

117

the first experiment, a growth curve was produced. Briefly, three glass beakers 5

ACS Paragon Plus Environment

Environmental Science & Technology

118

containing 1.5 L aerated tap-water were used, and each beaker contained 75 neonates

119

(< 24 h) to ensure 20 mL per individual. The experiment was performed using a

120

well-characterized single clone of D. magna, maintained indefinitely as pure

121

parthenogenetic cultures as reported before.28 The culture media were renewed

122

completely every two days and body length (from the apex of the helmet to the base

123

of the tail spine) was measured every two or four days depending on the

124

developmental stages of D. magna. This step was repeated until all the animals

125

expired. Body length data were used to produce the growth curve. In the second

126

experiment, reproductive and survival curves were obtained. Another three glass

127

beakers containing 1 L culture medium and 50 neonates (< 24 h) were used to ensure

128

20 mL per individual. During this experiment period, the number of offspring

129

produced in each beaker was counted every day to generate the reproductive curve.

130

Additionally, the number of dead D. magna in each beaker was recorded and a

131

corresponding survival curve was generated. Similarly, this step was repeated until all

132

the individuals expired.

133

Transcriptomic sequencing. Based on growth, reproductive and survival curves,

134

animals were sampled at day 0 (< 24 h), day 6, day 32 and day 62 for transcriptomic

135

sequencing. The four time points (day 0, day 6, day 32 and day 62) were considered

136

to be representatives of the stages of birth, growth, reproduction and death,

137

respectively. Briefly, three glass beakers containing 2 L aerated tap-water were used,

138

and each beaker contained 100 neonates (< 24 h) to ensure 20 mL per individual. D.

139

magna were cultured as above and 50, 30, 10 and 10 individuals from each beaker 6

ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

Environmental Science & Technology

140

were collected and pooled to generate one sample at day 0, day 6, day 32 and day 62,

141

respectively, for transcriptomic sequencing. Transcriptomic sequencing was

142

performed by Novogene Bioinformatics Technology Co., Ltd (Beijing, China). Briefly,

143

total RNA was isolated by use of TRIzol reagent and the purity, concentration and

144

integrity of RNA were determined. Once RNA samples were qualified, Oligo-dT

145

beads were used for the isolation of mRNA. After that, the isolated mRNA was

146

fragmented into short sections and used for the synthesis of first-strand cDNA. The

147

synthesis of second-strand cDNA was conducted by use of dNTPs, DNA polymerase I

148

and RNAase H. The purification of double-stranded cDNA was performed using

149

AMPure XP beads, and then the purified double-stranded cDNA were used for end

150

reparation, “A” base addition and ligation of sequencing adapters. After size selection

151

with AMPure XP beads, PCR amplification was performed and the production was

152

purified using AMPure XP beads for library construction. Finally, the library was

153

sequenced as 2 × 150 bp paired-end reads on Illumina HiSeqTM P150 sequencer

154

(Illumina, San Diego, CA, USA). Raw data were cleaned by removing reads with

155

adaptors, low quality (> 50 %) or high-proportion unknown bases (> 10 %) and then

156

were assembled by use of Trinity methodology.29,30 The annotation of transcriptome

157

sequences were performed using the following seven databases: Nr (NCBI

158

non-redundant protein sequences), Nt (NCBI nucleotide sequences), Pfam (Protein

159

family), KOG/COG (eukaryotic Ortholog Groups and Clusters of Orthologous Groups

160

of proteins), Swiss-Prot (A manually annotated and reviewed protein sequence

161

database), KEGG (Kyoto Encyclopedia of Genes and Genomes) and GO (Gene 7

ACS Paragon Plus Environment

Environmental Science & Technology

162

Ontology). Gene expression was calculated using FPKM (expected number of

163

fragments per kilobase of transcript sequence per millions base pairs sequenced), and

164

differentially expressed genes between day 0 and day 6, day 6 and day 32, and day 32

165

and day 62 were selected based on padj < 0.05. KEGG pathway analysis was

166

conducted and corrected P value (FDR) cut-off was set at 0.05.

167

Development of stage-specific PCR arrays. The life process of D. magna was

168

artificially divided into four stages: birth stage, growth stage, reproduction stage and

169

death stage. It was hypothesized that, compared with the last stage, when entering the

170

next stage, expressions of some genes would be up-regulated to fulfill the biological

171

function of the next stage. Therefore, in this study, up-regulated genes between two

172

continuous life stages were screened for the development of stage-specific PCR arrays.

173

Similar hypotheses were frequently used in previous studies to screen stage-specific

174

functional genes.31,32 For the development of the present stage-specific PCR arrays,

175

three steps were included. Firstly, differently expressed genes between day 0 and day

176

6, day 6 and day 32, and day 32 and day 62 were obtained, and those up-regulated

177

genes screened in each stage were used for KEGG pathway analysis; Secondly, those

178

significantly enriched pathways obtained in KEGG pathway analysis were picked out

179

and ranged based on the value of rich factor (RF); Thirdly, the top 5 up-regulated

180

pathways in each stage were selected, and in each pathway the top 5 genes (readcount

181

≥ 10) based on expression of fold change were used for the development of

182

stage-specific PCR arrays. The gene list for each of the stage-specific arrays was

183

provided in Table 1. 8

ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

Environmental Science & Technology

184

TDCIPP Exposure Experiments

185

TDCIPP exposure protocol. D. magna neonates (< 24 h) were exposed to 0, 300 or

186

3000 ng/L TDCIPP in glass beakers until 90 days when survival rate of control group

187

was < 50 %, and each beaker contained 1 L exposure solution and 50 animals and

188

each concentration contained 3 replicate beakers. The exposure concentrations (300

189

and 3000 ng/L) were selected based on a previous study.9 The 300 ng/L was an

190

environmentally relevant concentration and 3000 ng/L was used a greater

191

concentration. During the exposure period, the exposure media were renewed

192

completely with the corresponding concentrations of solutions every two days and

193

animals were fed everyday with the mixture of Chlorella pyrenoidosa and

194

Scenedesmus obliquus at a concentration of 2.5 × 106 cells/mL. Fecundity and

195

survival were monitored and recorded daily for calculating accumulated production of

196

offspring and survival rates, respectively. At day 6, day 32 and day 62, animals were

197

sampled for the measurement of body length and isolation of total RNA. For the

198

isolation of total RNA, there were three replicates for each concentration, and each

199

replicate included 15, 10 and 5 individuals at day 6, day 32 and day 62, respectively.

200

Control and exposure groups received 0.01 % aqueous DMSO.

201

Quantification of TDCIPP in exposure solutions. TDCIPP in the exposure solutions

202

was extracted and quantified according to a previous study.33 Briefly, an aliquot of 3

203

mL of exposure solution was transferred into a centrifuge tube, and then 4 mL hexane

204

and 10 µL internal standard solution (10 µg/mL in methanol), d15-TDCIPP, were

205

added. The mixture was blended and then ultrasonicated. After that, the supernatant 9

ACS Paragon Plus Environment

Environmental Science & Technology

206

was concentrated and then transferred into a clean glass tube. The extraction and

207

supernatant collection was repeated one more time. Combined with the first

208

supernatant and evaporated. Three hundred µL of methanol was added into the tube,

209

and then vortexed. The concentration multiple of the water samples was tenfold in the

210

present study. The resulting methanol solution was transferred into a liquid

211

chromatography (LC) vial for TDCIPP analysis. Quantification of TDCIPP was

212

performed using LC-tandem mass spectrometer system (LC-MS/MS) consisting of

213

5500 QTRAP MS/MS system (AB SCIEX, Singapore) and Waters ACQUITY™

214

system (Waters Corp., Milford, MA, USA). Chromatographic separation was

215

conducted using ACQUITY UPLC ® HSS T3 (1.8 µm, 2.1 mm × 100 mm, Waters),

216

and Water (A) and methanol (B) were used as mobile phases. The ESI source in

217

positive mode was as follows: ion spray (IS) voltage: 5500 V; curtain gas (CUR): 20

218

psi; nitrogen collision gas (CAD): 8 psi; nebulizer gas (GS1): 40 psi; auxiliary gas

219

(GS2): 40 psi; source temperature: 400 ℃. The instrumental and method detection

220

limit of TDCIPP was 0.3 µg/L and 0.03 µg/L, respectively.

221

Isolation of Total RNA and Quantitative Real-time PCR

222

Total RNA was extracted with TRIzol reagent (Takara, Dalian, Liaoning, China)

223

following the manufacturer’s instructions. After examination of the purity and

224

concentration of RNA, 500 ng of total RNA was used for reverse transcription using

225

the Prime ScriptTM RT reagent kits (Takara, Dalian, Liaoning, China). Real-time

226

quantitative PCR were performed with SYBR Green Premix Ex TaqII kits (Takara,

227

Dalian, Liaoning, China). In order to monitor the quantity assurance and quality 10

ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

Environmental Science & Technology

228

control (QA & QC) of PCR reactions, melting curves and standard curves were

229

utilized to check the purity of PCR productions and amplification efficiency of PCR,

230

respectively.

231

DAPPUDRAFT_311135) kept unchanged between different developmental stages or

232

after exposure to different concentrations of TDCIPP, and thus this gene was used as

233

an internal control or housekeeping reference gene. Primer sequences were designed

234

with Primer Premier 6 software (Premier Company, Canada), and all the primers used

235

in this study were tested for efficiency, and their efficiencies (95%-103%) were

236

comparable. Real-time PCR thermal cycling was set as follow: 95 ℃ for 3 min, and

237

then 40 cycles of 95 ℃ for 15 s and followed by 60 ℃ for 1 min. Expressions of genes

238

were calculated by 2−∆∆Ct method, and were presented as fold change compared with

239

the control.

240

Statistical Analyses

241

Statistical analyses were completed with Kyplot Demo 3.0 software (Tokyo, Japan).

242

Normality of variances for parameters was determined using Kolmogorov−Smirnow

243

and the homogeneity was checked using Levene’s test. One-way analysis of variance

244

and tukey’s multiple range test were used to determine significant differences of

245

parameters tested between control and exposure groups. The level of significance for

246

all statistical analyses was set at p < 0.05.

247

RESULTS AND DISCUSSION

248

Normal Culture of D. magna

249

Growth,

Expression

reproduction

of

and

gene

c9685_g1

survival

curves.

11

ACS Paragon Plus Environment

(hypothetical

The

protein

whole-life-stage

Environmental Science & Technology

250

characterization in growth, reproduction and survival of D. magna was investigated in

251

this study, and the corresponding curves are shown in Figure 1. D. magna has been a

252

model species in the study of environmental toxicology;34-37 however, its basic

253

biological information remains limited. To the best of our knowledge, we report for

254

the first time on the characterization of growth, reproduction and survival of D.

255

magna from birth to natural death. For D. magna, there was a rapid growing stage

256

from day 0 to day 12 after birth, and during this period the body length of D. magna

257

increased rapidly from 934.27 ± 44.98 µm at day 0 to 2733.47 ± 92.02 µm at day 12

258

(Figure 1a). After day 12, the body length of D. magna remained essentially constant

259

until death (Figure 1a). Therefore, the growth curve suggested that the days from 0 to

260

12 were a key period for the growth of D. magna, and thus the period might be a

261

sensitive window for the evaluation of effects of chemicals on growth. For

262

reproduction, D. magna started to produce offspring at day 10, and since then the

263

individuals kept reproducing until death. The accumulated number of offspring was

264

353.04 ± 7.85 at day 94 for each individual (Figure 1b). Previous studies reported that

265

D. magna has two reproductive modes: cyclic parthenogenesis and sexual

266

reproduction, and under favorable rearing conditions, cyclic parthenogenesis is the

267

only reproductive mode.38,39 In this study, sexual reproduction was not observed since

268

no resting eggs were produced throughout the whole-life-stage experiment, and thus

269

all the offspring produced were considered to be from cyclic parthenogenesis.

270

Additionally, our results demonstrated that the offspring of D. magna were initially

271

produced at day 10, and were consistent or comparable with results from previous 12

ACS Paragon Plus Environment

Page 12 of 34

Page 13 of 34

Environmental Science & Technology

272

studies.4,40 Furthermore, for the first time, we found that under favorable rearing

273

conditions, D. magna could maintain high reproductive capacity until the end of life.

274

This might be a unique reproductive strategy for the organisms with parthenogenesis.

275

For survival, from day 0 to day 53, no mortality was observed (Figure 1c). However,

276

survival rates of D. magna decreased from 100 % at day 0 to day 53 to 0 % at day 94

277

(Figure 1c). The survival curve results suggested that D. magna had a relatively broad

278

death period compared with its life span, and this phenomenon has also been reported

279

for other organisms such as mice.41

280

Transcriptomic

281

Transcriptomic profiles between two continuous developmental stages were compared

282

and stage-specific PCR arrays were developed. More than 48,890 genes were detected

283

at day 0, day 6, day 32 and day 62 and the numbers of differently expressed genes

284

between day 0 and day 6, day 6 and day 32, and day 32 and day 62 were 7479, 2971,

285

and 886, respectively. In order to validate the results of transcriptomic sequencing,

286

seventeen genes were randomly selected for quantitative real-time PCR validation,

287

and the results demonstrated that expressions of genes selected were comparable in

288

the two parts: transcriptomic sequencing and quantitative real-time PCR (Supporting

289

Information, Figure S1).

290

Development of stage-specific PCR arrays. Compared with the birth stage (at day

291

0), twenty-one pathways were significantly up-regulated in the growth stage (at day 6)

292

(Figure 2; Table S1, see Supporting Information) and the top 5 up-regulated pathways

293

were all related to DNA replication and DNA repair. Previously, it has been reported

profiles

and

quantitative

real-time

13

ACS Paragon Plus Environment

PCR

validation.

Environmental Science & Technology

294

that during the growth process of multi-cellular organisms, cell proliferation is the

295

most essential biological behavior, and the core part of cell proliferation is DNA

296

replication.42,43 Therefore, our results were consistent with the findings in these

297

studies. Furthermore, compared with growth stage, five pathways were up-regulated

298

in reproduction stage (Figure 2; Table S1, see Supporting Information), and three

299

pathways were related to DNA replication and RNA transcription. In this study,

300

parthenogenesis was the only reproductive mode for D. magna and where oocytes are

301

produced by stem cells in the ovary and undergo only the equational meiotic division

302

in accordance with previous study.44 During the process of parthenogenesis, DNA

303

replication and RNA transcription play important roles. Therefore, the up-regulation

304

of the three pathways was considered to be responsible for the reproduction of D.

305

magna. There were 6 up-regulated pathways in the death stage (at day 62) compared

306

with the reproduction stage (at day 32) (Figure 2; Table S1, see Supporting

307

Information). In the 6 up-regulated pathways, three were related to cardiopathy, and

308

thus we hypothesized that D. magna died from functional loss of the heart. For the

309

stage-specific PCR arrays, 15 pathways and 73 genes were included according to the

310

protocol of development of stage-specific PCR arrays above. The detailed information,

311

including names of pathways and genes involved in each pathway of the arrays, were

312

shown in Table 1. Previously, PCR arrays have been documented as a useful tool to

313

evaluate the effects of chemicals on expressions of genes involved in some certain

314

pathways and explore possible molecular mechanisms.45-47 Therefore, by the use of

315

the method of comparative transcriptome, stage-specific PCR arrays were developed 14

ACS Paragon Plus Environment

Page 14 of 34

Page 15 of 34

Environmental Science & Technology

316

for evaluation of the effects of TDCIPP on the expressions of genes which were

317

responsible for growth, reproduction and death in this study.

318

TDCIPP Exposure Experiments

319

Concentrations of TDCIPP in exposure solutions. Before water renewal, the actual

320

measured TDCIPP concentrations of the nominal 300 and 3000 ng/L exposure groups

321

were 191 ± 9.1 and 3138 ± 463 ng/L, respectively. After water renewal, the measured

322

TDCIPP concentrations were 238 ± 5.7 and 3930 ± 205 ng/L for the nominal 300 and

323

3000 ng/L exposure groups, respectively. Mean concentrations of TDCIPP for

324

samples taken before and after water renewing were 215 ± 7.2 and 3534 ± 298 ng

325

TDCIPP/L. No TDCIPP was detected in the control group samples. Therefore, these

326

data suggested that there were no significant depletion of TDCIPP in this D. magna

327

exposure system, and the results were consistent with a previous study.9

328

Exposure to TDCIPP inhibited growth and reproduction and decreased

329

survival rate. A whole-life-stage exposure protocol was developed using TDCIPP as

330

a case study chemical, and the effects on growth, reproduction and survival were

331

simultaneously evaluated within one life span of D. magna. In this study, exposure to

332

environmentally relevant or greater concentrations of TDCIPP significantly decreased

333

body length of D. magna at day 32 and day 62, while no significant effect was

334

observed at day 6 (Figure 3a). At day 32, the body length of control group was

335

3262.81 ± 57.03 µm, and after TDCIPP exposure the body length was significantly

336

decreased to 3081.29 ± 25.89 µm in 3534 ng/L group (Figure 3a). No significant

337

effect on body length was observed after exposure to 215 ng/L TDCIPP when 15

ACS Paragon Plus Environment

Environmental Science & Technology

338

compared with the control. When the exposure time was prolonged to day 62, the

339

body length was significantly decreased in both the 215 and 3534 ng/L groups (Figure

340

3a). Previous studies demonstrated that treatment with environmentally relevant or

341

greater concentrations of TDCIPP caused developmental toxicity in various

342

organisms.21-27 For example, a recent study reported that environmentally relevant or

343

greater concentrations of TDCIPP significant inhibited body length of D. magna after

344

exposure for 21 and 28 days.9 Therefore, our results were consistent with this study.

345

Additionally, the rapid growth stage of D. magna was from day 0 to day 12, but the

346

decrease in body length was only observed at day 32 and day 62 after TDCIPP

347

exposure. No significant effect on body length was observed at day 6. The possible

348

reason was that the decrease in body length was actually caused from day 6 to day 12,

349

which resulted in the observed decrease at day 32.

350

Besides the inhibition of growth, reproduction of D. magna was also inhibited by

351

TDCIPP (Figure 3b). At day 90, accumulated number of offspring in the control group

352

was 316.62 ± 6.70, but exposure to 215 or 3534 ng/L significantly decreased the

353

accumulated number of offspring to 274.22 ± 6.85 and 263.45 ± 6.19, respectively

354

(Figure 3b). The effects of chemicals on reproduction are usually examined in the

355

evaluation of the environmental risk of chemicals.48,49 In the present study, exposure

356

to environmentally relevant or greater concentrations of TDCIPP for 21 days did not

357

change the number of accumulated offspring of D. magna, and this result was

358

comparable with a recent study.9 However, when the exposure time was prolonged to

359

90 days, a significant inhibitory effect on the production of accumulated offspring was 16

ACS Paragon Plus Environment

Page 16 of 34

Page 17 of 34

Environmental Science & Technology

360

observed. The results demonstrated that a 21-day reproductive toxicity test might

361

underestimate environmental hazards of chemicals, especially for those organisms

362

with reproductive capacity throughout their lives. For survival, exposure to 215 or

363

3534 ng/L TDCIPP did not change the survival rates at day 6 and day 32, but when

364

the exposure time was prolonged to 62 days, the survival rate was significantly

365

decreased by 14.76 % compared with the control (Figure 3c). It should be noted that

366

the decrease in survival rate might be resulted from both natural death and TDCIPP

367

exposure. Previously, most relatively low concentrations in chemical exposure

368

experiments were used in a partial life stage, and thus lethal effects were infrequently

369

observed. However, in the present study, we found that the whole-life-stage exposure

370

of TDCIPP could shorten the life span of D. magna. Therefore, our results indicated a

371

potential risk of environmental chemical pollutants for survival of wildlife, especially

372

for those organisms that were continuously exposed to chemical pollutants during the

373

whole life stage.

374

Transcriptional responses to TDCIPP. Based on the developed PCR arrays,

375

transcriptional responses to TDCIPP were examined at day 6, day 32 and day 62. At

376

day 6, twenty-four of 25 genes involved in the 5 up-regulated pathways selected at the

377

growth stage were significantly down-regulated upon TDCIPP exposure, and among

378

them expressions of the 15 genes were down-regulated in both the 215 and 3534 ng/L

379

groups (Figure 4). The five enriched pathways, including homologous recombination,

380

DNA replication, mismatch repair, nucleotide excision repair and cytosolic

381

DNA-sensing pathway, were all related to DNA replication and DNA repair, which 17

ACS Paragon Plus Environment

Environmental Science & Technology

382

was a core part for growth of multi-cellular organisms.44 Therefore, expressions of

383

genes involved in the five pathways should be up-regulated in fast growth stage

384

compared with the birth stage. The hypothesis was consistent with our findings in

385

comparative transcriptomic study in this study, where expressions of all the genes

386

involved in the five pathways were up-regulated. Therefore the down-regulation due

387

to TDCIPP exposure was considered to be responsible for the inhibition of growth

388

observed in this study. Here, it should be noted that no significant effect on body

389

length was observed after TDCIPP exposure at day 6, although the expressions of

390

those genes screened were down-regulated. The possible explanation was that there

391

might be a delayed effect on body length. At day 32, exposure to 215 or 3534 ng/L

392

TDCIPP significantly altered the expressions of 16 genes involved in the top 5

393

up-regulated pathways of reproduction stage, and expressions of the other 9 genes

394

involved in these pathways were not changed. Among the 16 genes, expressions of 13

395

genes were down-regulated and expressions of 3 genes were up-regulated compared

396

with the control (Figure 5). For the parthenogenesis of D. magna, oocytes are

397

produced by stem cells in the ovary and undergo only the equational meiotic

398

division.50 Therefore, changes in the expressions of genes involved in cell cycle,

399

splicesome and transcriptional misregulation in cancer pathways after TDCIPP

400

exposure were considered as a possible reason for inhibition of reproduction observed

401

in this study. At day 62, expressions of all the genes involved in dilated

402

cardiomyopathy and hypertrophic cardiomyopathy pathways were significantly

403

up-regulated after TDCIPP exposure compared with the control ((Figure 6). In the 18

ACS Paragon Plus Environment

Page 18 of 34

Page 19 of 34

Environmental Science & Technology

404

comparative transcriptomic study, we found that the dilated cardiomyopathy and

405

hypertrophic cardiomyopathy pathways were up-regulated at day 62 compared with

406

day 32, and it was calculated that D. magna might die from functional loss of the

407

heart. Therefore, the up-regulation of genes involved in the two pathways due to

408

TDCIPP might promote functional loss of heart and accelerate death of D. magna. It

409

needs to be emphasized that the exact relationships between the up-regulation of

410

dilated cardiomyopathy and hypertrophic cardiomyopathy pathways, functional loss

411

of heart and decrease of survival rate were not firmly established in this study, and

412

further studies were needed to explore these relationships. Additionally, it should be

413

noted that the changes in the expressions of genes might have resulted from both

414

natural death and TDCIPP exposure. The expressions of genes involved in

415

GABAergic synapse, serotonergic synapse and neuroactive ligand-receptor interaction

416

pathways were also significantly up-regulated. However, the linkage between the

417

up-regulation of these genes and decreased survival rates needs to be explored further.

418

SUPPORTING INFORMATION

419

The Supporting Information is available free of change on the ACS Publications

420

websites. The significantly up-regulated pathways at the stage of growth, reproduction

421

and death for D. magna (Table S1); Sequences of primers for the genes tested in D.

422

magna (Table S2); Expressions of genes in D. magna randomly selected in

423

transcriptomic sequencing and real-time quantitative PCR (Figure S1).

424

ACKNOWLEDGEMENT

425

This work was supported by the Fundamental Research Funds for the Central 19

ACS Paragon Plus Environment

Environmental Science & Technology

426

Universities (2662015PY036) and National Natural Science Foundation of China

427

(21622702). Dr. G. Su was supported by Natural Science Foundation of Jiangsu

428

Province (Grant No. 400 BK20170830), and “the Fundamental Research Funds for

429

the Central Universities” (Grant No. 30917011305).

430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 20

ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34

Environmental Science & Technology

448

REFERENCES:

449

1 Gwinn, M. R.; Axelrad, D. A.; Bahadori, T.; Bussard, D.; Cascio, W. E.; Deener, K.; Dix, D.;

450

Thomas, R. S.; Kavlock, R. J.; Burke, T. A. Chemical Risk Assessment: Traditional vs Public

451

Health

452

10.2105/AJPH.2017.303771.

Perspectives.

AJPH.

Risk.

Assessment.

2017,

107

(7),

1032-1039;

DOI

453

2 Birnbaum, L. S.; Burke, T. A.; Jones, J. J. Informing 21st-century risk assessments with

454

21st-century science. Environ. Health. Perspect. 2016, 124 (4), A60-A63; DOI

455

10.1289/ehp.1511135.

456

3 National Research Council. Toxicity Testing in the 21st Century: A Vision and a Strategy; National Academies Press: Washington, DC, 2007.

457 458

4 Blewett, T. A.; Delompre, P. L. M.; He, Y.; Folkerts, E. J.; Flynn, S. L.; Alessi, D. S.; Goss, G.

459

G. Sublethal and Reproductive Effects of Acute and Chronic Exposure to Flowback and

460

Produced Water from Hydraulic Fracturing on the Water Flea Daphnia magna. Environ. Sci.

461

Technol. 2017, 51, 3032-3039; DOI 10.1021/acs.est.6b05179.

462

5 Manach, S. L.; Khenfech, N.; Huet, H.; Qiao, Q.; Duval, C.; Marie, A.; Bolbach, G; Clodic, G.;

463

Djediat, C.; Bernard, C.; Edery, M.; Marie, B. Gender-Specific Toxicological Effects of Chronic

464

Exposure to Pure Microcystin-LR or Complex Microcystis aeruginosa Extracts on Adult

465

Medaka Fish. Environ. Sci. Technol. 2016, 50, 8324-8334; DOI 10.1021/acs.est.6b01903.

466

6 Wang, G. W.; Du, Z. K.; Chen, H. Y.; Su, Y.; Gao, S. X.; Mao, L. Tissue-Specific Accumulation,

467

Depuration, and Transformation of Triphenyl Phosphate (TPHP) in Adult Zebrafish (Danio

468

rerio), Environ. Sci. Technol. 2016, 50, 13555-13564; DOI 10.1021/acs.est.6b04697.

469

7 Stollewerk, A. “The water flea Daphnia—a “new” model system for ecology and evolution?” J. Biol. 2010, 9, 21; DOI 10.1186/jbiol212.

470 471

8 Otte, K. A.; Fröhlich, T.; Arnold, G. J.; Laforsch, C. Proteomic analysis of Daphnia magna hints

472

at molecular pathways involved in defensive plastic responses. BMC. Genomics. 2014, 15, 306;

473

DOI 10.1186/1471-2164-15-306.

474

9 Li, H.; Su, G. Y.; Zou, M.; Yu, L. Q.; Letcher, R. J.; Yu, H. X.; Giesy, J. P.; Zhou, B. S.; Liu, C.

475

S. Effects of Tris(1,3-dichloro-2-propyl) Phosphate on Growth, Reproduction, and Gene

476

Transcription of Daphnia magna at Environmentally Relevant Concentrations. Environ. Sci.

477

Technol. 2015, 49 (21), 12975-12983; DOI 10.1021/acs.est.5b03294.

478

10 Betts, K. S. Exposure to TDCPP Appears Widespread. Environ. Health. Perspect. 2013, 121 (5), a150; DOI 10.1289/ehp.121-a150.

479 480

11

Hu, M.; Li, J.; Zhang, B.; Cui, G.; Wei, S.; Yu, H. Regional distribution of halogenated

481

organsphosphate flame retardants in seawater samples from three coastal cities in China. Mar.

482

Pollut. Bull. 2014, 86 (1-2), 569-574; DOI 10.1016/j.marpolbul.2014.06.009.

483

12 Sundkvist, A. M.; Olofsson, U.; Haglund, P. Organophosphorus flame retardants and

484

plasticizers in marine and fresh water biota and in human milk. J. Environ. Monit. 2010, 12 21

ACS Paragon Plus Environment

Environmental Science & Technology

485

(4), 943-951; DOI 10.1039/B921910B.

486

13 Carignan, C.C.; McClean, M.D.; Cooper, E.M.; Watkins, D.J.; Fraser, A.J.; Heiger-Bernays,

487

W.; Stapleton, H. M.; Webster, T. F. Predictors of tris (1,3-dichloro-2-propyl) phosphate

488

metabolite in the urine of office workers. Environ. Int. 2013, 55, 56-61; DOI

489

10.1016/jenvint.2013.02.004.

490

14 Meeker, J. D.; Cooper, E. M.; Stapleton, H. M.; Hauser, R. Urinary metabolites of

491

organophosphate flame retardants: temporal variability and correlations with house dust

492

concentrations. Environ. Health Persp. 2013, 121 (5), 585-592; DOI 10.1289/ehp.1205907.

493

15 Ike, van der Veen.; Jacob, de Boer. Phosphorus flame retardants: Properties, production,

494

environmental occurrence, toxicity and analysis. Chemosphere. 2012, 88 (10), 1119-1153;

495

DOI 10.1016/j.chemosphere.2012.03.067.

496

16 Hu, M. Y.; Li, J.; Zhang, B. B.; Cui, Q. L.; Si, W.; Yu, H. X. Regional distribution of

497

halogenated organsphosphate flame retardants in seawater samples from three coastal cities

498

in China. Mar. Pollut. Bull. 2014, 86 (1-2), 569-574; DOI 10.1016/j.marpolbul.2014.06.009.

499

17 Dishaw, L. V.; Powers, C. M.; Ryde, L. T.; Roberts, S. C.; Seidler, F. J.; Slotkin, T. A.;

500

Stapleton, H. M. Is the pentaBDE replacement, tris(1,3-dichloro-2-propyl) phosphate

501

(TDCIPP), a developmental neurotoxicant? Studies in PC12 cells. Toxicol. Appl. Pharmacol.

502

2011, 256, 281-289; DOI 10.1016/j.tapp.2011.01.005.

503

18 Dishaw, L. V.; Hunter, D. L.; Padnos, B.; Padilla, S.; Stapleton, H. M. Developmental

504

exposure to organophosphate flame retardants elicits overt toxicity and alters behavior in

505

early life stage zebrafish (Danio rerio). Toxicol. Sci. 2014, 142 (2), 445-454; DOI

506

10.1093/toxsci/kfu194.

507

19 Ta, N.; Li, C. N.; Fang, Y. J.; Liu, H. L.; Lin, B. C.; Jin, H.; Tian, L.; Zhang, H.; Zhang, W.;

508

Xi, Z. Toxicity of TDCPP and TCEP on PC12 cell: changes in AMKII, GAP43, tubulin and

509

NF-H gene and protein expression. Toxicol. Lett. 2014, 227 (3), 164-171; DOI

510

10.1016/j.toxlet.2014.03.023.

511

20 Wang, Q. W.; Lam, J. C. W.; Man, Y. C.; Lai, N. L. S.; Kwok, K. Y.; Guo, Y. Y.; Lam, P. W. S.;

512

Zhou, B. S. Bioconcentration, metabolism and neurotoxicity of the organophorous flame

513

retardant 1,3-dichloro-2-propyl phosphate (TDCIPP) to zebrafish. Aquat. Toxicol. 2015, 158,

514

108-115; DOI 10.1016/j.aquatox.2014.11.001.

515

21 Wang, Q. W.; Lam, J. C. W.; Han, J.; Wang, X. F.; Guo, Y. Y.; Lam, P. K. S.; Zhou, B.

516

Developmental exposure to the organophosphorus flame retardant tris(1,3-dichloro-2-propyl)

517

phosphate: estrogenic activity, endocrine disruption and reproductive effects on zebrafish. 22

ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34

Environmental Science & Technology

518

Aquat. Toxicol. 2015, 160, 163-171; DOI 10.1016/j.aquatox.2015.01.014.

519

22 Farhat, A.; Crump, D.; Chiu, S.; Williams, K. L.; Letcher, R. J.; Gauthier, L. T.; Kennedy, S.

520

W. In ovo effects of two organophosphate flame retardants-TCPP and TDCPP-on pipping

521

success, development, mRNA expression, and thyroid hormone levels in chicken embryos.

522

Toxicol. Sci. 2013, 134 (1), 92-104; DOI 10.1093/toxsci/kft100.

523

23 McGee, S. P.; Cooper, E.; Stapleton, H. M.; Volz, D. C. Early zebrafish embryogenesis is

524

susceptible to developmental TDCPP exposure. Environ. Health Perspect. 2012, 120 (11),

525

1585-1591; DOI 10.1289/ehp.1205316.

526

24 Fu, J.; Han, J.; Zhou, B. S.; Gong, Z. Y.; Santos, E. M.; Huo, X. Zheng, W.; Liu, H.; Yu, H.;

527

Liu, C. Toxicogenomic responses of zebrafish embryos/larvae to tris (1,3-dichloro-2-propyl)

528

phosphate (TDCPP) reveal possible molecular mechanisms of developmental toxicity.

529

Environ. Sci. Technol. 2013, 47 (18), 10574-10582; DOI 10.1021/es401265q.

530

25 Wang, Q. W.; Liang, K.; Liu, J. F.; Yang, L. H.; Guo, Y. Y.; Liu, C.; Zhou, B. Exposure of

531

zebrafish embryos/larvae to TDCPP alters concentrations of thyroid hormones and

532

transcriptions of genes involved in the hypothalamic-pituitary-thyroid axis. Aquat. Toxicol.

533

2013, 126, 207-213; DOI 10.1016/j.aquatox.2012.11.009.

534

26 Liu, X.; Ji, K.; Choi, K. Endocrine disruption potentials of organophosphate flame retardants

535

and related mechanisms in H295R and MVLN cell lines and in zebrafish. Aquat. Toxicol.

536

2012, 114-115, 173-181; DOI 10.1016/j.aquatox.2012.02.019.

537

27 Liu,X.S.;Ji,K.;Jo,A.;Moon,H.B.;Choi,K.Effects of TDCPP or TPP on gene transcriptions and

538

hormones of HPG axis, and their consequences on reproduction in adult zebrafish (Danio

539

rerio). Aquat. Toxicol. 2013, 134-135, 104-111; DOI 10.1016/j.aquatox.2013.03.013.

540

28 Campos, B.; Piña, B.; Barata, C. Mechanisms of action of selective serotonin reuptake

541

inhibitors in Daphnia magna. Environ. Sci. Technol. 2012, 46, 2943-2950; DOI

542

org/10.1021/es203157f.

543

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

544

X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; Chen, Z.; Mauceli, E.; Hacohen, N.; Gnirke, A.;

545

Rhind, N.; di Palma, F.; Birren, B. W.; Nusbaum, C.; Lindblad-Toh, K.; Friedman, N.; Regev,

546

A. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat.

547

Biotechnol. 2011, 29, 644-652; DOI 10.1038/nbt.1883. 23

ACS Paragon Plus Environment

Environmental Science & Technology

Page 24 of 34

548

30 Haas, B. J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P. D.; Bowden, J.; Couger,

549

M. B.; Eccles, D.; Li, B.; Lieber, M.; Macmanes, M. D.; Ott, M.; Orvis, J.; Pochet, N.;

550

Strozzi, F.; Weeks, N.; Westerman, R.; William, T.; Dewey, C. N.; Henschel, R.; Leduc, R. D.;

551

Friedman, N.; Regev, A. De novo transcript sequence reconstruction from RNA-seq using the

552

Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8 (8), 1494-1512;

553

DOI 10.1038/nprot.2013.084.

554

31 Mathavan, S.; Lee, S. G. P.; Mak, A.; Miller, L. D.; Murthy, K. R. K.; Govindarajan, K. R.;

555

Tong, Y.; Wu, Y. L.; Lam, S. H.; Yang, H.; Ruan, Y. J.; Korzh, V.; Gong, Z. Y.; Liu, E. T.;

556

Lufkin, T. Transcriptome Analysis of Zebrafish Embryogenesis Using Microarrays. PLoS.

557

Genet. 2005, 1 (2), e29; DOI:10.1371/journal.pgen.0010029.

558

32 Viñas, J.; Piferrer, F. Stage-Specific Gene Expression During Fish Spermatogenesis as

559

Determined by Laser-Capture Microdissection and Quantitative-PCR in Sea Bass

560

(Dicentrarchus

561

DOI:10.1095/biolreprod.108.069708.

labrax)

Gonads.

Biol.

Reprod.

2008,

79

(4),

738-747;

562

33 Li, J.; Ma, X. F.; Su, J. Y.; Giesy, J. P.; Xiao, Y.; Zhou, B. S,; Letcher, R, J.; Liu, C. S.

563

Multigenerational effects of tris(1,3-dichloro-2-propyl) phosphate on the free-living ciliate

564

protozoa Tetrahymena thermophila exposed to environmentally relevant concentrations and

565

after

566

10.1016/j.envpol.2016.08.034.

subsequent

recovery.

Environ.

Pollut.

2016,

218,

50e58;

DOI

567

34 Khan, F. R.; Pual, K. B.; Dybowska, A. D.; Valsami-Jones, E.; Lead, J. R.; Stone, V.;

568

Fernandes, T. F. Accumulation Dynamics and Acute Toxicity of Silver Nanoparticles to

569

Daphnia magna and Lumbriculus variegatus: Implications for Metal Modeling Approaches.

570

Environ. Sci. Technol. 2015, 49 (7),4389-4397; DOI 10.1021/es506124x.

571

35 Li, L. Z.; Wu, H. F.; Ji, C. L.; Van Gestel, C. A. M.; Allen, H. E.; Peijnenburg, W. J. G. M.

572

A metabolomic study on the responses of Daphnia magna exposed to silver nitrate and

573

coated

574

10.1016/j.ecoenv.2015.05.005.

silver

nanoparticles.

Ecotox.

Environ.

Safe.

2015,

119,

66-73;

DOI

575

36 Simpson, A. M.; Jeyasingh, P. D.; Belden, J. B. Variation in toxicity of a current-use

576

insecticide among resurrected Daphnia pulicaria genotypes. Ecotoxicology. 2015, 24,

577

488-496; DOI 10.1007/s10646-014-1397-1. 24

ACS Paragon Plus Environment

Page 25 of 34

Environmental Science & Technology

578

37 Zhang, L. L.; Liu, J. L.; Liu, H. Y.; Wan, G. S.; Zhang, S. W. The occurrence and ecological

579

risk assessment of phthalate esters (PAEs) in urban aquatic environments of China.

580

Ecotoxicology. 2015, 24, 967-984; DOI 10.1007/s10646-015-1446-4.

581

38 Smith, A. S.; Acharya, K.; Jack, J. Overcrowding, food and phosphorus limitation effects on

582

ephippia production and population dynamics in the invasive species Daphnia lumholtzi.

583

Hydrobiologia. 2009, 618, 47-56; DOI 10.1007/s10750-008-9546-2.

584

39 Acharya, K.; Kyle, M.; Elser, J. J. Biological stoichiometry of Daphnia growth: an

585

ecophysiological test of the growth rate hypothesis. Limnol. Oceanogr. 2004, 49, 656-665;

586

DOI 10.4319/lo.2004.49.3.0656.

587

40 Holderbaum, D. F.; Cuhra, M.; Wickson, F.; Orth, A. I.; Nodari, R. O.; Bøhn, T. Chronic

588

Responses of Daphnia magna Under Dietary Exposure to Leaves of a Transgenic (Event

589

Mon810) Bt–maized Hybrid and Its Conventional Near-isoline. J. Toxicol. Env Health. 2015,

590

78, 993-1007; DOI 10.1080/15287394.2015.1037877.

591

41 Bitto, A.; Ito, T. K.; Pineda, V. V.; LeTexier, N. J.; Huang, H. Z.; Sutlief, E.; Tung, H.; Vizzini,

592

N.; Chen, B.; Smith, K.; Meza, D.; Yajima, M.; Beyer, R. P.; Kerr, K. F.; Davis, D. J.;

593

Gillespie, C. H.; Snyder, J. M.; Treuting, P. M.; Kaeberlein, M. Transient rapamycin

594

treatment can increase lifespan and healthspan in middle-aged mice. ELife. 2016, 5, e16351;

595

DOI 10.7554/eLife.16351.

596

42 Murray, H.; Koh, A. Multiple Regulatory Systems Coordinate DNA Replication with Cell

597

Growth

in

Bacillus

subtilis.

598

10.1371/journal.pgen.1004731.

PLOS.

Genet.

2014,

10

(10),

e1004731;

DOI

599

43 Piunti, A.; Rossi1, A.; Cerutti, A.; Albert, M.; Jammula, S,; Scelfo, A.; Cedrone, L.; Fragola,

600

G.; Olsson, L.; Koseki, H.; Testa1, G.; Casola, D.; Helin, K.; Fagagna, F. A.; Pasini, D.

601

Polycomb proteins control proliferation and transformation independently of cell cycle

602

checkpoints by regulating DNA replication. Nat. Commun. 2014, 5:3649; DOI

603

10.1038/ncomms4649.

604

44 Kami D. M. Harris, Nicholas J. Bartlett, and Vett K. Lloyd. Daphnia as an Emerging

605

Epigenetic

606

10.1155/2012/147892.

607

Model

Organism.

Genet.

Res.

Int.

2012,

2012:147892;

DOI

45 Pripuzova1, N.; Wang, R.; Tsai, S.; Li, B. J.; Hung, G.; Ptak, R, G.; Lo, S. Development of 25

ACS Paragon Plus Environment

Environmental Science & Technology

608

Real-Time PCR Array for Simultaneous Detection of Eight Human Blood-Borne Viral

609

Pathogens. PLoS. ONE. 7 (8), e43246; DOI 10.1371/journal.pone.0043246.

610

46 Boone, D. R.; Micci, M.; Taglialatela, I. G.; Hellmich, J. L.; Weisz, H. A.; Bi, M.; Prough, D.

611

S.; DeWitt, D. S.; Hellmich, H. L. Pathway-Focused PCR Array Profiling of Enriched

612

Populations of Laser Capture Microdissected Hippocampal Cells after Traumatic Brain Injury.

613

PLoS. ONE. 10 (5), e0127287; DOI 10.1371/journal.pone.0127287.

614

47 Veselenak, R. L.; Miller, A. L.; Milligan, G. N.; Bourne, N.; Pyles, R. B. Development and

615

Utilization of a Custom PCR Array Workflow: Analysis of Gene Expression in Mycoplasma

616

genitalium and Guinea Pig (Cavia porcellus). Mol. Biotechnol. 2015, 57, 172-183; DOI

617

10.1007/s12033-014-9813-6.

618

48 King-Heiden, T. C.; Mehta, V.; Xiong, K. M.; Lanham, K. A.; Antkiewicz, D. S.; Ganser, A.;

619

Heideman, W.; Peterson, R. E. Reproductive and Developmental Toxicity of Dioxin in Fish.

620

Mol. Cell. Endocrinol. 2012, 354 (1-2), 121-138; DOI 10.1016/j.mce.2011.09.027.

621

49 Duonga, A.; Steinmausa, C.; McHalea, C. M.; Vaughanc, C. P.; Zhang, L. P. Reproductive

622

and Developmental Toxicity of Formaldehyde: A Systematic Review. Mutat. Res. 2011, 728

623

(3), 118-138; DOI 10.1016/j.mrrev.2011.07.003.

624

50 Riera, A.; Barbon, M.; Noguchi, Y.; Reuter, L. M.; Schneider, S.; Speck, C. From structure to

625

mechanism—understanding initiation of DNA replication. Gene. Dev. 2017, 31, 1073-1088;

626

DOI 10.1101/gad.298232.

627 628 629 630 631 632 633 634 635 636 637 26

ACS Paragon Plus Environment

Page 26 of 34

Page 27 of 34

Environmental Science & Technology

638 639 640

641 642

Figure 1. The whole-life-stage growth (a), reproduction (b) and survival (c) curves of D. magna.

643

Values are the mean ± SD (n=3).

644 645 646 647 648 649 650 651

27

ACS Paragon Plus Environment

Environmental Science & Technology

652 653

Figure 2. Numbers of significantly up-regulated gene pathways between different developmental

654

stages in D. magna.

655 656 657 658 659

28

ACS Paragon Plus Environment

Page 28 of 34

Page 29 of 34

Environmental Science & Technology

660 661 662

Figure 3. The body length (a), accumulated offspring (b) and survival rate (c) of D. magna after

663

exposure to 215 or 3534 ng/L TDCIPP for 6, 32 and 62 days. Values are the mean ± SD (n=3).

664

One-way analysis of variance and tukey’s multiple range test were used to determine

665

significant differences of parameters tested between control and exposure groups. * p

666

< 0.05.

667 668 669 670 671 672

29

ACS Paragon Plus Environment

Environmental Science & Technology

673

674 675

Figure 4. Expression of genes involved in the top 5 up-regulated pathways in the growth stage of

676

D. magna were significantly down-regulated after exposure to 215 or 3534 ng/L TDCIPP for 6

677

days. Values are the mean ± SD (n=3). One-way analysis of variance and tukey’s multiple

678

range test were used to determine significant differences of parameters tested between

679

control and exposure groups. * p < 0.05.

680 681 682 683 684 685 686

30

ACS Paragon Plus Environment

Page 30 of 34

Page 31 of 34

Environmental Science & Technology

687 688 689

Figure 5. Expressions of genes involved in the top 5 up-regulated pathways of D. Magna in the

690

reproduction stage were significantly changed after exposure to 215 or 3534 ng/L TDCIPP for 32

691

days. Values are the mean ± SD (n=3). One-way analysis of variance and tukey’s multiple

692

range test were used to determine significant differences of parameters tested between

693

control and exposure groups. * p < 0.05.

694 695 696 697 698 699 700

31

ACS Paragon Plus Environment

Environmental Science & Technology

701 702 703 704

Figure 6. Expressions of genes involved in the top 5 up-regulated pathways in the death stage D.

705

magna were significantly altered after exposure to 215 or 3534 ng/L TDCIPP for 62 days. Values

706

are the mean ± SD (n=3). One-way analysis of variance and tukey’s multiple range test

707

were used to determine significant differences of parameters tested between control.

708

and exposure groups. * p < 0.05.

709 710 711 712 713 714 715 716 32

ACS Paragon Plus Environment

Page 32 of 34

Page 33 of 34

Environmental Science & Technology

717

Table 1. A list of the top 5 genes and the corresponding pathways in the developed

718

stage-specific PCR arrays Life stage

Stage of growth

Stage of reproduction

Stage of death

Up-regulated pathway

Genes

Homologous recombination

rpa1, pold1, rad54l, rad52, mre11

DNA replication

rnha, rnaseh1, rfc3_5, pold1, mcm7

Mismatch repair

rpa1, pms2, podl1, msh6, rfc3_5

Nucleotide excision repair

pms2, pold1, pcna, msh6, rfc3_5

Cytosolic DNA-sensing pathway

rpc1_1, rpc1_2, rpc3, rpc1_3, rpc1_4

Chronic myeloid leukemia

runx1, shc1, abl1, cdk6, pik3c

Cell cycle

mcm2, ccnb, ccne, cdc25c, cdk6

Transcriptional misregulation in cancer

etv6_7, utx, flt1, prom1, yan

Spliceosome

hnrnpk, hspa1_8_1, hspa1_8_2, snrnp200, hspa1_8_3

Axon guidance

ablim, sgrap, plxna, smo, ptch1

GABAergic synapse

slc6a1_1, gnao, slc6a1_2, plcl1, slc32a

Dilated cardiomyopathy

sgcd, tpm3, actb_g1, adcy5_1

Neuroactive ligand-receptor interaction

grik2, sstr5, glra1, takr2, crhr1

Hypertrophic cardiomyopathy

sgcd, tpm3, actb_g1

Serotonergic synapse

ddc, gnao, cacna1b, slc18a1_2, adcy5_1

Internal control:

c9685_g1

719 720 721 722 723 33

ACS Paragon Plus Environment

Environmental Science & Technology

724

TOC art

725

34

ACS Paragon Plus Environment

Page 34 of 34