The cytochrome P450 (CYP) superfamily is a diverse and important enzyme family, playing a central role in chemical defense and in synthesis and metabolism of major biological signaling molecules. The CYPomes of four cnidarian genomes (Hydra vulgaris, Acropora digitifera, Aurelia aurita, Nematostella vectensis) were annotated; phylogenetic analyses determined the evolutionary relationships amongst the sequences and with existing metazoan CYPs. 155 functional CYPs were identified and 90 fragments. Genes were from 24 new CYP families and several new subfamilies; genes were in 9 of the 12 established metazoan CYP clans. All species had large expansions of clan 2 diversity, with H. vulgaris having reduced diversity for both clan 3 and mitochondrial clan. We identified potential candidates for xenobiotic metabolism and steroidogenesis. That each genome contained multiple, novel CYP families may reflect the large evolutionary distance within the cnidarians, unique physiology in the cnidarian classes, and/or different ecology of the individual species.

Background

Many persistent organic pollutants (POPs) accumulate readily in polar bears because of their position as apex predators in Arctic food webs. The pregnane X receptor (PXR, formally NR1I2, here proposed to be named promiscuous xenobiotic receptor) is a xenobiotic sensor that is directly involved in metabolizing pathways of a wide range of environmental contaminants.

Objectives

In the present study, we comparably assess the ability of 51 selected pharmaceuticals, pesticides and emerging contaminants to activate PXRs from polar bears and humans using an in vitro luciferase reporter gene assay.

Results

We found that polar bear PXR is activated by a wide range of our test compounds (68%) but has a slightly more narrow ligand specificity than human PXR that was activated by 86% of the 51 test compounds. The majority of the agonists identified (70%) produces a stronger induction of the reporter gene via human PXR than via polar bear PXR, however with some notable and environmentally relevant exceptions.

Conclusions

Due to the observed differences in activation of polar bear and human PXRs, exposure of each species to environmental agents is likely to induce biotransformation differently in the two species. Bioinformatics analyses and structural modeling studies suggest that amino acids that are not part of the ligand-binding domain and do not interact with the ligand can modulate receptor activation.

Regulators of biotransformation are of particular interest in pharmacology and toxicology, determining in part the metabolism, disposition, and toxicity of chemicals. The nuclear receptor NR1I2 (pregnane X receptor, PXR) is a prominent xenosensor that regulates the expression of biotransformation enzymes governing elimination of many exogenous as well as endogenous compounds. Zebrafish (Danio rerio) has only one gene locus for pxr, but different genetic variants have been identified in zebrafish. However, the prevalence and significance of these variants are unknown. We hypothesize that sequence variation occurring in the Pxr gene of zebrafish may affect the action and fate of many chemicals in this species, a key model organism in various fields of research, including environmental toxicology. Here, we examine variation in Pxr sequences from four different strains of zebrafish and assess the responses of each Pxr to clotrimazole and butyl-4-aminobenzoate. The Pxr variants differed in both their ability to bind these structurally different ligands and to regulate reporter gene expression in vitro. We infer that the observed sequence variations in zebrafish Pxrs likely affect the response to putative Pxr agonists in vivo and potentially cause strain-specific biotransformation of xenobiotics in zebrafish. Thus, the choice of zebrafish strain could affect the outcome of downstream toxicological studies.

Atlantic killifish (Fundulus heteroclitus) thrive in New Bedford Harbor (NBH), MA, highly contaminated with polychlorinated biphenyls (PCBs). Resident killifish have evolved tolerance to dioxin-like (DL) PCBs, whose toxic effects through the aryl hydrocarbon receptor (AhR) are well studied. In NBH, non-dioxin like PCBs (NDL PCBs), which lack activity toward the AhR, vastly exceed levels of DL congeners yet how killifish counter NDL toxic effects has not been explored. In mammals and fish, NDL PCBs are potent activators of ryanodine receptors (RyR), Ca(2+) release channels necessary for a vast array of physiological processes. In the current study we compared the expression and function of RyR related pathways in NBH killifish with killifish from the reference site at Scorton Creek (SC, MA). Relative to the SC fish, adults from NBH displayed increased levels of skeletal muscle RyR1 protein, and increased levels of FK506-binding protein 12 kDa (FKBP12) an accessory protein essential for NDL PCB-triggered changes in RyR channel function. In accordance with increased RyR1 levels, NBH killifish displayed increased maximal ligand binding, increased maximal response to Ca(2+) activation and increased maximal response to activation by the NDL PCB congener PCB 95. Compared to SC, NBH embryos and larvae had increased levels of mtor and ryr2 transcripts at multiple stages of development, and generations, while levels of serca2 were decreased at 9 days post-fertilization in the F1 and F2 generations. These findings suggest that there are compensatory and heritable changes in RyR mediated Ca(2+) signaling proteins or potential signaling partners in NBH killifish.

