Genome inventory and analysis of nuclear hormone receptors in Tetraodon nigroviridis

Agonistic and potentiating effects of perfluoroalkyl substances (PFAS) on the Atlantic cod (Gadus morhua) peroxisome proliferator-activated receptors (Ppars)

Toxicity mediated by per- and polyfluoroalkyl substances (PFAS), and especially perfluoroalkyl acids (PFAAs), has been linked to activation of peroxisome proliferator-activated receptors (Ppar) in many vertebrates. Here, we present the primary structures, phylogeny, and tissue-specific distributions of the Atlantic cod (Gadus morhua) gmPpara1, gmPpara2, gmPparb, and gmPparg, and demonstrate that the carboxylic acids PFHxA, PFOA, PFNA, as well as the sulfonic acid PFHxS, activate gmPpara1 in vitro, which was also supported by in silico analyses. Intriguingly, a binary mixture of PFOA and the non-activating PFOS produced a higher activation of gmPpara1 compared to PFOA alone, suggesting that PFOS has a potentiating effect on receptor activation. Supporting the experimental data, docking and molecular dynamics simulations of single and double-ligand complexes led to the identification of a putative allosteric binding site, which upon binding of PFOS stabilizes an active conformation of gmPpara1. Notably, binary exposures of gmPpara1, gmPpara2, and gmPparb to model-agonists and PFAAs produced similar potentiating effects. This study provides novel mechanistic insights into how PFAAs may modulate the Ppar signaling pathway by either binding the canonical ligand-binding pocket or by interacting with an allosteric binding site. Thus, individual PFAAs, or mixtures, could potentially modulate the Ppar-signaling pathway in Atlantic cod by interfering with at least one gmPpar subtype.

Genome-wide identification, evolution of DNA methyltransferases and their expression during gonadal development in Nile tilapia

DNA methyltransferases (dnmts) are responsible for DNA methylation and play important roles in organism development. In this study, seven dnmts genes (dnmt1, dnmt2, dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1, dnmt3bb.2) were identified in Nile tilapia. Comprehensive analyses of dnmts were performed using available genome databases from representative animal species. Phylogenetic analysis revealed that the dnmts family were highly conserved in teleosts. Based on transcriptome data from eight adult tilapia tissues, the dnmts were found to be dominantly expressed in the head kidney, testis and ovary. Analyses of the gonadal transcriptome data in different developmental stages revealed that all dnmts were expressed in both ovary and testis, and four de novo dnmts (dnmt3aa, dnmt3ab, dnmt3bb.1, dnmt3bb.2) showed higher expression in the testis than in the ovary. Furthermore, during sex reversal induced by Fadrozole, the expression of these four de novo dnmts increased significantly in treated group compared to female control group. By in situ hybridization, the seven dnmts were found to be expressed mainly in phase I and II oocytes of the ovary and spermatocytes of the testis. When gonads were incubated with a methyltransferase inhibitor (5-AzaCdR) in vitro, the expression of dnmts genes were down-regulated significantly, while the expression of cyp19a1a (a key gene in female pathway) and dmrt1 (a key gene in male pathway) increased significantly. Our results revealed the conservation of dnmts during evolution and indicated a potential role of dnmts in epigenetic regulation of gonadal development.

Annotation of the Nuclear Receptors in an Estuarine Fish species, Fundulus heteroclitus

The nuclear receptors (NRs) are ligand-dependent transcription factors that respond to various internal as well as external cues such as nutrients, pheromones, and steroid hormones that play crucial roles in regulation and maintenance of homeostasis and orchestrating the physiological and stress responses of an organism. We annotated the Fundulus heteroclitus (mummichog; Atlantic killifish) nuclear receptors. Mummichog are a non-migratory, estuarine fish with a limited home range often used in environmental research as a field model for studying ecological and evolutionary responses to variable environmental conditions such as salinity, oxygen, temperature, pH, and toxic compounds because of their hardiness. F. heteroclitus have at least 74 NRs spanning all seven gene subfamilies. F. heteroclitus is unique in that no RXRα member was found within the genome. Interestingly, some of the NRs are highly conserved between species, while others show a higher degree of divergence such as PXR, SF1, and ARα. Fundulus like other fish species show expansion of the RAR (NR1B), Rev-erb (NR1D), ROR (NR1F), COUPTF (NR2F), ERR (NR3B), RXR (NR2B), and to a lesser extent the NGF (NR4A), and NR3C steroid receptors (GR/AR). Of particular interest is the co-expansion of opposing NRs, Reverb-ROR, and RAR/RXR-COUPTF.

Genome-wide identification, evolution of chromobox family genes and their expression in Nile tilapia

