The nuclear receptor gene family in the Pacific oyster, Crassostrea gigas, contains a novel subfamily group

De novo transcriptome assembly of the marine gastropod Reishia clavigera for supporting toxic mechanism studies

The intertidal whelk Reishia clavigera is commonly used as a biomonitor of chemical contamination in the marine environment along Western Pacific region, and as a model for mechanistic studies of organotin-mediated imposex development. However, limited genomic resources of R. clavigera have restricted its role for the investigation of molecular mechanisms of such endocrine disruptions. This study, therefore, aimed to establish tissue-specific transcriptomes of the digestive gland, gonad, head ganglia, penis and the remaining body part of the male and female R. clavigera. By combining the results, a global transcriptome was obtained. A total of 578,134,720 high-quality filtered reads were obtained using Illumina sequencing. The R. clavigera transcriptome comprised of 38,466 transcripts and 32,798 unigenes with predicted open reading frames. The average length of transcripts was 1,709bp with N50 of 2,236bp. Based on sequence similarity searches against public databases, 28,657 transcripts and 24,403 unigenes had at least one BLAST hit. There were 17,530 transcripts and 14,897 unigenes annotated with at least one Gene Ontology (GO) term. Moreover, 5,776 transcripts and 5,137 unigenes were associated with 333 Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways. The numbers of unigenes were similar among the five target tissues and between sexes, but tissue-specific expression profiles were revealed by multivariate analyses. Based on the functional annotation, putative steroid hormone-associated unigenes were identified. In particular, we highlighted the presence of steroid hormone receptor homologues that could be the targets for mechanistic studies of the organotin-mediated imposex development in marine gastropods. This newly generated transcriptome assembly of R. clavigera provides a valuable molecular resource for ecotoxicological and environmental genomic studies.

Steroid Androgen Exposure during Development Has No Effect on Reproductive Physiology of Biomphalaria glabrata

Gastropod mollusks have been proposed as alternative models for male reproductive toxicity testing, due to similarities in their reproductive anatomy compared to mammals, together with evidence that endocrine disrupting chemicals can cause effects in some mollusks analogous to those seen in mammals. To test this hypothesis, we used the freshwater pulmonate snail, Biomphalaria glabrata, for which various genetic tools and a draft genome have recently become available, to investigate the effects of two steroid androgens on the development of mollusk secondary sexual organs. Here we present the results of exposures to two potent androgens, the vertebrate steroid; 5α-dihydrotestosterone (DHT) and the pharmaceutical anabolic steroid; 17α-methyltestosterone (MT), under continuous flow-through conditions throughout embryonic development and up to sexual maturity. Secondary sexual gland morphology, histopathology and differential gene expression analysis were used to determine whether steroid androgens stimulated or inhibited organ development. No significant differences between tissues from control and exposed snails were identified, suggesting that these androgens elicited no biologically detectable response normally associated with exposure to androgens in vertebrate model systems. Identifying no effect of androgens in this mollusk is significant, not only in the context of the suitability of mollusks as alternative model organisms for testing vertebrate androgen receptor agonists but also, if applicable to other similar mollusks, in terms of the likely impacts of androgens and anti-androgenic pollutants present in the aquatic environment.

A mollusk VDR/PXR/CAR-like (NR1J) nuclear receptor provides insight into ancient detoxification mechanisms

The origin and diversification of the metazoan endocrine systems represents a fundamental research issue in biology. Nuclear receptors are critical components of these systems. A particular group named VDR/PXR/CAR (NR1I/J) is central in the mediation of detoxification responses. While orthologues have been thoroughly characterized in vertebrates, a sparse representation is currently available for invertebrates. Here, we provide the first isolation and characterization of a lophotrochozoan protostome VDR/PXR/CAR nuclear receptor (NR1J), in the estuarine bivalve the peppery furrow shell (Scrobicularia plana). Using a reporter gene assay, we evaluated the xenobiotic receptor plasticity comparing the human PXR with the S. plana NR1Jβ. Our results show that the molluscan receptor responds to a natural toxin (okadaic acid) in a similar fashion to that reported for other invertebrates. In contrast, the pesticide esfenvalerate displayed a unique response, since it down regulated transactivation at higher concentrations, while for triclosan no response was observed. Additionally, we uncovered lineage specific gene duplications and gene loss in the gene group encoding NRs in protostomes with likely impacts on the complexity of detoxification mechanisms across different phyla. Our findings pave the way for the development of multi-specific sensor tools to screen xenobiotic compounds acting via the NR1I/J group.

