Of Retinoids and Organotins: The Evolution of the Retinoid X Receptor in Metazoa.

Nuclear receptors (NRs) are transcription factors accomplishing a multiplicity of functions, essential for organismal homeostasis. Among their numerous members, the retinoid X receptor (RXR) is a central player of the endocrine system, with a singular ability to operate as a homodimer or a heterodimer with other NRs. Additionally, RXR has been found to be a critical actor in various processes of endocrine disruption resulting from the exposure to a known class of xenobiotics termed organotins (e.g., tributyltin (TBT)), including imposex in gastropod molluscs and lipid perturbation across different metazoan lineages. Thus, given its prominent physiological and endocrine role, RXR is present in the genomes of most extant metazoan species examined to date. Here, we expand on the phylogenetic distribution of RXR across the metazoan tree of life by exploring multiple next-generation sequencing projects of protostome lineages. By addressing amino acid residue conservation in combination with cell-based functional assays, we show that RXR induction by 9-cis retinoic acid (9cisRA) and TBT is conserved in more phyla than previously described. Yet, our results highlight distinct activation efficacies and alternative modes of RXR exploitation by the organotin TBT, emphasizing the need for broader species sampling to clarify the mechanistic activation of RXR.

The ongoing explosion in the assembly of full genome sequences from multiple species/lineages is key to understand NR gene diversification. RXR, for instance, has been found in the genomes of most extant metazoan species examined to date [15,16]. Nevertheless, the lack of data, especially from invertebrate protostomes, and the poor taxonomic sampling of most studies can lead to

Sampling and RNA Extraction
Bugula neritina specimens were collected at the Marina da Póvoa de Varzim, Portugal; Phoronopsis californica, and Megathiris detruncata specimens were collected at Madeira Island, Portugal; Bonellia viridis was sampled at Peniche, Portugal; Priapulus caudatus and Xenoturbella bocki specimens were collected at the Gullmarn fjord, Sweden. The specimens were preserved on RNAlater (Invitrogen, Carlsbad, CA, USA) immediately after their collection for further total RNA extraction. A small portion of each specimen was homogenized with PureZOL RNA Isolation Reagent ® (Bio-Rad, Hercules, CA, USA), and chloroform was used to extract the nucleic acids according to the manufacturer's instructions. The illustra RNAspin Mini RNA Isolation (GE Healthcare, Chicago, IL, USA) kit was used to isolate the total RNA from the aqueous phase obtained in the first step. and genomic contamination was prevented by an on-column DNAse I digestion step. RNA was eluted in 30 µL of RNase-free water. The synthesis of first-strand cDNA (1 µg) and of 5'and 3´cDNA for rapid amplification of cDNA ends (RACE) was performed with the iScript™cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA) and the SMARTer™ RACE cDNA Amplification kit (Clontech, Mountain View, CA, USA), respectively, following to the manufacturer's guidelines.

RXR Gene Isolation
The open reading frames (ORF) of RXR genes were obtained by a combination of RACE, nested, hemi-nested, and/or degenerated polymerase chain reactions (PCRs), using Phusion Flash High-Fidelity PCR Master Mix (Thermo Fisher Scientific Co., Waltham, MA, USA) according to the manufacturer's protocol. Partial gene sequence of B. viridis RXR and B. neritina RXR1 and RXR2 and two partial RXR-like sequences from X. bocki were retrieved through a BLAST approach using Sequence Read Archive (SRA) files. A partial gene sequence of P. californica RXR and M. detruncata RXR were obtained by PCR with degenerated primers (Supplementary Material Table S1) designed with the CODEHOP software [45]. The sequences obtained with degenerated PCRs and BLAST searches were used to design specific primers (Supplementary Material Table S1) with the Primer3 (v.0.4.0) software [46,47]. All PCR products were purified using NZYGelpure (Nzytech, Lisbon, Portugal), cloned into Nzy5α competent cells (Nzytech, Lisbon, Portugal) using the pGEM-T Easy Vector System (Promega, Madison, WI, USA), and the sequences were confirmed by sequencing (Eurofins GATC, Constance, Germany). The sequences have been deposited in GenBank (Accession numbers: P. californica RXR MT264997, M. detruncata RXR MT264998, B. viridis RXR MT264999, B. neritina RXR1 MT265000, B. neritina RXR2 MT265001, and X. bocki RXR MT265002). Trichoplax adhaerens and Aurelia aurita RXR ORFs (Table 1) were chemically synthetized by NZYTech (Lisbon, Portugal).

