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