Arsenic Reduces Gene Expression Response to ... - ACS Publications

Arsenic Impairs Salinity Acclimation in Killifish [Running Head]. 1. Thomas H. Hampton ... 25. Page 2 of 39. ACS Paragon Plus Environment. Environment...
0 downloads 0 Views 2MB Size
Subscriber access provided by Technical University of Munich University Library

Ecotoxicology and Human Environmental Health

Arsenic Reduces Gene Expression Response to Changing Salinity in Killifish Thomas Hampton, Craig Jackson, Dawoon Jung, Celia Y. Chen, Stephen P. Glaholt, Bruce A. Stanton, John K Colbourne, and Joseph R. Shaw Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01550 • Publication Date (Web): 06 Jul 2018 Downloaded from http://pubs.acs.org on July 7, 2018

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

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

Page 1 of 39

Environmental Science & Technology

1

Arsenic Impairs Salinity Acclimation in Killifish [Running Head]

2

Thomas H. Hampton

3

Geisel School of Medicine at Dartmouth

4

Department of Microbiology and Immunology

5

HB 7550, Remsen, Room 517

6

Hanover, NH 03755

7

USA

8

Tel: 603-650-1184

9

Fax: 603-650-1130

10

email: [email protected]

ACS Paragon Plus Environment

Environmental Science & Technology

11

Arsenic Reduces Gene Expression Response to Changing

12

Salinity in Killifish

13

Thomas H. Hampton*1,2, Craig Jackson3, Dawoon Jung,2,4, Celia Y. Chen5, Stephen P.

14

Glaholt3, Bruce A. Stanton2, John K. Colbourne1, Joseph R. Shaw1,3

15 16 17 18 19 20

1. Environmental Genomics Group, School of Biosciences, University of Birmingham, Birmingham, United Kingdom 2. Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH 3. The School of Public and Environmental Affairs, Indiana University, Bloomington

21

4. Korea Environment Institute, Republic of Korea

22

5. Department of Biological Sciences, Dartmouth College, Hanover, NH

23 24 25

ACS Paragon Plus Environment

Page 2 of 39

Page 3 of 39

Environmental Science & Technology

26 27

ACS Paragon Plus Environment

Environmental Science & Technology

28 29 30

Abstract

31

Toxicogenomic approaches can detect and classify adverse interactions between

32

environmental toxicants and other environmental stressors but require more complex

33

experimental designs and analytical approaches. Here we use novel toxicogenomic

34

techniques to analyze the effect of arsenic exposure in wild killifish populations

35

acclimating to changing salinity. Fish from three populations were acclimated to full

36

strength seawater and transferred to fresh water for 1 h or 24 h. Linear models of gene

37

expression in gill tissue identified 31 genes that responded to osmotic shock at 1 h and

38

178 genes that responded at 24 h. Arsenic exposure (100 µg/l) diminished the responses

39

(reaction norms) of these genes by 22% at 1h (p = 1.0 e-6) and by 10% at 24 h (p =

40

3.0 e-10). Arsenic also significantly reduced gene co-regulation in gene regulatory

41

networks (p = 0.002, paired Levene’s test), and interactions between arsenic and salinity

42

acclimation were uniformly antagonistic at the biological pathway level (p < 0.05,

43

binomial test). Arsenic’s systematic interference with gene expression reaction norms

44

was validated in a mouse multi-stressor experiment, demonstrating the ability of these

45

toxicogenomic approaches to identify biologically relevant adverse interactions between

46

environmental toxicants and other environmental stressors.

47 48

Table of Contents (TOC)/Abstract Art.

ACS Paragon Plus Environment

Page 4 of 39

Page 5 of 39

Environmental Science & Technology

49 Osmotic Shock

Gene Expression Response 3.5

***

***

Delta Effect Estimate (log2)

2.5

1.5

0.5

−0.5

−1.5

NaAs02

Pathway Activation Regulatory Networks

*** ** *** * * * *** ** *** *** * * *** ** * ***

***

**

* * **

** *** ** **

−2.5

50

** * *** * *

*

*

−3.5 1h

24h

* *

* *

51 52

Introduction

53

Arsenic is naturally present in bedrock and is released into aquatic systems, particularly

54

in reducing conditions at high pH 1. Dissolved inorganic arsenic enters cells through

