Skip to main content

Nuclear receptor complement of the cnidarian Nematostella vectensis: phylogenetic relationships and developmental expression patterns



Nuclear receptors are a superfamily of metazoan transcription factors that regulate diverse developmental and physiological processes. Sequenced genomes from an increasing number of bilaterians have provided a more complete picture of duplication and loss of nuclear receptors in protostomes and deuterostomes but have left open the question of which nuclear receptors were present in the cnidarian-bilaterian ancestor. In addition, nuclear receptor expression and function are largely uncharacterized within cnidarians, preventing determination of conserved and novel nuclear receptor functions in the context of animal evolution.


Here we report the first complete set of nuclear receptors from a cnidarian, the starlet sea anemone Nematostella vectensis. Genomic searches using conserved DNA- and ligand-binding domains revealed seventeen nuclear receptors in N. vectensis. Phylogenetic analyses support N. vectensis orthologs of bilaterian nuclear receptors in four nuclear receptor subfamilies within nuclear receptor family 2 (COUP-TF, TLL, HNF4, TR2/4) and one putative ortholog of GCNF (nuclear receptor family 6). Other N. vectensis genes grouped well with nuclear receptor family 2 but represented lineage-specific duplications somewhere within the cnidarian lineage and were not clear orthologs of bilaterian genes. Three nuclear receptors were not well-supported within any particular nuclear receptor family. The seventeen nuclear receptors exhibited distinct developmental expression patterns, with expression of several nuclear receptors limited to a subset of developmental stages.


N. vectensis contains a diverse complement of nuclear receptors including orthologs of several bilaterian nuclear receptors. Novel nuclear receptors in N. vectensis may be ancient genes lost from triploblastic lineages or may represent cnidarian-specific radiations. Nuclear receptors exhibited distinct developmental expression patterns, which are consistent with diverse regulatory roles for these genes. Understanding the evolutionary relationships and developmental expression of the N. vectensis nuclear receptor complement provides insight into the evolution of the nuclear receptor superfamily and a foundation for mechanistic characterization of cnidarian nuclear receptor function.


Nuclear receptors (NRs) are one of the most abundant classes of transcription factors in metazoans and coordinate diverse processes ranging from developmental patterning to physiological regulation. NRs appear to be restricted to animals because extensive work has shown nuclear receptors are present in all metazoan phyla (e.g., sponges [1, 2], cnidarians [3, 4], and bilaterians [5]) but not in plants, fungi or choanoflagellates [6]. The signature motifs of NRs are the DNA-binding domain (DBD), which includes a Cys-Cys zinc coordinating region, and the ligand-binding domain (LBD), a carboxy-terminal domain that binds ligands and facilitates receptor dimerization and coactivator recruitment (Figure 1). High sequence conservation within these domains, particularly the DBD, has permitted identification of NRs from a number of animals and phylogenetic analysis of the NR superfamily [5].

Figure 1

Nuclear receptor gene structure. Nuclear receptors are composed of two well-conserved regions, the DNA binding domain (DBD) and the ligand binding domain (LBD).

Genes in the NR superfamily have been classified into six major NR families, many of which are further divided into subfamilies of orthologs and paralogs. All six families are represented in both the protostome and deuterostome lineages, supporting the hypothesis that NRs had undergone extensive radiation prior to the divergence of triploblastic animals. An analysis by Bertrand et al. [5] of NRs from nine bilaterian genomes did much to resolve the evolutionary history of this superfamily. They inferred 25 NRs that likely existed in the Urbilaterian with at least one gene from each of the six families. After this divergence, many NRs underwent extensive duplication in the vertebrates, particularly fishes, and one gene underwent extensive duplication in nematodes [7]. With the increasing availability of sequenced genomes and PCR-based surveys of NR diversity, particularly from non-ecdysozoan protostomes [810], it has become apparent that a number of NRs have also been lost since the Urbilaterian ancestor, particularly in insects and nematodes.

Because the NR superfamily had already diversified prior the divergence of the protostomes and deuterostomes, it is essential to characterize NRs in "basal" metazoans to understand not only when the NR families evolved but also the evolutionary relationships among these families. NRs have been reported from a handful of cnidarians from three of the four classes (Anthozoa [4], Cubozoa [11], Hydrozoa [3]) and two sponges [1, 2]. Sequence and phylogenetic analyses of identified NRs have led these authors to suggest that many or all of these genes belong to NR family 2 with homologs of TLL (tailless, NR2E) in corals [4], HNF4 (hepatocyte nuclear factor 4, NR2A) in coral [4] and sponge [1], RXR (retinoic × receptor, NR2B) in the box jelly [11], and COUP-TF (chicken ovalbumin upstream promoter transcription factor, NR2F) in coral and hydra ([3, 4]).

Despite the importance of understanding how NR regulation of development and physiology has evolved, little has been reported regarding the expression and function of NRs in non-bilaterian animals. Some cnidarian NRs have been shown to share conserved features with their vertebrate orthologs. For example, the RXR homolog in the box jelly Tripedalia cystophora specifically binds 9-cis retinoic acid with an affinity similar to that of vertebrate RXRs, supporting a conserved function of RXR in these evolutionary distant taxa [11]. The COUP-TF homolog from hydra is expressed in neuronal cells, suggesting an evolutionary origin in nerve cell specification [3]. Many of the identified NRs in cnidarians are orphan receptors (lacking known ligands). Cnidarian tissues contain diverse lipids including prostaglandins, fatty acids, sterols and steroids, and it is plausible that some of these compounds act as ligands for nuclear receptors [1218]. In addition, the developmental expression of most cnidarian NRs remains uncharacterized. Temporal patterns of expression can provide significant insight into the potential role for NRs during development and metamorphosis as well as the opportunity for co-expressed NRs to interact with one another [19].

Here we present the first genomic analysis of nuclear receptors in the starlet sea anemone, Nematostella vectensis, and quantify expression of all identified NRs throughout the life cycle. We identified seventeen nuclear receptors in the genome of N. vectensis. Many of the NRs did not show clear orthology to bilaterian NRs, and some N. vectensis NRs appear to have resulted from lineage specific gene duplications. Expression of most NRs varied significantly during development, with clusters of NRs exhibiting similar expression patterns.


Phylogenetic relationships of N. vectensisnuclear receptors

We identified 17 NRs from N. vectensis with a combination of bioinformatic searches of the assembled N. vectensis genome [20]. We have designated these as NvNR1-17 (Accession Numbers Table 1, numbering is arbitrary and does not indicate phylogenetic position within the NR superfamily; for amino acid sequences see Additional file 1). All of these genes were amplified from cDNA preparations and confirmed as expressed transcripts.

Table 1 Accession numbers of N. vectensis nuclear receptors from NCBI GenBank, JGI Nematostella Database, and ESTs (where available).

When DBD plus LBD were used for phylogenetic analyses of N. vectensis and bilaterian NRs, the topology of the tree with the highest likelihood was consistent with monophyletic relationships of the recognized NR families. Nodes for NR families 1-4 were each supported by high bootstrap scores (Figure 2, for a more detailed tree with nodes expanded see Additional file 2). Although NR families 5 and 6 were each weakly supported (BS = 40 and 37, respectively), a clade grouping these two families was fairly well supported (BS = 79) as previously shown by Thornton [21]. Phylogenetic analyses using just the DBD or portion of the LBD resulted in more poorly supported nodes throughout the tree and in some cases failed to recover monophyletic NR families (see Additional files 3 and 4). In addition, a complementary analysis using neighbor joining with the DBD and LBD also resulted in weakly supported nodes for the NR families (see Additional file 5). Therefore, we inferred the phylogenetic relationships of the N. vectensis NRs from the likelihood analysis based on the DBD plus LBD.

Figure 2