Sequence and Phylogenetic Analysis
The amino acid sequences from the isolated RXR genes were aligned with human RAR, HNF4, COUP and RXR amino acid sequences and RXR amino acid sequences from other metazoans retrieved from the GenBank database and blast researches on SRA files ( Table 1) The following parameters were used: generation number = 10,000,000, rate matrix for aa = mixed (Jones), nruns = 2, nchains = 4, temp = 0.20, sampling set to 1000, and burnin to 0.25. The statistical support for each branch was expressed as Bayesian posterior probabilities [50] and indicated at the nodes. FigTree (v.1.3.1) was used to visualize the tree. The previous alignment was visualized and edited in Geneious ® (v7.1.7) to identify the residues known to interact with 9cisRA and TBT, based on previous studies [29,42,[51][52][53].

Statistical Analysis
The differences between the mean of the technical replicates were tested for significance using Student's t-test and one-way analysis of variance (ANOVA), followed by the Holm-Sidak test. Data were transformed whenever the homogeneity of variances and/or the normality failed. The level of significance (P-value) was set to 0.05. Statistical analysis was performed using SigmaPlot software (v.11.0).

Phylogenetic and Sequence Analyses of RXR in Metazoan Lineages
In the present study, we investigated a set of currently available genomes and transcriptomes from multiple metazoan phyla to retrieve partial RXR amino acid sequences. By combining database mining with PCR approaches, we were able to successfully isolate full-length RXR ortholog genes from several species, validating the occurrence of RXRs within the selected Metazoa lineages (Table 1). Phylogenetic analysis robustly corroborated the orthology of the new isolated RXR genes ( Figure 1). For most newly examined species, a single RXR ortholog gene was retrieved. A notable exception were bryozoans (B. neritina and M. membranacea), for which two RXR genes were identified ( Figure 1). The amino acid sequence alignment of the LBD of RXRs revealed that the 18 amino acid residues located in the α-helices H3, H5, H7, and H11 and in the β-turn of the LBD, which confer the ability of binding to 9cisRA and induce a transcriptional response in human RXRα [51,52], are overall conserved in the new isolated RXRs (Figure 2), suggesting that these receptors could be able to interact with 9cisRA. The exceptions were RXR2 from B. neritina (Bryozoa), exhibiting one nonconserved residue, and the RXRs from P. caudatus (Priapulida), T. adhaerens (Placozoa), and A. aurita (Cnidaria), with one, four, and three non-conserved residues, respectively ( Figure 2). On the other hand, the cysteine residue in helix H11, crucial for transactivation of human RXRα by organotins (TBT and TPT) [42,56], is substituted in bryozoan (B. neritina and M. membranacea) RXR1s, cnidarian (A. aurita), and placozoan (T. adhaerens) RXRs, suggesting a possible inability to The amino acid sequence alignment of the LBD of RXRs revealed that the 18 amino acid residues located in the α-helices H3, H5, H7, and H11 and in the β-turn of the LBD, which confer the ability of binding to 9cisRA and induce a transcriptional response in human RXRα [51,52], are overall conserved in the new isolated RXRs (Figure 2), suggesting that these receptors could be able to interact with 9cisRA. The exceptions were RXR2 from B. neritina (Bryozoa), exhibiting one non-conserved residue, and the RXRs from P. caudatus (Priapulida), T. adhaerens (Placozoa), and A. aurita (Cnidaria), with one, four, and three non-conserved residues, respectively ( Figure 2). On the other hand, the cysteine residue in helix H11, crucial for transactivation of human RXRα by organotins (TBT and TPT) [42,56], is substituted in bryozoan (B. neritina and M. membranacea) RXR1s, cnidarian (A. aurita), and placozoan (T. adhaerens) RXRs, suggesting a possible inability to respond to such compounds. In addition, and unlike other invertebrate RXRs, 17 to 21 amino acid insertions were also found in B. neritina and M. membranacea RXR2s, respectively. respond to such compounds. In addition, and unlike other invertebrate RXRs, 17 to 21 amino acid insertions were also found in B. neritina and M. membranacea RXR2s, respectively.

In vitro Interaction of RXRs with 9cisRA
Since the main residues of the human RXRα pocket known to interact with 9cisRA are conserved in most of the isolated RXRs, we tested their ability to activate gene expression upon binding to 9cisRA, using a luciferase reporter gene assay ( Figure 3). Our results showed a significant activation (P < 0.05) of the chimeric RXRs in the presence of 1 μM 9cisRA, yielding weaker activations than the human RXRα, with the exceptions of the phoronid P. californica (P = 0.073) and the placozoan T. adhaerens (P = 0.176) RXRs. Unexpectedly, the brachiopod M. detruncata RXR elicited an activation stronger than the human control (P = 0.002). Also, among the analyzed RXRs, low affinity towards 9cisRA was verified with respect to the bryozoan B. neritina RXR2, the priapulid P. caudatus, the xenoturbellid X. bocki, and the cnidarian A. aurita RXRs. Interestingly, the efficacy/effectiveness of 9cisRA was significantly different (P < 0.001) for the bryozoan RXRs.

In vitro Interaction of RXRs with 9cisRA
Since the main residues of the human RXRα pocket known to interact with 9cisRA are conserved in most of the isolated RXRs, we tested their ability to activate gene expression upon binding to 9cisRA, using a luciferase reporter gene assay (Figure 3). Our results showed a significant activation (P < 0.05) of the chimeric RXRs in the presence of 1 µM 9cisRA, yielding weaker activations than the human RXRα, with the exceptions of the phoronid P. californica (P = 0.073) and the placozoan T. adhaerens (P = 0.176) RXRs. Unexpectedly, the brachiopod M. detruncata RXR elicited an activation stronger than the human control (P = 0.002). Also, among the analyzed RXRs, low affinity towards 9cisRA was verified with respect to the bryozoan B. neritina RXR2, the priapulid P. caudatus, the xenoturbellid X. bocki, and the cnidarian A. aurita RXRs. Interestingly, the efficacy/effectiveness of 9cisRA was significantly different (P < 0.001) for the bryozoan RXRs.

In vitro Interaction of RXRs with TBT
As classical targets of endocrine disruption by organotins [24,26,38,39], we tested these RXRs for their ability to be activated by TBT and to promote the expression of the luciferase reporter gene (Figure 4). With the exception of the placozoan T. adhaerens RXR, all the tested RXRs were able to significantly (P < 0.05) activate gene expression, at least at the higher tested concentrations of TBT (100 and 250 nM). Once again, we found a stronger response for the brachiopod M. detruncata RXR when compared to the human RXRα, and weaker responses for the remaining receptors. Even without conservation of the previously described cysteine residue, suggested to be key for TBT binding and mediation of transactivation [42,56], the cnidarian A. aurita RXR and the bryozoan B. neritina RXR1 were able to promote gene transcription in the presence of TBT, similarly to the rotifer B. koreanus RXR [40]. Curiously, and despite the conservation of the key cysteine residue in the bryozoan B. neritina RXR2, both subtypes responded similarly to TBT (P > 0.05), indicating that TBT possibly interacts with the pocket environment of both receptors in a similar fashion. Data represent means ± SEM from three independent experiments (n = 3). The results were normalized to the control condition (DMSO without ligand). Significant differences between 9cisRA and the solvent control were inferred with Student's t-test. Asterisks denote significant differences (* P ≤ 0.05, ** P ≤ 0.01, *** P ≤ 0.001). Significant differences among the different RXRs in the test condition were inferred with one-way ANOVA (source of variation: F = 61.9, P < 0.001). The different letters denote significant differences. Hsa, H. sapiens; Mde, M. detruncata; Bvi, B. viridis; Pcal, P. californica; Bne, B. neritina, Pcau, P. caudatus; Xbo, X. bocki; Tad, T. adhaerens; Aau, A. aurita.

In vitro Interaction of RXRs with TBT
As classical targets of endocrine disruption by organotins [24,26,38,39], we tested these RXRs for their ability to be activated by TBT and to promote the expression of the luciferase reporter gene (Figure 4). With the exception of the placozoan T. adhaerens RXR, all the tested RXRs were able to significantly (P < 0.05) activate gene expression, at least at the higher tested concentrations of TBT (100 and 250 nM). Once again, we found a stronger response for the brachiopod M. detruncata RXR when compared to the human RXRα, and weaker responses for the remaining receptors. Even without conservation of the previously described cysteine residue, suggested to be key for TBT binding and mediation of transactivation [42,56], the cnidarian A. aurita RXR and the bryozoan B. neritina RXR1 were able to promote gene transcription in the presence of TBT, similarly to the rotifer B. koreanus RXR [40]. Curiously, and despite the conservation of the key cysteine residue in the bryozoan B. neritina RXR2, both subtypes responded similarly to TBT (P > 0.05), indicating that TBT possibly interacts with the pocket environment of both receptors in a similar fashion. . The results were normalized to the control condition (DMSO without ligand). Significant differences between the tested concentrations of TBT and the solvent control were inferred with one-way ANOVA. The F and P values of the source of variation are indicated below each species name abbreviation. Asterisks denote significant differences (* P ≤ 0.05, ** P ≤ 0.01, *** P ≤ 0.001). Significant differences between BneRXR1 and BneRXR2 with TBT at the same concentration were inferred with Student's t-test. The green rectangles indicate the RXRs which have the cysteine residue previously described as crucial for binding to TBT. Hsa, H. sapiens

Mutagenesis of Bryozoan Specific RXR Orthologs
To further clarify the transcriptional responses obtained with the bryozoan B. neritina RXRs, with both 9cisRA and TBT, site-directed and deletion mutants were tested ( Figure 5). First, the bryozoan RXR1 Thr388Cys mutant (Bne1 T > C) was assayed for its ability to elicit a stronger response to TBT, when compared to the wild-type receptor (Bne1 WT) ( Figure 5A). In effect, the replacement  3). The results were normalized to the control condition (DMSO without ligand). Significant differences between the tested concentrations of TBT and the solvent control were inferred with one-way ANOVA. The F and P values of the source of variation are indicated below each species name abbreviation. Asterisks denote significant differences (* P ≤ 0.05, ** P ≤ 0.01, *** P ≤ 0.001). Significant differences between BneRXR1 and BneRXR2 with TBT at the same concentration were inferred with Student's t-test. The green rectangles indicate the RXRs which have the cysteine residue previously described as crucial for binding to TBT. Hsa, H. sapiens; Mde, M. detruncata; Bvi, B. viridis; Pcal, P. californica; Bne, B. neritina, Pcau, P. caudatus; Xbo, X. bocki; Tad, T. adhaerens; Aau, A. aurita.