Chromobox (Cbx) family proteins are transcriptional repressors that involved in epigenetic and developmental processes. In this study, comprehensive analyses of Cbxs were performed using available genome databases from representative animal species. The Cbx family were originated from one Polycomb (Pc) gene like the yeast Pc, which duplicated into two and gave rise to the Pc and the Heterochromatin protein 1 (Hp1) identified in invertebrates from protozoon to lancelet. Rapid expansion of Cbx family members was observed in vertebrates as ~8 (5 Pc and 3 Hp1) were identified in spotted gar, coelacanth and tetrapods. Further expansion of the members to ~14 (9 Pc and 5 Hp1) was observed in teleosts due to the third round genome duplication (3R). Based on transcriptome data from eight adult tilapia tissues, most of the Cbxs were found to be dominantly expressed in the brain, testis, ovary and heart. Analyses of the gonadal transcriptome data from four developmental stages revealed that all Cbxs were expressed in both ovary and testis except Cbx7b, with significant increase of the total and average RPKM from 5 to 90dah (days after hatching). By in situ hybridization, the three most highly and sexual dimorphically expressed Cbx genes in gonads, Cbx1b, Cbx3a and Cbx5, were found to be expressed in phase I and II oocytes of the ovary, and in secondary spermatocytes (Cbx1b and Cbx3a) and spermatids (Cbx5) of the testis. Our results revealed the evolution of Cbx genes and indicated a potential role of Cbxs in epigenetic regulation of gametogenesis.

Genome-wide identification, evolution and expression analysis of nuclear receptor superfamily in Nile tilapia, Oreochromis niloticus

The nuclear receptor (NR) superfamily, which is divided into 7 subfamilies, constitutes one of the largest classes of transcription factors. In this study, through comprehensive database search, we identified all NRs (including 4 novel members) from the tilapia (75), common carp (137), zebrafish (73), fugu (73), tetraodon (72), stickleback (70), medaka (69), coelacanth (55), spotted gar (51) and elephant shark (50). For 21 NRs, two duplicates were found in teleosts, while only one in tetrapods. These duplicates, except those of DAX1, SHP and GCNF found in the elephant shark, were derived from 3R (third round of genome duplication). The linkage duplication of 5 syntenic blocks (comprising 14 duplicated NR couples) in teleosts further supported their 3R origin. Based on transcriptome data from adult tilapia, 53 NRs were found to be expressed in more than one tissue (brain, head kidney, heart, liver, kidney, muscle, ovary and testis), and 4 were tissue-specific, indicating their essential roles in the corresponding tissue. Based on the XX and XY gonadal transcriptome data from four developmental stages, 65 NRs were detected in gonads, with 21, 31, 11 and 29 expressed sexual dimorphically at 5, 30, 90 and 180days after hatching, respectively. The expression of four selected genes was examined by in situ hybridization (ISH) and quantitative PCR (qPCR) to validate the spatial and temporal expression profiles of NRs. Comparative analyses of the expression profiles of duplicated NRs revealed divergence in gene expression as well as gene function. Our results demonstrated that NRs may play important roles in sex determination and gonadal development in teleosts.

Molecular cloning, functional characterization, and evolutionary analysis of vitamin D receptors isolated from basal vertebrates

The vertebrate genome is a result of two rapid and successive rounds of whole genome duplication, referred to as 1R and 2R. Furthermore, teleost fish have undergone a third whole genome duplication (3R) specific to their lineage, resulting in the retention of multiple gene paralogs. The more recent 3R event in teleosts provides a unique opportunity to gain insight into how genes evolve through specific evolutionary processes. In this study we compare molecular activities of vitamin D receptors (VDR) from basal species that diverged at key points in vertebrate evolution in order to infer derived and ancestral VDR functions of teleost paralogs. Species include the sea lamprey (Petromyzon marinus), a 1R jawless fish; the little skate (Leucoraja erinacea), a cartilaginous fish that diverged after the 2R event; and the Senegal bichir (Polypterus senegalus), a primitive 2R ray-finned fish. Saturation binding assays and gel mobility shift assays demonstrate high affinity ligand binding and classic DNA binding characteristics of VDR has been conserved across vertebrate evolution. Concentration response curves in transient transfection assays reveal EC50 values in the low nanomolar range, however maximum transactivational efficacy varies significantly between receptor orthologs. Protein-protein interactions were investigated using co-transfection, mammalian 2-hybrid assays, and mutations of coregulator activation domains. We then combined these results with our previous study of VDR paralogs from 3R teleosts into a bioinformatics analysis. Our results suggest that 1, 25D3 acts as a partial agonist in basal species. Furthermore, our bioinformatics analysis suggests that functional differences between VDR orthologs and paralogs are influenced by differential protein interactions with essential coregulator proteins. We speculate that we may be observing a change in the pharmacodynamics relationship between VDR and 1, 25D3 throughout vertebrate evolution that may have been driven by changes in protein-protein interactions between VDR and essential coregulators.