Phylogenetic relationships of N. vectensis NRs (NvNRs) depicted by a maximum-likelihood tree with percent bootstrap support of 1000 bootstraps for nodes of evolutionary significance. The alignment was constructed using the DBD and a portion of the LBD (additional details in text). Monophyletic relationships for bilaterian sequences in particular families or subfamilies were recovered in all cases but one (NR2E) and are depicted as horizontal triangle. A full tree with all taxa depicted as individual branches is presented in Additional File 2. The tree was rooted with NR family 2 because this family is monophyletic, represents a defined nuclear receptor family, and includes what is likely the original nuclear receptor, represented by an HNF4-like homolog from the sponge Amphimedon queenslandica [1]. Most N. vectensis NRs group within NR family 2 and many (e.g., NvNR11- NvNR14) are supported as independent radiations of subfamilies within this family. N. vectensis also has a NR related to GCNF (NR family 6), but with low support. Three NRs did not group with any previously described family (NvNR1-3). These genes may represent ancestral genes that later diversified into one or more of the NR 1, 3, and 4 families.

NvNR1, 2, and 3 cannot be unambiguously assigned to a family and nest as equally related to NR families 1 and 4. A majority of N. vectensis NRs (NvNR4-16) are strongly supported as homologs of genes in NR family 2. These N. vectensis NRs are associated with four subfamilies within family 2. NvNR4 is an ortholog of HNF4 (subfamily 2A), which has also been identified in vertebrates, insects, and nematodes [22, 23]. NvNR5-9 cluster with subfamily 2E (TLL, FAX, PNR). NvNR5 and 6 are most closely related to TLL with moderate bootstrap support. NvNR7-9 do not have a clear orthologs with the bilaterian genes from this subfamily and thus may represent a cnidarian-specific radiation. NvNR10-14 group with subfamily 2F (COUP-TFs), and among these, NvNR10 is most closely related to the bilaterian COUP-TFs. NvNR15-16 are supported as homologs of the testicular orphan receptors (TR2/TR4, subfamily 2C/D). An RXR homolog (jRXR, subfamily 2B) has been identified in the box jelly Tripedalia cystophora [11], but an RXR homolog is not apparent in the N. vectensis genome. NvNR17 groups with NR family 6 (GCNF, germ cell nuclear factor) but with relatively weak support. In the neighbor joining analysis, NvNR17 was positioned as an outgroup to NR families 5 and 6 but with low bootstrap support. A GCNF-like gene has not been previously reported from a cnidarian.

Gene structure

N. vectensis NRs are composed of between 3 and 7 exons and span 2059 to 11741 bp of genomic sequence (see Additional file 6). Because we currently have complete transcript sequence with 5- and 3- prime untranslated region (UTR) for fewer than half of the NRs (6 of 17), additional exons may be identified when complete transcripts are confirmed. In some cases (e.g., NvNR10), exons were composed completely of UTR sequence; thus it is likely that additional exons will be identified for some of the NRs with uncharacterized UTR sequence.

Phylogenetic analyses indicated that several NRs in N. vectensis are most closely related to one another suggesting potential lineage-specific duplication. With the exception of two genes (NvNR11 and NvNR12), the NRs in N. vectensis are located on different scaffolds. Thus, the NR diversity in N. vectensis does not appear to have resulted from recent tandem duplications, but instead may represent more ancient duplications within the cnidarian lineage. Some closely related genes also shared conserved intron-exon structure (see Additional file 7). For example, NvNR8 and NvNR9 group together with 100% bootstrap support and share two of three splice sites.

Developmental expression of nuclear receptors

One-way ANOVA indicated that 13 of the 17 NRs varied significantly in their level of expression during development (Figure 3). Even with relatively coarse developmental staging, we observed clear temporal patterns of gene expression (Figure 4a). NvNR8, 11, 12 and 13 were primarily expressed during embryonic and/or early larval stages with little expression in the juvenile and adult stages. NvNR1, 7, 14 and 15 were most highly expressed during larval and juvenile stages. NvNR4, 5, 6, 9, 16, and 17 showed increased expression during larval and juvenile stages and continued high expression in adults. Finally, NvNR2, 3, and 10 were primarily expressed in adults, and NvNR10 was only minimally expressed in earlier stages. 18S expression did not vary significantly across the stages and is apparently a good housekeeping gene with respect to N. vectensis development. We also tested GAPDH as a potential control gene to standardize gene expression. GAPDH expression was highly variable across these developmental stages and was not used for comparisons of NR expression among stages.

Figure 3

Developmental expression of Nematostella vectensis nuclear receptors (NR1-NR17) and putative housekeeping genes (18S and GAPDH). Expression measured using qPCR, as described in text. Expression in molecules per μl cDNA indicated on the y-axis. Bar indicate the mean ± standard error of 3-5 biological replicates. The first bars ("cleavage") represents a single sample, which was not included in statistical analysis. Expression patterns were analyzed with a one-way ANOVA followed by pairwise comparisons with Tukey's Honestly Significant Difference Test. Letters indicate groups that were different at a level of p = 0.05. GAPDH showed variable expression over the sampled stages and was not used as a standard for comparing gene expression of NRs among developmental stages.

Figure 4

Heat maps indicating developmental expression of Nematostella vectensis nuclear receptors. Cladogram indicates the topology from Figure 2 with other taxa removed. (A) For each gene, expression is shown on a linear scale as a percentage of the stage with the highest expression. (B) Expression is shown on a log scale as molecules per microliter, calculated in relation to standard curves constructed from serially diluted plasmids.

Absolute levels of gene expression, normalized to plasmid standards, indicated large differences in number of expressed transcripts among genes (Figure 4b). Although its expression was developmentally variable (Figure 3, 4a), NvNR4 was the most highly expressed NR over all stages (Figure 4b). NvNR3, 9 and 17 also showed high levels of overall expression.

In several cases, phylogenetically closely related genes showed divergent temporal expression patterns during development. For example, NR11 was expressed at very low levels (Figure 4b), primarily during embryonic and early larval stages (Figure 4a). In contrast, NR10 was expressed at moderate levels, mostly in adults.

Principal components analysis (PCA) is a statistical technique for reducing the dimensionality of a dataset, in this case facilitating identification of NRs with similar expression patterns. PCA was conducted for all 17 genes from all developmental stages (Figure 5). The first principal component (PC1) accounted for 39% of the overall variance in gene expression and second principal component (PC2) accounted for 36% of the variance. When the values for PC1 and PC2 are shown on a scatter plot (Figure 5B), three clusters emerged, NvNR11, 12 and 13 clustered to the exclusion of other genes along PC1 and were each highly expressed in early developmental stages (embryo, early larva). A second cluster contained NvNR2, 3, and 10, which are most strongly expressed in adults. The remaining genes constituted a third cluster and were not strongly distinguished from one another by this analysis.

Figure 5

Principal components analysis of nuclear receptor relative expression data. Analysis was conducted on gene expression normalized to the maximum expression for each gene, as in Figure 4A. (A) Eigen vectors showing relative contribution of each life stage to each principal component. Longer bars indicate stronger relative contribution. The first and second principal components explained 39% and 36% of the variance, respectively. (B) Scatter plot of first and second principal component scores for all seventeen nuclear receptors. Colors indicate three distinct clusters of NRs showing similar expression patterns. Numbers 1-17 indicate N. vectensis nuclear receptors 1-17, respectively.


Cnidarians both occupy a key evolutionary position and are ecologically important as predators, prey and structure-builders in marine and freshwater environments. As "basal" metazoans, cnidarians form an outgroup to the bilaterian animals and are intermediate in complexity between sponges and bilaterians. In spite of the acknowledged importance of understanding cnidarians from both an ecological and an evolutionary perspective, cnidarian physiology is poorly understood, particularly at the molecular level. Thus, characterizing cnidarian NRs provides insight both into evolution of NR signaling and bioregulatory processes in a major group of aquatic animals.