Oct-GnRH, the first protostomian gonadotropin-releasing hormone-like peptide and a critical mini-review of the presence of vertebrate sex steroids in molluscs

In protostome and deuterosome invertebrates, neurosecretory cells play major roles in the endocrine system. The optic glands of cephalopods are indicators of sexual maturation. In mature octopuses, optic glands enlarge and secrete a gonadotropic hormone. A peptide with structural features similar to that of vertebrate gonadotropin-releasing hormone (GnRH) was isolated from the octopus, Octopus vulgaris, and was named oct-GnRH. The discovery of oct-GnRH has triggered structural determinations and predictions of other mollusc GnRH-like peptides in biochemical and in silico studies. Interestingly, cephalopods studied so far are characterized by a single molecular form of oct-GnRH with a C-terminal -Pro-Gly-NH2 sequence, which is critical for gonadotropin-releasing activity in vertebrates. Other molluscan GnRH-like peptides lack the C-terminal -Pro-Gly-NH2 sequence but have -X-NH2 or -Pro-Gly although all protostome GnRH-like peptides have yet to be sequenced. In marine molluscs, relationships between GnRH-like peptides and sex steroids have been studied to verify the hypothesis that molluscs have vertebrate-type sex steroid system. However, it is currently questionable whether such sex steroids are present and whether they play endogenous roles in the reproductive system of molluscs. Because molluscs uptake and store steroids from the environment and fishes release sex steroids into the external environment, it is impossible to rule out the contamination of vertebrate sex steroids in molluscs. The function of key enzymes of steroidogenesis within molluscs remains unclear. Thus, evidence to deny the existence of the vertebrate-type sex steroid system in molluscs has been accumulated. The elucidation of substances, which regulate the maturation and maintenance of gonads and other reproductive functions in molluscs will require rigorous and progressive scientific study.

Pharmacological evaluation of the mechanisms involved in increased adiposity in zebrafish triggered by the environmental contaminant tributyltin

One proposed contributing factor to the rise in overweight and obesity is exposure to endocrine disrupting chemicals. Tributyltin chloride (TBT), an organotin, induces adipogenesis in cell culture models and may increases adipose mass in vivo in vertebrate model organisms. It has been hypothesized that TBT acts via the peroxisome proliferator activated receptor (PPAR)γ-dependent pathway. However, the mechanisms involved in the effects of TBT exposure on in vivo adipose tissue metabolism remain unexplored. Semitransparent zebrafish larvae, with their well-developed white adipose tissue, offer a unique opportunity for studying the effects of toxicant chemicals and pharmaceuticals on adipocyte biology and whole-organism adiposity in a vertebrate model. Within hours, zebrafish larvae, treated at environmentally-relevant nanomolar concentrations of TBT, exhibited a remarkable increase in adiposity linked to adipocyte hypertrophy. Under the experimental conditions used, we also demonstrated that zebrafish larvae adipose tissue proved to be highly responsive to selected human nuclear receptor agonists and antagonists. Retinoid X receptor (RXR) homodimers and RXR/liver X receptor heterodimers were suggested to be in vivo effectors of the obesogenic effect of TBT on zebrafish white adipose tissue. RXR/PPARγ heterodimers may be recruited to modulate adiposity in zebrafish but were not a necessary requirement for the short term in vivo TBT obesogenic effect. Together, the present results suggest that TBT may induce the promotion of triacylglycerol storage in adipocytes via RXR-dependent pathways without necessary using PPAR isoforms.

Genetic and structural analyses of cytochrome P450 hydroxylases in sex hormone biosynthesis: Sequential origin and subsequent coevolution