55

sodium/phosphate cotransporters and water channels 2, where it may take the place of

56

phosphorous in biochemical reactions or interact with thiol groups in proteins and

57

peptides 3 leading to toxic effects. Many organisms have evolved mechanisms to detoxify

58

arsenic, e.g., arsenic methyltransferases, that are highly conserved in nature 4

59

demonstrating that arsenic is a broadly relevant environmental toxicant. Nonetheless,

60

environmental arsenic is best known as a threat to human health, affecting over 100

61

million people 5 by increasing their risk of dozens of diseases 6.

62

Many organisms have been used to understand arsenic’s diverse effects at low doses. For

63

example, studies in chronically exposed mice reveal reduced immune signaling 7, an

64

effect also observed in zebrafish embryos 8. Xenopus tadpoles chronically exposed to

65

arsenic show impaired metamorphosis 9, chick embryos exposed to arsenic respond less

66

to estrogen 10, and mice chronically exposed to arsenic respond less effectively to

ACS Paragon Plus Environment

Environmental Science & Technology

67

influenza infection 11. Collectively, these studies suggest that arsenic exposure diminishes

68

the ability of organisms to respond effectively to stress, and our overarching hypothesis is

69

therefore that arsenic’s ability to diminish stress responses is associated with reduced

70

gene expression responses.

71

Therefore, we tested the hypothesis that exposure to arsenic during a response to a second

72

stressor reduces responses at the gene expression level using killifish as a model system.

73

Like most marine fish, killifish living in salt water maintain internal sodium chloride

74

concentrations that are about 60% less than salt water by excreting excess chloride

75

through their gills 12. Unlike most fish, killifish that are acclimated to salt water can

76

survive the hypo-osmotic shock of being placed in fresh water. They undergo a process of

77

gill tissue remodeling over several days to produce a gill phenotype appropriate to fresh

78

water. Fresh water gills differ from salt water gills primarily in terms of the density and

79

disposition of ionocytes, ion channel activity, epithelial tight junctions and blood supply

80

13

81

studied at various time points out to 10 days, and killifish from different populations are

82

known to vary in their ability to respond to salinity 14-18. We assessed gene expression in

83

three populations of killifish, at two distinctly different time points, to identify conserved

84

responses.

85

Killifish survive in changing salinity because they can alter their gill phenotype to suit

86

environmental conditions, in other words, they exhibit phenotypic plasticity 19 20. Altered

87

gene expression has been reported to facilitate phenotypic plasticity in extreme

88

temperatures 21,22 or variability in water availability 23. The ability of an organism to

89

regulate gene expression can be quantified as gene expression reaction norms 24. Here we

. Changes in killifish gill gene expression during gill remodeling have been extensively

ACS Paragon Plus Environment

Page 6 of 39

Page 7 of 39

Environmental Science & Technology

90

applied the concept of gene expression reaction norms to quantify the impact of arsenic

91

on salinity responses in killifish, hypothesizing that reduced gene expression reaction

92

norms during acclimation to salinity in the killifish gill will reduce phenotypic plasticity.

93

Gene regulatory networks coordinate gene expression and have been shown to facilitate

94

phenotypic plasticity 15,25. We used a simple statistic based upon Pearson correlation 26 to

95

quantify the level of gene regulatory phenotypic plasticity as connectivity in our system

96

and the effect of arsenic on that network connectivity. We hypothesized that exposure to

97

arsenic would reduce network connectivity and interfere with the ability of genes to

98

respond in a coordinated fashion.

99

Specific gene regulatory networks have been annotated as biological pathways that

100

perform specific functions, e.g., those belonging to the Kyoto Encyclopedia of Genes and

101

Genomes (KEGG) 27. We hypothesized that if arsenic interference with gene expression

102

responses leads to reduced phenotypic plasticity, arsenic interference might target

103

pathways with functions that facilitate the early or late phases of tissue remodeling in the

104

killifish gill. To assess this, we used Pathway Activation Analysis (PAA) to score

105

whether genes in a given pathway are systematically turned on or off by arsenic or other

106

treatment effects. The general principle is similar to Differential Expression for Pathways

107

(DEAP) 28 but PAA does not require explicit knowledge of regulatory relationships. In

108

addition, PAA does not rely on arbitrary significance cutoffs, which are a substantial