Homology and expression of NvNRs

Analysis of the phylogenetic relationships among NRs from the starlet sea anemone, Nematostella vectensis, provides strong support for the diversification of NR family 2 prior to the divergence of the cnidarian and bilaterian lineages. Phylogenetic analyses of N. vectensis and bilaterian NRs revealed cnidarian orthologs of HNF4, COUP, TLL, and TR2/4. Grasso et al. [4] identified 10 NRs in the coral Acropora millepora, including several members of NR family 2. Our results support their conclusion that family 2 was well diversified prior to the split between the cnidarian and bilaterian lineages.

NvNR4 is an apparent homolog of hepatocyte nuclear factor 4 (HNF4, subfamily 2A), which has also been identified in vertebrates, insects, nematodes, and corals [4, 22, 23]. In mammals, HNF4α binds endogenous fatty acid ligands [24, 25] and regulates hepatocyte differentiation, energetic metabolism and xenobiotic detoxification. HNF4 also regulates insect gut development [26, 27]. HNF4 underwent extensive duplications in C. elegans [28] such that inferring homologous functions for HNF4 from the nematode is difficult. One C. elegans paralog, nhr-49, is involved in energetic metabolism as a regulator of fat storage [29]. In addition, an HNF4 homolog has been cloned from sponges and may be similar to the ancestral NR [1]. Sponge HNF4 is expressed throughout development in ciliated column cells of the outer epithelium [1]. In N. vectensis, NvNR4 was expressed at all developmental time points with higher mean expression from larval to adult stages. We predict that NvNR4 will be expressed in endodermal tissue, which would be consistent with conserved roles in development of the digestive epithelium, energetic metabolism, and/or detoxification.

NvNR5- NvNR9 cluster with subfamily 2E (TLL/TLX, FAX, PNR); members of this subfamily have previously been reported in corals [4, 30]. In our phylogenetic analyses, we were unable to fully resolve the evolutionary relationships of N. vectensis and bilaterian NRs within subfamily 2E (Figure 2). In other animals, members of subfamily 2E are involved in neural differentiation. Tailless homologs in mammals (TLX) and insects (TLL) are particularly important for eye and forebrain development [3134] and embryonic anterior-posterior patterning [35]. FAX-1 regulates neural patterning in C. elegans [36, 37], and PNR plays a more specialized role by regulating retinal development in vertebrates [38]. Several N. vectensis members of subfamily 2E (i.e., NvNR5, 7, 8) were highly expressed during early developmental stages coinciding with neurogenesis and embryonic patterning [39]. The stage of maximal expression varied among members of this subfamily, indicating potential diversification of function within the cnidarian lineage.

NvNR10-NvNR14 group with subfamily 2F (COUP-TFs); COUP-TFs have also been identified in corals, hydra, flatworms, sea urchin, and lancelets [3, 4, 4042]. NvNR10 was most closely related to bilaterian COUP-TFs, while NvNR11- NvNR14 are supported as an independent radiation with this subfamily. COUP-TFs generally act as transcriptional repressors and regulate development of muscles, the heart and the nervous system, particularly differentiation of the hindbrain and photoreceptors [4346]. In hydra, COUP-TF is expressed in nematoblasts (nematocyte precursors) and in a subset of neuronal cells, consistent with a conserved role of COUP-TF in regulating neural differentiation [3]. COUP-TFs can also affect reproductive processes through cross-talk with estrogen receptors and ecdysteroid receptors in vertebrates and insects, respectively [4750]. COUP-TF homologs from N. vectensis showed diverse expression patterns during development. For example, NvNR12 and13 were primarily expressed during embryonic and/or early larval stages with little expression in the juvenile and adult stages. The expression of these COUP-TF-like genes coincides with neurogenesis and embryo patterning [39]. In contrast, we observed relatively high expression of NvNR10 in adults but lower expression during embryogenesis and larval development.

We identified two members of the subfamily 2C/D (TR2/4, NvNR15-16). TR2/4 homologs have been identified in a range of animals including vertebrates, sea urchins, ascidians, nematodes, insects, flatworms, and corals [5154]. Protostomes and most deuterostomes have a single TR2/4 homolog. Interestingly, both N. vectensis and the coral Acropora millepora [4] contain two TR2/4 homologs, suggesting a cnidarian-specific duplication. In general, TR2/4 homologs act as transcriptional repressors through several mechanisms including competition with other nuclear receptors for binding sites and co-factors [55, 56]. TR2/4 homologs are broadly expressed in vertebrate tissues [57] and throughout development in mammals [57], ascidians [51], and flatworms [58]. We observed similar ubiquitous expression in N. vectensis where both NvNR15 and 16 were expressed throughout all developmental stages.

NvNR17 forms a monophyletic grouping with NR family 6 (GCNFs, germ cell nuclear factors). Although bootstrap support for this assignment is relatively weak, likelihood analyses of the DBD plus LBD and each of these domains individually group NvNR17 as a homolog to the GCNFs. While neighbor joining analysis placed NvNR17 as an outgroup to NR families 5 and 6, bootstrap support for deep nodes of this tree were low relative to the maximum likelihood analyses, potentially obscuring evolutionary relationships. GCNF has been reported from a number of bilaterian phyla. This is the first report of a putative cnidarian homolog from this NR family. Expression of GCNF homologs varies greatly among taxa. In D. melanogaster, DmHR4 is expressed following pulses of 20-hydroxyecdysone during discrete stages of embryogenesis, the second larval instar and prepupal development [59, 60]. In C. elegans, CeNHR-91 is expressed in response to ATP-binding cassette protein E in embryos, larvae and several adult tissues [61]. In vertebrates, GCNF is expressed primarily expressed in gonadal tissues and during embryogenesis [62, 63]. In vertebrates, GCNF is involved in neuronal differentiation, germ cell development, and axial patterning, and gametogenesis [6365]. In our study, NvNR17 expression increased throughout the sampled developmental stages. Highest expression was in adult tissues, which may reflect a role in gametogenesis.

NvNR1-3 did not group well with any previously defined NR family. These genes are most closely related to NR Family 1 and 4. With the available data, we cannot conclude whether these represent the descendents of ancient NRs that later diversified into recognized NR families or if these genes represent an independent radiation somewhere within the cnidarian lineage. Further sampling within the Cnidaria will be necessary to help resolve the evolutionary relationship of these NRs.

We did not identify a member of NR family 3 in N. vectensis. NR family 3 is represented in protostomes and deuterostomes [5, 21, 66, 67], and has been more recently identified in the placozoan Trichoplax adhaerens [[68], Reitzel and Tarrant, unpublished data] but has not yet been reported from any cnidarian. Thus, NR 3 homologs have apparently been lost early in the cnidarian lineage. NR family 3 contains the vertebrate-type steroid (i.e., non-ecdysteroid) receptors. While NR3 family members were not identified in N. vectensis, estrogens and other steroids have been detected in cnidarian tissues [12, 15, 69], are apparently released during spawning events [12, 69], and can experimentally affect coral growth and reproduction [70]. The mechanism for steroid action in cnidarians is currently unknown and may be mediated through nuclear receptors that are not orthologs of the vertebrate steroid receptors or through alternative mechanisms.

RXR loss in N. vectensis