Mutagenesis of Bryozoan Specific RXR Orthologs
To further clarify the transcriptional responses obtained with the bryozoan B. neritina RXRs, with both 9cisRA and TBT, site-directed and deletion mutants were tested ( Figure 5). First, the bryozoan RXR1 Thr388Cys mutant (Bne1 T > C) was assayed for its ability to elicit a stronger response to TBT, when compared to the wild-type receptor (Bne1 WT) ( Figure 5A). In effect, the replacement of the threonine residue (Thr) by the classical key cysteine residue (Cys) prompted a significant higher response to all the tested concentrations of TBT, without affecting the response to 9cisRA. Next, using the bryozoan RXR2, three mutants were produced: the RXR2 Met231Ala (Bne2 M > A) and RXR2 Val294Ala (Bne2 V > A) mutants were tested expecting higher responses to 9cisRA, since the stereochemistry of Bne2M > A and Bne2V > A mutants was expected to be more similar to the one of RXR1, while the RXR2 ∆Ser211-Asp227 (Bne2 del) mutant was tested expecting a favorable conformational change in the receptor binding pocket, affecting interactions to both 9cisRA and TBT ( Figure 5B). A slightly higher sensitivity to 9cisRA was detected for Bne2 M > A, Bne2 V > A, and Bne2 del mutants, however, the increase was not significant (P > 0.05). Regarding the interaction with TBT, the Bne2 del mutant presented the highest levels of activation, being indicative that, apart from the cysteine residue, the whole structural conformation of the LDB can affect the responsiveness to TBT. of the threonine residue (Thr) by the classical key cysteine residue (Cys) prompted a significant higher response to all the tested concentrations of TBT, without affecting the response to 9cisRA. Next, using the bryozoan RXR2, three mutants were produced: the RXR2 Met231Ala (Bne2 M > A) and RXR2 Val294Ala (Bne2 V > A) mutants were tested expecting higher responses to 9cisRA, since the stereochemistry of Bne2M > A and Bne2V > A mutants was expected to be more similar to the one of RXR1, while the RXR2 ∆Ser211-Asp227 (Bne2 del) mutant was tested expecting a favorable conformational change in the receptor binding pocket, affecting interactions to both 9cisRA and TBT ( Figure 5B). A slightly higher sensitivity to 9cisRA was detected for Bne2 M > A, Bne2 V > A, and Bne2 del mutants, however, the increase was not significant (P > 0.05). Regarding the interaction with TBT, the Bne2 del mutant presented the highest levels of activation, being indicative that, apart from the cysteine residue, the whole structural conformation of the LDB can affect the responsiveness to TBT. Data represent means ± SEM from three separate experiments (n = 3). The results were normalized to the control condition (DMSO without ligand). In each condition, the significant differences between B. neritina RXR1 wild-type and mutant and between B. neritina RXR2 wild-type and mutants were inferred with Student's t-test and one-way ANOVA, respectively. The F and P values of the source of variation for each ANOVA analysis are indicated below each group condition. Asterisks denote significant differences (* P ≤ 0.05, ** P ≤ 0.01). Data represent means ± SEM from three separate experiments (n = 3). The results were normalized to the control condition (DMSO without ligand). In each condition, the significant differences between B. neritina RXR1 wild-type and mutant and between B. neritina RXR2 wild-type and mutants were inferred with Student's t-test and one-way ANOVA, respectively. The F and P values of the source of variation for each ANOVA analysis are indicated below each group condition. Asterisks denote significant differences (* P ≤ 0.05, ** P ≤ 0.01).

