The evolutionary exploitation of vertebrate peroxisome proliferator

2 days ago - A striking example highlights the role of a group of compounds known as obesogens. In mammals, most examples involve the modulation of ...
1 downloads 0 Views 2MB Size
Ecotoxicology and Human is published by the Environmental Health

American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Subscriber access provided by UNIV OF Chemical Society. LOUISIANA However, no copyright

The evolutionary exploitation of vertebrate peroxisome

proliferatoractivatedisreceptor published by the American Chemical Society. 1155 Sixteenth # by organotins Street N.W., Washington,

Ana Capitão, Mónica DC 20036 LopesPublished by American Marques, Yoichiro Ishii, Chemical Society. © American Raquel Ruivo, Copyright Elza Fonseca,

Subscriber access provided by UNIV OF Chemical Society. LOUISIANA However, no copyright

Inês Páscoa, Rodolfo Jorge, Mélanie Barbosa, Youhei is published by the Hiromori, Takayuki American Miyagi, Chemical Society. 1155Miguel Sixteenth Tsuyoshi NAKANISHI, Street N.W., Washington, Santos, and L. Filipe DC 20036C. Castro Published by American

Environ. Sci. Technol., Just Chemical Society. Accepted Manuscript DOI: Copyright ©• American

Subscriber access provided by UNIV OF Chemical Society. LOUISIANA However, no copyright

10.1021/acs.est.8b04399 • Publication Date (Web): 06 Nov 2018 published by the Downloadedisfrom http:// American Chemical pubs.acs.org on November 2018 Society. 1155 7, Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Just Subscriber Accepted access provided by UNIV OF Chemical Society. LOUISIANA However, no copyright

“Just Accepted” manuscripts have been pee online prior to technical editing, is published by the formatting fo American Chemical as a service Society provides “Just Accepted” Society. 1155 Sixteenth of scientific material as soon as possible a Street N.W., Washington, full in PDF formatDCaccompanied by an HTM 20036 Published by American peer reviewed, but should not be considere Chemical Society. Digital Object Identifier Copyright(DOI®). © American“Just Accep Subscriber access provided by UNIV OFnot inclu Chemical the “Just Accepted” WebSociety. site may LOUISIANA However, no copyright

a manuscript is technically edited and form site and published as an ASAP article. No by the graphics wh to the manuscriptis published text and/or American Chemical ethical guidelinesSociety. that 1155 apply to the journal Sixteenth Street N.W., consequences arising fromWashington, the use of infor DC 20036 Published by American Chemical Society. Copyright © American Subscriber access provided by UNIV OF Chemical Society. LOUISIANA However, no copyright

TBT

Page 1γ ofEnvironmental 25 Science & Technology PPAR

ACS Paragon Plus Environment TBT agonist affinity towards PPARg : gain; loss.

Environmental Science & Technology

Page 2 of 25

1

The evolutionary exploitation of vertebrate peroxisome proliferator-activated receptor γ by

2

organotins

3 4

Ana M. F. Capitão1,2●, Mónica S. Lopes-Marques1●, Yoichiro Ishii3, Raquel Ruivo1, Elza S. S.

5

Fonseca1,2, Inês Páscoa1, Rodolfo P. Jorge1, Mélanie A. G. Barbosa1,2, Youhei Hiromori3,4,

6

Takayuki Miyagi3, Tsuyoshi Nakanishi3, Miguel M. Santos1,2* and L. Filipe C. Castro1,2*

7 8 9

1Interdisciplinary

Centre of Marine and Environmental Research (CIIMAR/CIMAR), University

10

of Porto (U.Porto), 4450-208 Matosinhos, Portugal

11

2Faculty

12

Porto, Portugal

13

3Laboratory

14

1-25-4 Daigaku-nishi, Gifu, 501-1196, Japan

15

4Faculty

16

Minamitamagaki, Suzuka, Mie 513-8670, Japan.

of Sciences (FCUP), Department of Biology, University of Porto (U.Porto), 4169-007

of Hygienic Chemistry and Molecular Toxicology, Gifu Pharmaceutical University,

of Pharmaceutical Sciences, Suzuka University of Medical Science 3500-3,

17 18

●Equal

19

*Address

contribution correspondence to:

20

Luís Filipe Costa Castro and Miguel Santos

21

Interdisciplinary Centre for Marine and Environmental Research (CIIMAR)

22

Av. General Norton de Matos s/n, 4450-208 Matosinhos, Portugal

23

Tel.: +351223401800

24

Fax: +351223390608

25

E-mail: [email protected]; [email protected]

26

Keywords: PPAR; chordates; evolutionary toxicology; obesogens

27 28 29

ACS Paragon Plus Environment

1

Page 3 of 25

Environmental Science & Technology

30

Abstract

31

Globally persistent man-made chemicals display ever-growing ecosystemic consequences, a

32

hallmark of the Anthropocene. In this context, the assessment of how lineage-specific gene

33

repertoires influence organism sensitivity towards endocrine disruptors is a central question in

34