N. vectensis appears to have lost the ortholog of RXR because this gene is present in another cnidarian, the cubozoan Tripedalia cystophora [11]. The RXR homolog (jRXR) from T. cystophora bound 9-cis retinoic acid, and heterodimerized with vertebrate thyroid hormone receptor [11]. This evidence suggests that this cnidarian RXR functions in a similar way as vertebrate RXRs. In contrast to the box jelly, an RXR ortholog was not identified in a PCR-based survey of NRs from the coral Acropora millepora [4]. However, RXR immunoreactivity has been reported in epithelial nerve tissue of the coral A. millepora and the sea pansy Renilla koellikeri, but the genes were not identified [71]. Indeed the authors of the immunological study indicated that the epitope for the antibody has highest sequence similarity to A. millepora nuclear receptor 8, but phylogenetic analyses indicate that AmNR8 is most similar to TR2/4 [4]. The absence of an RXR ortholog in N. vectensis and in corals is particularly interesting from a functional perspective because vertebrate and insect RXRs form heterodimers with many other nuclear receptors to regulate gene expression [48, 7274]. In a parallel example, C. elegans independently lost an RXR ortholog, but another NR (nhr-49, HNF4-like) can dimerize with several other NRs [75]. Thus, it has been suggested that nhr-49 may have assumed an RXR-like function as a common heterodimerization partner [76]. Similarly, N. vectensis and other cnidarians may have one or more receptors that "promiscuously" form heterodimers with other NRs in a manner similar to bilaterian RXRs. Alternatively, it is possible that RXR has been lost and not replaced with a functional equivalent.

Gene duplication and divergence

Within NR family 2, our data from N. vectensis support independent radiations of several subfamilies. Phylogenetic analyses support two TR2/4 genes and multiple COUP-TF-like NRs in N. vectensis. Similarly, we observed 5 N. vectensis NRs that group with bilaterian NRs in the subfamily NR 2E. These genes are likely cnidarian-specific radiations within these NR 2 subfamilies and may represent interesting cases for better understanding the diversification of NR function within an early diverging metazoan.

After gene duplications, paralogs may be selected for divergent functions that may be reflected by different expression patterns [77, 78]. Previous reports of N. vectensis paralogs have suggested that duplicated genes have varying degrees of distinct spatial and temporal gene expression patterns (e.g., snail [79], Otx [80]). In some cases, NvNRs supported as likely lineage specific duplications showed divergent temporal gene expression patterns consistent with subfunctionalization. For example, NvNR1 was highly expressed in the larval and juvenile stages, and the closely related NvNR2 was most highly expressed during cleavage and in adults. In situ hybridization studies of N. vectensis will reveal spatial expression patterns (e.g., anterior-posterior gradients, neural expression), providing further insight into the potential subfuctionalization of duplicated NRs.


In this study, we have identified 17 nuclear receptors in the sea anemone Nematostella vectensis. Phylogenetic analyses indicate that a majority of these NRs are members of NR family 2 and are consistent with previous studies indicating that this family had diversified prior to the cnidarian-bilaterian split. We also identified a putative homolog to GCNF in NR family 6 as well as three NRs that do not group with any previously identified NR family. Expression of N. vectensis NRs varied during development, consistent with stage-specific functions related to development, metamorphosis and adult physiology. Future studies of NR function in N. vectensis will provide a critical understanding of the evolution the nuclear receptor superfamily in animals.



The starlet sea anemone, Nematostella vectensis, is a small burrowing anemone found in estuaries along the Atlantic Coast of the United States, as well as in populations in England and along the Pacific Coast of the United States. We collected adults from Great Sippewissett Marsh, MA USA and cultured them under standard conditions (0.45 μm filtered seawater diluted to 13 ppt, room temperature, fed twice weekly brine shrimp and once mussel ovary; similar to previously described conditions [81, 82]). Under these conditions, N. vectensis undergoes a predictable reproductive cycle that includes weekly spawning. Within 2-3 hours after spawning, embryos were gathered from dishes, rinsed with diluted seawater, and placed in sterile dishes. At specific stages (i.e., cleavage, early planula, late planula [pyramid shape], and juvenile, [81]) we collected individuals with a pipette, lightly centrifuged, decanted excess seawater, and snap froze with liquid nitrogen for later analysis. Adults were left unfed for at least one week prior to freezing to minimize contamination from food sources.

RNA extraction and cDNA synthesis

For cloning of nuclear receptors, total RNA was extracted from N. vectensis of varying developmental stages, ranging from embryos to adult. RNA was extracted by homogenization of whole individuals in STAT-60 (Tel-Test, Inc). RNA for cloning was then pooled from developmental stages. RNA purity and yield were quantified using a ND-1000 spectrophotometer (Nanodrop). RNA quality was visualized for representative samples on a denaturing agarose gel. Polyadenylated-enriched RNA (polyA-RNA) was prepared from total RNA using the MicroPoly(A)Purist Kit (Ambion). Complementary DNA (cDNA) was synthesized from polyA-RNA using the Iscript cDNA Synthesis Kit (Bio-Rad) using 1 μg of polyA-RNA per 20 μl reaction.

For characterization of developmental expression patterns, RNA was extracted using STAT-60 and the Aurum Total RNA Mini Kit with on-column DNAse digestion, as described previously [83]. Each extraction included multiple individuals from a single spawning event or in the case of adults, multiple individuals from the breeding stock. From the total RNA, cDNA was synthesized with the Iscript cDNA Synthesis Kit using 2 μg of RNA per 30 μl reaction. cDNA was prepared from the following developmental stages: cleavage (0.5 days post fertilization [dpf], n = 1), embryos (1 dpf, n = 3), early larvae (3-4 dpf, n = 4), late larvae (7-11 dpf, not metamorphosed, n = 4), juvenile (8-23 dpf, metamorphosed, n = 4), and adult (>100 dpf, n = 4).

Database searching and cloning

N. vectensis nuclear receptors were identified through a combination of BLASTp and domain searches through the Joint Genome Institute N. vectensis assembly. Genes were additionally queried against the EST database at JGI in order to annotate more complete transcripts. To confirm the predicted sequences we amplified, cloned, and sequenced large portions of each predicted NR transcript (primers designed with Primer3, primer sequences listed in Table S1). The PCR mixture consisted of 11.8 μl of molecular biology grade distilled water, 0.1 μl of AmpliTaq Gold (Applied Biosystems), 2 μl of the accompanying 10× Buffer, 3.75 mM MgCl2, 0.8 mM dNTPs, 0.5 μl of 10 μM gene-specific primers, and 0.5 μl of 1:5 diluted polyA cDNA. PCR conditions were as follows: 95°C for 5 min; 40 cycles of 95°C for 15 s, 60°C for 30 s, and 72°C for 60 s, followed by a 10 min extension at 72°C. Products were gel purified, ligated into pGEM-T Easy vector (Promega), and sequenced to confirm the targeted amplicon.

Phylogenetic analysis

We used a likelihood based approach to determine evolutionary relationships of the N. vectensis NRs with bilaterian NRs. NRs from Homo sapiens, Danio rerio, Xenopus laevis, Drosophila melanogaster, and Caenorhabditis elegans from Bertrand et al. [5] were used as representative bilaterians. We included sequences for these species from each of the defined families and subfamilies in order to fully represent the broad-scale diversity of NRs. We also included NR sequences from Ciona intestinalis reported in [51] and an estrogen receptor homolog from the mollusc Aplysia californica [84]. Additional reports of NR diversity from other animals [9, 67, 85, 86] have shown that these taxa fully represent the diversity of NRs to the sub-family level, despite lineage-specific losses in some ecdysozoans and gene-specific duplications in a variety of taxa. Full length sequences for all taxa were aligned with Muscle 3.6 [87] and edited manually in the case of clear errors. Maximum likelihood analyses were run using RAxML (version 7.0.4, [88]) with a JTT+G matrix (model determined by AIC criteria with ProtTest v1.4, [89]). Separate analyses were conducted using the DBD only, a portion of the LDB, and the DBD plus a portion of LBD. We also performed a complementary neighbor joining analysis of the DBD and LBD alignment (PHYLIP v3.6). Trees were visualized and illustrated with FigTree v1.1.2

Gene structure

