Subscriber access provided by University of Winnipeg Library
Environmental Measurements Methods
BLOOD TRANSCRIPTOMICS ANALYSIS OF FISH EXPOSED TO PERFLUORO ALKYLS SUBSTANCES: ASSESSMENT OF A NON-DESTRUCTIVE SAMPLING TECHNIQUE FOR ADVANCING AQUATIC TOXICOLOGY RESEARCH. Ignacio Alejandro Rodriguez-Jorquera, Reyna Cristina Colli-Dula, Kevin J Kroll, B. Sumith Jayasinghe, Maria Virginia Parachu Marco, Cecilia Silva-Sanchez, Gurpal S. Toor, and Nancy D. Denslow Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b03603 • Publication Date (Web): 20 Dec 2018 Downloaded from http://pubs.acs.org on December 24, 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 37
Environmental Science & Technology
1 2 3
BLOOD TRANSCRIPTOMICS ANALYSIS OF FISH EXPOSED TO PERFLUORO ALKYLS SUBSTANCES: ASSESSMENT OF A NON-LETHAL SAMPLING TECHNIQUE FOR ADVANCING AQUATIC TOXICOLOGY RESEARCH.
4 5 6
Ignacio A. Rodríguez-Jorquera, 1,2,3* R. Cristina Colli-Dula, 3,4 Kevin Kroll, 3 B. Sumith Jayasinghe, 3 Maria V. Parachu Marco, 5,6,7 Cecilia Silva-Sanchez, 3 Gurpal S. Toor, 8 and Nancy D. Denslow 3
7 8
1 Centro
de Humedales Río Cruces (CEHUM), Universidad Austral de Chile, Independencia 641,Valdivia, Región de los Ríos, Chile.
9 10
2 Interdisciplinary
11 12
3 Department
13 14
4
15 16 17 18
5
19 20
6
21 22 23
7 Proyecto
24 25
8
26
*Corresponding Author:
[email protected] Ecology Program, School of Natural Resources and Environment, Soil and Water Science Department, University of Florida, Gainesville, FL, USA. of Physiological Sciences & Center for Environmental and Human Toxicology, University of Florida, Gainesville, FL, USA. Departamento de Recursos el Mar, Centro de Investigación y Estudios Avanzados del Instituto Politécnico Nacional, Unidad Mérida, México. Laboratorio de Biología Celular y Molecular Aplicada (LBCMA), Instituto de Ciencias Veterinarias del Litoral (ICiVet - Litoral), Universidad Nacional del Litoral (UNL) / Consejo Nacional de Investigaciones Científicas y Tecnológicas (CONICET), CP 3080, Esperanza, Santa Fe, Argentina. Facultad de Humanidades y Ciencias - Universidad Nacional del Litoral, Paraje El Pozo s/n, CP 3000, Santa Fe, Argentina. Yacaré. Laboratorio de Zoología Aplicada: anexo Vertebrados (Facultad de Humanidades y Ciencias - Universidad de Nacional del Litoral/ MASPyMA), Aristóbulo del Valle 8700, 3000 Santa Fe (Santa Fe) Argentina. Department of Environmental Science and Technology, University of Maryland, College Park, MD, USA.
ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 2 of 37
27
Abstract
28
In contrast to mammals, the blood from other vertebrates such as fish contains nucleated red
29
cells. Using a fathead minnow (Pimephales promelas) oligonucleotide microarray, we compared
30
altered transcripts in the liver and whole blood after exposure to environmentally relevant
31
concentrations of perfluorooctanesulfonic acid (PFOS) and a mixture of seven types of perfluoro
32
alkyl substances (PFAS) including perflurooctanoic acid (PFOA). We used qPCR and cell based
33
assays to confirm the main effects and found that blood responded with a greater number of
34
altered genes than the liver. The exposure to PFAS altered similar genes with central roles in a
35
cellular pathway in both tissues, including estrogen receptor alpha and peroxisome proliferator
36
activator beta, indicating that the genes previously associated with PFAS exposures are
37
differentially expressed in blood and liver. The altered transcripts are involved with cholesterol
38
metabolism and mitochondrial function. Our data confirmed that PFAS are weak xenoestrogens
39
and exert effects on DNA integrity. Gene expression profiling from blood samples not related
40
with the immune system including very low-density lipoprotein, vitellogenin, estrogen receptor,
41
and thyroid hormone receptor demonstrated that blood is a useful tissue to assess endocrine
42
disruption in non-mammalian vertebrates. We conclude that the use of blood as non-lethal
43
sampling in genomics studies is informative and particularly useful to assess effects of pollution
44
in endangered species. Lastly, using blood will reduce animal use and widen the experimental
45
design options to study effects of contaminants exposure on wildlife.
46 47
Keywords: Perfluoro Alkyl Substances, Non-lethal Sampling, DNA damage, Xenoestrogens,
48
Cholesterol Metabolism, Endocrine Disruption, Re-sampling.
ACS Paragon Plus Environment
2
Page 3 of 37
49
Environmental Science & Technology
TABLE OF CONTENTS (TOC) ART
50
51 52
ACS Paragon Plus Environment
3
Environmental Science & Technology
53 54
Page 4 of 37
INTRODUCTION Aquatic pollution affects biodiversity at local, regional, and global scales 1. In
55
response to this issue, several techniques have been developed to research the impact of
56
pollution on biota. Recently, genomic tools have improved our understanding of how
57
chemicals impact the environment, revolutionizing both the fields of toxicology and
58
ecotoxicology. 2 Specifically, transcriptomic analyses have been used to understand the
59
effects of aquatic pollution in biota.3–5 Transcriptomics profiling allows researchers to
60
perform analysis of thousands of genes at a time, to identify expression changes in response
61
to chemical exposures, which can help to better understand the mode of action associated
62
with biological responses to these exposures.2
63
Despite the major advances in techniques in environmental toxicology and
64
ecological risk assessment, improvements in data collection are necessary to understand the
65
effects of pollution at the population level. For instance, the duration and consistency of
66
pollutant effects, adaption and recovery periods collected in laboratory micro- and macro-
67
mesocosm experiments should be evaluated over time to ensure a realistic response.6 In this
68
context, the development of non-lethal techniques to assess pollution effects on the same
69
individual will help to advance the understanding and action of pollution. Therefore, the
70
use of blood emerges as an ideal tissue to assess changes in physiology in vertebrates
71
caused by external stressors, particularly the non-mammalian vertebrates that have
72
nucleated red cells and functional mitochondria.7
73
Advances in genomics approaches, especially for dealing with specific issues
74
affecting human subjects, allow the use of blood as a non-lethal sampling method to
75
analyze gene expression signatures for toxicity determination.8,9 Moreover, blood is a
ACS Paragon Plus Environment
4
Page 5 of 37
Environmental Science & Technology
76
logical choice for assessing immunotoxicity and is a much easier and cheaper tissue that
77
can be sampled repeatedly from the same individual. This approach allows a refinement of
78
sampling methodologies of experimental designs for toxicology/ecotoxicology and also
79
reduces the number of animals required since the same individuals can be used for multiple
80
comparisons (repeated measurements).10 In the case of genomics, the use of blood will
81
eliminate fish exposure history differences as well as public concern about animal welfare
82
and their use in experimental designs. The application of this approach could avoid the
83
necessity to euthanize animals to assess toxicant effects on target organs, which is a major
84
drawback of these applications with regards to monitoring wildlife in ecotoxicology
85
studies, particularly for endangered species.
86
Fully fluorinated alkyls, also known as per and polyfluoro alkyl substances (PFAS),
87
are considered one of the most ubiquitous types of endocrine disrupting compounds
88
(EDCs)11. Their wide use, high persistence and bioaccumulative properties have led to an
89
omnipresent occurrence in the environment. Particularly, the presence of PFAS in the
90
aquatic environment has been commonly documented.12,13 The accumulation of PFAS has
91
been widely detected in freshwater and marine fish14 with the liver and blood serving as the
92
main body compartments for the bioaccumulation of PFAS in fish including the fathead
93
minnow.15,16 Toxic effects of PFAS on biota have been studied for almost four decades in
94
animals17 and are largely documented in the scientific literature.18–20
95
Reported health-related changes in fish range from general effects such as reduced
96
growth, survival and reproductive organ impairment to alteration of specific cellular
97
pathways like lipid metabolism, oxidative stress, DNA damage, immunotoxicity, and
98
estrogenic effects.19 The induction of the peroxisome proliferator receptors (PPAR) seem to
ACS Paragon Plus Environment
5
Environmental Science & Technology
99
Page 6 of 37
have a pivotal role in triggering several deleterious effects in mammals such as DNA
100
damage, production of reactive oxygen species (ROS) and alteration of lipid metabolism.21
101
In some fish species, such as rainbow trout, exposure to PFAS increases the expression of
102
the estrogen receptor (ER), a response related with carcinogenesis, more so than alterations
103
of expression of PPAR isoforms. The effect of PFAS on the immune system of fish has also
104
been previously observed,22 including effects on the pro-inflammatory cytokines (TNF-α,
105
IL-1β, and IL-8), which have been highlighted as indicators of immune-toxicity produced
106
by perfluoroctane sulfonic acid (PFOS) exposure.23,24
107
Our hypothesis was that blood, an alternative tissue, would respond to low levels of
108
exposure to PFAS. This response would be different from liver (which is the primary
109
target-tissue for toxicological studies) but as informative to assess toxicity of the aquatic
110
environment using fish. We used a transcriptomics approach to investigate the differential
111
expression of gene responses of liver and blood in order to study whole blood as a non-
112
lethal sampling option.
113
MATERIALS AND METHODS
114
Exposure Water Preparation
115
Four different treatments were prepared for a static exposure of the fish (Table S1,
116
supporting information). A control (0 g/L); two PFOS (Perfluorooctanesulfonic acid,
117
CAS: 2795-39-3) concentrations: PFOS High (25 g/L), PFOS Low (0.5 g/L); and a
118
PFAS Mix, consisting of 7 types of PFAS at concentrations similar to those found
119
previously downstream from a wastewater treatment plant 25 to mimic the mixtures usually
120
found in waters with wastewater influence. The mixture treatment contained the following
ACS Paragon Plus Environment
6
Page 7 of 37
Environmental Science & Technology
121
seven types of PFAS: PFBA (Perfluorobutanoate, CAS: 375-22-4); PFHxA
122
(Perfluorohexanoic acid, CAS: 307-24-4); PFHpA (Perfluoroheptanoic acid, CAS: 375-85-
123
9); PFOA (Perfluorooctanoic acid, CAS: 335-67-1); PFOS (Perfluorooctanesulfonic acid,
124
CAS: 2795-39-3); PFNA (Perfluorononanoic acid, CAS: 375-95-1); PFDA
125
(Perfluorodecanoic acid, CAS: 335-76-2). We added 50 mL PFAS treatment preparations
126
into a pre-cleaned fiberglass distribution cylinder containing 38 L of carbon filtered and de-
127
chlorinated municipal water (City of Gainesville, Florida, USA) to give the desired final
128
concentrations. For the control group, Milli-Q water and a carrier Triethylene glycol (TEG)
129
were mixed in equivalent volumes for all treatments. All PFAS were purchased from
130
Wellington labs (Ontario, Canada) with purity over 99%. At the end of the exposures, 1 L
131
of water samples were collected from one tank of each treatment in polypropylene bottles
132
for analytical chemistry analysis of seven types of PFAS using EPA Method 537.26
133
Water Analyses
134
To corroborate the concentrations used in the exposures, the seven types of PFAS
135
were analyzed. The detection of PFAS was performed using a high performance liquid
136
chromatography system coupled with tandem mass spectrometry (HPLC/ESI-MS/MS)
137
allowing for a limit of detection (LOD) of 1.0 ng/L, and a limit of quantification (LOQ) of
138
4 ng/L, as previously reported by Rodriguez-Jorquera et al. 25 Quality control and standards
139
used were as previously reported.27,28 More details about water chemistry analysis are
140
included in the supplementary information section.
141
Fish Exposure and Tissue Collection
ACS Paragon Plus Environment
7
Environmental Science & Technology
142
Page 8 of 37
Reproductively mature pond-reared fathead minnow (FHM) were purchased from
143
Andersen Minnow Farm (Lonoke, Arkansas) approximately 6-months before the
144
exposures. Thirty-two males were separated from the common tank 2-weeks before the
145
experiment and placed in the treatment aquaria for 48 h. Each exposure was conducted in
146
quadruplicate and each aquarium contained two male FHM in 4 L of treatment water. All
147
exposures were static. Positions of the treatment tanks were randomized and test initiation
148
times were staggered to ensure an exposure/sampling interval of 48 h. The fish were not fed
149
during the experiment. The temperature range of the water was 24–26 °C and the
150
photoperiod of 16 h light: 8 h dark was used. At the conclusion of the exposures, fish were
151
anesthetized with MS-222 and weighed to the nearest 0.1 g. Average fish weights were
152
3.68 g (control); 3.19 g (PFOS Low 0.5 ug/L); 3.99 g (PFOS High 25 ug/L); and 3.47 g
153
(PFAS Mix). The testes were removed and preserved for histological analysis to confirm
154
sex and sexual maturity stage. Liver tissue and whole blood was flash frozen using liquid
155
nitrogen and stored at -80 °C until RNA extraction. Samples of liver and blood (a minimum
156
of 40 uL) from the same individuals (four for each treatment) were used for RNA
157
extraction and microarray processing (Figure 1). All procedures involving live fish were
158
reviewed and approved by the University of Florida Institutional Animal Care and Use
159
Committee (IACUC).
160
Microarray Analysis
161
Total RNA was extracted following the RNA STAT-60 reagent protocol (Tel-Test,
162
Friendswood, TX, USA). RNA was reconstituted in RNAsecure (Ambion; New York,
163
USA), and DNase-treated with Turbo DNA-free (Ambion; New York, USA). In order to
164
extract RNA from blood, column purification was used (RNeasy® Mini Kit Qiagen
ACS Paragon Plus Environment
8
Page 9 of 37
Environmental Science & Technology
165
®;Limburg, Germany). The RNA quantity for microarray analysis was measured using the
166
NanoDrop ND-1000 (Nanodrop Technologies, Wilmington, DE) and RNA quality was
167
evaluated using the Agilent 2100 BioAnalyzer with the RNA 6000 Nanochip. RNA
168
integrity numbers (RIN) ranges were 8.1 to 9.3 for blood samples and 7.5 to 9.4 for liver
169
samples.
170
An FHM 8 X 15K oligonucleotide microarray manufactured by Agilent (Palo Alto,
171
CA) and designed in our laboratory (GEO: GPL9248) was used in this study.29 The array
172
consisted of 15,744 oligonucleotide probes based on publicly available EST sequences at
173
the time plus genes that had been identified through suppressive subtractive hybridization
174
in the laboratory and also contained 536 control probes. Four biological replicates were
175
used for RNA isolation from FHM livers and blood using the same individuals to eliminate
176
exposure history as a variable when comparing the two gene expression profiles.
177
Microarray hybridizations were performed according to the Agilent "One-color
178
microarray”- (document no. G4140-90040 v6.5) using Cyanine 3 (Cy3) (Agilent, Palo Alto,
179
CA). For blood and liver samples, 25 ng and 50 ng of total RNA per sample respectively
180
were used to produce cDNA by reverse transcription using the poly (A) RNA present in the
181
starting total RNA sample. Each sample contained a specific activity >8.15 pmol Cy3/ug,
182
and amounts were adjusted to a final mass of 600 ng per sample for hybridization. A final
183
volume of 25 μL containing fragmented cRNA was added to the microarrays, and then
184
hybridization proceeded for 17 h at 65 ºC. Then, microarrays were washed according to the
185
Agilent protocol and kept in the dark until scanning on an Agilent G2505B microarray
186
scanner (same day). Data extraction was performed using Agilent Feature Extraction
187
software (v9.5).
ACS Paragon Plus Environment
9
Environmental Science & Technology
Page 10 of 37
188 189 190
Figure 1. General procedures of blood and liver gene expression analysis in fish using cDNA microarrays. In this study, the fish were euthanized.
ACS Paragon Plus Environment
10
Page 11 of 37
191
Environmental Science & Technology
Bioinformatics
192
Raw expression data (gProcessedSignal) were imported into JMP Genomics v5
193
(SAS, Cary, NC) and log 2-transformed and normalized by LOESS before performing
194
ANOVA to identify differentially regulated transcripts. Differentially regulated transcripts
195
p < 0.01; fold change greater than ± 1.5 were subjected to hierarchical clustering. Distance
196
calculations were performed with the program Cluster 3.0 30 using Euclidean distance as a
197
similarity metric and average linkage as a clustering method and visualized using the Java
198
Tree View software.31
199
To evaluate gene expression changes through pathway analysis, human orthologs
200
were found for FHM genes. Gene Ontology (GO) analysis was used to determine the
201
biological processes containing over represented genes using Fisher Exact Test (p