toxicology. A striking example highlights the role of a group of compounds known as obesogens.

35

In mammals, most examples involve the modulation of the nuclear receptor, peroxisome

36

proliferator-activated receptor  (PPAR). To address the structural and biological determinants

37

of PPARγ exploitation by a model obesogen, tributyltin (TBT), in chordates, we employed

38

comparative genomics, transactivation and ligand binding assays, homology modeling and site-

39

directed-mutagenesis. We show that the emergence of multiple PPARs (,  and ) in vertebrate

40

ancestry coincides with the acquisition of TBT agonist affinity, as can be deduced from the

41

conserved transactivation and binding affinity of the chondrichthyan and mammalian PPAR. The

42

amphioxus single copy PPAR is irresponsive to TBT; as well as the investigated teleosts, a

43

probable consequence of a specific mutational remodeling of the ligand binding pocket. Our

44

findings endorse the modulatory ability of man-made chemicals and suggest an evolutionary

45

diverse setting, with impacts for environmental risk assessment.

ACS Paragon Plus Environment

2

Environmental Science & Technology

Page 4 of 25

47

Introduction

48

Over the past decades, various examples from disparate phylogenetic lineages have clearly

49

established the central role of nuclear receptor (NR) exploitation by xenobiotics in the generation

50

of endocrine imbalances. A critical question is to decipher how variations in NR gene function

51

and repertoire underscore variable endocrine responses to xenobiotic exposure1. One of these

52

examples has hypothesized that the exposure to an ample group of chemicals (termed obesogens)

53

parallels the ongoing obesity epidemic

54

extensive and includes a wide variety of compounds including phthalates, bisphenol A and

55

organotins 5. Importantly, a string of studies using mice has documented that the obesity

56

phenotype induced by obesogens such as tributyltin (TBT), results from the activation of a master

57

gene regulator of adipogenesis: the peroxisome proliferator–activated receptor γ (PPARγ)

58

In addition to mammals, this class of endocrine disruptors has also been shown to induce lipid

59

alterations in other animal phyla 7-13, with several lipid metabolic modifications described upon

60

obesogen exposure in teleosts and amphibians 7, 9, 10, 14-16. Yet, while the toxico-metabolic impact

61

has been described in various taxa, is it mostly unknown whether these effects underpin conserved

62

(or divergent) mechanisms of action, exploited by compounds such as TBT and acting via PPARs

63

8, 17, 18.

64

partial)

65

revealed that TBT is a potent repressor of PPARα and β but inactive on PPARγ 22. In effect,

66

animals exhibit conserved and divergent sensitivities to external stimuli yet, the role of PPARγ as

67

a direct target of obesogens has been addressed in a minute number of species 3, 7, 23, 24. Thus, the

68

extent to which genetic variation underlies such dissimilar responses is far from fully assessed.

69

Therefore, we have carried out an extensive examination of PPAR gene diversity in the main

70

chordate lineages, including cephalochordates (amphioxus), chondrichthyans (sharks, skates and

71

chimaeras), ray-finned fish and amphibians. We surveyed a phylogenetic blueprint to address the

72

impact of a model obesogen, TBT, acting via PPARγ, a critical regulator of adipogenesis. By

73

investigating representative species of the chordate phylum (Figure S1), we provide a detailed

2-4.

The list of confirmed and suspected obesogens is

3, 6, 7.