We assembled the most complete transcripts for each of N. vectensis' NRs with a combination of cloned sequence, ESTs, and, where necessary, gene prediction models from the Joint Genome Institute (JGI). Intron-exon structure was determined by aligning the assembled transcripts to the most recent genome scaffolds in the JGI N. vectensis database. Gene structure was illustrated with GenePalette v1.21 [90].

Quantitative real-time RT-PCR (qPCR)

Oligonucleotide primers (see Additional file 8) were designed to amplify each N. vectensis NR, as well as glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and 18S ribosomal RNA (18S). Primers were 20-21 nt, with a GC content of 40-60%, either overlapped predicted exon-exon boundaries by 3-4 bp or spanned a large intron, and produced predicted amplicons of 55-146 bp with minimal predicted secondary structure (m-fold, [91]). A standard curve was constructed from serially-diluted plasmids containing the amplicon of interest. The standard curve was used in qPCR reactions to quantify amplification efficiency and to calculate the number of molecules per reaction (as in [92]). qPCR was performed using iQ SYBR Green Supermix (Bio-Rad), and reactions were run a MyCycler Real-Time PCR detection system (Bio-Rad). For each gene, standards were run in triplicate wells and experimental samples were run in duplicate wells (technical replicates) on a single plate. The PCR mixture consisted of 11.5 μl of molecular biology grade distilled water, 12.5 μl of IQ SYBR Green Supermix, 0.5 μl of 10 μM gene-specific primers, and 0.5 μl of cDNA. PCR conditions were as follows: 95°C for 3 min; 40 cycles of 95°C for 15 s and 64°C for 45 s. After 40 cycles, the PCR products from each reaction were subjected to melt curve analysis to ensure that only a single product was amplified. Selected reactions for each gene were visualized on 15% TBE gels (Bio-Rad) and consistently yielded single bands of the predicted size. The number of molecules per μl for each gene was calculated by comparing the threshold cycle (Ct) from the sample with the standard curve. Expression was compared among developmental stages using one-way analysis of variance (ANOVA) with Tukey's Honestly Significant Difference Test as a posthoc test (SYSTAT 12). Relative gene expression patterns were compared with principal component analysis (PCA, calculated with SYSTAT 12). This statistical approach is similar to that previously described by Tarrant and colleagues [83].


  1. 1.

    Larroux C, Fahey B, Liubicich D, Hinman V, Gauthier M, Gongora M, Green K, Wörheide G, Leys S, Degnan B: Developmental expression of transcription factor genes in a demosponge: insights into the origin of metazoan multicellularity. Evolution and Development. 2006, 8: 150-173. 10.1111/j.1525-142X.2006.00086.x.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Wiens M, Batel R, Korzhev M, Müller W: Retinoid × receptor and retinoic acid response in the marine sponge Suberites domuncula. Journal of Experimental Biology. 2003, 206: 3261-71. 10.1242/jeb.00541.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Gauchat D, Escriva H, Miljkovic-Licina M, Chera S, Langlois M-C, Begue A, Laudet V, Galliot B: The orphan COUP-TF nuclear receptors are markers for neurogenesis from cnidarians to vertebrates. Developmental Biology. 2004, 275: 104-23. 10.1016/j.ydbio.2004.07.037.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Grasso L, Hayward D, Trueman J, Hardie K, Janssens P, Ball E: The evolution of nuclear receptors: evidence from the coral Acropora. Molecular Phylogenetics and Evolution. 2001, 21: 93-102. 10.1006/mpev.2001.0994.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Bertrand S, Brunet F, Escriva H, Parmentier G, Laudet V, Robinson-Rechavi M: Evolutionary genomics of nuclear receptors: from twenty-five ancestral genes to derived endocrine systems. Molecular Biology and Evolution. 2004, 21: 1923-37. 10.1093/molbev/msh200.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    King N, Westbrook MJ, Young SL, Kuo A, Abedin M, Chapman J, Fairclough S, Hellsten U, Isogai Y, Letunic I, et al: The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans. Nature. 2008, 451: 783-788. 10.1038/nature06617.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  7. 7.

    Sluder A, Matthews S, Hough D, Yin V, Maina C: The nuclear receptor superfamily has undergone extensive proliferation and diversification in nematodes. Genome Research. 1999, 9: 103-20.

    CAS  PubMed  Google Scholar 

  8. 8.

    Keay J, Bridgham J, Thornton J: The Octopus vulgaris estrogen receptor is a constitutive transcriptional activator: evolutionary and functional implications. Endocrinology. 2006, 147: 3861-9. 10.1210/en.2006-0363.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Wu W, Niles EG, El-Sayed N, Berriman M, LoVerde PT: Schistosoma mansoni (Platyhelminthes, Trematoda) nuclear receptors: Sixteen new members and a novel subfamily. Gene. 2006, 366: 303-315. 10.1016/j.gene.2005.09.013.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Wu W, Niles EG, Hirai H, LoVerde PT: Evolution of a novel subfamily of nuclear receptors with members that each contain two DNA binding domains. BMC Evolutionary Biology. 2007, 7: 27-10.1186/1471-2148-7-27.

    PubMed Central  Article  PubMed  Google Scholar 

  11. 11.

    Kostrouch Z, Kostrouchova M, Love W, Jannini E, Piatigorsky J, Rall J: Retinoic acid × receptor in the diploblast, Tripedalia cystophora. Proceedings of the National Academy of Science, USA. 1998, 95: 13442-7. 10.1073/pnas.95.23.13442.

    CAS  Article  Google Scholar 

  12. 12.

    Tarrant A, Atkinson S, Atkinson M: Estrone and estradiol-17β concentration in tissue of the scleractinian coral, Montipora verrucosa. Comparative Biochemistry and Physiology A. 1999, 122:

    Google Scholar 

  13. 13.

    Twan W-H, Hwang J-S, Chang C-F: Sex steroids in scleractinian coral, Euphyllia ancora: implication in mass spawning. Biol Reprod. 2003, 68: 2255-60. 10.1095/biolreprod.102.012450.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Slattery M, Hines G, Starmer J, Paul V: Chemical signals in gametogenesis, spawning, and larval settlement and defense of the soft coral Sinularia polydactyla. Coral Reefs. 1999, 18: 75-84. 10.1007/s003380050158.

    Article  Google Scholar 

  15. 15.

    Pernet V, Anctil M: Annual variations and sex-related differences of estradiol-17beta levels in the anthozoan Renilla koellikeri. General and Comparative Endocrinology. 2002, 129: 63-8. 10.1016/S0016-6480(02)00513-0.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Eterovic V, Hann R, Ferchmin P, Rodriques A, Li L, Lee Y, McNamee M: Diterpenoids from Caribbean gorgonians act as noncompetitive inhibitors of the nicotinic acetylcholine receptor. Cell Mol Neurobiol. 1993, 13: 99-110. 10.1007/BF00735367.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Fleury B, Lages B, Barbosa J, Kaiser C, Pinto A: New hemiketal steroid from the introduced soft coral Chromonephthea braziliensis is a chemical defense against predatory fishes. Journal of Chemical Ecology. 2008, 34: 987-93. 10.1007/s10886-008-9499-y.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Koljak R, Järving I, Kurg R, Boeglin W, Varvas K, Valmsen K, Ustav M, Brash A, Samel N: The basis of prostaglandin synthesis in coral. Journal of Biological Chemistry. 2001, 276: 7033-40. 10.1074/jbc.M009803200.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Sullivan AA, Thummel CS: Temporal profiles of nuclear receptor gene expression reveal coordinate transcriptional responses during Drosophila development. Mol Endocrinol. 2003, 17: 2125-2137. 10.1210/me.2002-0430.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Putnam N, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov V, et al: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317: 86-94. 10.1126/science.1139158.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Thornton JW: Nuclear receptor diversity: phylogeny, evolution and endocrine disruption. Pure and Applied Chemistry. 2003, 75: 1827-1839. 10.1351/pac200375111827.

    CAS  Article  Google Scholar 

  22. 22.

    Maglich J, Sluder A, Guan X, Shi Y, McKee D, Carrick K, Kamdar K, Willson T, Moore J: Comparison of the complete nuclear receptor sets from the human, Caenorhabditis elegans, and Drosophila genomes. Genome Biology. 2001, 2: 1-7. 10.1186/gb-2001-2-8-research0029.

    Article  Google Scholar 

  23. 23.

    Sladek F, Zhong W, Lai E, Darnell J: Liver-enriched transcription factors HNF-4 is a novel member of the steroid hormone receptor superfamily. Genes Dev. 1990, 4: 2353-65. 10.1101/gad.4.12b.2353.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Dhe-Paganon S, Duda K, Iwamoto M, Chi Y, Shoelson S: Crystal structure of the HNF4 alpha ligand binding domain in complex with endogenous fatty acid ligand. Journal of Biological Chemistry. 2002, 277: 37973-6. 10.1074/jbc.C200420200.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Wisely G, Miller A, Davis R, Thornquest AJ, Johnson R, Spitzer T, Sefler A, Shearer B, Moore J, Miller A, et al: Hepatocyte nuclear factor 4 is a transcription factor that constitutively binds fatty acids. Structure. 2002, 10: 1225-34. 10.1016/S0969-2126(02)00829-8.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Watt A, Garrison W, Duncan S: HNF4: a central regulator of hepatocyte differentiation and function. Hepatology. 2003, 37: 1249-53. 10.1053/jhep.2003.50273.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Zhong W, Sladek F, Darnell J: The expression patttern of a Drosophila homolog to the mouse transcription factor HNF-4 suggests a determinative role in gut formation. The EMBO Journal. 1993, 12: 537-44.

    PubMed Central  CAS  PubMed  Google Scholar 

  28. 28.

    Robinson-Rechavi M, Maina C, Gissendanne C, Laudet V, Sluder A: Explosive lineage-specific expansion of the orphan nuclear receptor HNF4 in nematodes. Journal of Molecular Evolution. 2004, 60: 577-86. 10.1007/s00239-004-0175-8.

    Article  Google Scholar 

  29. 29.

    Van Gilst MR, Hadjivassiliou H, Jolly A, Yamamoto KR: Nuclear hormone receptor NHR-49 controls fat consumption, storage, and composition in C. elegans. PLoS Biol. 2004, 3: e53-10.1371/journal.pbio.0030053.

    Article  Google Scholar 

  30. 30.

    Tarrant A, Cortés J, Atkinson M, Atkinson S, Johanning K, Chiang T-C, Vargas J, McLachlan J: Three orphan nuclear receptors in the scleractinian coral Pocillopora damicornis. Revista de Biologia Tropical. 56 (Suppl 4): 39-48.

  31. 31.

    Stenman J, Yu R, Evans R, Campbell K: Tlx and Pax6 co-operate genetically to establish the pallio-subpallial boundary in the embryonic mouse telencephalon. Development. 2003, 130: 1113-22. 10.1242/dev.00328.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Yu R, Chiang M-Y, Tanabe T, Kobayashi M, Yasuda K, Evans R, Umesono K: The orphan nuclear receptors Tlx regulates Pax2 and is essential for vision. Proceedings of the National Academy of Science, USA. 2000, 97: 2621-5. 10.1073/pnas.050566897.

    CAS  Article  Google Scholar 

  33. 33.

    Kobayashi M, Takezawa S, Hara K, Yu R, Umesono K, Agata K, Taniwaki M, Yasuda K, Umesono K: Identification of a photoreceptor cell-specific nuclear receptor. Proceedings of the National Academy of Science, USA. 1999, 96: 4814-9. 10.1073/pnas.96.9.4814.

    CAS  Article  Google Scholar 

  34. 34.

    Roy K, Thiels E, Monaghan A: Loss of the tailless gene affects forebrain development and emotional behavior. Physiol Behav. 2002, 77: 595-600. 10.1016/S0031-9384(02)00902-2.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  35. 35.

    Steingrimsson E, Pignoni F, Liaw G, Lengyel J: Dual role of the Drosophila pattern gene tailless in embryonic termini. Science. 1991, 254: 418-21. 10.1126/science.1925599.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Much JW, Slade DJ, Klampert K, Garriga G, Wightman B: The fax-1 nuclear hormone receptor regulates axon pathfinding and neurotransmitter expression. Development. 2000, 127: 703-712.

    CAS  PubMed  Google Scholar 

  37. 37.

    Wightman B, Ebert B, Carmean N, Weber K, Clever S: The C. elegans nuclear receptor gene fax-1 and homeobox gene unc-42 coordinate interneuron identity by regulating the expression of glutamate receptor subunits and other neuron-specific genes. Developmental Biology. 2005, 287: 74-85. 10.1016/j.ydbio.2005.08.032.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Kobayashi M, Hara K, Yu RT, Yasuda K: Expression and functional analysis of Nr2e3, a photoreceptor-specific nuclear receptor, suggest common mechanisms in retinal development between avians and mammals. Development Genes and Evolution. 2008, 218: 10.1007/s00427-008-0232-1.

    Google Scholar 

  39. 39.

    Matus D, Thomsen G, Martindale M: FGF signaling in gastrulation and neural development in Nematostella vectensis, an anthozoan cnidarian. Development Genes and Evolution. 2007, 217: 137-148. 10.1007/s00427-006-0122-3.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  40. 40.

    Escriva H, Safi R, Hanni C, Langlois M-C, Saumitou-Laprade P, Stehelin D, Capron A, Pierce R, Laudet V: Ligand binding was acquired during the evolution of nuclear receptors. Proceedings of the National Academy of Science, USA. 1997, 94: 6803-8. 10.1073/pnas.94.13.6803.

    CAS  Article  Google Scholar 

  41. 41.

    Langlois M-C, Vanacker J-M, Holland N, Escriva H, Queva C, Laudet V, Holland L: Amphicoup-TF, a nuclear orphan receptor of the lancelet Branchiostoma floridae, is implicated in retinoic acid signalling pathways. Development, Genes and Evolution. 2000, 210: 471-82. 10.1007/s004270000087.

    CAS  Article  Google Scholar 

  42. 42.

    Chan S-M, Xu N, Niemeyer C, Bone J, Flytzanis C: SpCOUP-TF: a sea urchin member of the steroid/thyroid hormone receptor family. Proceedings of the National Academy of Science, USA. 1992, 89: 10568-72. 10.1073/pnas.89.22.10568.

    CAS  Article  Google Scholar 

  43. 43.

    Anderson GW, Larson RJ, Oas DR, Sandhofer CR, Schwartz HL, Mariash CN, Oppenheimer JH: Chicken ovalbumin upstream promoter-transcription factor (COUP-TF) modulates expression of the Purkinje cell protein-2 gene. A potential role for COUP-TF in repressing premature thyroid hormone action in the developing brain. J Biol Chem. 1998, 273: 16391-16399. 10.1074/jbc.273.26.16391.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Yamamoto M, Fujinuma M, Tanaka M, Dräger UC, McCaffery P: Sagittal band expression of COUP-TF2 gene in the developing cerebellum. Mechanisms of Development. 1999, 84: 143-146. 10.1016/S0925-4773(99)00054-4.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Kerber B, Fellert S, Hoch M: Seven-up, the Drosophila homolog of the COUP-TF orphan receptors, controls cell proliferation in the insect kidney. Genes Dev. 1998, 12: 1781-1786. 10.1101/gad.12.12.1781.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  46. 46.

    Fjose A, Nornes S, Weber U, Mlodzik M: Functional conservation of vertebrate seven-up related genes in neurogenesis and eye development. The EMBO Journal. 1993, 12: 1403-14.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. 47.

    Liu Y, Yang N, Teng CT: COUP-TF acts as a competitive repressor for estrogen receptor-mediated activation of the mouse lactoferrin gene. Mol Cell Biol. 1993, 13: 1836-1846.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  48. 48.

    Zhu J, Miuma K, Chen L, Raikhel A: Cyclicity of mosquito vitellogenic ecdysteroid-mediated signaling is modulated by alternative dimerization of the RXR homolog Ultraspiraclei. Proceedings of the National Academy of Science, USA. 2003, 100: 544-9. 10.1073/pnas.0235695100.

    CAS  Article  Google Scholar 

  49. 49.

    Miura K, Zhu J, Dittmer N, Chen L, Raikhel A: A COUP-TF/Svp homolog is highly expressed during vitellogenesis in the mosquito Aedes aegypti. Journal of Molecular Endocrinology. 2002, 29: 223-38. 10.1677/jme.0.0290223.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Klinge CM: Role of estrogen receptor ligand and estrogen response element sequence on interaction with chicken ovalbumin upstream promoter transcription factor (COUP-TF). The Journal of Steroid Biochemistry and Molecular Biology. 1999, 71: 1-19. 10.1016/S0960-0760(99)00124-7.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Devine C, Hinman VF, Degnan BM: Evolution and developmental expression of nuclear receptor genes in the ascidian Herdmania. International Journal of Developmental Biology. 2002, 46: 687-692.

    CAS  PubMed  Google Scholar 

  52. 52.

    Lee Y-F, Lee H-J, Chang C: Recent advances in the TR2 and TR4 orphan receptors of the nuclear receptor superfamily. The Journal of Steroid Biochemistry and Molecular Biology. 2002, 81: 291-308. 10.1016/S0960-0760(02)00118-8.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Hu Y, Shyr C, Che W, Mu X, Kim E, Chang C: Suppression of estrogen receptor-mediated transcription and cell growth by interaction with TR2 orphan receptor. Journal of Biological Chemistry. 2002, 277: 33571-9. 10.1074/jbc.M203531200.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Yagi K, Satou Y, Mazet F, Shimeld S, Degnan B, Rokhsar D, Levine M, Kohara Y, Satoh N: A genome-wide survey of developmental genes in Ciona intestinalis III. Genes for Fox, ETS, nuclear receptors and NF-kB. Development, Genes and Evolution. 2003, 213: 235-44. 10.1007/s00427-003-0322-z.

    CAS  Article  Google Scholar 

  55. 55.

    Tanabe O, Shen Y, Liu Q, Campbell AD, Kuroha T, Yamamoto M, Engel JD: The TR2 and TR4 orphan nuclear receptors repress Gata1 transcription. Genes Dev. 2007, 21: 2832-2844. 10.1101/gad.1593307.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  56. 56.

    Zhang Y, Dufau ML, Gerald L: Gene Silencing by Nuclear Orphan Receptors. Vitamins & Hormones. 2004, Academic Press, 68: 1-48.

    Google Scholar 

  57. 57.

    Lee Y-F, Lee H-J, Chang C: Recent advances in the TR2 and TR4 orphan receptors of the nuclear receptor superfamily. Journal of Steroid Biochemistry and Molecular Biology. 2002, 81: 291-308. 10.1016/S0960-0760(02)00118-8.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Hu R, Wu W, Niles EG, LoVerde PT: SmTR2/4, a Schistosoma mansoni homologue of TR2/TR4 orphan nuclear receptor. International Journal for Parasitology. 2006, 36: 1113-1122. 10.1016/j.ijpara.2006.06.003.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Sullivan A, Thummel C: Temporal profiles of nuclear receptor gene expression reveal coordinate transcriptional responses during Drosophila development. Molecular Endocrinology. 2003, 17: 2125-37. 10.1210/me.2002-0430.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    King-Jones K, Charles J-P, Lam G, Thummel C: The ecdysone-induced DHR4 orphan nuclear receptor coordinates growth and maturation in Drosophila. Cell. 2005, 121: 773-84. 10.1016/j.cell.2005.03.030.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Zhao Z, Fang L, Johnsen R, Baillie D: ATP-binding cassette protein E is involved in gene transcription and translation in Caenorhabditis elegans. Biochemical and Biophysical Research Communications. 2004, 323: 104-111. 10.1016/j.bbrc.2004.08.068.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Lan Z-J, Gu P, Xu X, Cooney A: Expression of the orphan nuclear receptor, germ cell nuclear factor in mouse gonads and preimplantation embryos. Biology of Reproduction. 2003, 68: 282-9. 10.1095/biolreprod.102.008151.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Chung ACK, Katz D, Pereira FA, Jackson KJ, DeMayo FJ, Cooney AJ, O'Malley BW: Loss of orphan receptor Germ Cell Nuclear Factor function results in ectopic development of the tail bud and a novel posterior truncation. Mol Cell Biol. 2001, 21: 663-677. 10.1128/MCB.21.2.663-677.2001.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  64. 64.

    Borgmeyer U: Dimeric binding of the mouse germ cell nuclear factor. European Journal of Biochemistry. 1997, 244: 120-7. 10.1111/j.1432-1033.1997.00120.x.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Cooney A, Hummelke G, Herman T, Chen F, Jackson K: Germ cell nuclear factor is a response element-specfic repressor of transcription. Biochemical and Biophysical Research Communications. 1998, 245: 94-100. 10.1006/bbrc.1998.8391.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Keay J, Bridgham JT, Thornton JW: The Octopus vulgaris estrogen receptor is a constitutive transcriptional activator: evolutionary and functional implications. Endocrinology. 2006, 147: 3861-3869. 10.1210/en.2006-0363.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Keay J, Thornton JW: Hormone-Activated Estrogen Receptors in Annelid Invertebrates: Implications for Evolution and Endocrine Disruption. Endocrinology. 2009, 150: 1731-1738. 10.1210/en.2008-1338.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  68. 68.

    Baker ME: Trichoplax, the simplest known animal, contains an estrogen-related receptor but no estrogen receptor: Implications for estrogen receptor evolution. Biochemical and Biophysical Research Communications. 2008, 375: 623-627. 10.1016/j.bbrc.2008.08.047.

    CAS  Article  PubMed  Google Scholar 

  69. 69.

    Atkinson S, Atkinson M: Detection of estradiol-17β during a mass coral spawn. Coral Reefs. 1992, 11: 33-5. 10.1007/BF00291932.

    Article  Google Scholar 

  70. 70.

    Tarrant A, Atkinson M, Atkinson S: Effects of steroidal estrogens on coral growth and reproduction. Marine Ecology Progress Series. 2004, 269: 121-9. 10.3354/meps269121.

    CAS  Article  Google Scholar 

  71. 71.

    Bouzaiene M, Angers A, Anctil M: Immunohistochemical localization of a retinoic acid-like receptor in nerve cells of two colonial anthozoans (Cnidaria). Tissue and Cell. 2007, 39: 123-130. 10.1016/j.tice.2007.02.001.

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    Chawla A, Repa J, Evans R, Mangelsdorf D: Nuclear receptors and lipid physiology: opening the X-files. Science. 2001, 294: 1866-70. 10.1126/science.294.5548.1866.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Mangelsdorf D, Evans R: The RXR heterodimers and orphan receptors. Cell. 1995, 83: 841-50. 10.1016/0092-8674(95)90200-7.

    CAS  Article  PubMed  Google Scholar 

  74. 74.

    Amoutzias G, Pichler E, Mian N, De Graaf D, Imsiridou A, Robinson-Rechavi M, Bornberg-Bauer E, Robertson D, Oliver S: A protein interaction atlas for the nuclear receptors: properties and quality of a hub-based dimerisation network. BMC Systems Biology. 2007, 1: 34-10.1186/1752-0509-1-34.

    PubMed Central  Article  PubMed  Google Scholar 

  75. 75.

    Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, Vidalain P-O, Han J-DJ, Chesneau A, Hao T, et al: A map of the interactome network of the metazoan C. elegans. Science. 2004, 303: 540-543. 10.1126/science.1091403.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  76. 76.

    Antebi A: Nuclear hormone receptors in C. elegans. 2006, Wormbook, The C. elegans Research Community, doi/10.1895/wormbook.1.64.1, []

    Google Scholar 

  77. 77.

    Ohno S: Evolution by Gene Duplication. 1970, Heidelberg Springer-Verlag

    Chapter  Google Scholar 

  78. 78.

    Zhang J: Evolution by gene duplication: an update. Trends in Ecology & Evolution. 2003, 18: 292-298.

    Article  Google Scholar 

  79. 79.

    Martindale MQ, Pang K, Finnerty JR: Investigating the origins of triploblasty: 'mesodermal' gene expression in a diploblastic animal, the sea anemone Nematostella vectensis (phylum, Cnidaria; class, Anthozoa). Development. 2004, 131: 2463-2474.

    CAS  Article  PubMed  Google Scholar 

  80. 80.

    Mazza ME, Pang K, Martindale MQ, Finnerty JR: Genomic organization, gene structure, and developmental expression of three Clustered otx genes in the sea anemone Nematostella vectensis. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution. 2007, 308B: 494-506.

    CAS  Article  Google Scholar 

  81. 81.

    Reitzel AM, Burton PM, Krone C, Finnerty JR: Comparison of developmental trajectories in the starlet sea anemone Nematostella vectensis: embryogenesis, regeneration, and two forms of asexual fission. Invertebrate Biology. 2007, 126: 99-112.

    Article  Google Scholar 

  82. 82.

    Fritzenwanker J, Technau U: Induction of gametogenesis in the basal cnidarian Nematostella vectensis (Anthozoa). Development, Genes and Evolution. 2002, 212: 99-103.

    Article  Google Scholar 

  83. 83.

    Tarrant A, Baumgartner M, Verslycke T, Johnson C: Differential gene expression in diapausing and active Calanus finmarchicus (Copepoda). Marine Ecology Progress Series. 2008, 355: 193-207.

    CAS  Article  Google Scholar 

  84. 84.

    Thornton JW, Need E, Crews D: Resurrecting the ancestral steroid receptor: ancient origin of estrogen signaling. Science. 2003, 301: 1714-1717.

    CAS  Article  PubMed  Google Scholar 

  85. 85.

    Maglich JM, Caravella JA, Lambert MH, Willson TM, Moore JT, Ramamurthy L: The first completed genome sequence from a teleost fish (Fugu rubripes) adds significant diversity to the nuclear receptor superfamily. Nucl Acids Res. 2003, 31: 4051-4058.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  86. 86.

    Schubert M, Brunet F, Paris M, Bertrand S, Benoit G, Laudet V: Nuclear hormone receptor signaling in amphioxus. Development Genes and Evolution. 2008, 218: 651-665.

    CAS  Article  PubMed  Google Scholar 

  87. 87.

    Edgar R: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32: 1792-7.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  88. 88.

    Stamatakis A: RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-90.

    CAS  Article  PubMed  Google Scholar 

  89. 89.

    Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-5.

    CAS  Article  PubMed  Google Scholar 

  90. 90.

    Rebeiz M, Posakony J: GenePalette: a universal software tool for genome sequence visualization and analysis. Developmental Biology. 2004, 271: 431-8.

    CAS  Article  PubMed  Google Scholar 

  91. 91.

    Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31: 3406-3415.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  92. 92.

    Tarrant A, Greytak S, Callard G, Hahn M: Estrogen receptor-related receptors in the killifish Fundulus heteroclitus: diversity, expression, and estrogen responsiveness. Journal of Molecular Endocrinology. 2006, 37: 105-20.

    CAS  Article  PubMed  Google Scholar 