Biosynthesis of steroid hormones in vertebrates involves three cytochrome P450 hydroxylases, CYP11A1, CYP17A1 and CYP19A1, which catalyze sequential steps in steroidogenesis. These enzymes are conserved in the vertebrates, but their origin and existence in other chordate subphyla (Tunicata and Cephalochordata) have not been clearly established. In this study, selected protein sequences of CYP11A1, CYP17A1 and CYP19A1 were compiled and analyzed using multiple sequence alignment and phylogenetic analysis. Our analyses show that cephalochordates have sequences orthologous to vertebrate CYP11A1, CYP17A1 or CYP19A1, and that echinoderms and hemichordates possess CYP11-like but not CYP19 genes. While the cephalochordate sequences have low identity with the vertebrate sequences, reflecting evolutionary distance, the data show apparent origin of CYP11 prior to the evolution of CYP19 and possibly CYP17, thus indicating a sequential origin of these functionally related steroidogenic CYPs. Co-occurrence of the three CYPs in early chordates suggests that the three genes may have coevolved thereafter, and that functional conservation should be reflected in functionally important residues in the proteins. CYP19A1 has the largest number of conserved residues while CYP11A1 sequences are less conserved. Structural analyses of human CYP11A1, CYP17A1 and CYP19A1 show that critical substrate binding site residues are highly conserved in each enzyme family. The results emphasize that the steroidogenic pathways producing glucocorticoids and reproductive steroids are several hundred million years old and that the catalytic structural elements of the enzymes have been conserved over the same period of time. Analysis of these elements may help to identify when precursor functions linked to these enzymes first arose.

Identification of Thyroid Hormones and Functional Characterization of Thyroid Hormone Receptor in the Pacific Oyster Crassostrea gigas Provide Insight into Evolution of the Thyroid Hormone System

Thyroid hormones (THs) play important roles in development, metamorphosis, and metabolism in vertebrates. During the past century, TH functions were regarded as a synapomorphy of vertebrates. More recently, accumulating evidence has gradually convinced us that TH functions also occur in invertebrate chordates. To date, however, TH-related studies in non-chordate invertebrates have been limited. In this study, THs were qualitatively detected by two reliable methods (HPLC and LC/MS) in a well-studied molluscan species, the Pacific oyster Crassostrea gigas. Quantitative measurement of THs during the development of C. gigas showed high TH contents during embryogenesis and that oyster embryos may synthesize THs endogenously. As a first step in elucidating the TH signaling cascade, an ortholog of vertebrate TH receptor (TR), the most critical gene mediating TH effects, was cloned in C. gigas. The sequence of CgTR has conserved DNA-binding and ligand-binding domains that normally characterize these receptors. Experimental results demonstrated that CgTR can repress gene expression through binding to promoters of target genes and can interact with oyster retinoid X receptor. Moreover, CgTR mRNA expression was activated by T4 and the transcriptional activity of CgTR promoter was repressed by unliganded CgTR protein. An atypical thyroid hormone response element (CgDR5) was found in the promoter of CgTR, which was verified by electrophoretic mobility shift assay (EMSA). These results indicated that some of the CgTR function is conserved. However, the EMSA assay showed that DNA binding specificity of CgTR was different from that of the vertebrate TR and experiments with two dual-luciferase reporter systems indicated that l-thyroxine, 3,3′,5-triiodothyronine, and triiodothyroacetic acid failed to activate the transcriptional activity of CgTR. This is the first study to functionally characterize TR in mollusks. The presence of THs and the functions of CgTR in mollusks contribute to better understanding of the evolution of the TH system.

A Label-Free Proteomic Analysis on Competent Larvae and Juveniles of the Pacific Oyster Crassostrea gigas