109

drawback of first generation pathway analysis approaches including using Fisher’s exact

110

test to identify pathways enriched in genes found significant in a statistical test29.

111

Novel methods presented here are relevant to any toxicogenomic study involving

112

multiple stressors, and demonstrate that they can detect significant shifts in gene

ACS Paragon Plus Environment

Environmental Science & Technology

113

expression that increase our understanding of interactions between biotic and abiotic

114

stressors 30 31. Specifically, we found that arsenic reduced gene expression reaction norms

115

in killifish gill tissue responding to osmotic shock, and that this effect generalized to a

116

second animal model (i.e., mouse). We also report that arsenic diminished network

117

connectivity, and limited activation of biological pathways. An analysis of publicly

118

available data showed that arsenic also reduces reaction norms in response to

119

dexamethasone, a synthetic corticosteroid stress hormone. Taken together, our findings

120

suggest that arsenic toxicity may increase as a function of environmental stress.

121

Materials and Methods

122

Killifish, Arsenic Exposures, and Tissue Collection

123

Gill tissue was obtained from male wild killifish sampled from two locations in Maine

124

and one in Virginia. We used fish from Northeast Creek, Maine because we have used

125

these fish extensively in previous experiments 2,15,32-39. Fish from Horseshoe Cove Maine,

126

and King’s Creek Virginia were included in the study to ensure that our results would

127

focus on responses shared by multiple killifish populations. We used a nearby Maine

128

population (Horseshoe Cove) because they were easy to collect and expected to be

129

genetically similar to the Northeast Creek population. The population from King’s Creek,

130

Virginia was supplied by colleagues from the Duke Ecotoxicology Lab, and was included

131

because it is representative of the southern killifish clade, which includes all killifish

132

south of the Hudson River in New York 40. Fish (Northeast Creek, 1.62±0.47g,

133

Horseshoe Cove, 1.68±0.32g, King’s Creek, 4.74±0.62g) were pre-acclimated to

134

common garden conditions, housed, maintained and exposed to arsenic as previously

ACS Paragon Plus Environment

Page 8 of 39

Page 9 of 39

Environmental Science & Technology

135

described 15 in accordance with IACUC # MDIBL 13-01 from Mount Desert Island

136

Biological Laboratory, Salisbury Cove, Maine. Briefly, 72 killifish were maintained for at

137

least two weeks in seawater (pH 8.1 ± 0.4; salinity 33 ± 0.5‰, 15±1 °C) and exposed to

138

natural light cycle (16:8 h light:dark), to assure all were fully acclimated. Experimental

139

factors included population (Northeast Creek, ME, Horseshoe Cove, ME, King’s Cove,

140

VA), exposure to 100 µg/l arsenic (yes/no), and time exposed to fresh water (0 h, 1 h, 24

141

h) to yield a balanced design of 18 groups with 4 fish per group. Freshwater and seawater

142

conditions were established as previously described 41 15. Fish were fed commercial flake

143

food (48% protein, 9% fat; Tetracichlid, Tetra, Blacksburg, VA) once a day that

144

contained no detectable inorganic arsenic, or monomethyl- or dimethyl arsenic (Shaw et

145

al. 2010). Arsenic concentrations in exposed samples were assayed at 100.2 µg/l ± 1.4

146

µg/l (SD), total arsenic, 98.9% inorganic As III using a 7700x ICP-MS (Agilent, Santa

147

Clara, CA) mass spectrometer. Arsenic concentrations in unexposed samples were below

148

the detection limit of 1.0 µg/l arsenic.

149

Fish were anesthetized, pithed, and gills were removed and stored in RNAlater (Qiagen,

150

Valencia, CA) according to manufacturer’s recommendations.

151

RNA Isolation, Hybridization and Normalization

152

Tissue samples were removed from RNAlater and rinsed before being homogenized

153

using a Tissue-Tearor (Biospec Products, Bartlesville, OK). The homogenized tissue was

154

then processed using RNeasy kits with DNase treatment (Qiagen, Valencia, CA)

155

according to manufacturer’s protocols to extract RNA. RNA quality and concentration

156

were determined by an Agilent 2100 bioanalyzer. All samples used in genome expression

ACS Paragon Plus Environment

Environmental Science & Technology

157