Download references


We are grateful for financial support from the Woods Hole Oceanographic Institution (WHOI) through the Tropical Research Initiative, the Ocean Life Institute (AMT), the Academic Programs Office, and to the Beacon Institute for Rivers and Estuaries (AMR). We thank Natasha Rabinowitz for assistance with animal care. We would also like to thank Mark Hahn, Jed Goldstone, John McLachlan, and Tung-chin Chiang for assistance in formulating ideas and approaches during this study.

Author information



Corresponding author

Correspondence to Ann M Tarrant.

Additional information

Authors' contributions

AMR conducted bioinformatic, gene structure, and phylogenetic analysis, designed primers, cloned nuclear receptor fragments, and drafted the manuscript. AMT conceived of the study, conducted qPCR assays and statistical analysis, and contributed to the draft. The authors together interpreted the data and revised the manuscript. Both authors read and approved of the final manuscript.

Electronic supplementary material


Additional file 1: N. vectensisdeduced amino acid sequences used for phylogenetic analyses. List of amino acid sequences. (DOC 28 KB)

Maximum likelihood tree using the DBD and portion of the LBD

Additional file 2: . NRs were taken from N. vectensis, Homo sapiens, Danio rerio, Xenopus laevis, Drosophila melanogaster, Caenorhabditis elegans, Aplysia californica, and Ciona intestinalis. Tree was rooted with the NR 2 family. Values above branches indicate percent of 1000 bootstraps. (EPS 2 MB)