Current understandings on the molecular mechanisms underlying bivalve metamorphosis are still fragmentary, and a comprehensive description is required. In this study, using a large-scale label-free proteomic approach, we described and compared the proteomes of competent larvae (CL) and juveniles (JU) of the Pacific oyster, Crassostrea gigas. A total of 788 proteins were identified: 392 in the CL proteome and 636 in the JU proteome. Gene Ontology analysis of the proteome from each sample revealed active metabolic processes in both stages. Further quantitative analyses revealed 117 proteins that were differentially expressed between the two samples. These proteins were divided into eight groups: cytoskeleton and cell adhesion, protein synthesis and degradation, immunity and stress response, development of particular tissues, signal regulation, metabolism and energy supply, transport, and other proteins. A certification experiment using real-time PCR assay confirmed 20 of 30 examined genes exhibited the same trends at the mRNA and protein levels. The differentially expressed proteins may play roles in tissue remodeling, signal transduction, and organ development during and after metamorphosis. Novel roles were proposed for some differentially expressed proteins, such as chymotrypsin. The results of this work provide an overview of metamorphosis and post-metamorphosis development of C. gigas at the protein level. Future studies on the functions of the differentially expressed proteins will help to obtain a more in-depth understanding of bivalve metamorphosis.

In Vivo Administration of Scallop GnRH-Like Peptide Influences on Gonad Development in the Yesso Scallop, Patinopecten yessoensis

Existing research on the role of gonadotropin-releasing hormone (GnRH) in bivalve reproduction is inadequate, even though a few bivalve GnRH orthologs have been cloned. The objective of this paper was to elucidate the in vivo effect of GnRH administration in Yesso scallop reproduction. We performed in vivo administration of scallop GnRH (py-GnRH) synthetic peptide into the developing gonad, and analyzed its effect on gonad development for 6 weeks during the reproductive season. The resulting sex ratio in the GnRH administered (GnRH(+)) group might be male biased, whereas the control (GnRH(-)) group had an equal sex ratio throughout the experiment. The gonad index (GI) of males in the GnRH(+) group increased from week 2 to 24.8% at week 6. By contrast the GI of the GnRH(-) group peaked in week 4 at 16.6%. No significant difference was seen in female GI between the GnRH(+) and GnRH(-) groups at any sampling point. Oocyte diameter in the GnRH(+) group remained constant (about 42–45 μm) throughout the experiment, while in the GnRH(-) group it increased from 45 to 68 μm i.e. normal oocyte growth. The number of spermatogonia in the germinal acini of males in the GnRH(+) group increased from week 4 to 6. Hermaphrodites appeared in the GnRH(+) group in weeks 2 and 4. Their gonads contained many apoptotic cells including oocytes. In conclusion, this study suggests that py-GnRH administration could have a potential to accelerate spermatogenesis and cause an inhibitory effect on oocyte growth in scallops.

Evolution of a novel nuclear receptor subfamily with emphasis on the member from the Pacific oyster Crassostrea gigas

Nuclear receptors (NRs) belong to the transcription factor superfamily that regulates development, homeostasis, differentiation, and reproduction in metazoans via control of gene expression. Recently, rapid advances in genome projects on various metazoans have provided new opportunities for studying the evolution and function of NRs. Typically structured NRs are divided into six subfamilies. Here, the gene for a typically structured NR (CgNR8A1) was cloned from the Pacific oyster Crassostrea gigas. However, this novel receptor could not be assigned to a known NR subfamily. By data mining, nine other CgNR8A1 gene homologs were identified in metazoans such as cnidarians, mollusks, annelids, echinoderms, hemichordates, and cephalochordates. Phylogenetic analysis showed that these receptors belonged to a novel NR subfamily, hereafter designated as NR8. Evolutionary analysis revealed that the NR8 subfamily was phylogenetically the third-oldest NR subfamily, and it originated from a common ancestor of Eumetazoa; several gene loss events occurred independently in ancestors of vertebrates, ecdysozoans, and platyhelminths, which do not have NR8 members. Furthermore, the function of CgNR8A1 was investigated to provide an insight into the functions of this novel NR subfamily. A nuclear localization signal peptide, GKHRNKKPRLD, was identified in CgNR8A1, and a recombinant full-length protein of CgNR8A1 was localized in the nuclei of HeLa cells. The mRNA expression profile of CgNR8A1 suggested that it might be involved in the embryogenesis of C. gigas. The electrophoretic mobility shift assay showed that CgNR8A1 binds strongly to conserved DNA core motifs DR0, DR2, and DR4 and weakly to DR1, DR3, DR5, Half, and Pal0. In summary, the novel NR8 subfamily identified in this study improves our understanding of NR evolution, and the functional analysis of CgNR8A1 provided further insights into the functions of NR8A1s.