Discussion
In this work, we successfully isolated full-length RXR gene orthologs in all of the examined metazoan species. These covered an ample set of phyla, configuring the conclusion that RXR is largely conserved throughout Metazoa and likely represents an essential gene for animal homeostasis and development [10,11]. Curiously, we were also able to deduce the occurrence of two RXR genes in bryozoan species (B. neritina and M. membranacea). Multiple RXR isoforms were previously identified in other invertebrate species such as gastropod molluscs [25,43], the platyhelminth parasite Schistosoma mansoni [57,58], and the bivalve Chlamys farreri [59]; yet, the physiological consequences of this novelty remain unknown and should be investigated in the future. The inspection of the LBD amino acid alignment from the isolated RXRs puts into perspective the possibility of these RXRs to be activated upon binding to 9cisRA, despite the substitution of some residues reported as important for the interaction with 9cisRA [51][52][53]. In fact, some studies demonstrated that the strict conservation of the residues lining the ligand-binding pocket is not sufficient to guarantee ligand-receptor interaction, as demonstrated by the insect Tribolium castaneum USP, which is strongly conserved but does not bind to 9cisRA [31]. On the other hand, amino acid substitutions within the RXR ligand-binding pocket from the placozoan T. adhaerens or the cnidarian T. cystophora, closely related to A. aurita, did not impair 9cisRA binding [14,47,51]. In the analyzed species, 9cisRA was able to induce RXR-mediated transcription, albeit with different effectiveness. With the exception of the brachiopod M. detruncata RXR, all examined RXRs yielded weaker responses when compared with the human homologue. These results are in accordance with previous studies of RXRs from the cephalochordate B. floridae [29], the molluscs Biomphalaria glabrata [25] and Reishia clavigera [60], and the annelid Platynereis dumerilii [27]. Specifically, in bryozoans both RXRs presented different sensitivities towards 9cisRA, similarly to the alternative splicing isoforms found in gastropod molluscs [43,60]. This finding might indicate that the two RXRs of this lineage should have different endogenous ligands and/or different physiological functions. Overall, these results further clarify the previously proposed evolutionary history of RXR, in which an ancestral function of RXR able to bind 9cisRA was improved in vertebrates and secondarily lost at the base of the insects [29]. Yet, regarding insects, and given that the earlier branching orthopteran L. migratoria [61] USP was suggested to bind with high affinity to 9cisRA and ATRA [32], unlike other analyzed insect USPs [30,31], a more detailed phylogenetic functional mapping should be performed to resolve episodes of loss and gain of function towards 9cisRA. Still, the present results reinforce the hypothesis of an ancestral RXR with low affinity for 9cisRA and, in addition to vertebrates, we propose an increase of affinity for 9cisRA in Brachiopoda, Phoronida, Bryozoa, Annelida, and Placozoa ( Figure 6).
In addition to binding to 9cisRA, RXR has been described as a target of endocrine disruption mediated by organotins (TBT, TPT). The induction of imposex in marine snails is a classic example of the RXR-dependent disruption by organotins at nanomolar levels [24,26,60,[62][63][64]. Furthermore, the evolutionary conservation of RXR-dependent modulation by organotins was suggested among lophotrochozoans, vertebrates, and crustaceans [27,33,65]. Recently, an RXR-dependent signaling pathway was suggested to promote the changes in the fatty acid profiles of the rotifer B. koreanus exposed to TBT [40]. The Organisation for Economic Co-operation and Development (OECD) suggests in vitro assays as a level 2 approach to test and assess endocrine disruption by environmental contaminants and denotes NRs as key players in this framework [17]. As follows, we addressed if TBT was able to activate the isolated RXRs in order to verify if RXR exploitation by TBT was also conserved during metazoan evolution. Our results demonstrated that the interaction of TBT with RXR is conserved among Metazoa, with the exception of Placozoa (T. adhaerens) ( Figure 6). Moreover, we found that the condition described in rotifers [40] is not unique, since the cnidarian A. aurita RXR lacks the key cysteine and was nonetheless activated by TBT. The bryozoan B. neritina RXRs, on the other hand, yielded puzzling results, as similar responses were found independently of the conservation status of the cysteine residue. However, the substitution of the threonine residue by a cysteine led to a drastic increase of B. neritina RXR1 responsiveness. Regarding B. neritina RXR2, the deletion of the amino acid insertion found within the ligand-binding pocket also improved the receptor's ability to bind TBT, albeit with lower potency than the RXR1 Thr388Cys mutant. These observations suggest that other ligand-pocket interactions could accommodate TBT binding, given that some pocket environments seem favorable to TBT binding even in the absence of the key cysteine residue. Yet, although not always strictly required, in more structurally conserved pockets, the substitution by a cysteine residue seems to enhance this interaction. Altogether, the analyses presented here indicate that the structural remodeling of RXR LBD during evolution favored the interaction of TBT with the cysteine residue, promoting a stronger sensitivity to this organotin. In addition to binding to 9cisRA, RXR has been described as a target of endocrine disruption mediated by organotins (TBT, TPT). The induction of imposex in marine snails is a classic example of the RXR-dependent disruption by organotins at nanomolar levels [24,26,60,62,63,64]. Furthermore, the evolutionary conservation of RXR-dependent modulation by organotins was suggested among lophotrochozoans, vertebrates, and crustaceans [27,33,65]. Recently, an RXR-dependent signaling pathway was suggested to promote the changes in the fatty acid profiles of the rotifer B. koreanus exposed to TBT [40]. The Organisation for Economic Co-operation and Development (OECD) suggests in vitro assays as a level 2 approach to test and assess endocrine disruption by environmental contaminants and denotes NRs as key players in this framework [17]. As follows, we addressed if TBT was able to activate the isolated RXRs in order to verify if RXR exploitation by TBT was also conserved during metazoan evolution. Our results demonstrated that the interaction of TBT with RXR is conserved among Metazoa, with the exception of Placozoa (T. adhaerens) ( Figure 6). Moreover, we found that the condition described in rotifers [40] is not unique, since the cnidarian A. aurita RXR lacks the key cysteine and was nonetheless activated by TBT. The bryozoan B. neritina RXRs, on the other hand, yielded puzzling results, as similar responses were found independently of the conservation status of the cysteine residue. However, the substitution of the threonine residue by a cysteine led to a drastic increase of B. neritina RXR1 responsiveness. Regarding B. neritina RXR2, Figure 6. Evolutionary conservation of RXR modulation by 9cisRA and tributyltin. * The ultraspiracle (USP) ortholog from the earlier branching orthopteran Locusta migratoria retains high affinity towards 9cisRA and all-trans retinoic acid.

Conclusions
In this study, we screened multiple metazoan genomes for the occurrence of RXR. This NR is present in all of the examined species, although with some important differences. We found that the ability of 9cisRA and TBT to transactivate RXR is conserved in more phyla than previously described. Moreover, we demonstrated that receptors with highly conserved sequences have different responses to ligand binding, highlighting the importance of studies with a broader species sampling. Yet, we were able to identify more cases of unconventional modulation of RXR by the organotin TBT, emphasizing the conservation of a common mechanism of endocrine disruption by TBT in Metazoa.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2218-273X/10/4/594/s1, Table S1: List of primers used to isolate RXR genes, Table S2: List of primers used to amplify RXR hinge and LBD regions to be cloned into pBIND expression vector, Table S3: Ratio between firefly and Renilla luciferase luminescent values.