Additional file 3: Maximum likelihood tree using only the DBD portion of the NR alignment. Tree was rooted with the NR 2 family. Values above branches indicate percent of 1000 bootstraps. (EPS 2 MB)


Additional file 4: Maximum likelihood tree using only the alignable portions of the LBD of the NR alignment. Tree was rooted with the NR 2 family. Values above branches indicate percent of 1000 bootstraps. (EPS 2 MB)

Neighbor joining tree constructed using the DBD and portion of LDB of representative NRs

Additional file 5: . NRs from N. vectensis, Homo sapiens, Danio rerio, Xenopus laevis, Drosophila melanogaster, Caenorhabditis elegans, Aplysia californica, and Ciona intestinalis. Tree was rooted with the NR 2 family. Values above branches indicate percent of 1000 bootstraps. (EPS 2 MB)


Additional file 6: Intron-exon structure for seventeen N. vectensisNRs. White boxes indicate exons with coding sequence, gray boxes are exons or portions of exons containing the untranslated regions, and lines between boxes indicate introns. (EPS 993 KB)


Additional file 7: Exon boundaries (vertical bars) for two sets of inferred cnidarian-specific NR radiations. (A) NvNR8 and 9, TLL-like genes, share 2 of the 3 exon positions. We could not compare the first boundary in NvNR9 with NvNR8 because we do not have corresponding sequence data for this portion of the protein. (B) Four NvNRs inferred to be duplications related to COUP-TF. NvNR12, 13 and 14 each has 3 exons with two conserved splice sites in all three genes. These two positions are also conserved in NvNR11, but this gene has two additional exon boundaries without conservation in the other three NRs from this subfamily. (EPS 2 MB)


Additional file 8: Primer sequences for amplifying pieces of NvNRs for cloning and for qPCR. Table of primer sequences. (DOC 36 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Reitzel, A.M., Tarrant, A.M. Nuclear receptor complement of the cnidarian Nematostella vectensis: phylogenetic relationships and developmental expression patterns. BMC Evol Biol 9, 230 (2009).

Download citation


  • Nuclear Receptor
  • Nematostella Vectensis
  • Nuclear Receptor Family
  • Temporal Gene Expression Pattern
  • Germ Cell Nuclear Factor