The nuclear receptors of Biomphalaria glabrata and Lottia gigantea: implications for developing new model organisms

Nuclear receptors (NRs) are transcription regulators involved in an array of diverse physiological functions including key roles in endocrine and metabolic function. The aim of this study was to identify nuclear receptors in the fully sequenced genome of the gastropod snail, Biomphalaria glabrata, intermediate host for Schistosoma mansoni and compare these to known vertebrate NRs, with a view to assessing the snail's potential as a invertebrate model organism for endocrine function, both as a prospective new test organism and to elucidate the fundamental genetic and mechanistic causes of disease. For comparative purposes, the genome of a second gastropod, the owl limpet, Lottia gigantea was also investigated for nuclear receptors. Thirty-nine and thirty-three putative NRs were identified from the B. glabrata and L. gigantea genomes respectively, based on the presence of a conserved DNA-binding domain and/or ligand-binding domain. Nuclear receptor transcript expression was confirmed and sequences were subjected to a comparative phylogenetic analysis, which demonstrated that these molluscs have representatives of all the major NR subfamilies (1-6). Many of the identified NRs are conserved between vertebrates and invertebrates, however differences exist, most notably, the absence of receptors of Group 3C, which includes some of the vertebrate endocrine hormone targets. The mollusc genomes also contain NR homologues that are present in insects and nematodes but not in vertebrates, such as Group 1J (HR48/DAF12/HR96). The identification of many shared receptors between humans and molluscs indicates the potential for molluscs as model organisms; however the absence of several steroid hormone receptors indicates snail endocrine systems are fundamentally different.

Families of nuclear receptors in vertebrate models: characteristic and comparative toxicological perspective

Various synthetic chemicals are ligands for nuclear receptors (NRs) and can cause adverse effects in vertebrates mediated by NRs. While several model vertebrates, such as mouse, chicken, western clawed frog and zebrafish, are widely used in toxicity testing, few NRs have been well described for most of these classes. In this report, NRs in genomes of 12 vertebrates are characterized via bioinformatics approaches. Although numbers of NRs varied among species, with 40-42 genes in birds to 66-74 genes in teleost fishes, all NRs had clear homologs in human and could be categorized into seven subfamilies defined as NR0B-NR6A. Phylogenetic analysis revealed conservative evolutionary relationships for most NRs, which were consistent with traditional morphology-based systematics, except for some exceptions in Dolphin (Tursiops truncatus). Evolution of PXR and CAR exhibited unexpected multiple patterns and the existence of CAR possibly being traced back to ancient lobe-finned fishes and tetrapods (Sarcopterygii). Compared to the more conservative DBD of NRs, sequences of LBD were less conserved: Sequences of THRs, RARs and RXRs were ≥90% similar to those of the human, ERs, AR, GR, ERRs and PPARs were more variable with similarities of 60%-100% and PXR, CAR, DAX1 and SHP were least conserved among species.

Marine invertebrate xenobiotic-activated nuclear receptors: their application as sensor elements in high-throughput bioassays for marine bioactive compounds