studies achieved an RNA integrity number (RIN) score >7, as required by the in-house

158

protocols of the IU Center for Genomics and Informatics microarray processing core, and

159

widely suggested in the literature42. RNA was amplified using MessageAmp II kits

160

(Ambion) and hybridized to a custom NimbleGen array that interrogates 135,000 probes

161

associated with 69,426 unique contigs, that map to 16,104 unique genes 15. Raw

162

fluorescence values were quantile normalized across arrays using RMA43. Complete gene

163

expression data has been deposited at the NCBI Gene Expression Omnibus 44 and are

164

accessible through GEO Series accession number

165

GSE104218 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104218).

166

Identifying Genes that Significantly Responded to Hypo-osmotic Shock

167

As certain killifish genes were associated with more than one probe on the microarray,

168

we selected the probe with the highest median expression across all conditions to

169

represent each unique killifish gene. The Robust MultiArray Average (RMA) normalized

170

log2 expression values for each gene were then analyzed using various linear models in R.

171

False discovery rates (FDR) were calculated from linear model p values using the method

172

of Benjamini and Yekutieli 45, and genes with an FDR less than 0.05 were deemed

173

significant. First, ordinary linear models were used to assess main effects of population,

174

arsenic, acclimation time and interactions between arsenic and acclimation. The

175

expression of about 60% of all genes differed significantly (FDR < 0.05) between

176

populations, and Maine populations were as different from each other as they were from

177

the Virginia population. Neither arsenic as a factor, nor interactions with arsenic

178

identified specific genes that were differentially expressed (FDR < 0.05) in this model,

179

though genes were identified as differentially expressed at 1 h and 24 h. To increase

ACS Paragon Plus Environment

Page 10 of 39

Page 11 of 39

Environmental Science & Technology

180

sensitivity to arsenic and its interactions with changing salinity, we modeled gene

181

expression as a function of arsenic, acclimation time and their interactions, treating

182

population as a random effect using the nlme package 46. Fixed effect linear models and

183

mixed effects linear models on data subsets were also used to explore the data.

ACS Paragon Plus Environment

Environmental Science & Technology

184

Visualizing Differences between of 1 h and 24 h Responses

185

We used Venn diagrams from the R gplots package (https://CRAN.R-

186

project.org/package=gplots) and pairs plots from the R GGally package

187