Functional diversification of vitamin D receptor paralogs in teleost fish after a whole genome duplication event

The diversity and success of teleost fishes (Actinopterygii) has been attributed to three successive rounds of whole-genome duplication (WGD). WGDs provide a source of raw genetic material for evolutionary forces to act upon, resulting in the divergence of genes with altered or novel functions. The retention of multiple gene pairs (paralogs) in teleosts provides a unique opportunity to study how genes diversify and evolve after a WGD. This study examines the hypothesis that vitamin D receptor (VDR) paralogs (VDRα and VDRβ) from two distantly related teleost orders have undergone functional divergence subsequent to the teleost-specific WGD. VDRα and VDRβ paralogs were cloned from the Japanese medaka (Beloniformes) and the zebrafish (Cypriniformes). Initial transactivation studies using 1α, 25-dihydroxyvitamin D3 revealed that although VDRα and VDRβ maintain similar ligand potency, the maximum efficacy of VDRβ was significantly attenuated compared with VDRα in both species. Subsequent analyses revealed that VDRα and VDRβ maintain highly similar ligand affinities; however, VDRα demonstrated preferential DNA binding compared with VDRβ. Protein-protein interactions between the VDR paralogs and essential nuclear receptor coactivators were investigated using transactivation and mammalian two-hybrid assays. Our results imply that functional differences between VDRα and VDRβ occurred early in teleost evolution because they are conserved between distantly related species. Our results further suggest that the observed differences may be associated with differential protein-protein interactions between the VDR paralogs and coactivators. We speculate that the observed functional differences are due to subtle ligand-induced conformational differences between the two paralogs, leading to divergent downstream functions.

Evolution of promiscuous nuclear hormone receptors: LXR, FXR, VDR, PXR, and CAR

Nuclear hormone receptors (NHRs) are transcription factors that work in concert with co-activators and co-repressors to regulate gene expression. Some examples of ligands for NHRs include endogenous compounds such as bile acids, retinoids, steroid hormones, thyroid hormone, and vitamin D. This review describes the evolution of liver X receptors α and β (NR1H3 and 1H2, respectively), farnesoid X receptor (NR1H4), vitamin D receptor (NR1I1), pregnane X receptor (NR1I2), and constitutive androstane receptor (NR1I3). These NHRs participate in complex, overlapping transcriptional regulation networks involving cholesterol homeostasis and energy metabolism. Some of these receptors, particularly PXR and CAR, are promiscuous with respect to the structurally wide range of ligands that act as agonists. A combination of functional and computational analyses has shed light on the evolutionary changes of NR1H and NR1I receptors across vertebrates, and how these receptors may have diverged from ancestral receptors that first appeared in invertebrates.

Two farnesoid X receptor alpha isoforms in Japanese medaka (Oryzias latipes) are differentially activated in vitro

The nuclear receptor farnesoid X receptor alpha (FXRalpha, NR1H4) is activated by bile acids in multiple species including mouse, rat, and human and in this study we have identified two isoforms of Fxralpha in Japanese medaka (Oryzias latipes), a small freshwater teleost. Both isoforms share a high amino acid sequence identity to mammalian FXRalpha (approximately 70% in the ligand-binding domain). Fxralpha1 and Fxralpha2 differ within the AF1 domain due to alternative splicing at the fourth intron-exon boundary. This process results in Fxralpha1 having an extended N-terminus compared to Fxralpha2. A Gal4DBD-FxralphaLBD fusion construct was activated by chenodeoxycholic, cholic, deoxycholic and lithocholic acids, and the synthetic agonist GW4064 in transient transactivation assays. Activation of the Gal4DBD-FxralphaLBD fusion construct was enhanced by addition of PGC-1alpha, as demonstrated through titration assays. Surprisingly, when the full-length versions of the two Fxralpha isoforms were compared in transient transfection assays, Fxralpha2 was activated by C(24) bile acids and GW4064, while Fxralpha1 was not significantly activated by any of the compounds tested. Since the only significant difference between the full-length constructs was sequence in the AF1 domain, these experiments highlight a key functional region in the Fxralpha AF1 domain. Furthermore, mammalian two-hybrid studies demonstrated the ability of Fxralpha2, but not Fxralpha1, to interact with PGC-1alpha and SRC-1, and supported our results from the transient transfection reporter gene activation assays. These data demonstrate that both mammalian and teleost FXR (Fxralpha2 isoform) are activated by primary and secondary bile acids.

Cloning and expression pattern of peroxisome proliferator-activated receptors, estrogen receptor alpha and retinoid X receptor alpha in the thicklip grey mullet Chelon labrosus