Developing high-throughput assays to screen marine extracts for bioactive compounds presents both conceptual and technical challenges. One major challenge is to develop assays that have well-grounded ecological and evolutionary rationales. In this review we propose that a specific group of ligand-activated transcription factors are particularly well-suited to act as sensors in such bioassays. More specifically, xenobiotic-activated nuclear receptors (XANRs) regulate transcription of genes involved in xenobiotic detoxification. XANR ligand-binding domains (LBDs) may adaptively evolve to bind those bioactive, and potentially toxic, compounds to which organisms are normally exposed to through their specific diets. A brief overview of the function and taxonomic distribution of both vertebrate and invertebrate XANRs is first provided. Proof-of-concept experiments are then described which confirm that a filter-feeding marine invertebrate XANR LBD is activated by marine bioactive compounds. We speculate that increasing access to marine invertebrate genome sequence data, in combination with the expression of functional recombinant marine invertebrate XANR LBDs, will facilitate the generation of high-throughput bioassays/biosensors of widely differing specificities, but all based on activation of XANR LBDs. Such assays may find application in screening marine extracts for bioactive compounds that could act as drug lead compounds.

Genome-wide identification of nuclear receptor (NR) superfamily genes in the copepod Tigriopus japonicus

Background

Nuclear receptors (NRs) are a large superfamily of proteins defined by a DNA-binding domain (DBD) and a ligand-binding domain (LBD). They function as transcriptional regulators to control expression of genes involved in development, homeostasis, and metabolism. The number of NRs differs from species to species, because of gene duplications and/or lineage-specific gene losses during metazoan evolution. Many NRs in arthropods interact with the ecdysteroid hormone and are involved in ecdysone-mediated signaling in arthropods. The nuclear receptor superfamily complement has been reported in several arthropods, including crustaceans, but not in copepods. We identified the entire NR repertoire of the copepod Tigriopus japonicus, which is an important marine model species for ecotoxicology and environmental genomics.

Results

Using whole genome and transcriptome sequences, we identified a total of 31 nuclear receptors in the genome of T. japonicus. Nomenclature of the nuclear receptors was determined based on the sequence similarities of the DNA-binding domain (DBD) and ligand-binding domain (LBD). The 7 subfamilies of NRs separate into five major clades (subfamilies NR1, NR2, NR3, NR4, and NR5/6). Although the repertoire of NR members in, T. japonicus was similar to that reported for other arthropods, there was an expansion of the NR1 subfamily in Tigriopus japonicus. The twelve unique nuclear receptors identified in T. japonicus are members of NR1L. This expansion may be a unique lineage-specific feature of crustaceans. Interestingly, E78 and HR83, which are present in other arthropods, were absent from the genomes of T. japonicus and two congeneric copepod species (T. japonicus and Tigriopus californicus), suggesting copepod lineage-specific gene loss.

Conclusions

We identified all NR receptors present in the copepod, T. japonicus. Knowledge of the copepod nuclear receptor repertoire will contribute to a better understanding of copepod- and crustacean-specific NR evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-993) contains supplementary material, which is available to authorized users.

Retinoid metabolism in invertebrates: when evolution meets endocrine disruption

Recent genomic and biochemical evidence in invertebrate species pushes back the origin of the retinoid metabolic and signaling modules to the last common ancestor of all bilaterians. However, the evolution of retinoid pathways are far from fully understood. In the majority of non-chordate invertebrate lineages, the ongoing functional characterization of retinoid-related genes (metabolism and signaling pathways), as well as the characterization of the endogenous retinoid content (precursors and active retinoids), is still incomplete. Despite limited, the available data supports the presence of biologically active retinoid pathways in invertebrates. Yet, the mechanisms controlling the spatial and temporal distribution of retinoids as well as their physiological significance share similarities and differences with vertebrates. For instance, retinol storage in the form of retinyl esters, a key feature for the maintenance of retinoid homeostatic balance in vertebrates, was only recently demonstrated in some mollusk species, suggesting that such ability is older than previously anticipated. In contrast, the enzymatic repertoire involved in this process is probably unlike that of vertebrates. The suggested ancestry of active retinoid pathways implies that many more metazoan species might be potential targets for endocrine disrupting chemicals. Here, we review the current knowledge about the occurrence and functionality of retinoid metabolic and signaling pathways in invertebrate lineages, paying special attention to the evolutionary origin of retinoid storage mechanisms. Additionally, we summarize existing information on the endocrine disruption of invertebrate retinoid modules by environmental chemicals. Research priorities in the field are highlighted.