Some species of butterflyfish have had preyed upon corals for millions of years, yet the mechanism of butterflyfish specialized coral feeding strategy remains poorly understood. Certain butterflyfish have the ability to feed on allelochemically rich soft corals, e.g. Sinularia maxima. Cytochrome P450 (CYP) is the predominant enzyme system responsible for the detoxification of dietary allelochemicals. CYP2-like and CYP3A-like content have been associated with butterflyfish that preferentially consumes allelochemically rich soft corals. To investigate the role of butterflyfish CYP2 and CYP3A enzymes in dietary preference, we conducted oral feeding experiments using homogenates of S. maxima and a toxin isolated from the coral in four species of butterflyfish with different feeding strategies. After oral exposure to the S. maxima toxin 5-episinulaptolide (5ESL), which is not normally encountered in the Hawaiian butterflyfish diet, an endemic specialist, Chaetodon multicinctus experienced 100% mortality compared to a generalist, Chaetodon auriga, which had significantly more (3-6 fold higher) CYP3A-like basal content and catalytic activity. The specialist, Chaetodon unimaculatus, which preferentially feed on S. maxima in Guam, but not in Hawaii, had 100% survival, a significant induction of 8-12 fold CYP3A-like content, and an increased ability (2-fold) to metabolize 5ESL over other species. Computer modeling data of CYP3A4 with 5ESL were consistent with microsomal transformation of 5ESL to a C15-16 epoxide from livers of C. unimaculatus. Epoxide formation correlated with CYP3A-like content, catalytic activity, induction, and NADPH-dependent metabolism of 5ESL. These results suggest a potentially important role for the CYP3A family in butterflyfish-coral diet selection through allelochemical detoxification.

Some species of butterflyfish have had preyed upon corals for millions of years, yet the mechanism of butterflyfish specialized coral feeding strategy remains poorly understood. Certain butterflyfish have the ability to feed on allelochemically rich soft corals, e.g. Sinularia maxima. Cytochrome P450 (CYP) is the predominant enzyme system responsible for the detoxification of dietary allelochemicals. CYP2-like and CYP3A-like content have been associated with butterflyfish that preferentially consumes allelochemically rich soft corals. To investigate the role of butterflyfish CYP2 and CYP3A enzymes in dietary preference, we conducted oral feeding experiments using homogenates of S. maxima and a toxin isolated from the coral in four species of butterflyfish with different feeding strategies. After oral exposure to the S. maxima toxin 5-episinulaptolide (5ESL), which is not normally encountered in the Hawaiian butterflyfish diet, an endemic specialist, Chaetodon multicinctus experienced 100% mortality compared to a generalist, Chaetodon auriga, which had significantly more (3-6 fold higher) CYP3A-like basal content and catalytic activity. The specialist, Chaetodon unimaculatus, which preferentially feed on S. maxima in Guam, but not in Hawaii, had 100% survival, a significant induction of 8-12 fold CYP3A-like content, and an increased ability (2-fold) to metabolize 5ESL over other species. Computer modeling data of CYP3A4 with 5ESL were consistent with microsomal transformation of 5ESL to a C15-16 epoxide from livers of C. unimaculatus. Epoxide formation correlated with CYP3A-like content, catalytic activity, induction, and NADPH-dependent metabolism of 5ESL. These results suggest a potentially important role for the CYP3A family in butterflyfish-coral diet selection through allelochemical detoxification.

Non-dioxin-like polychlorinated biphenyls (NDL PCBs) activate ryanodine receptors (RyR), microsomal Ca²⁺ channels of broad significance. Teleost fish may be important models for NDL PCB neurotoxicity, and we used sequencing databases to characterize teleost RyR and FK506 binding protein 12 or 12.6kDa (genes FKBP1A; FKBP1B), which promote NDL PCB-triggered Ca²⁺ dysregulation. Particular focus was placed on describing genes in the Atlantic killifish (Fundulus heteroclitus) genome and searching available RNA-sequencing datasets for single nucleotide variants (SNV) between PCB tolerant killifish from New Bedford Harbor (NBH) versus sensitive killifish from Scorton Creek (SC), MA. Consistent with the teleost whole genome duplication (tWGD), killifish have six RyR genes, corresponding to a and b paralogs of mammalian RyR1, 2 and 3. The presence of six RyR genes was consistent in all teleosts investigated including zebrafish. Killifish have four FKBP1; one FKBP1b and three FKBP1a named FKBP1aa, FKBP1ab, likely from the tWGD and a single gene duplicate FKBP1a3 suggested to have arisen in Atherinomorphae. The RyR and FKBP1 genes displayed tissue and developmental stage-specific mRNA expression, and the previously uncharacterized RyR3, herein named RyR3b, and all FKBP1 genes were prominent in brain. We identified a SNV in RyR3b encoding missense mutation E1458D. In NBH killifish, 57% were heterozygous and 28% were homozygous for this SNV, whereas almost all SC killifish (94%) lacked the variant (n≥39 per population). The outlined sequence differences between mammalian and teleost RyR and FKBP1 together with outlined population differences in SNV frequency may contribute to our understanding of NDL PCB neurotoxicity.