For instance, while in mammals TBT was shown to act as a PPARγ agonist 7, 19 (although 20, 21,

in the teleost fish Pleuronectes platessa (European plaice) a reporter gene assay

ACS Paragon Plus Environment

3

Page 5 of 25

Environmental Science & Technology

74

map of TBT activation/repression of PPARγ, as a proxy for potential obesogenic disorder. Our

75

study illustrates the central importance of comparative approaches to unravel the full ecosystem

76

impact of endocrine disruptors and contributes for the informed selection of model animal systems

77

in human health risk assessments.

ACS Paragon Plus Environment

4

Environmental Science & Technology

78

Materials and Methods

79

Sequence investigation

Page 6 of 25

80

Tblastn and blastp searches using human PPARα (Q07869), PPARβ (Q03181) and

81

PPARγ (P37231) as query were performed in the available databases NCBI Ensembl and JGI.

82

Ppar sequences were retrieved for the major vertebrate lineages: mammals (Homo sapiens, Mus

83

musculus); birds (Gallus gallus) and reptiles (Anolis carolinesis); amphibians (Xenopus

84

tropicalis); coelacanths (Latimeria chalumnae); Lepisosteiformes (Lepisosteus oculatus),

85

Osteoglossomorpha (Pantodon buchholzi) and Cypriniforme (Danio rerio), Chondrichthyes

86

(Callorhinchus milii and Leucoraja erinacea) and for the following invertebrates, tunicates

87

(Ciona intestinalis) cephalochordate (Branchiostoma floridae, Branchiostoma lanceolatum,

88

Branchiostoma belcheri) (Figure S1 and Table S1 for accession numbers).

89

Phylogenetic analysis

90

Phylogenetic analysis was performed using the set of sequences through database search

91

as well as the isolated sequences. A total of 36 sequences were aligned in MAFFT 25, 26 with L-

92

INS-I method. The resulting sequence alignment was stripped of all columns containing 90% gaps

93

leaving a total of 515 positions for phylogenetic analysis. Maximum likelihood phylogenetic

94

analysis was performed in PhyML V3.0 27 and the evolutionary model was determined using the

95

smart model selection (SMS) option resulting in a JTT +G. The branch support was calculated

96

using aBayes. The resulting tree was analyzed in Fig Tree V1.3.1 available at

97

http://tree.bio.ed.ac.uk/software/figtree/ and rooted with the invertebrate sequences.

98

Synteny maps

99

Comparative synteny maps were constructed using NCBI Gene database, using as

100

reference the latest available genome assemblies for the following species: H. sapiens

101

(GCF_000001405.33), L. oculatus (GCF_000242695.1), C. milii (GCF_000165045.1), B.

102

floridae (GG666565.1) and B. belcheri (AYSS01004225.)1. For each species, with the exception

ACS Paragon Plus Environment

5

Page 7 of 25

Environmental Science & Technology

103

of the cephalochordates that present a sole PPAR, we analysed the genomic location of PPARα,

104

PPARβ and PPARγ gene. The genomic locus of each PPAR gene was retrieved, as well as, the

105

five contiguous flanking genes to each side, when possible. Following the assembly of the synteny

106

maps, we proceeded to identify and localize the corresponding human orthologues of non-

107

conserved neighboring genes in the B. floridae and B belcheri. Orthology was determined through

108

the Ensembl orthologue-paralogue pipeline and our own phylogenetic analysis (not shown).

109

Finally, synteny maps and annotated orthologues were then used to infer the localization of the

110

ancestral PPAR gene in the reconstructed genome of the vertebrate ancestor using as reference

111

the reconstruction presented by Nakatani and colleagues

112

Putnam and colleagues 29.

113

Gene isolation and cloning

28

and the reconstruction presented by

114

Hinge and ligand binding domain (LBD) of PPARγ genes from Homo sapiens, X.

115

tropicalis, P. buchholzi, D. rerio, L. oculatus, L. erinacea and PPAR from B. lanceolatum (Table

116

S1) were predicted using NCBI's conserved domain database30 and isolated using a polymerase

117

chain reaction (PCR) approach (Table S2 for primer sequence). The hinge and LBD were then

118

cloned into the pBIND vector to produce a nuclear receptor (NR) LBD-Gal4 hybrid protein. This

119

hybrid protein, contain the DNA binding domain (DBD) of Gal4 and acts on an upstream

120

activation sequence (UAS) response element. The hinge and LBD of PPARγ from L. erinacea

121

was then cloned into the pCold-TF vector (Takara bio) to produce a PPARγ LBD-His6-tagged

122

trigger factor hybrid protein (HisTF-LePPARγ). Plasmid sequences were confirmed using Sanger

123

sequencing (GATC).

124

Homology modeling

125

Homology models were calculated in SWISS-MODEL homology modelling workspace

126

in alignment mode

127

submitted to SWISS-MODEL for homology modelling. Homology models were calculated using

128

ProMod3 and human PPARγ crystal structure 3WJ4 was selected as template. All models were

31, 32.

For each species PPARγ sequence corresponding to the LBD was

ACS Paragon Plus Environment

6

Environmental Science & Technology

Page 8 of 25

129

evaluated for the following parameters: sequence identity, GQME (Global model Quality

130

estimation), QMEAN 33, 34 and rendered as reliable (Table S3). Finally, models were visualized,

131

inspected and aligned to human crystal structure 3WJ4 in PyMOL v1.74 35.

132

Chemicals and solutions

133

All chemicals and reagents were obtained from Sigma-Aldrich unless stated otherwise in

134

the text. Three TBT-Cl (tributyltin chloride) and three TPT-Cl (triphenyltin chloride) (250μM,

135

100μM and 10μM), three ARA (Arachidonic acid), three EPA (cis-5,8,11,14,17-

136

Eicosapentaenoic acid) and three erucic acid (200mM, 100mM and 50mM), one Rosiglitazone

137

(Cayman) (1mM) and one T0070907 (Cayman) (1mM) solutions were prepared in DMSO

138

(dimethyl sulfoxide).

139

Cell culture and in vitro assays

140

Cell culture and transactivation assays were performed as described in André et al. 36.

141

Briefly, Cos-1 cells were seeded in 24-well culture plates. After 24 hours cells

142

were transfected using lipofectamine 2000 reagent (Invitrogen), in Opti-MEM reduced

143

serum medium (Gibco, Thermo Fisher), according manufacturer’s indications and 0.5 μg

144

of pBIND PPAR LBD-Gal4 or RXR LBD-Gal4 and 1μg of pGL4.31 luciferase reporter

145

vector, containing five UAS elements upstream the firefly luciferase reporter gene. For

146

heterodimer transfection assays 0.5 μg of pBind PPAR LBD-Gal4, 0.5 μg of pACT and

147

1μg of pGL4.31, were used. After 5 hours of incubation the medium was replaced with medium

148

containing the test compounds (TBT – 10nM, 100nM and 250nM; TPT– 10nM, 100nM and

149

250nM; ARA - 50µM, 100µM and 200µM; EPA - 50µM, 100µM and 200µM; Erucic acid -

150

50µM, 100µM and 200µM; Rosiglitazone (Rosi) - 1µM and T0070907 - 1µM) dissolved in

151

DMSO (0.1%). DMSO (0.1%) was used as solvent control. 24 hours later cells were lysed and

152

the luminescent activities of Firefly luciferase (pGL4.31) and Renilla luciferase (pBIND) were

153

determined. All assays were performed, independently, five times and each time two technical

ACS Paragon Plus Environment

7

Page 9 of 25

Environmental Science & Technology

154

replicates per condition were performed to validate the results. Two different batches of Cos-1

155

cells were used; three replicates per batch were included per assay.

156

Ligand Binding Assay

157

The L. erinacea PPARγ LBD-His6-tagged trigger factor hybrid protein were expressed

158

in Escherichia coli BL21 (DE3) cells and purified by using His-select nickel affinity gel (Sigma

159

Aldrich). Ligand binding assay was assessed as previously described (39-41). In brief, the purified

160

protein (100, 200, 400 or 800 µg/ml) were incubated with 200 nM of tri[U-14C]phenyltin

161

hydroxide ([14C]TPT; radiochemical purity > 96.6 %, 2.04 GBq/mmol; Amersham Biosciences)

162

or 50 nM of [3H]rosiglitazone ([3H]rosiglitazone; radiochemical purity > 97 %, 2.168 TBq/mmol;

163

PerkinElmer). After incubation at 4 ◦C for 1 h, hydroxyapatite was added to precipitate the

164

receptor protein and bound radioactive compounds. The hydroxyapatite pellets was washed and

165

then radioactivity in the pellet was determined by liquid scintillation. Binding in the presence of

166

a 100-fold molar excess of unlabeled compound was defined as nonspecific binding; specific

167

binding was defined as total binding minus nonspecific binding. For competition binding assay,

168

the purified protein (800 µg/ml) were incubated with 200 nM of [14C]TPT and increasing

169

concentrations of test compounds.

170

Statistical analysis

171

Transactivation data was normalized dividing by the solvent control (DMSO) and the

172

mean of the technical replicates was used in the statistical analysis. One-way ANOVA followed

173

by Tukey post hoc test was used to analyze differences when parametric criteria was achieved,

174

otherwise a non-parametric Kruskal-wallis ANOVA followed by a Games-Howel test was used.

175

All analyses were performed using IBM SPSS Statistics 24. P Thr and Phe363>Met. Additionally, Cys285 is replaced by a tyrosine in

323

butterfly fish, though conserved in the spotted gar (Figure 5A). The butterfly fish homology model

324

reveals that the aromatic ring of tyrosine projects into the LBP, altering the pocket organization

325

as observed in human PPARγ (Figures 5A, D and E). Also, the replacement of Ile281 by a Phe in

326

Helix 3 in both spotted gar and butterfly fish results in a second aromatic ring protruding into the

ACS Paragon Plus Environment

17

Page 19 of 25

Environmental Science & Technology

327

LBP (Figures 5A, 5D and 5E). We next constructed an homology model of the LBP of zebrafish

328

PPARγ, a popular comparative model in lipid endocrine disruption studies

329

interacting residues of the PPARγ-LBP in zebrafish are conserved with those of the butterfly fish,

330

in line with the incapacity of TBT to bind and transactivate PPARγ (Figures 5A and 5E). Finally,

331

B. lanceolatum presents the lowest number of conserved TBT-interacting residues; aside from the

332

replacement of the Cys285 by a leucine, we find six additional substitutions, including amino

333

acids that were found to be totally conserved in vertebrates such as Phe282>Leu and Phe360>Leu

334

(Figure 5G). Overall, our analysis strongly suggests that, Cys285 is essential for PPARγ

335

transactivation in the presence of TBT.

9, 10.

The TBT-

336 337

Conserved and derived molecular features in PPARγ exploitation by organotins

338

Animal endocrine systems encompass an intricate network of signaling cascades, that

339

when perturbated lead to episodes of disruption with developmental and physiological

340

impairment 48. In this context, a critical and complex issue is that of how man-made chemicals

341

differently exploit the genetic and physiological make-up of animal species and influence

342

homeostasis. This question assumes a key importance, since animal model-derived studies are

343

essential to evaluate the risk of man-made chemicals in human populations. In recent years, the

344

concept of obesogens has gathered increased attention as the obesity epidemic gained pace in

345

western societies. The overall feature of these compounds is their capacity to activate the

346

mammalian adipogenesis master control gene PPAR; yet, PPAR is largely uncharacterized in

347

most vertebrate lineages and in some species seems irresponsive to a model obesogen, TBT 22, 23,

348

47.

349

heterodimer emerged in the ancestor of vertebrates concomitant with the appearance of PPAR

350

(Fig.7). This might result in potential obesogenic impacts at environmental conditions which

351

require further investigations. The conversion of activation to a neutral/negative profile observed

352

in the tested teleost represents a secondary event dictated by changes in the ligand pocket of

353

PPAR. These observations raise a paradoxical scenario. While both mammals and teleosts

Considering this set-up, we show that the TBT agonist capacity towards PPARγ/RXR

ACS Paragon Plus Environment

18

Environmental Science & Technology

Page 20 of 25

354

experience obesogenic outcomes upon organotin exposure

355

mechanisms of such disruption are possibly different, since the tested teleost PPAR does not

356

activate transcription in the presence of TBT

357

apparent absence of obesogenic impacts in fish. Yet, TBT was suggested to inhibit teleost

358

PPAR/

359

accumulation 9. Other hypothesis to consider includes the activation of the heterodimer

360

PPARγ/RXR by TBT through the RXR monomer 10, although our in vitro results are contradictory

361

with this scenario. Overall, our results pave the way for future functional and in vivo studies to

362

investigate the diversity of physiological landscapes resulting from the exposure to obesogens.

22,

22

7, 9, 10, 49

(Fig. 7), the molecular

(this work) (Fig. 7). This would involve an

eventually suppressing lipid oxidation in vivo, which could potentially cause fat

363

364 365 366

Fig. 7 – The chordate taxonomic scope of PPARγ exploitation by TBT. PPARs transactivation results in the presence of TBT.

367 368

Acknowledgments

369

This work was supported by Norte2020 and FEDER (Coral - Sustainable Ocean Exploitation -

370

Norte-01-0145-FEDER-000036), Project No. 031342 co-financed by COMPETE 2020, Portugal

371

2020 and the European Union through the ERDF, and by FCT through national funds. We

372

acknowledge the Branchiostoma lanceolatum genome consortium for providing access to the

373

genome sequence and Fundação para a Ciência e a Tecnologia for the support to AC

374

(SFRH/BD/90664/2012), RR (SFRH/BPD/72519/2010) and to EF (SFRH/BD/100262/2014). We

375

acknowledge Ana André and Ricardo Capela for their help in the sampling process.

376

The authors declare no competing financial interest.

ACS Paragon Plus Environment

19

Page 21 of 25

Environmental Science & Technology

377 378 379

Supporting Information.

380

Additional results and experimental details regarding gene isolation, homology modeling,

381

transactivation assays and synteny maps.

382

ACS Paragon Plus Environment

20

Environmental Science & Technology

Page 22 of 25

383

References

384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435

1. Eide, M.; Rydbeck, H.; Tørresen, O. K.; Lille-Langøy, R.; Puntervoll, P.; Goldstone, J. V.; Jakobsen, K. S.; Stegeman, J.; Goksøyr, A.; Karlsen, O. A., Independent losses of a xenobiotic receptor across teleost evolution. Scientific reports 2018, 8, (1), 10404. 2. Baillie-Hamilton, P. F., Chemical toxins: a hypothesis to explain the global obesity epidemic. Journal of alternative and complementary medicine (New York, N.Y.) 2002, 8, (2), 185-92. 3. Grun, F.; Blumberg, B., Environmental obesogens: organotins and endocrine disruption via nuclear receptor signaling. Endocrinology 2006, 147, (6 Suppl), S50-5. 4. Kanayama, T.; Kobayashi, N.; Mamiya, S.; Nakanishi, T.; Nishikawa, J., Organotin compounds promote adipocyte differentiation as agonists of the peroxisome proliferatoractivated receptor gamma/retinoid X receptor pathway. Molecular pharmacology 2005, 67, (3), 766-74. 5. Holtcamp, W., Obesogens: an environmental link to obesity. Environmental health perspectives 2012, 120, (2), a62-8. 6. Evans, R. M.; Barish, G. D.; Wang, Y. X., PPARs and the complex journey to obesity. Nature medicine 2004, 10, (4), 355-61. 7. Grun, F.; Watanabe, H.; Zamanian, Z.; Maeda, L.; Arima, K.; Cubacha, R.; Gardiner, D. M.; Kanno, J.; Iguchi, T.; Blumberg, B., Endocrine-disrupting organotin compounds are potent inducers of adipogenesis in vertebrates. Molecular endocrinology (Baltimore, Md.) 2006, 20, (9), 2141-55. 8. Capitao, A.; Lyssimachou, A.; Castro, L. F. C.; Santos, M. M., Obesogens in the aquatic environment: an evolutionary and toxicological perspective. Environment international 2017, 106, 153-169. 9. Lyssimachou, A.; Santos, J. G.; Andre, A.; Soares, J.; Lima, D.; Guimaraes, L.; Almeida, C. M.; Teixeira, C.; Castro, L. F.; Santos, M. M., The Mammalian "Obesogen" Tributyltin Targets Hepatic Triglyceride Accumulation and the Transcriptional Regulation of Lipid Metabolism in the Liver and Brain of Zebrafish. PLoS One 2015, 10, (12), e0143911. 10. Ouadah-Boussouf, N.; Babin, P. J., Pharmacological evaluation of the mechanisms involved in increased adiposity in zebrafish triggered by the environmental contaminant tributyltin. Toxicology and applied pharmacology 2016, 294, 32-42. 11. Janer, G.; Navarro, J. C.; Porte, C., Exposure to TBT increases accumulation of lipids and alters fatty acid homeostasis in the ramshorn snail Marisa cornuarietis. Comparative biochemistry and physiology. Toxicology & pharmacology : CBP 2007, 146, (3), 368-74. 12. Jordao, R.; Casas, J.; Fabrias, G.; Campos, B.; Pina, B.; Lemos, M. F.; Soares, A. M.; Tauler, R.; Barata, C., Obesogens beyond Vertebrates: Lipid Perturbation by Tributyltin in the Crustacean Daphnia magna. Environmental health perspectives 2015, 123, (8), 813-9. 13. Jordao, R.; Campos, B.; Pina, B.; Tauler, R.; Soares, A. M.; Barata, C., Mechanisms of Action of Compounds That Enhance Storage Lipid Accumulation in Daphnia magna. Environmental science & technology 2016, 50, (24), 13565-13573. 14. Riu, A.; McCollum, C. W.; Pinto, C. L.; Grimaldi, M.; Hillenweck, A.; Perdu, E.; Zalko, D.; Bernard, L.; Laudet, V.; Balaguer, P.; Bondesson, M.; Gustafsson, J. A., Halogenated bisphenol-A analogs act as obesogens in zebrafish larvae (Danio rerio). Toxicological sciences : an official journal of the Society of Toxicology 2014, 139, (1), 48-58. 15. Zhang, J.; Zuo, Z.; He, C.; Cai, J.; Wang, Y.; Chen, Y.; Wang, C., Effect of tributyltin on testicular development in Sebastiscus marmoratus and the mechanism involved. Environmental toxicology and chemistry 2009, 28, (7), 1528-35. 16. Zhang, J.; Zuo, Z.; Xiong, J.; Sun, P.; Chen, Y.; Wang, C., Tributyltin exposure causes lipotoxicity responses in the ovaries of rockfish, Sebastiscus marmoratus. Chemosphere 2013, 90, (3), 1294-9. 17. Castro, L. F.; Santos, M. M., To bind or not to bind: the taxonomic scope of nuclear receptor mediated endocrine disruption in invertebrate phyla. Environmental science & technology 2014, 48, (10), 5361-3.

ACS Paragon Plus Environment

21

Page 23 of 25

436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490

Environmental Science & Technology

18. Santos, M.; Reis-Henriques, M.; Castro, L., Lipid homeostasis perturbation by organotins : Effects on vertebrates and invertebrates. . In Biochemical and Physiological Effects of Organotins, 2012; pp 83-96. 19. Yanik, S. C.; Baker, A. H.; Mann, K. K.; Schlezinger, J. J., Organotins are potent activators of PPARgamma and adipocyte differentiation in bone marrow multipotent mesenchymal stromal cells. Toxicological sciences : an official journal of the Society of Toxicology 2011, 122, (2), 476-88. 20. le Maire, A.; Grimaldi, M.; Roecklin, D.; Dagnino, S.; Vivat-Hannah, V.; Balaguer, P.; Bourguet, W., Activation of RXR-PPAR heterodimers by organotin environmental endocrine disruptors. EMBO reports 2009, 10, (4), 367-73. 21. Harada, S.; Hiromori, Y.; Nakamura, S.; Kawahara, K.; Fukakusa, S.; Maruno, T.; Noda, M.; Uchiyama, S.; Fukui, K.; Nishikawa, J.-i.; Nagase, H.; Kobayashi, Y.; Yoshida, T.; Ohkubo, T.; Nakanishi, T., Structural basis for PPARγ transactivation by endocrine-disrupting organotin compounds. Scientific reports 2015, 5, 8520. 22. Colliar, L.; Sturm, A.; Leaver, M. J., Tributyltin is a potent inhibitor of piscine peroxisome proliferator-activated receptor alpha and beta. Comparative biochemistry and physiology. Toxicology & pharmacology : CBP 2011, 153, (1), 168-73. 23. Janesick, A.; Blumberg, B., Minireview: PPARgamma as the target of obesogens. The Journal of steroid biochemistry and molecular biology 2011, 127, (1-2), 4-8. 24. Pereira-Fernandes, A.; Vanparys, C.; Hectors, T. L.; Vergauwen, L.; Knapen, D.; Jorens, P. G.; Blust, R., Unraveling the mode of action of an obesogen: mechanistic analysis of the model obesogen tributyltin in the 3T3-L1 cell line. Molecular and cellular endocrinology 2013, 370, (1-2), 52-64. 25. Katoh, K.; Kuma, K.; Toh, H.; Miyata, T., MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic acids research 2005, 33, (2), 511-8. 26. Katoh, K.; Toh, H., Recent developments in the MAFFT multiple sequence alignment program. Briefings in bioinformatics 2008, 9, (4), 286-98. 27. Guindon, S.; Dufayard, J. F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O., New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic biology 2010, 59, (3), 307-21. 28. Nakatani, Y.; Takeda, H.; Kohara, Y.; Morishita, S., Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates. Genome research 2007, 17, (9), 1254-65. 29. Putnam, N. H.; Butts, T.; Ferrier, D. E.; Furlong, R. F.; Hellsten, U.; Kawashima, T.; Robinson-Rechavi, M.; Shoguchi, E.; Terry, A.; Yu, J. K.; Benito-Gutierrez, E. L.; Dubchak, I.; Garcia-Fernandez, J.; Gibson-Brown, J. J.; Grigoriev, I. V.; Horton, A. C.; de Jong, P. J.; Jurka, J.; Kapitonov, V. V.; Kohara, Y.; Kuroki, Y.; Lindquist, E.; Lucas, S.; Osoegawa, K.; Pennacchio, L. A.; Salamov, A. A.; Satou, Y.; Sauka-Spengler, T.; Schmutz, J.; Shin, I. T.; Toyoda, A.; Bronner-Fraser, M.; Fujiyama, A.; Holland, L. Z.; Holland, P. W.; Satoh, N.; Rokhsar, D. S., The amphioxus genome and the evolution of the chordate karyotype. Nature 2008, 453, (7198), 1064-71. 30. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C. J.; Lu, S.; Chitsaz, F.; Derbyshire, M. K.; Geer, R. C.; Gonzales, N. R.; Gwadz, M.; Hurwitz, D. I.; Lu, F.; Marchler, G. H.; Song, J. S.; Thanki, N.; Wang, Z.; Yamashita, R. A.; Zhang, D.; Zheng, C.; Geer, L. Y.; Bryant, S. H., CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic acids research 2017, 45, (D1), D200-D203. 31. Arnold, K.; Bordoli, L.; Kopp, J.; Schwede, T., The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics (Oxford, England) 2006, 22, (2), 195-201. 32. Biasini, M.; Bienert, S.; Waterhouse, A.; Arnold, K.; Studer, G.; Schmidt, T.; Kiefer, F.; Gallo Cassarino, T.; Bertoni, M.; Bordoli, L.; Schwede, T., SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic acids research 2014, 42, (Web Server issue), W252-8. 33. Benkert, P.; Biasini, M.; Schwede, T., Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics (Oxford, England) 2011, 27, (3), 343-50.

ACS Paragon Plus Environment

22

Environmental Science & Technology

491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545

Page 24 of 25

34. Benkert, P.; Kunzli, M.; Schwede, T., QMEAN server for protein model quality estimation. Nucleic acids research 2009, 37, (Web Server issue), W510-4. 35. Schrodinger, L. The PyMOL Molecular Graphics System, Version 1.3. , 2010. 36. Andre, A.; Ruivo, R.; Capitao, A.; Froufe, E.; Pascoa, I.; Costa Castro, L. F.; Santos, M. M., Cloning and functional characterization of a retinoid X receptor orthologue in Platynereis dumerilii: An evolutionary and toxicological perspective. Chemosphere 2017, 182, 753-761. 37. Bertrand, S.; Thisse, B.; Tavares, R.; Sachs, L.; Chaumot, A.; Bardet, P. L.; Escriva, H.; Duffraisse, M.; Marchand, O.; Safi, R.; Thisse, C.; Laudet, V., Unexpected novel relational links uncovered by extensive developmental profiling of nuclear receptor expression. PLoS genetics 2007, 3, (11), e188. 38. Madureira, T. V.; Pinheiro, I.; de Paula Freire, R.; Rocha, E.; Castro, L. F.; Urbatzka, R., Genome specific PPARalphaB duplicates in salmonids and insights into estrogenic regulation in brown trout. Comparative biochemistry and physiology. Part B, Biochemistry & molecular biology 2017, 208-209, 94-101. 39. Leaver, M. J.; Boukouvala, E.; Antonopoulou, E.; Diez, A.; Favre-Krey, L.; Ezaz, M. T.; Bautista, J. M.; Tocher, D. R.; Krey, G., Three peroxisome proliferator-activated receptor isotypes from each of two species of marine fish. Endocrinology 2005, 146, (7), 3150-62. 40. Kaur, S.; Jobling, S.; Jones, C. S.; Noble, L. R.; Routledge, E. J.; Lockyer, A. E., The Nuclear Receptors of Biomphalaria glabrata and Lottia gigantea: Implications for Developing New Model Organisms. PLOS ONE 2015, 10, (4), e0121259. 41. Vogeler, S.; Galloway, T. S.; Isupov, M.; Bean, T. P., Cloning retinoid and peroxisome proliferator-activated nuclear receptors of the Pacific oyster and in silico binding to environmental chemicals. PLoS One 2017, 12, (4), e0176024. 42. Fonseca, E.; Ruivo, R.; Lopes-Marques, M.; Zhang, H.; Santos, M. M.; Venkatesh, B.; Castro, L. F., LXRalpha and LXRbeta Nuclear Receptors Evolved in the Common Ancestor of Gnathostomes. Genome biology and evolution 2017, 9, (1), 222-230. 43. Jaillon, O.; Aury, J. M.; Brunet, F.; Petit, J. L.; Stange-Thomann, N.; Mauceli, E.; Bouneau, L.; Fischer, C.; Ozouf-Costaz, C.; Bernot, A.; Nicaud, S.; Jaffe, D.; Fisher, S.; Lutfalla, G.; Dossat, C.; Segurens, B.; Dasilva, C.; Salanoubat, M.; Levy, M.; Boudet, N.; Castellano, S.; Anthouard, V.; Jubin, C.; Castelli, V.; Katinka, M.; Vacherie, B.; Biemont, C.; Skalli, Z.; Cattolico, L.; Poulain, J.; De Berardinis, V.; Cruaud, C.; Duprat, S.; Brottier, P.; Coutanceau, J. P.; Gouzy, J.; Parra, G.; Lardier, G.; Chapple, C.; McKernan, K. J.; McEwan, P.; Bosak, S.; Kellis, M.; Volff, J. N.; Guigo, R.; Zody, M. C.; Mesirov, J.; Lindblad-Toh, K.; Birren, B.; Nusbaum, C.; Kahn, D.; Robinson-Rechavi, M.; Laudet, V.; Schachter, V.; Quetier, F.; Saurin, W.; Scarpelli, C.; Wincker, P.; Lander, E. S.; Weissenbach, J.; Roest Crollius, H., Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate protokaryotype. Nature 2004, 431, (7011), 946-57. 44. Braasch, I.; Gehrke, A. R.; Smith, J. J.; Kawasaki, K.; Manousaki, T.; Pasquier, J.; Amores, A.; Desvignes, T.; Batzel, P.; Catchen, J.; Berlin, A. M.; Campbell, M. S.; Barrell, D.; Martin, K. J.; Mulley, J. F.; Ravi, V.; Lee, A. P.; Nakamura, T.; Chalopin, D.; Fan, S.; Wcisel, D.; Canestro, C.; Sydes, J.; Beaudry, F. E.; Sun, Y.; Hertel, J.; Beam, M. J.; Fasold, M.; Ishiyama, M.; Johnson, J.; Kehr, S.; Lara, M.; Letaw, J. H.; Litman, G. W.; Litman, R. T.; Mikami, M.; Ota, T.; Saha, N. R.; Williams, L.; Stadler, P. F.; Wang, H.; Taylor, J. S.; Fontenot, Q.; Ferrara, A.; Searle, S. M.; Aken, B.; Yandell, M.; Schneider, I.; Yoder, J. A.; Volff, J. N.; Meyer, A.; Amemiya, C. T.; Venkatesh, B.; Holland, P. W.; Guiguen, Y.; Bobe, J.; Shubin, N. H.; Di Palma, F.; Alfoldi, J.; Lindblad-Toh, K.; Postlethwait, J. H., The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons. Nature genetics 2016, 48, (4), 427-37. 45. Lammi, J.; Perlmann, T.; Aarnisalo, P., Corepressor interaction differentiates the permissive and non-permissive retinoid X receptor heterodimers. Archives of biochemistry and biophysics 2008, 472, (2), 105-14. 46. Gampe, R. T., Jr.; Montana, V. G.; Lambert, M. H.; Miller, A. B.; Bledsoe, R. K.; Milburn, M. V.; Kliewer, S. A.; Willson, T. M.; Xu, H. E., Asymmetry in the PPARg/RXRa; Crystal Structure Reveals the Molecular Basis of Heterodimerization among Nuclear Receptors. Molecular cell 2000, 5, (3), 545-555.

ACS Paragon Plus Environment

23

Page 25 of 25

546 547 548 549 550 551 552 553 554

Environmental Science & Technology

47. Harada, S.; Hiromori, Y.; Nakamura, S.; Kawahara, K.; Fukakusa, S.; Maruno, T.; Noda, M.; Uchiyama, S.; Fukui, K.; Nishikawa, J.; Nagase, H.; Kobayashi, Y.; Yoshida, T.; Ohkubo, T.; Nakanishi, T., Structural basis for PPARgamma transactivation by endocrinedisrupting organotin compounds. Scientific reports 2015, 5, 8520. 48. Kabir, E. R.; Rahman, M. S.; Rahman, I., A review on endocrine disruptors and their possible impacts on human health. Environmental toxicology and pharmacology 2015, 40, (1), 241-58. 49. Tingaud-Sequeira, A.; Ouadah, N.; Babin, P. J., Zebrafish obesogenic test: a tool for screening molecules that target adiposity. Journal of lipid research 2011, 52, (9), 1765-72.

555

ACS Paragon Plus Environment

24