(https://CRAN.R-project.org/package=GGally) to establish that genes responding at 1 h

188

are distinct from each other and respond differently at the two time points.

189

Identifying Genes that Significantly Responded to Arsenic

190

Procedures that control false discovery rates to correct for multiple hypothesis testing set

191

a higher bar for significance and therefore decrease sensitivity, that is, procedures that

192

decrease the false positive rate increase the false negative rate. In the toxicogenomic

193

setting, the capacity to detect true effects may be more important than the ability to avoid

194

detecting spurious effects, so we used a graphical analysis of p value distributions to

195

assess whether FDR corrected statistics might underestimate differential gene expression

196

responses to arsenic and interactions between arsenic and exposure to changing salinity.

197

We used ggplot2 (http://ggplot2.org) to visualize that the number of observed p values

198

that were less than 0.05 exceeded the fraction that would be expected based on the null

199

hypothesis and gene expression independence.

200

Systematic Assessment of Gene Expression Reaction Norms

201

The number of genes that respond to treatment is a useful proxy for biological response,

202

and therefore relevant to toxicogenomics 47,48. Therefore, we calculated the impact of

203

arsenic on gene expression reaction norms in salinity responsive genes. First, we used

204

mixed effect linear models with population as a random effect to identify genes that

205

respond significantly (FDR < 0.05) to salinity. Second, we calculated the average

ACS Paragon Plus Environment

Page 12 of 39

Page 13 of 39

Environmental Science & Technology

206

absolute log2 response (gene expression reaction norm) in these genes at 1 h and 24 h

207

compared to 0 h, in fish that were not experimentally exposed to arsenic. Analogous

208

reaction norm calculations were made for arsenic exposed fish, and the impact of arsenic

209

on hypo-osmotic shock gene expression reaction norms was calculated by subtracting

210

arsenic unexposed reaction norms from arsenic exposed reaction norms. The impact of

211

arsenic on reaction norms was visualized for genes that significantly responded to hypo-

212

osmotic shock using pirate plots in yarrr 49, and significant differences were identified by

213

a one sample t test where the mean difference attributable to arsenic under the null

214

hypothesis was predicted to be zero.

215

Quantifying Changes in Gene Regulatory Network Structure

216

Genes regulate other genes, and if arsenic interferes with the ability of killifish to respond

217

to changes in salinity, we hypothesized that it may do so by interfering with regulatory

218

relationships among genes. Regulatory relationships among genes may be visualized as

219

networks and inferred from the level of correlation among genes 50. The impact of arsenic

220

on regulatory network structure was quantified using correlation matrices as follows.

221

Separate correlation matrices were constructed for fish exposed to either 100 µg/l or 0

222

µg/l arsenic to explore regulatory relationships between genes that respond to hypo-

223

osmotic shock. Pairwise Pearson correlations between all genes significantly responding

224

to salinity at a given time point (1 h, 24 h) were calculated for all samples using the

225

corr.test package in the R psych library 51, resulting in four matrices. Correlation values

226

from each time point (plus and minus arsenic) were then compared in terms of their

227

variability using levene.var.test in the R package PairedData 52 which implements a

228

paired Levene’s test using Wilcox’s approach 53. Smaller variability in correlation

ACS Paragon Plus Environment

Environmental Science & Technology

229

indicates a larger proportion of genes with near zero correlation, and hence less

230

connectivity in the gene regulatory network.

231

Pathway Activation Analysis

232

While it is useful to quantify the impact of arsenic on gene reaction norms and regulatory

233

networks, these assessments do not identify specific biological functions affected by

234

arsenic. To address this shortcoming, we mapped our killifish genes to orthologous genes

235

in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways based on zebrafish

236

orthologs using OrthoMCL54 and reciprocal blastp as described in the publication of the

237

killifish reference genome 55. Effect estimates from linear models for all genes (regardless

238

of statistical significance) in 166 KEGG paths were tested for significant pathway level

239

biases toward induction or repression using binomial tests with a null hypothesis that

240

50% of the genes in any pathway are induced. Significance levels of binomial tests were

241

FDR-corrected for multiple hypothesis testing.

242

Validation of Arsenic Impact on Gene Expression Reaction Norms

243

Finally, we used publicly available data to validate that the patterns of arsenic

244

interference observed in killifish during hypo-osmotic shock generalize to other systems

245

exposed to multiple stressors. Validation used normalized gene expression data from two

246

studies in the Gene Expression Omnibus (GEO) 56. First, we used GSE47035, a study of

247

killifish exposed to arsenic during acclimation to hyper-osmotic shock 15. Second, we

248

used GSE11056, a study of gene expression responses in mouse lung tissue 8 h after IP

249

injection with arsenic (1 mg/kg), or dexamethasone (1 mg/kg) or both. We chose this

250

study because we are familiar with the experimental details of the data, and because

ACS Paragon Plus Environment

Page 14 of 39

Page 15 of 39

Environmental Science & Technology

251

arsenic was combined with dexamethasone, a synthetic steroid hormone similar to

252

cortisol which evokes a broad stress response 57.

253

Results and Discussion

254

Arsenic’s diverse toxic effects in human populations range from cancer, heart disease,

255

and diabetes to reduced IQ 6, suggesting that environmental exposure to arsenic might

256

affect organisms living in dynamic environments in many different ways. Our

257

toxicogenomic analysis in killifish reveals that arsenic interferes with gene expression

258

responses to osmotic shock at three levels of biological organization: individual genes,

259

regulatory networks, and biological pathways. Publicly available data, analyzed with our

260

techniques, suggest that arsenic’s interference with gene expression responses extends to

261

other species.

262

We had previously observed that acute exposure to high doses of arsenic (12 mg/l) during

263

acclimation from fresh water to salt water led to high rates of mortality 39, although this

264

arsenic concentration failed to increase mortality in fish maintained in stable conditions,

265

suggesting adverse interactions between arsenic and osmotic shock. We have also shown

266

that 100 µg/l arsenic interferes with the induction of serum and glucocorticoid-regulated

267

kinase 1 (SGK1) that normally occurs when fresh water acclimated killifish are exposed

268

to salt water36. Finally, in a recent whole transcriptome experiment in a single killifish

269

population we showed that arsenic’s impact on gene expression during salinity

270

acclimation extends far beyond SGK1 by identifying roughly 400 killifish genes with a

271

significant interactions between exposure to arsenic (100 µg/l) and changed salinity 15. In

272

the present study, we assessed the impact of arsenic on gene expression in three

ACS Paragon Plus Environment

Environmental Science & Technology

273

populations of sea water acclimated killifish during acclimation to fresh water, looking

274

for arsenic-mediated differences in gene expression.

275

Identifying Differentially Expressed Genes

276

This report presents methods to identify biologically relevant patterns in gene expression

277

responses. As patterns require many genes to define, we began by estimating how many

278

genes respond to each treatment condition. Assuming that the roughly 16,000 genes

279

measured on our array respond independently, about 800 (5%) of these genes are

280

expected to achieve a p value of 0.05 or less in any statistical test under the null

281

hypothesis. Recognizing this, toxicogenomic studies often use false discovery rate (FDR)

282

calculations 45 to reduce the number of false positive calls, but this approach can also

283

increase false negative calls, eliminating genes whose expression truly varies and making

284

downstream analysis problematic 58. FDR calculations may therefore underestimate the

285

number of genes differentially expressed in a toxicogenomic experiment.

286

Our design included 24 fish from each of three separate populations, making our results

287

more likely to represent killifish in general that an experimental design using a single

288

population. Although this study was not intended to explore population differences in

289

detail, fixed effect linear models were used to independently assess each of the three

290

populations (Northeast Creek, Maine (NEC), King’s Creek Virginia (KC), and Horseshoe

291

Cove, Maine (HC)). These results were compared to mixed effect linear models treating

292

population as a random effect.

293

Mixed effects linear models were more sensitive to differential gene expression than

294

independent, fixed effect models, identifying 178 genes that responded significantly

ACS Paragon Plus Environment

Page 16 of 39

Page 17 of 39

Environmental Science & Technology

295

(FDR < 0.05) to fresh water exposure at 24 h (Figure 1A, red bar) compared to a total of

296

105 unique genes identified by fixed effect models in either Northeast Creek (green bar)

297

King’s Creek (blue bar) or Horseshoe Cove (purple bar). Neither fixed nor mixed effect

298

linear models identified large numbers of arsenic responsive genes or genes with arsenic

299

interactions. This either suggests that 100 µg/l arsenic affects very few genes in the

300

killifish gill, or that our design and statistics simply failed to detect arsenic responsive

301

genes.

302 303

Figure 1. A Number of significant genes (FDR < 0.05) identified by mixed effect

304

linear models (Mixed) or by fixed effect linear models of each killifish population

305

(NEC = Northeast Creek, KC = King’s Creek, HC = Horseshoe Cove). Counts

306

shown for each model term including time of exposure to fresh water, exposure to

ACS Paragon Plus Environment

Environmental Science & Technology

307

arsenic or interactions between arsenic and time (As : 1h, As : 24 h). B Distribution

308

of mixed effect p values for all 16,104 unique genes associated with each model

309

term. Horizontal line shows expected number of genes (100) in each bin assuming

310

uniform distribution. Vertical line demarks a nominal p value of 0.05.

311

Complete mixed effect linear model output is available in Supplemental Table 1.

312

An analysis of p value distributions suggests that using FDR correction to identify

313

treatment effects underestimates the number of arsenic-responsive genes, particularly

314

with respect to interactions between arsenic exposure and exposure to fresh water for 1 h.

315

Figure 1B shows the observed number of p values within a specified range and the

316

theoretical expectation (red horizontal line) based on the null hypothesis and gene

317

expression independence. P value bin ranges in Figure 1B were selected to contain

318

exactly 100 genes in the case that the null hypothesis is universally true, that is, in the

319

case where there are no true treatment effects. The area under the curve in Figure 1B that

320

lies above the red line therefore represents the number of genes that achieved significance

321

for some other reason than simultaneously testing many genes. For all model terms, the

322

total area under the curve, to the left of the vertical bars (p < 0.05) and above the

323

horizontal line identifies 1,793 instances of significance that are not attributable to the

324

multiple hypotheses testing burden -- about 8 times as many significantly differentially

325

expressed genes as were identified in Figure 1A by mixed effect models with an FDR