Transcriptomic and Functional Analyses on the Effects of Dioxin on

Sep 7, 2017 - In this study, transcriptomic and Ingenuity Pathway Analysis (IPA) underlined that an ex-vivo TCDD treatment (0.1 nM) stimulated insulin...
0 downloads 10 Views 2MB Size
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