Subscriber access provided by GRIFFITH UNIVERSITY
Article
Transcriptomic and Functional Analyses on the Effects of Dioxin on Insulin Secretion of Pancreatic Islets and #-cells Keng Po Lai, Hin Ting Wan, Alice Hoi-Man Ng, Jing Woei Li, Ting-Fung Chan, and Chris Kong-Chu Wong Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b02830 • Publication Date (Web): 07 Sep 2017 Downloaded from http://pubs.acs.org on September 8, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 32
Environmental Science & Technology
1
2 3 4
6
Transcriptomic and Functional Analyses on the Effects of Dioxin on Insulin Secretion of Pancreatic Islets and β-cells.
7
Short title: Transcriptomic and Functional Analysis of TCDD-treated Islets.
5
8 Keng Po Lai2#, Hin Ting Wan1#, Alice Hoi-Man Ng1, Jing Woei Li 2,3, Ting Fung Chan3, Chris Kong-Chu Wong1*
9 10 11 12 13
1
Croucher Institute for Environmental Sciences, Partner State Key Laboratory of Environmental and Biological Analysis, Department of Biology, Hong Kong Baptist University, 2Department of Chemistry, City University of Hong Kong, 3 School of Life Sciences, Hong Kong Bioinformatics Centre, The Chinese University of Hong Kong, Hong Kong SAR, China.
14 15 16 17 18 19 20 21 22 23 24 25 26 27
#
Equal contribution *Corresponding author, Dr Chris KC Wong Croucher Institute for Environmental Sciences, Partner State Key Laboratory of Environmental and Biological Analysis, Department of Biology, Hong Kong Baptist University, Email address:
[email protected]. Phone: 852-34117053
Keywords: Akt-mTOR, AMPK, ATP, GSIS, pyruvate dehydrogenase kinase.
28 29
The authors declare they have no actual or potential competing financial interests.
30 31
1 ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 32
32
ABSTRACT
33
In this study, transcriptomic and Ingenuity Pathway Analysis (IPA) underlined that an ex-vivo
34
TCDD treatment (0.1 nM) stimulated insulin-release in mouse pancreatic islets via the effect on the
35
Akt-mTOR-p70S6K, AMPK and ERK1/2 pathways. Functional studies using both ex-vivo islets
36
and the mouse β-cell-line (Min-6) validated the stimulatory effects of TCDD (0.1 and 1 nM) on
37
basal-insulin secretion. At 0.1 nM TCDD treatment on Min-6, western blot analysis showed
38
activation of ERK1/2 and decreased expression of pyruvate dehydrogenase kinase (PDK). A
39
reduction of PDK expression is associated with an increase of pyruvate dehydrogenase flux. This
40
observation was supported by the detection of significantly higher cellular ATP levels, an
41
increase of glucose-stimulated-insulin-secretion (GSIS), and an inhibition of the AMPK pathway.
42
At 1 nM TCDD treatment on Min-6, significant inhibitions of the Akt-mTOR pathway, cellular ATP
43
production, and GSIS were evident. The experimental studies in Min-6 supported the IPA of
44
transcriptomic data in pancreatic islets. Collectively, TCDD treatment caused an elevated basal-
45
insulin release in both islets and β-cell cultures. Moreover, our data revealed that the modulation of
46
the Akt-mTOR-p70S6K, AMPK and ERK1/2 pathways might be an important component of the
47
mechanism for the TCDD-perturbing effects on ATP production in β-cells in affecting insulin
48
secretion.
49
2 ACS Paragon Plus Environment
Page 3 of 32
Environmental Science & Technology
50
INTRODUCTION
51
Type-2 diabetes (T2D) is a pandemic and constitutes 90% of people with diabetes. Genetic and
52
nutritional factors are known to be major risk factors for this disease. However, these factors
53
alone are still not able to account for the high prevalence of metabolic disease worldwide in
54
recent years1,2. Therefore, additional efforts are required to reveal the pathogenesis of the disease.
55
Recently, scientific researches have focused on characterizing other potential risk factors that
56
may disrupt body energy homeostasis. The role of environmental chemical contaminants has
57
drawn significant attention as potential contributors, independent of diet and physical activity3.
58
Possible roles of chemical contaminants in metabolic diseases have been suggested, based on the
59
observation of positive correlations between the increase in incidence of metabolic disorders and
60
the use of industrial chemicals in the past 50 years4. Moreover, epidemiological studies in
61
chemical-spill incidents and laboratory exposure studies hint at the possible link between
62
pollutant exposure to disturbed glucose homeostasis and insulin resistance5.
63
The well-known anthropogenic pollutant dioxin (i.e. 2,3,7,8-tetrachlorodibenzo-p-dioxin,
64
TCDD) is proposed to be a putative diabetogenic factor6, based on the analysis of numerous
65
human epidemiological and laboratory animal data. During the Vietnam War in the early 60s,
66
dioxin exposure (via Agent Orange) was reported to be associated with an increase in the
67
prevalence of diabetes mellitus, with 111,726 identified cases in Korean Vietnam veterans7.
68
From the residential exposure to dioxins at the chemical explosion in Seveso in 1976, an inverse
69
association between maternal serum dioxins with birth weight was recognized in 1,211 cases of
70
post-explosion pregnancies8. A connection between size at birth to the risk of developing insulin
71
resistance and cardiovascular diseases is known9. Moreover, a positive association between
72
serum dioxin concentrations and metabolic syndrome was found in women who were ≤ 12 years
3 ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 32
73
of age at the time of exposure in Seveso10. An increased risk for diabetes was reported in other
74
incidences of dioxin exposure, like in phenoxy-herbicide and chlorophenol production workers11,
75
and cohorts in the Yucheng accident12. Recent studies showed positive associations between
76
serum dioxin concentrations and insulin resistance13,14, hypertension and glucose intolerance15,16.
77
However, a prospective association analysis in elderly Swedish adults showed no direct
78
association between serum dioxin levels and T2D17. These observations suggest that the risk of
79
dioxins may be greater when exposed at younger age. Nevertheless, these studies do not establish
80
causality, because the cross-sectional study design makes it difficult to assess the exposure-
81
disease association. On the other hand, laboratory data provided some mechanistic data to
82
support the notion that dioxin exposure is associated with dysglycemia. An earlier study of
83
dioxin exposure in rats reported a pathological outcome of wasting syndrome, where diminishing
84
availability of glucose for metabolism was noted18,19. The hallmark of wasting syndrome, a
85
decrease of hepatic gluconeogenesis, was found to be mediated by the activation of the cytosolic
86
receptor, aryl hydrocarbon receptor (AhR), which stimulates the expression of TCDD-inducible
87
poly(ADP-ribose) polymerase to suppress hepatic glucose synthesis20,21. In addition, inhibitory
88
effects of dioxins on glucose uptake in the brain, adipose and pancreas were demonstrated in
89
guinea pigs20,22, via modulation of glucose transporting activities23.
90
In fact, establishing an association between TCDD exposure and the development of diabetes is
91
inherently difficult. A comparison of the gene expression pattern of pancreatic islets between
92
normal and T2D people revealed a significant reduction in the levels of AhR nuclear translocator
93
(ARNT) in patients with the disease24. ARNT is a member of the bHLH/PAS family of
94
transcriptional factors, and is a heterodimer partner of AhR or hypoxia-inducible factor-1α25.
95
Depletion of ARNT expression in pancreatic β-cells of ARNT-knockout mice caused an 4 ACS Paragon Plus Environment
Page 5 of 32
Environmental Science & Technology
96
impairment of glucose-stimulated insulin release (GSIS), with no significant effects on insulin
97
content and the total number of isolated islets24. The actions of ARNT were found to regulate
98
changes in islet gene expression. Since the dysfunction of pancreatic β-cells is known to be a key
99
component in the pathogenesis of T2D, and the AhR/ARNT heterodimer regulates many
100
pancreatic β-cell genes24, it warrants further analysis of the underlying mechanistic effects of
101
dioxin on pancreatic islet functions.
102 103
5 ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 32
104
MATERIALS AND METHODS
105
Experimental Animals and Chemicals
106
CD-1 mice were kept and maintained according to the “Guidelines and Regulations” of the
107
Department of Health, Hong Kong Special Administrative Region. The animal handling
108
procedure was approved by the Animal Committee of the Hong Kong Baptist University (Permit
109
no. 261812). Male CD-1 mice (6-8 weeks old) were received from the “Animal Unit” of the
110
University of Hong Kong. The animals, housed in polypropylene cages with sterilized bedding,
111
were kept under controlled temperature (22 °C) and 12-h/12-h light-dark cycle. TCDD was
112
purchased from Sigma Aldrich (CA, USA). A stock solution of TCDD (1 µM in dimethyl
113
sulfoxide, DMSO) was made. Further dilutions were prepared in DMSO and phenol red-free
114
DMEM-high glucose medium with charcoal-stripped serum. The final DMSO concentration was
115
less than 0.1%. The corresponding percentage of DMSO was added as the vehicle control.
116 117
Pancreatic islets isolation
118
We followed the method of pancreatic islet isolation described previously by Nunemaker and
119
coworkers26. Briefly, male mice were euthanized by cervical dislocation. For each mouse, the
120
animal was dissected and the common bile duct attached to the duodenum at the papilla was
121
identified. A 30-gauge 0.5-inch needle secured to a 5 ml syringe was cannulated into the bile
122
duct and tied with a suture. A 3 ml aliquot of collagenase solution [1.4 mg/ml Collagenase-P
123
(Roche) in HBSS/1% BSA] was injected slowly to completely inflate the pancreas. The inflated
124
pancreas was then dissected and placed into a 50 ml conical tube containing 3 ml of collagenase
125
solution, and then incubated at 37 oC for 11 min, without agitation. The tissue suspension was
126
then shaken vigorously by hand for approximately 5 s. The enzymatic digestion was then slowed
6 ACS Paragon Plus Environment
Page 7 of 32
Environmental Science & Technology
127
down by adding 15 ml of ice-cold HBSS/1% BSA. The suspension was centrifuged for 2 min at
128
1200 rpm. The pellet was washed with 10 ml of HBSS/1% BSA twice and the final suspension
129
was filtered through a 420 µm sized sieve into a new 50 ml conical tube. The filtered suspension
130
was filled up to 20 ml with HBSS/1% BSA, followed by centrifugation for 2 min at 1200 rpm.
131
The pellet was resuspended in 10 ml of Histopaque 1100 solution (Sigma), and then centrifuged
132
for 20 min at 1200 rpm. The supernatant was collected and diluted with 25 ml of HBSS/1% BSA
133
in a new 50 ml conical tube, and then centrifuged for 4 min at 1500 rpm. The pellet was washed
134
in 10 ml of HBSS/1% BSA and finally resuspended in RPMI 1640 (10% FBS, 1% Pen/Strep).
135
After overnight incubation, islets from the preparation of the individual pancreas were picked
136
under a light microscope. Islets from two mice were pooled as one biological sample (~150 islets
137
per replicate), and three replicates (6 animals per group) were used for the control and TCDD
138
treatment. The islets were then treated with either DMSO (solvent control) or 0.1 nM TCDD for
139
24 h.
140 141
Transcriptomic Analysis of Ex-vivo Islets
142
Total RNA was isolated from the control and TCDD-treated islets using the mirVanaTM RNA
143
isolation kit (Applied Biosystems). The RNA quality was determined (Agilent 2100 Bioanalyzer
144
system), RNA with a RNA Integrity Number > 8 was used for library construction and
145
sequencing (pair-end reads, 150 bp read-length) using the Illumina MiSeq sequencer. Sequence-
146
reads were trimmed according to BWA's – q algorithm as described previously27. Quality-
147
trimmed sequence reads were mapped onto the mouse genome reference GRCm38 using STAR
148
v2.5.2b. Read-count data were used for differential expression analysis using the edgeR
149
package28. Genes with a FDR < 0.05 were considered as differentially expressed genes (DEGs).
7 ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 32
150
The DEGs were analysed using the Database for Annotation, Visualization and Integrated
151
Discovery (DAVID) v6.829,30 for functional annotation and to determine the biological processes.
152
Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was conducted to identify the
153
alteration of pathways caused by TCDD treatment. Moreover, Ingenuity Pathway Analysis
154
(IPA®, QIAGEN Redwood City, www.qiagen.com/ingenuity) was applied to underline the
155
molecular interaction networks. Canonical pathways with p < 0.05 were set as statistically
156
significant.
157 158
Real-time PCR analysis
159
Total RNA with an A260/A280 > 1.8 were used for quantitative real-time PCR analyses. For each
160
sample, complementary DNA (cDNA) was prepared from 150 ng of total RNA using the VILO
161
cDNA synthesis Kit (Invitrogen, ThermoFisher). Gene-specific primers were designed according
162
to the published sequences (Supplementary Table 1). Real-time PCR was implemented with a
163
program of 10 min at 95 oC, followed by 40 cycles of 98 °C for 20 s, 60 °C for 1 s and 72 °C for
164
60 s. The copy numbers of standards and cDNAs of the samples were measured using the
165
StepOne Real-time PCR system with the KAPA HiFi Hotstart ReadyMix PCR kit (KAPA
166
Biosystems). The data were normalized with the transcript levels of the housekeeping gene,
167
mouse actin. Control amplifications were performed either without the step of reverse
168
transcription or RNA.
169 170
Cell culture of mouse pancreatic β-cell-line
171
Min-6 cells were maintained in a humidified atmosphere with 5% CO2 in a complete DMEM-
172
high glucose medium (Gibco, ThermoFisher) supplemented with 10% heat-inactivated fetal
8 ACS Paragon Plus Environment
Page 9 of 32
Environmental Science & Technology
173
bovine serum, 50 µM 2-mercaptoethanol, 100 U/ml penicillin and 100 µg/ml streptomycin
174
(Invitrogen). Cell cytotoxicity was assessed by the MTT assay [3-(4,5-Dimethylthiazol-2-yl)-2,5-
175
Diphenyltetrazolium Bromide, Invitrogen] following the manufacturer’s instructions. In brief,
176
cells were incubated in MTT solution at 37 °C for 4 h, followed by the addition of 100 µl of
177
DMSO and mixed thoroughly. After incubating at 37 °C for 10 min, the absorbance at 540 nm
178
was read using a microplate reader.
179 180
Basal and Glucose-stimulated Insulin Secretion
181
Min-6 cells were seeded at a density of 5 x 104 cells/cm2 in a 24-well plate for 48 h before drug
182
treatment. Isolated pancreatic islets were evenly distributed in a 24-well-plate (60 islets per well).
183
Each culture was treated with either DMSO, 0.1 nM or 1 nM TCDD for 24 h. The Min-6 cells
184
and islets were then washed twice with Krebs-Ringer bicarbonate buffer (KRBB, containing 2.4
185
mmol/L CaCl2, 5.4 mmol/L KCl, 1.2 mmol/L KH2PO4, 1.2 mmol/L MgSO4, 120 mmol/L NaCl,
186
20 mmol/L HEPES, 0.1% BSA, pH 7.0) to remove leftover glucose in the medium. Cells were
187
then pre-incubated in KRBB containing 1.67 mmol/L glucose for 1 h in the incubator, and
188
afterwards incubated with 1.67 mmol/L, followed by 16.7 mmol/L glucose, 1 h for each. The
189
supernatants were collected as basal and stimulated insulin secretions respectively. The levels of
190
insulin were measured by using the Insulin AlphaLISA detection kit and the fluorescent signals
191
were measured by the EnSight Multimode plate reader (PerkinElmer). Data were normalized by
192
measuring protein concentrations using the Bio-Rad DC protein assay.
193 194
ATP Determination
9 ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 32
195
Min-6 cells were seeded at a density of 5 x 104 cells/cm2 in a 24-well plate for 48 h before drug
196
treatment. After 24 h of TCDD (0.1 and 1 nM) incubation, the cells were lysed in 1X passive
197
lysis buffer (Promega) followed by centrifugation at 10,000 x g for 2 min. Aliquots of 10 µl of
198
samples were subjected to ATP determination using the ATP determination kit (Invitrogen) and
199
read by the VICTOR X4 multilabel plate reader (PerkinElmer). Data were normalized by
200
measuring the protein concentrations.
201 202
Western Blotting and Immunofluorescence Analysis (IF)
203
Min-6 cells were seeded at a density of 2 x 105 cells/cm2 in a 6-well plate for 48 h before TCDD
204
treatment. The protein samples were lysed in the radioimmunoprecipitation assay buffer (RIPA:
205
150 mM sodium chloride, 1% NP-40, 0.1% SDS, 0.5% sodium deoxycholate, and 50 mM Tris-
206
Cl, pH 8.0) mixed with a cocktail of protease inhibitors (Sigma-Aldrich) at 1:100 dilution.
207
Concentrations of protein were measured and 15-20 µg of samples were subjected to
208
electrophoresis in 7.5% or 10% polyacrylamide gels. The gels were then blotted onto a PVDF
209
membrane (Bio-Rad) and blocked in 5% non-fat dry milk (wt/vol) in PBS-0.1% Tween-20 buffer.
210
Primary and secondary antibodies (Supplementary Table 2) were prepared using 0.1% BSA
211
(wt/vol) or 3% non-fat dry milk in PBS-Tris buffer (10 mM Tris, pH 7.4). The protein blots
212
were incubated with the primary antibody overnight at 4oC, then secondary antibody for 1 h at
213
room temperature. The blots were incubated with WESTSAVE Up (AbFrontier) and protein
214
bands were analyzed using ImageJ software.
215
For immunofluorescent staining, the cells were seeded at a density of 1 x 105 cells/cm2 on a
216
cover glass in each well of a 12-well plate for 48 h before TCDD treatments. After 24 h of the
217
treatment, cells were fixed with 4% paraformaldehyde in PBS (pH 7.4) at room temperature for
10 ACS Paragon Plus Environment
Page 11 of 32
Environmental Science & Technology
218
10 min. The cells were permeabilized using 0.1% Triton X-100 in PBS for 10 min. The cells
219
were then incubated with 1% (wt/vol) BSA in a PBS-Tris buffer, followed by the primary
220
antibody of guinea pig anti-insulin (DAKO) or rabbit anti-pdx-1 (Millipore) overnight, and a 1 h
221
incubation with Alexa-Fluor-conjugated secondary antibodies (Invitrogen, Thermofisher). Nuclei
222
were visualized with Vectashield antifade mounting medium with DAPI (4’,6’-diamidino-2-
223
phenylindole, Vector laboratories). Images in TIFF format were examined using Adobe
224
PhotoShop in Creative Suite 6. The imaging data were the representative results from at least
225
four independent experiments.
226 227
Statistical analysis
228
Statistical evaluations were implemented using SigmaStat version 3.5. Data were analyzed by the
229
Student’s t-test or one-way analysis of variance (ANOVA) followed by Duncan’s multiple range
230
test. All data are presented as statistical mean ± SD. A p value < 0.05 was used as the cutoff for
231
statistical significance.
232 233 234
11 ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 32
235
RESULTS AND DISCUSSION
236
Comparative transcriptomic analysis revealed the activation of the insulin signalling pathway
237
upon TCDD treatment.
238
In an attempt to understand the molecular mechanisms underlying the effects of TCDD on pancreatic
239
islets, comparative transcriptomic analysis on 0.1 nM TCDD-treated and DMSO (solvent control)-
240
treated islets were conducted. Quality-trimmed Illumina reads of 20.65 M and 23.14 M were
241
acquired from the TCDD- and DMSO -treated islets, respectively. A total of 3.3 Gb of clean bases
242
was retrieved from the RNA sequencing (Supplementary Table S3). A total of 5484 differentially
243
expressed genes (DEGs), including 5179 up- and 305 down-regulated genes, were identified from
244
comparing the data between TCDD- and DMSO-treated islets (Supplementary Table S4). The DEGs
245
were subjected to the Database for Annotation, Visualization and Integrated Discovery (DAVID)
246
v6.7 and Ingenuity Pathway Analysis (IPA) to identify the altered gene networks and canonical
247
pathways. The result of KEGG from DAVID analysis showed that TCDD treatment perturbed
248
molecular and biochemical pathways that affected islet functions (Fig. 1A and Supplementary Table
249
S5). The pathways directly related to the islet function and insulin secretion include the (i)
250
phosphatidylinositol signaling system, (ii) calcium signaling pathway, (iii) mammalian Target of
251
Rapamycin (mTOR) signaling pathway, (iv) Wnt signaling pathway, (v) insulin signaling pathway
252
and (vi) Type II diabetes mellitus. Analysis of the Ingenuity Canonical Pathway (IPA) revealed that
253
the TCDD treatment significantly modulated 275 pathways (p < 0.05), in which 154 pathways were
254
stimulated and 13 pathways were inhibited. The canonical pathways directly associated with islet
255
function and insulin secretion are (i) PI3K/AKT, (ii) mTOR signaling, (iii) p70S6K signaling, (iv)
256
AMP-activated protein kinase (AMPK) signaling, (v) calcium signaling, (vi) insulin receptor
257
signaling, (vii) Type I diabetes mellitus signaling, and (viii) Type II diabetes mellitus signaling
12 ACS Paragon Plus Environment
Page 13 of 32
Environmental Science & Technology
258
(Supplementary Table S6). The canonical pathways associated with TCDD binding include (i) aryl
259
hydrocarbon receptor signaling and (ii) estrogen receptor signaling. Based on the pathway analyses,
260
several candidate genes were selected for validation by real-time PCR assays using an independent
261
cohort of isolated pancreatic islets (Fig. 1B). A highlighted IPA canonical pathway – insulin receptor
262
signaling was schemed (Fig. 1C). Perturbation of the Akt-mTOR-p70S6K and AMPK pathways
263
was observed, which predicted a stimulation of insulin release upon exposure to 0.1 nM TCDD.
264
Intriguingly, the mTOR and AMPK-pathways are known to be important in nutrient and energy
265
sensing in pancreatic β-cells31. This observation prompted us to conduct experiments using ex vivo
266
pancreatic islet culture and the mouse islet β-cells, Min-6, to elucidate the effects of TCDD on insulin
267
secretion.
268 269
TCDD stimulated the basal insulin secretion
270
To verify the predicted effects of TCDD on islet functions, mouse pancreatic islets were isolated to
271
measure insulin release upon TCDD treatment (Fig. 2A). Ex-vivo islet culture incubated with 0.1 nM
272
or 1 nM TCDD for 24 h showed increases of basal insulin secretion (Fig. 2B). No significant
273
differences were observed in the glucose-stimulated-insulin-release (GSIS) in both the control and
274
TCDD treated groups. However, the ratio of GSIS/basal insulin secretion was significantly lower in
275
the TCDD-treated islets (Fig. 2C). To validate the effects and elucidate the mechanistic action of
276
TCDD, the mouse β-cell line, Min-6 was used in the subsequent experiments. No noticeable effect
277
on cell viability was observed from increasing the concentrations of TCDD treatment (1 pM – 4 nM)
278
(Fig. 3A). Effects of low (1.67 mM) and high glucose (16.7 mM) levels on cellular ATP levels of the
279
control and TCDD treated cells were measured (Fig. 3B). TCDD treatments (0.1 and 1 nM) caused
280
significant increases of ATP levels in cells incubated at the low glucose level. Only 0.1 nM TCDD
13 ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 32
281
treatment increased ATP production in cells incubated at the high glucose medium. The increases of
282
cellular ATP levels were positively correlated with the levels of insulin release at both low and high
283
glucose incubation (Fig. 3C). At the low glucose level, the basal insulin release in TCDD treatment
284
was significantly elevated. At the high glucose incubation, the level of insulin release (GSIS) was
285
significantly induced in 0.1 nM TCDD treatment as compared with the control. Similar to the results
286
of the ex-vivo pancreatic culture, the ratios of GSIS/basal insulin secretion were significantly lower
287
in the TCDD-treated cells (Fig. 3D). We then tested if the TCDD treatment affected the expression
288
levels of AhR-target genes (i.e. CYP1A1) and the transcriptional factors for β-cell differentiation
289
and maturation [i.e. paired Box 4 (Pax4), pancreatic and duodenal homeobox 1 (Pdx1)]. A dose-
290
dependent induction of CYP1A1 protein was observed (Fig. 3E). These results implied that the
291
AhR-ARNT pathway responded to TCDD treatment in the β-cells. No significant changes in the
292
expression levels of Pax4 and Pdx1 were observed among the control and treatment groups (Fig.
293
3E). No noticeable changes in the cellular levels of Pdx1 and insulin were observed in the
294
immunofluorescent staining (Fig. 4). Taken together, both ex vivo and cell-line culture models
295
showed elevations of basal insulin secretion upon TCDD treatment. The functional analysis
296
supported the results of the transcriptomic and bioinformatics analyses.
297 298
TCDD altered cell signaling pathways in β-cells
299
With reference to the canonical pathways identified in the transcriptomic data, western blot
300
analyses of the AKT-mTOR-p70S6K and AMPK signaling pathways were implemented to
301
elucidate the modulations of these signaling pathways upon TCDD treatment in Min-6. At 0.1 nM
302
TCDD treatment of Min-6, the western blot analysis showed activation of ERK1/2 but decreased
303
levels of pyruvate dehydrogenase kinase (PDK) (Fig. 5A). Functionally, the stimulation of
14 ACS Paragon Plus Environment
Page 15 of 32
Environmental Science & Technology
304
ERK1/2 pathway could positively regulate pyruvate dehydrogenase (PDH) flux via the reduction
305
of PDK expression32, resulting in an increase of mitochondria-derived ATP. Consistently,
306
significant higher cellular ATP levels were detected in the cells treated with 0.1 nM TCDD at
307
both low and high glucose conditions. In addition, the TCDD treatment significantly reduced the
308
levels of pAMPKα (Fig. 5A). Since AMPK is a key regulator of cellular energy status in
309
response to changes in cellular ATP levels33, these results supported the data of greater cellular
310
ATP contents in the 0.1 nM TCDD treated cells. The treatment however showed no noticeable
311
effect on the phosphorylation of the protein kinases Akt, mTOR and p70S6K as compared with the
312
control. In fact, multiple factors control basal insulin secretion and GSIS. The conventional model
313
establishes that an increase of cellular ATP leads to the closure of ATP-sensitive K+-channels,
314
followed by the opening of voltage Ca2+-channels, which triggers insulin exocytosis34. Intriguingly, a
315
previous study using the insulin-secreting β-cell line (INS-1E) demonstrated that TCDD treatment
316
at high doses (50-100 nM) could increase the cellular Ca2+-levels to trigger insulin exocytosis,
317
independent of glucose stimulation35. Taken together, our data and the previous study supported
318
our observation that 0.1 nM TCDD treatment elicited greater levels of basal insulin release and GSIS
319
in the culture.
320
On the other hand, inhibition of the Akt/mTOR signaling pathway at 1 nM TCDD treatment was
321
observed. The treatment showed significant reductions in the levels of p-Akt (ser 773) and p-mTOR
322
(ser 2448) (Fig. 5B). No significant changes in the levels of pAMPKα, pERK1/2 and PDK were
323
found. The inhibition of the Akt-mTOR pathway supported the observation of significantly lower
324
cellular ATP levels and GSIS in the cells at the high glucose condition, as compared to the 0.1 nM
325
TCDD treatment. Even though the 1 nM TCDD treatment did not cause cytotoxic effects on the
326
cells as demonstrated by the MTT assay, suppression on the mTOR (the sensor for cellular
15 ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 32
327
nutrient levels)36 for the coordination of energy metabolism was evident. Indeed, another study
328
showed the cytotoxic effects of high concentrations of TCDD treatment (12.5 and 25 nM) on
329
INS-1E, mediated by increased rates of Ca2+ influx and mitochondria depolarization37. While the
330
treatment of INS-1E with 1 nM TCDD reduced GSIS. Similar findings on TCDD-impaired GSIS
331
were also reported in other animal studies. In TCDD-treated rats (1 µg/kg b.w. i.p.), a decrease
332
of glucose uptake in isolated pancreatic islets was found to be associated with an inhibitory
333
effect on GSIS38. Moreover, an impairment of GSIS was shown using TCDD pre-treated mice
334
(10 µg/kg b.w. i.p.) and the mechanistic action was found to be AhR-dependent39. Although our
335
data of using 1 nM TCDD treatment on Min-6 did not demonstrate an inhibition of GSIS, treatment
336
at this concentration showed a significant decrease of GSIS when compared with the cells incubated
337
at 0.1 nM TCDD. The dissimilar observation in GSIS at the two different doses of TCDD treatments
338
might be due to the differential effects on the energy and nutrient sensors, AMPK-signaling (0.1 nM
339
TCDD) and AKT-mTOR-signaling (1 nM TCDD), resulting in the perturbation of cellular ATP
340
levels and insulin secretion.
341
Collectively, the present study characterized the transcriptional and functional effects of TCDD on
342
mouse pancreatic islets and the β-cells (Min-6). The transcriptome and IPA of the pancreatic islets
343
underlined the effects of TCDD on the perturbations of the Akt-mTOR-p70S6K, AMPK and
344
ERK1/2 pathways, predicting the outcome of stimulating insulin release. Subsequent experimental
345
studies confirmed the stimulatory effects of TCDD on basal insulin secretion in the two culture
346
systems. With the benefit of hindsight and the data of the present study, the stimulation of insulin
347
release was suggested to be mediated by an elevation of cytosolic Ca2+ and ATP levels. Moreover,
348
western blot analysis of the mouse β-cell-line supported the data of IPA and identified the effects of
349
TCDD on AMPK and mTOR signaling, the energy and nutrient sensors in the β-cells31. The
16 ACS Paragon Plus Environment
Page 17 of 32
Environmental Science & Technology
350
modulation of these signaling pathways might be an important component of the mechanism for
351
TCDD-perturbing effects on ATP production in β-cells in affecting basal and glucose-stimulated
352
insulin release.
353 354
ACKNOWLEDGEMENTS
355
This work was supported by the Partner State Key Laboratory of Environmental and Biological
356
Analysis (SKLP-16-17-P01, HKBU) and Department of Biology (Strategy Development Fund,
357
40-49-038, HKBU) to Dr Chris KC Wong (Hong Kong Baptist University).
358 359 360 361 362 363 364
17 ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 32
365
FIGURE LEGENDS
366
Figure 1. Comparative transcriptomic analysis of mouse pancreatic islets upon TCDD
367
treatment. Pancreatic islets were prepared separately from individual mouse. Islets from two
368
mice were pooled as one biological sample (~150 islets per replicate), and three replicates (6
369
animals per group) were used for the control and TCDD treatment. (A) KEGG pathway analysis
370
using DAVID on TCDD-dysregulated genes. (B) Validation of transcriptomic results on
371
independent cohorts of islet-culture using qPCR analysis (experimental replicates n = 4). (C)
372
Canonical pathway analysis of IPA shows the activation of insulin receptor signaling. Grey shade
373
represents overexpression; dot represents suppression. The intensity of color shade or dot represents
374
the degree of overexpression or suppression, respectively.
375 376
Figure 2. Effects of TCDD exposure for 24 h on ex-vivo pancreatic islets. (A) A
377
photomicrograph of isolated mouse islets. Effects of TCDD on (B) the levels of the basal insulin and
378
glucose-stimulated insulin release (n = 4, experimental replicates). The levels of basal insulin
379
secretion were significant greater in the TCDD-treated islets as compared with the control group (p