Subscriber access provided by Purdue University Libraries
Article
Novel cadmium responsive microRNAs in Daphnia pulex Shuai Chen, Garrett Justin McKinney, Krista M. Nichols, John K Colbourne, and Maria S Sepulveda Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b03988 • Publication Date (Web): 09 Nov 2015 Downloaded from http://pubs.acs.org on November 12, 2015
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
Novel cadmium responsive microRNAs in Daphnia pulex
1 2 3
Shuai Chen1, Garrett J. McKinney2, Krista M. Nichols1,3, John K. Colbourne4 and Maria S. Sepúlveda1*
4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
1
Department of Forestry and Natural Resources, Purdue University, West Lafayette, Indiana, USA 2 School of Aquatic and Fishery Sciences, University of Washington, Seattle, Washington State, USA 3 Conservation Biology Division, Northwest Fisheries Science Center, National Marine Fisheries Service, National Ocean and Atmospheric Administration, Seattle, Washington State, USA 4 Environmental Genomics Group, School of Biosciences, University of Birmingham, Edgbaston, Birmingham, UK *Corresponding Author: Maria S. Sepúlveda Address: FORS Room 103 195 Marsteller St West Lafayette, IN 47907-2033 Phone: 765.496.3428 Fax: 765.496.2422 Email:
[email protected]
23 24 25 26 27 28 29 30 31 32
1
ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 34
33
Abstract
34
Daphnia pulex is a widely used toxicological model and is known for its sensitivity to
35
cadmium (Cd). Recent research suggests microRNAs (miRNAs) play a critical role
36
in animal responses to heavy metals.
37
miRNAs under Cd exposure, we analyzed the miRNA profiles of D. pulex after 48
38
hours using miRNA microarrays and validated our findings by q-PCR.
39
dpu-let-7 was identified as a stably expressed gene and used as a reference. We
40
identified 22 and 21 differentially expressed miRNAs under low (20 µg/L CdCl2) and
41
high (40) exposure concentrations compared to controls, respectively.
42
functions of predicted miRNA target Cd-responsive genes included oxidative stress,
43
ion transport, mitochondrial damage and DNA repair.
44
also identified in relation to several Cd-responsive miRNAs. The expression of three
45
predicted target genes for miR-71 and miR-210 were evaluated and expression of two
46
of them (SCN2A and SLC31A1) was negatively correlated with the expression of their
47
regulator miRNAs.
48
propose Cd and hypoxia induce miR-210 via a same HIF1α modulated pathway.
49
Collectively, this research advances our understanding on the role of miRNAs in
50
response to heavy metal exposure.
To investigate the functions of D. pulex
miRNA
Cellular
An insulin-related network was
We show miR-210 is hypoxia-responsive in D. pulex and
51 52 53
2
ACS Paragon Plus Environment
Page 3 of 34
Environmental Science & Technology
Introduction
54 55
Cadmium (Cd) is highly toxic to aquatic organisms1 and ranks seventh on the
56
substance priority list developed by the Agency for Toxic Substances and Disease
57
Registry (ATSDR)2. In aquatic animals, Cd ions are up-taken through passive or
58
facilitated diffusion3.
59
therefore perturb Ca+2 uptake, transport and homeostasis4,5,6. Furthermore, Cd ions
60
increase the production of reactive oxidative species (ROS) leading to oxidative
61
stress, DNA damage, enhanced lipid peroxidation and apoptosis7,8.
62
research also suggests ROS inducers such as hypoxia might also lead to cellular
63
stress via similar pathways9.
64
Within cells, they act as calcium (Ca+2) antagonists and
Previous
MicroRNAs (miRNAs) are a family of short non-coding endogenous RNAs that
65
play a critical role in response to heavy metals including Cd10,11.
66
post-transcriptional regulators, miRNAs are involved in a variety of cellular processes
67
via negative regulation of their target mRNA12,13.
68
gene expression by miRNAs requires the binding of miRNA to the target gene’s
69
3’-untranslated region (3’-UTR), resulting in cleavage or suppressed translation14. A
70
number of heavy metal responsive miRNAs have been identified in vertebrates; their
71
functions have been associated with regulatory transcription factors, metabolic
72
processes and stress responses15,16,17.
73
miRNAs and heavy metals makes miRNAs promising molecules for elucidating the
74
mechanisms involved in heavy metal response, given its role in gene regulation.
75
However, the role of miRNAs in response to heavy metal exposure in aquatic
As key
Generally, the inhibition of target
The rapidly growing evidence linking
3
ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 34
76
invertebrates is largely unknown, despite the fact that heavy metals are persistent in
77
aquatic environments and can bioaccumulate and biomagnify in the food chain18,19.
78
A better understanding on the regulatory role of miRNAs in these organisms after
79
heavy metal exposure is needed.
80
Daphnia pulex is a branchiopod crustacean that resides in continental waters and
81
are widely used as a toxicological model organism. This species is also the first
82
crustacean to have its genome sequenced and annotated20, thereby represents an
83
ideal candidate for aquatic invertebrate Cd-responsive miRNA research.
84
cornerstone species in freshwater ecosystems, Daphnia are important sentinels and
85
have much accumulated knowledge on their acute responsiveness and rapid
86
acclimation to Cd21,22,23. In daphiids, acclimation to Cd is thought to be mostly driven
87
by upregulation of metallothionein genes24.
88
miRNAs play critical roles during D. pulex development25, however, no studies have
89
examined the role of miRNAs on Cd acclimation in daphniids.
90
work are two-fold: 1) identify Cd-responsive D. pulex miRNAs; and 2) study the
91
potential regulatory mechanisms of these miRNAs under Cd exposure.
92
studies are the first to report on the role of miRNAs in aquatic invertebrates in
93
response to heavy metals.
As a
Our previous research revealed that
The objectives of this
These
94
Materials and Methods
95 96 97
Animal culture D. pulex were obtained from Carolina Biological (Burlington, NC, USA).
4
ACS Paragon Plus Environment
Page 5 of 34
Environmental Science & Technology
98
Cultures were maintained at a density of approximately 35 individuals/L in
99
reconstituted moderate hard water (RMHW)26 in 3 L plastic tanks in an environmental
100
chamber (20°C on a 16/8 light/dark cycle).
Culture water was changed twice a week
101
and D. pulex were fed 15 ml YCT mix (yeast, cereal, trout chow) after every water
102
change. In order to prevent potential maternal or grand maternal co-variance, D.
103
pulex were grown in the environmental chamber for four generations before the start
104
of exposures and only neonates (< 24 h) from the third through eighth brood were
105
used for exposures.
106
Daphnia tests and sample preparation
107
Cadmium exposure
Animals were not fed during the exposure test.
108
Cadmium exposures were conducted by exposing 70 neonates to 2 L of either 20
109
µl/L CdCl2 (Cd-low, ~48 h LC01 value), 40 µl/L CdCl2 (Cd-high, ~48 h LC10 value), or to
110
RMHW (control) for 48 h24.
111
were diluted from a stock solution of 100 mg/L CdCl2 (Sigma-Aldrich, St. Louis, MO,
112
USA) using RMHW.
113
(ThermoFinnigan, Bremen, Germany) mass spectrometer in the medium resolution
114
mode27 and all measured concentrations varied within 10% of the nominal
115
concentrations.
116
Hypoxia exposure
117
Each treatment included six replicates and Cd solutions
Cd concentrations were tested using an ELEMENT-2
Predicted targets of Cd-responsive miR-210 suggest it should also be involved in
118
D. pulex response to hypoxia.
Because both hypoxia and Cd can induce ROS
119
production28,29 and mammalian miR-210 is known to be up-regulated by
5
ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 34
120
hypoxia-induced ROS production30, we tested the hypothesis that miR-210 is
121
hypoxia-inducible in D. pulex.
122
our hypothesis.
A follow up hypoxia test was then carried out to test
123
Hypoxic exposures were carried out by introducing 70 neonates to 2 L of either
124
hypoxic RMHW (dissolved oxygen, DO, 1.0-1.5 mg/L), or normoxic RMHW (DO,
125
8.0-9.0 mg/L) for 48 h. Hypoxic or normoxic RMHW was produced by aerating N2 or
126
fresh air into RMHW31, respectively.
127
optical dissolved oxygen meter (YSI, Yellow Springs, Ohio, USA).
128
had three replicates.
DO was measured hourly using a YSI ProODO Each treatment
129
Immediately after exposure, 50 neonates from each replicate were randomly
130
selected and pooled for small RNA (including miRNAs) extraction using protocols
131
already published25.
132
sample was extracted using a Purelink Micro-to-Midi total RNA purification kit
133
(Invitrogen, Carlsbad, CA, USA) following the manufacture’s protocol. Small RNA
134
and total RNA extracts were treated with DNAse I (ThermoScientific, Waltham, MA,
135
USA) following manufacture’s protocols.
136
were quantified using NanoDrop (Invitrogen, Carlsbad, CA, USA) and stored at -80°C
137
until microarray or q-PCR analysis as described below.
138
miRNA microarray analysis
Total RNA (with exception of small RNAs) from the same
After DNAse treatment, all RNA samples
139
Affymetrix GeneChip 4.0 (Affymetrix, Santa Clara, CA, USA) miRNA microarrays
140
(catalog number: 902412) were used to determine differences in miRNA expression
141
upon exposure to Cd.
Microarray analyses were conducted at the Center for
6
ACS Paragon Plus Environment
Page 7 of 34
Environmental Science & Technology
142
Medical Genomics at Indiana University School of Medicine, Indianapolis, IN.
Small
143
RNA was extracted as already described and 300 ng from each of 3
144
samples/treatment was labeled using the Affymetrix FlashTag Bioton HSR kit.
145
Labeled samples were hybridized to the arrays for 17 hours and stained and washed
146
using standard protocols32.
147
used to scan the arrays and generate CEL files containing information on the intensity
148
values of the individual probes. Arrays were visually scanned for abnormalities or
149
defects.
150
USA) to generate RMA (robust multi-array average) expression levels using the RMA
151
normalization option32. Only miRNAs with a Log2 average expression level (Log2
152
AE) >2 were further investigated. Significant changes in miRNA expression were
153
identified by one-way ANOVA followed by Tukey’s HSD post-hoc test (Cd-low vs
154
Control, Cd-high vs Control) with Benjamini and Hochberg correction33 (FDR
167
4).
168
example, cte-miR-279 (log2AE = 13.34) and dwi-miR-279 (log2AE = 5.90) were
169
grouped in the miR-279 family.
170
sequence of miRNA with highest log2AE as a novel D. pulex miRNA sequence (e.g.
171
cte-miR-279 had the highest log2AE level in the miR-279 family, thus the sequence of
172
cte-miR-279 was selected as the miR-279 sequence in D. pulex. The probes on the
173
microarray differ from the annotated D. pulex locus for this gene by only 3 of 21
174
aligned nucleotides in the 3’ end.
175
miRNA cDNA preparation
To identify robustly expressed D. pulex miRNAs, we
Secondly, we grouped the robustly expressed miRNAs by family name. For
And lastly, from each miRNA family, we selected the
176
miRNA-specific stem-loop primers and forward primers (Table S2) were designed
177
following the method described by Chen et al.36. miRNAs were reversed-transcribed
178
by mixing 2 pmol of stem-loop primer, 0.5 mM dNTPs, 2 µg of RNA template and
179
incubated at 65°C for 5 min.
180
of RNase OUT and 2.5 units of SuperScript III were then added and incubated on ice
181
for 2 min.
182
30°C for 30 sec, 42°C for 30 sec and 50°C for 1 sec and then incubated at 85°C for 5
183
min.
184
-20°C until q-PCR analysis.
185
q-PCR reference gene selection and validation
Another mix of 1X first-strand buffer, 5 mM DTT, 2 units
The final mix was incubated at 16°C for 30 min followed by 60 cycles at
The reverse transcribed products were quantified by NanoDrop and stored at
8
ACS Paragon Plus Environment
Page 8 of 34
Page 9 of 34
Environmental Science & Technology
186
The protocol for miRNA q-PCR was adapted from Chen et al.25. Each q-PCR
187
reaction contained 4 µM SYBR Green I master mix (Bio-Rad, Hercules, CA, USA), 1
188
µM forward primer, 1 µM reverse primer and 2 µl RT product. q-PCR was conducted
189
on a Bio-Rad real-time PCR system (Bio-Rad, Hercules, CA, USA) using the following
190
program: 95°C for 5 min, 40 cycles of 95°C for 5 sec, 60°C for 10 sec and 72°C for 1
191
sec.
A melting curve analysis was performed to check for primer-dimers.
192
Although our previous research identified U6 ribosomal RNA as a stable
193
expressed reference gene during D. pulex development25, this gene is not stably
194
expressed under Cd exposure (data not shown).
195
gene selection process to identify stably expressed miRNAs under Cd exposure.
196
Fourteen highly expressed (log2AE > 4) but not differentially expressed miRNAs (let-7,
197
miR-1, miR-2, miR-7, miR-8, miR-9, miR-10, miR-12, miR-33, miR-34, miR-100,
198
miR-153, miR-283, and miR-317) identified by microarray were selected as candidate
199
reference genes.
200
Cd-high and control groups following the method described above.
201
included three biological replicates and each biological replicate included three
202
technical replicates. Two algorithms, geNorm37 and NormFinder38, were performed
203
to identify stably expressed reference gene under Cd exposure.
204
candidate genes were transformed to ∆Ct values and imported into geNorm.
205
geNorm calculates an M value for each gene, which represents the combined
206
variation within treatment group and between candidate genes.
207
expressing the smallest M values have the highest expression stability.
Thus, we conducted a reference
q-PCR was conducted using miRNA samples from the Cd-low,
9
ACS Paragon Plus Environment
Each treatment
The Ct values of
Transcripts The Ct
Environmental Science & Technology
Page 10 of 34
208
values were transformed into relative quantities and imported to NormFinder.
The
209
NormFinder algorithm calculates a stability (S) value for each gene based on the
210
inter- and intra-group variance.
211
expression stability.
212
of determination.
The gene with the lowest S value has the highest
The correlation between M and S was measured by coefficient
213
Seven differentially expressed miRNAs were selected for q-PCR validation, and
214
included six significantly upregulated miRNAs (miR-71, miR-92, miR-204, miR-210,
215
miR-252, and miR-279), and one significantly down-regulated miRNA (miR-7444).
216
q-PCR was performed using three miRNA samples from each treatment group with
217
three technical replicates for each sample following the method already described.
218
Fold-changes (Cd-low vs Control, Cd-high vs Control) were calculated based on the
219
expression of controls using ∆∆Ct method39 and normalized to the expression of
220
Dpu-let-7, the gene identified as the most stably expressed.
221
test, q-PCR was performed using three miRNA samples from hypoxia and normoxia
222
groups with three technical replicates.
223
calculated based on the expression of the normoxia group (control) using the ∆∆Ct
224
method39 and normalized to the expression of Dpu-let-7.
225
performed
226
Cd/hypoxia-exposed animals (p 4) in
291
D. pulex (Table S4), which included all previously predicted D. pulex miRNAs25 plus
292
the highly conserved dpu-let-7, which had not been predicted earlier.
293
miRNA clusters such as the miR-71/2 cluster were identified, supporting prior
294
research that the miR-71/2 cluster is Protostome-specific51.
295
Reference gene selection and q-PCR validation of microarray results
13
ACS Paragon Plus Environment
Several
Environmental Science & Technology
Page 14 of 34
296
Dpu-let-7 was identified as the most stably expressed reference gene by both
297
geNorm and NormFinder, having both the smallest M and S values, respectively
298
(Figure S2).
299
therefore considered as stably expressed), which further confirmed microarray results
300
that the expression of these miRNAs is no affected under Cd exposure. Let-7 has
301
also been identified as a stably expressed reference gene in other species52, which
302
further confirms its expression stability.
303
were able to validate microarray results by q-PCR; for all genes tested, the q-PCR
304
expression patterns matched the microarray results (Figure 2).
305
Predicted miRNA target genes and their relation to Cd toxicityf
306
All 14 candidate reference genes had M values of less than 1.5 (and
Using dpu-let-7 as a reference gene, we
We focused on highly conserved miRNAs related to Cd-induced cellular stress for
307
target gene prediction and validation of expression.
308
during computational target prediction, we compared miRNA targets predicted by both
309
the thermodynamic method and the seed complementary method. Only genes that
310
were 1) predicted by both methods; or 2) had multiple seed binding sites on its 3’UTR
311
for one specific miRNA were selected as promising miRNA targets.
312
inspection of predicted targets for the Cd-responsive miRNAs illustrates their function
313
is related with Cd-induced cellular stress that includes 1) alteration of ion hemostasis;
314
2) mitochondrial damage; and 3) DNA damage (Table S5).
315
To diminish false positives
A detailed
Alteration of ion hemostasis: Cd is a known Ca+2 antagonist that can block Ca+2
316
channels and modulate ion homeostasis5.
317
and miR-71 are likely involved in Cd-induced alteration of ion homeostasis by
Predicted target genes indicate miR-210
14
ACS Paragon Plus Environment
Page 15 of 34
Environmental Science & Technology
318
targeting
319
calcium-transporting ATPase type 2C (ATP2C2) and Ca+2-transporting ATPase
320
sarcoplasmic/endoplasmic reticulum type (Ca-P60A).
321
miR-71 should decrease ATP2C2 and Ca-P60A expression.
322
supported by previous research showing that Cd can inhibit the activity of Ca+2
323
ATPase53,54.
324
sodium
325
Ca+2-activated potassium channel protein (SK channel protein, KCNN1), sodium
326
channel
327
hyperpolarization-activated cyclic nucleotide-gated channel 2 (HCN2) are predicted to
328
be the targets of miR-210.
329
ion hemostasis.
330
Cd-induced miR-71, and miR-210 could down-regulate the sodium channel protein by
331
acting synergistically.
332
can inhibit sodium channel activity55.
333
channels, miR-210 and miR-71 might also co-regulate copper transport by targeting
334
copper-transporting ATPase 1 gene (ATP7A) and the high affinity copper transporter
335
gene (SLC31A1), respectively (Table S5).
336
a
and
number
of
key
ion
transport
proteins(Table
S5)
such
as
Cd up-regulated miR-210 and This finding is
Predicted targets also suggest that miR-210 and miR-71 regulate potassium
protein
ion
type
2
channels.
subunit
alpha
For
instance,
(SCN2A)
and
small
conductance
potassium/sodium
These ion channels play a critical role in maintenance of
SCN2A is also predicted as a target of miR-71, which suggests
This finding is in agreement with previous research that Cd In addition to Ca+2, sodium and potassium
Mitochondrial and DNA damage: Mitochondrial dysfunction is another important
337
aspect of Cd toxicity.
MiR-252 appears to be related to Cd-induced mitochondrial
338
damage by targeting two mitochondrial ribosomal proteins (MRPs)(Table S5).
339
Cd-induced miR-252 expression could decrease the expression of MRPs, which
15
ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 34
340
would cause inhibition of mitochondrial function. However, further investigation is
341
needed to test this hypothesis.
342
induce DNA damage by inhibiting mismatch repair56,57.
343
miR-210 target is DNA mismatch repair protein (MLH3) which is known to be inhibited
344
by Cd57.
345
Target gene expression profile
346
As a mutagen and carcinogen, Cd is also known to Another notable predicted
To examine the negative correlation of miRNA-target expression, we quantified
347
the expression of three predicted target genes: SCN2A, KCNN1 and SLC31A1.
As
348
expected, the expression of SCN2A was negatively correlated with that of its regulator
349
miR-210 (Figure 3). Similarly, the expression of SLC31A1 was negatively correlated
350
with the expression of its regulator miRNA miR-71, but only in animals from the
351
Cd-low treatment.
352
promising targets for miR-210 and miR-71.
353
significantly changed under either Cd treatment, which could be explained by
354
involvement of other regulatory mechanisms (e.g. transcription factors and DNA
355
methylation)58 or a false positive computational prediction.
356
validation such as luciferase assays are needed to prove these miRNA-target
357
relationships.
358
Cd-responsive miRNAs regulatory network
This result further confirms that SCN2A and SLC31A1 are The expression of KCNN1 was not
Further experimental
359
Our IPA analysis revealed only one significant network (p-score=22, p=10-22)
360
involving several Cd-responsive miRNAs in relation to insulin metabolism (Figure S3,
361
Table S6).
For instance, Cd-induced miR-92 (synonym of miR-25) and miR-204 can
16
ACS Paragon Plus Environment
Page 17 of 34
Environmental Science & Technology
362
block insulin production59,60 which results in an up-regulation of miR-210 and
363
miR-61561. In addition, down-regulation of miR-144 decreases the activity of PGTS2
364
(prostaglandin endoperoxide synthase 2), which can also reduce the secretion of
365
insulin.
366
growth factor, beta 1) which causes a down-regulation of miR-21662 and miR-21763.
367
The insulin pathway plays a critical role in predator-induced body size plasticity in D.
368
pulex64,65.
369
development66,67.
370
reductions of D. pulex body length after Cd exposure24,68. In addition, IPA revealed
371
that the “Top Diseases and Bio Function” related with differentially expressed miRNAs
372
is “organismal injury and abnormalities” which is related with Cd induced cellular
373
stress69.
374
Cd and hypoxia induce miRNA-210 via the HIFα pathway
375
Further, reduced insulin levels decrease the activity of TGFB1 (transforming
TGFB1 regulates cell growth and proliferation as well as skeleton Thus, a reduction of insulin and TGFB1 levels might explain
Other than the Cd-related genes mentioned above, a hypoxia up-regulated gene
376
HYOU1 was also identified as a target of miR-210.
This finding suggests that
377
miR-210 could also modulate cellular responses under hypoxia in D. pulex.
378
miR-210 was previously identified as the predominantly up-regulated miRNA under
379
hypoxia in mammalian cell lines30, we speculated that miR-210 should also be
380
up-regulated in D. pulex under hypoxia.
381
miR-210 expression under hypoxic and normoxic conditions and observed a
382
significant up-regulation of miR-210 under hypoxia (Figure 4).
383
that hypoxia induces miR-210 via a conserved pathway across taxa.
Because
To test this prediction, we profiled the
17
ACS Paragon Plus Environment
Our results indicate
Environmental Science & Technology
384
Hypoxia-induced
miR-210 1
expression α
(HIF1α)
occurs protein
Page 18 of 34
via
stabilization
which
binds
of
385
hypoxia-inducible-factor
to
386
hypoxia-responsive-elements (HRE) on the miR-210 promoter region70,71.
387
the HRE sequence ((A/G)CGTG) in the miR-210 promoter region is highly conserved
388
across mammalian species72.
389
HIF1α by increased ROS production28.
390
stabilization of HIF1α, which promotes miR-210 expression, we hypothesized that the
391
up-regulation of D. pulex miR-210 is due to Cd-induced HIF1α stabilization (Figure 4).
392
To address this hypothesis, we searched for the D. pulex homologs by BLAST73 of the
393
human HIF1α gene against the D. pulex genome and identified a homologous HIF1α
394
gene in D. pulex (E Value= 1e-44).
395
documented in OrthoDB74 (protein ID EFX84860).
396
miR-210 promoter region within the D. pulex genome and found a HRE motif ~400 bp
397
upstream from the coding region. We also analyzed the miR-210 promoter region of
398
D. magna - a closely related species to D. pulex.
399
motifs on the D. magna miR-210 promoter region, which implicate this HRE motif is
400
also conserved in Daphnia.
401
responsive to hypoxia and indicates a conserved regulatory pathway of HIF1α
402
induced miR-210 expression, which further support our hypothesis.
Further,
Previous research has shown Cd can also induce Since hypoxia and Cd can both induce
This homologous HIF1α gene has also been Secondly, we searched the
Our search identified two HRE
Our investigation revealed that miR-210 is also
403
In conclusion, our research provides the basis for D. pulex miRNA research,
404
including identification of a stably expressed reference gene under Cd exposure and a
405
list of Cd-responsive miRNAs with potential targets. To our knowledge, this is the
18
ACS Paragon Plus Environment
Page 19 of 34
Environmental Science & Technology
406
first study on miRNA responses in aquatic invertebrates after heavy metal exposure.
407
Finally, we have shown that miR-210 is also responsive to hypoxia via a proposed
408
HIFα modulated pathway.
409
environmental stress.
Our results indicate miR-210 is a promising marker of
410 411
Acknowledgements The authors would like to thank Dr. Don Gilbert, Indiana University, IN for
412 413
elaborations on the D. pulex genome annotations.
The authors gratefully
414
acknowledge funding from the Department of Forestry and Natural Resources,
415
Purdue University, IN in the form of an assistantship for Shuai Chen. Microarray
416
studies were carried out at the Center for Medical Genomics, Indiana University
417
School of Medicine, IN which is partially supported by the Indiana Genomic Initiative
418
at Indiana University (INGEN); INGEN is supported in part by the Lilly Endowment,
419
Inc. Our work benefits from and contributes to the Daphnia Genomics Consortium.
420 421
References
422 423 424
(1)
Borgmann, U.; Couillard, Y.; Doyle, P.; Dixon, D. G. Toxicity of sixty-three metals and metalloids to Hyalella azteca at two levels of water hardness. Environ. Toxicol. Chem. 2005, 24, 641–652.
425 426
(2)
Agency for Toxic Substances and Disease Registry Home Page. http://www.atsdr.cdc.gov/SPL/index.html (Accessed Dec 5, 2014).
427 428 429
(3)
Olsson, P.-E.; Hogstran, C. Subcellular distribution and binding of cadmium to metallothionein in tissues of rainbow trout after exposure to Cd in water. Environ. Toxicol. Chem. 1987, 6, 867–874.
19
ACS Paragon Plus Environment
Environmental Science & Technology
430 431 432
(4)
Shen, J. B.; Jiang, B.; Pappano, A. J. Comparison of L-type calcium channel blockade by nifedipine and/or cadmium in guinea pig ventricular myocytes. J. Pharmacol. Exp. Ther. 2000, 294, 562–570.
433 434 435
(5)
Verbost, P. M.; Flik, G.; Lock, R. A.; Wendelaar Bonga, S. E. Cadmium inhibition of Ca2+ uptake in rainbow trout gills. Am. J. Physiol. 1987, 253, R216–R221.
436 437 438
(6)
Verbost, P. M.; Flik, G.; Lock, R. A. C.; Wendelaar Bonga, S. E. Cadmium inhibits plasma membrane calcium transport. J Membr Biol. 1988, 102, 97– 104.
439 440 441 442
(7)
Risso-De Faverney, C.; Devaux, A.; Lafaurie, M.; Girard, J. P.; Bailly, B.; Rahmani, R. Cadmium induces apoptosis and genotoxicity in rainbow trout hepatocytes through generation of reactive oxygene species. Aquat. Toxicol. 2001, 53, 65–76.
443 444 445
(8)
Sarkar, S.; Yadav, P.; Trivedi, R.; Bansal, A. K.; Bhatnagar, D. Cadmium-induced lipid peroxidation and the status of the antioxidant system in rat tissues. J. Trace Elem. Med. Biol. 1995, 9, 144–149.
446 447 448 449
(9)
Onukwufor, J. O.; MacDonald, N.; Kibenge, F.; Stevens, D.; Kamunde, C. Hypoxia-cadmium interactions on rainbow trout (Oncorhynchus mykiss) mitochondrial bioenergetics: attenuation of hypoxia-induced proton leak by low doses of cadmium. J. Exp. Biol. 2014, 217, 831–840.
450 451
(10)
Ding, Y.-F.; Zhu, C. The role of microRNAs in copper and cadmium homeostasis. Biochem. Biophys. Res. Commun. 2009, 386, 6–10.
452 453
(11)
Hou, L.; Wang, D.; Baccarelli, A. Environmental chemicals and microRNAs. Mutat. Res. 2011, 714, 105–112.
454 455 456
(12)
Kagias, K.; Podolska, A.; Pocock, R. Reliable reference miRNAs for quantitative gene expression analysis of stress responses in C. elegans. BMC Genomics 2014, 15, 222.
457 458
(13)
Choudhuri, S. Small noncoding RNAs: Biogenesis, function, and emerging significance in toxicology. J. Biochem. Mol. Toxic. 2010, 24, 195–216.
459 460 461
(14)
Guo, H.-S.; Xie, Q.; Fei, J.-F.; Chua, N.-H. MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell 2005, 17, 1376–1386.
20
ACS Paragon Plus Environment
Page 20 of 34
Page 21 of 34
Environmental Science & Technology
462 463 464
(15)
Wang, L.; Bammler, T. K.; Beyer, R. P.; Gallagher, E. P. Copper-induced deregulation of microRNA expression in the zebrafish olfactory system. Environ. Sci. Technol. 2013, 47, 7466–7474.
465 466 467
(16)
Fabbri M., Urani C., Sacco M. G., Procaccianti C., Gribaldo L. Whole Genome Analysis and MicroRNAs Regulation in HepG2 Cells Exposed to Cadmium. ALTEX 2012, 29, 173–182.
468 469 470
(17)
Meng, X. Z.; Zheng, T. Sen; Chen, X.; Wang, J. B.; Zhang, W. H.; Pan, S. H.; Jiang, H. C.; Liu, L. X. microRNA expression alteration after arsenic trioxide treatment in HepG-2 cells. J. Gastroenterol. Hepatol. 2011, 26, 186–193.
471 472 473
(18)
Müller, M.; Anke, M. Distribution of cadmium in the food chain (soil-plant-human) of a cadmium exposed area and the health risks of the general population. Sci. Total Environ. 1994, 156, 151–158.
474 475
(19)
Rainbow, P. Trace metal concentrations in aquatic invertebrates: why and so what? Environ. Pollut. 2002, 120, 497–507.
476 477 478
(20)
Colbourne, J. K.; Pfrender, M. E.; Gilbert, D.; Thomas, W. K.; Tucker, A.; Oakley, T. H.; Tokishita, S.; Aerts, A.; Arnold, G. J.; Basu, M. K.; et al. The ecoresponsive genome of Daphnia pulex. Science. 2011, 331, 555–561.
479 480
(21)
Muyssen, B. T. a; Janssen, C. R. Multi-generation cadmium acclimation and tolerance in Daphnia magna Straus. Environ. Pollut. 2004, 130, 309–316.
481 482
(22)
Ward, T. J.; Robinson, W. E. Evolution of cadmium resistance in Daphnia magna. Environ. Toxicol. Chem. 2005, 24, 2341–2349.
483 484 485 486
(23)
Poynton, H. C.; Loguinov, A. V.; Varshavsky, J. R.; Chan, S.; Perkins, E. J.; Vulpe, C. D. Gene expression profiling in Daphnia magna Part I: concentration-dependent profiles provide support for the no observed transcriptional effect level. Environ. Sci. Technol. 2008, 42, 6250–6256.
487 488 489 490
(24)
Shaw, J. R.; Colbourne, J. K.; Davey, J. C.; Glaholt, S. P.; Hampton, T. H.; Chen, C. Y.; Folt, C. L.; Hamilton, J. W. Gene response profiles for Daphnia pulex exposed to the environmental stressor cadmium reveals novel crustacean metallothioneins. BMC Genomics 2007, 8, 477.
491 492 493
(25)
Chen, S.; McKinney, G. J.; Nichols, K. M.; Sepúlveda, M. S. In silico prediction and in vivo validation of Daphnia pulex microRNAs. PLoS One 2014, 9, e83708.
21
ACS Paragon Plus Environment
Environmental Science & Technology
494 495 496
(26)
EPA US: Methods for measuring the acute toxicity of effluents and receiving waters to freshwater and marine organisms, 5th addition. Washington, DC: U.S. EPA, Office of Water; 2002:275.
497 498 499
(27)
Wood, K. V.; Bonham, C. C.; Ng, J.; Hipskind, J.; Nicholson, R. Plasma desorption mass spectrometry of anthocyanidins. Rapid Commun. Mass Spectrom. 1993, 7, 400–403.
500 501 502 503 504
(28)
Jing, Y.; Liu, L. Z.; Jiang, Y.; Zhu, Y.; Guo, N. L.; Barnett, J.; Rojanasakul, Y.; Agani, F.; Jiang, B. H. Cadmium increases HIF-1 and VEGF expression through ROS, ERK, and AKT signaling pathways and induces malignant transformation of human bronchial epithelial cells. Toxicol. Sci. 2012, 125, 10– 19.
505 506 507
(29)
Kim, J. H.; Park, S. G.; Song, S.-Y.; Kim, J. K.; Sung, J.-H. Reactive oxygen species-responsive miR-210 regulates proliferation and migration of adipose-derived stem cells via PTPN2. Cell Death Dis. 2013, 4, e588.
508 509
(30)
Devlin, C.; Greco, S.; Martelli, F.; Ivan, M. MiR-210: More than a silent player in hypoxia. IUBMB Life, 2011, 63, 94–100.
510 511 512
(31)
Gorr, T. A.; Cahn, J. D.; Yamagata, H.; Bunn, H. F. Hypoxia-induced synthesis of hemoglobin in the crustacean Daphnia magna is hypoxia-inducible factor-dependent. J. Biol. Chem. 2004, 279, 36038–36047.
513 514 515
(32)
Irizarry, R. A.; Hobbs, B.; Collin, F.; Beazer-Barclay, Y. D.; Antonellis, K. J.; Scherf, U.; Speed, T. P. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4, 249–264.
516 517 518
(33)
Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. Roy. Statist. Soc. Ser. A, 1995, 57, 289–300.
519 520
(34)
Kozomara, A.; Griffiths-Jones, S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39, D152–D157.
521 522 523
(35)
Wheeler, B. M.; Heimberg, A. M.; Moy, V. N.; Sperling, E. a; Holstein, T. W.; Heber, S.; Peterson, K. J. The deep evolution of metazoan microRNAs. Evol. Dev. 2009, 11, 50–68.
524 525 526 527
(36)
Chen, C.; Ridzon, D. A.; Broomer, A. J.; Zhou, Z.; Lee, D. H.; Nguyen, J. T.; Barbisin, M.; Xu, N. L.; Mahuvakar, V. R.; Andersen, M. R.; et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33, e179.
22
ACS Paragon Plus Environment
Page 22 of 34
Page 23 of 34
Environmental Science & Technology
528 529 530 531
(37) Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.; Speleman, F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3, 1-12.
532 533 534 535
(38)
Andersen, C. L.; Jensen, J. L.; Ørntoft, T. F. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64, 5245–5250.
536 537 538
(39)
Livak, K. J.; Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408.
539 540 541
(40)
Lewis, B. P.; Burge, C. B.; Bartel, D. P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120, 15–20.
542 543
(41)
Gruber, A. R.; Lorenz, R.; Bernhart, S. H.; Neuböck, R.; Hofacker, I. L. The Vienna RNA websuite. Nucleic Acids Res. 2008, 36, 70–74.
544 545
(42)
Wang, X. Composition of seed sequence is a major determinant of microRNA targeting patterns. Bioinformatics 2014, 30, 1377–1383.
546 547 548 549
(43)
Rozen, S.; Skaletsky, H. J. Primer3 on the WWW for General Users and for Biologist Programmers. In Bioinformatics Methods and Protocols: Methods in Molecular Biology; Misener, S. &, Krawetz, S. A., Eds.; Humana Press Inc: Totowa (NJ), 2000; pp 365–386.
550 551
(44)
Colbourne, J.K.; V.R. Singan and D. Gilbert. wFleaBase: the Daphnia genome database. BMC Bioinformatics 2005, 6, 45.
553 554 555 556
(45)
Spanier, K. I.; Leese, F.; Mayer, C.; Colbourne, J. K.; Gilbert, D.; Pfrender, M. E.; Tollrian, R. Predator-induced defences in Daphnia pulex: selection and evaluation of internal reference genes for gene expression studies with real-time PCR. BMC Mol. Biol. 2010, 11, 50.
557 558 559 560 561 562
(46)
Martínez, C.; Vicario, M.; Ramos, L.; Lobo, B.; Mosquera, J. L.; Alonso, C.; Sánchez, A.; Guilarte, M.; Antolín, M.; de Torres, I.; et al. The jejunum of diarrhea-predominant irritable bowel syndrome shows molecular alterations in the tight junction signaling pathway that are associated with mucosal pathobiology and clinical manifestations. Am. J. Gastroenterol. 2012, 107, 736–746.
552
23
ACS Paragon Plus Environment
Environmental Science & Technology
563 564 565
(47)
Hsieh, Y.-W.; Chang, C.; Chuang, C.-F. The microRNA mir-71 inhibits calcium signaling by targeting the TIR-1/Sarm1 adaptor protein to control stochastic L/R neuronal asymmetry in C. elegans. PLoS Genet. 2012, 8, e1002864.
566 567 568
(48)
Boulias K, Horvitz HR. The C. elegans microRNA mir-71 acts in neurons to promote germline-mediated longevity through regulation of DAF-16/FOXO. Cell Metab. 2012,15:439–450.
569 570 571 572
(49)
He, M.; Lu, Y.; Xu, S.; Mao, L.; Zhang, L.; Duan, W.; Liu, C.; Pi, H.; Zhang, Y.; Zhong, M.; et al. MiRNA-210 modulates a nickel-induced cellular energy metabolism shift by repressing the iron-sulfur cluster assembly proteins ISCU1/2 in Neuro-2a cells. Cell Death Dis. 2014, 5, e1090.
573 574 575 576
(50)
Zhou, R.; Yuan, P.; Wang, Y.; Hunsberger, J. G.; Elkahloun, A.; Wei, Y.; Damschroder-Williams, P.; Du, J.; Chen, G.; Manji, H. K. Evidence for selective microRNAs and their effectors as common long-term targets for the actions of mood stabilizers. Neuropsychopharmacology 2009, 34, 1395–1405.
577 578 579
(51)
Sangokoya, C.; Telen, M. J.; Chi, J.-T. microRNA miR-144 modulates oxidative stress tolerance and associates with anemia severity in sickle cell disease. Blood 2010, 116, 4338–4348.
580 581 582 583
(52)
De Souza Gomes, M.; Donoghue, M. T. A.; Muniyappa, M.; Pereira, R. V.; Guerra-Sá, R.; Spillane, C. Computational identification and evolutionary relationships of the MicroRNA gene cluster miR-71/2 in protostomes. J. Mol. Evol. 2013, 76, 353–358.
584 585 586 587
(53)
Chen, X.; Liang, H.; Guan, D.; Wang, C.; Hu, X.; Cui, L.; Chen, S.; Zhang, C.; Zhang, J.; Zen, K.; et al. A combination of Let-7d, Let-7g and Let-7i serves as a stable reference for normalization of serum microRNAs. PLoS One 2013, 8(11), e79652.
588 589 590
(54)
Hechtenberg, S.; Beyersmann, D. Interference of cadmium with ATP-stimulated nuclear calcium uptake. Environ. Health Perspect. 1994, 102, 265–267.
591 592 593
(55)
Biagioli, M.; Pifferi, S.; Ragghianti, M.; Bucci, S.; Rizzuto, R.; Pinton, P. Endoplasmic reticulum stress and alteration in calcium homeostasis are involved in cadmium-induced apoptosis. Cell Calcium 2008, 43, 184–195.
594 595 596
(56)
Takeda, A. N.; Gautschi, I.; Van Bemmelen, M. X.; Schild, L. Cadmium trapping in an epithelial sodium channel pore mutant. J. Biol. Chem. 2007, 282, 31928– 31936.
24
ACS Paragon Plus Environment
Page 24 of 34
Page 25 of 34
Environmental Science & Technology
597 598 599
(57)
Jin, Y. H.; Clark, A. B.; Slebos, R. J. C.; Al-Refai, H.; Taylor, J. A.; Kunkel, T. A.; Resnick, M. A.; Gordenin, D. A. Cadmium is a mutagen that acts by inhibiting mismatch repair. Nat. Genet. 2003, 34, 326–329.
600 601 602 603
(58) Nazarov, P. V.; Reinsbach, S. E.; Muller, A.; Nicot, N.; Philippidou, D.; Vallar, L.; Kreis, S. Interplay of microRNAs, transcription factors and target genes: Linking dynamic expression changes to function. Nucleic Acids Res. 2013, 41, 2817–2831.
604 605 606
(59)
Banerjee, S.; Flores-Rozas, H. Cadmium inhibits mismatch repair by blocking the ATPase activity of the MSH2-MSH6 complex. Nucleic Acids Res. 2005, 33, 1410–1419.
607 608 609
(60)
Setyowati Karolina, D.; Sepramaniam, S.; Tan, H. Z.; Armugam, A.; Jeyaseelan, K. miR-25 and miR-92a regulate insulin I biosynthesis in rats. RNA Biol. 2013, 10, 1365–1378.
610 611
(61)
Xu, G.; Chen, J.; Jing, G.; Shalev, A. Thioredoxin-interacting protein regulates insulin transcription through microRNA-204. Nat. Med. 2013, 19, 1141–1146.
612 613 614 615 616 617
(62)
Granjon, A.; Gustin, M. P.; Rieusset, J.; Lefai, E.; Meugnier, E.; Güller, I.; Cerutti, C.; Paultre, C.; Disse, E.; Rabasa-Lhoret, R.; et al. The microRNA signature in response to insulin reveals its implication in the transcriptional action of insulin in human skeletal muscle and the role of a sterol regulatory element-binding protein-1c/myocyte enhancer factor 2C pathway. Diabetes 2009, 58, 2555–2564.
618 619 620 621
(63)
Kato, M.; Wang, L.; Putta, S.; Wang, M.; Yuan, H.; Sun, G.; Lanting, L.; Todorov, I.; Rossi, J. J.; Natarajan, R. Post-transcriptional up-regulation of Tsc-22 by Ybx1, a target of miR-216a, mediates TGF-??-induced collagen expression in kidney cells. J. Biol. Chem. 2010, 285, 34004–34015.
622 623
(64)
Pandit, K. V.; Milosevic, J.; Kaminski, N. MicroRNAs in idiopathic pulmonary fibrosis. Transl. Res. 2011, 157, 191–199.
624 625 626 627
(65)
Miyakawa, H.; Imai, M.; Sugimoto, N.; Ishikawa, Y.; Ishikawa, A.; Ishigaki, H.; Okada, Y.; Miyazaki, S.; Koshikawa, S.; Cornette, R.; et al. Gene up-regulation in response to predator kairomones in the water flea, Daphnia pulex. BMC Dev. Biol. 2010, 10, 45.
628 629 630
(66)
Boucher, P.; Ditlecadet, D.; Dubé, C.; Dufresne, F. Unusual duplication of the insulin-like receptor in the crustacean Daphnia pulex. BMC Evol. Biol. 2010, 10, 305.
25
ACS Paragon Plus Environment
Environmental Science & Technology
631 632
(67)
Kallapur, S.; Ormsby, I.; Doetschman, T. Strain dependency of TGFB1 function during embryogenesis. Mol. Reprod. Dev. 1999, 52, 341–349.
633 634 635
(68)
McLennan, I. S.; Poussart, Y.; Koishi, K. Development of skeletal muscles in transforming growth factor-beta 1 (TGF-beta1) null-mutant mice. Dev. Dyn. 2000, 217, 250–256.
636 637
(69)
Hyun, S. Body size regulation and insulin-like growth factor signaling. Cell Mol. Life Sci. 2013, 70, 2351–2365.
638 639 640
(70)
Bertin, G.; Averbeck, D. Cadmium: cellular effects, modifications of biomolecules, modulation of DNA repair and genotoxic consequences (a review). Biochimie 2006, 88, 1549–1559.
641 642 643 644
(71)
Zhang, Z.; Sun, H.; Dai, H.; Walsh, R. M.; Imakura, M.; Schelter, J.; Burchard, J.; Dai, X.; Chang, A. N.; Diaz, R. L.; et al. MicroRNA miR-210 modulates cellular response to hypoxia through the MYC antagonist MNT. Cell Cycle 2009, 8, 2756–2768.
645 646 647 648
(72)
Ortiz-Barahona, A.; Villar, D.; Pescador, N.; Amigo, J.; del Peso, L. Genome-wide identification of hypoxia-inducible factor binding sites and target genes by a probabilistic model integrating transcription-profiling data and in silico binding site prediction. Nucleic Acids Res. 2010, 38, 2332–2345.
649
(73)
Altschul, S. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410.
650 651 652 653
(74) Kriventseva, E. V.; Tegenfeldt, F.; Petty, T. J.; Waterhouse, R. M.; Simao, F. A.; Pozdnyakov, I. A.; Ioannidis, P.; Zdobnov, E. M. OrthoDB v8: update of the hierarchical catalog of orthologs and the underlying free software. Nucleic Acids Res. 2015, 43, D250–D256.
654 655 656 657 658 659 660 661 662 663 664 665 666
Figure and Table Legends Figure 1. (A) Microarray identified significantly altered miRNAs (p-value < 0.05 and fold change > 1.5 compared with control group) in Daphnia pulex neonates exposed to low (20 µg/L CdCl2) and high (40 µg/L CdCl2) cadmium (Cd) for 48 h. (B) Heatmap of significantly expressed miRNAs. Figure 2. Validation of miRNA microarray results by q-PCR. Bars plotted represent means and standard errors with black bars representing microarray results and gray bars representing q-PCR results. Asterisks denote significant differences from controls (p