Aquatic organisms are exposed to diverse xenobiotics that cause peroxisome proliferation and/or endocrine disruption, both modulated in vertebrates by transcription factors of the nuclear receptor (NR) superfamily. Peroxisome proliferators are agonists of peroxisome proliferator-activated receptors (PPARs) that heterodimerize with the retinoid X receptor (RXR). Many xenoestrogens activate the estrogen receptor (ER). Here, 1090 bp of PPARalpha, 1255 bp of PPARgamma, 278 bp of RXRalpha, and 578 bp of ERalpha of thicklip grey mullet Chelon labrosus were cloned. Sequences were highly conserved, although relevant changes with respect to mammalian homologs were identified in PPARgamma and ERalpha. Semi-quantitative RT-PCR was used to determine if these NRs were expressed in different tissues of male, female and undifferentiated mullets captured in January and June. Expression of PPARs was highest in liver and lowest in muscle. RXRalpha expression was homogeneous excepting a low expression in male and female gill in January and brain and heart of undifferentiated fish in January and June. ERalpha expression predominated in liver and female gonad in June. The expression level of PPARs and ERalpha was significantly higher in liver in January than in gills in January or June. The present results show tissue-dependent modulation of expression of NRs in mullets.

Paralogous vitamin D receptors in teleosts: transition of nuclear receptor function

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we describe the cloning and functional characterization of two novel vitamin D receptor (VDR) paralogs from the freshwater teleost medaka (Oryzias latipes). VDR sequences were identified through mining of the medaka genome database in which gene organization and structure was determined. Two distinct VDR genes were identified in the medaka genome and mapped to defined loci. Each VDR sequence exhibits unique intronic organization and dissimilar 5' untranslated regions, suggesting they are not isoforms of the same gene locus. Phylogenetic comparison with additional teleosts and mammalian VDR sequences illustrate that two distinct clusters are formed separating aquatic and terrestrial species. Nested within the teleost cluster are two separate clades for VDRalpha and VDRbeta. The topology of teleost VDR sequences is consistent with the notion of paralogous genes arising from a whole genome duplication event prior to teleost radiation. Functional characterization was conducted through the development of VDR expression vectors including Gal4 chimeras containing the yeast Gal4 DNA binding domain fused to the medaka VDR ligand binding domain and full-length protein. The common VDR ligand 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)(2)D(3)] resulted in significant transactivation activity with both the Gal4 and full-length constructs of medaka (m) VDRbeta. Comparatively, transactivation of mVDRalpha with 1alpha,25(OH)(2)D(3) was highly attenuated, suggesting a functional divergence between these two nuclear receptor paralogs. We additionally demonstrate through coactivator studies that mVDRalpha is still functional; however, it exhibits a different sensitivity to 1alpha,25(OH)(2)D(3), compared with VDRbeta. These results suggest that in mVDRalpha and VDRbeta have undergone a functional divergence through a process of sub- and/or neofunctionalization of VDR nuclear receptor gene pairs.

Evolution of pharmacologic specificity in the pregnane X receptor

Background

The pregnane X receptor (PXR) shows the highest degree of cross-species sequence diversity of any of the vertebrate nuclear hormone receptors. In this study, we determined the pharmacophores for activation of human, mouse, rat, rabbit, chicken, and zebrafish PXRs, using a common set of sixteen ligands. In addition, we compared in detail the selectivity of human and zebrafish PXRs for steroidal compounds and xenobiotics. The ligand activation properties of the Western clawed frog (Xenopus tropicalis) PXR and that of a putative vitamin D receptor (VDR)/PXR cloned in this study from the chordate invertebrate sea squirt (Ciona intestinalis) were also investigated.

Results

Using a common set of ligands, human, mouse, and rat PXRs share structurally similar pharmacophores consisting of hydrophobic features and widely spaced excluded volumes indicative of large binding pockets. Zebrafish PXR has the most sterically constrained pharmacophore of the PXRs analyzed, suggesting a smaller ligand-binding pocket than the other PXRs. Chicken PXR possesses a symmetrical pharmacophore with four hydrophobes, a hydrogen bond acceptor, as well as excluded volumes. Comparison of human and zebrafish PXRs for a wide range of possible activators revealed that zebrafish PXR is activated by a subset of human PXR agonists. The Ciona VDR/PXR showed low sequence identity to vertebrate VDRs and PXRs in the ligand-binding domain and was preferentially activated by planar xenobiotics including 6-formylindolo-[3,2-b]carbazole. Lastly, the Western clawed frog (Xenopus tropicalis) PXR was insensitive to vitamins and steroidal compounds and was activated only by benzoates.

Conclusion

In contrast to other nuclear hormone receptors, PXRs show significant differences in ligand specificity across species. By pharmacophore analysis, certain PXRs share similar features such as human, mouse, and rat PXRs, suggesting overlap of function and perhaps common evolutionary forces. The Western clawed frog PXR, like that described for African clawed frog PXRs, has diverged considerably in ligand selectivity from fish, bird, and mammalian PXRs.