TEMPERATURE AND SALINITY AFFECT DEVELOPMENT OF THE PARASITIC SEA ANEMONE EDWARDSIELLA LINEATA POTENTIALLY LIMITING ITS IMPACT AS A BIOLOGICAL CONTROL ON THE CTENOPHORE MNEMIOPSIS LEIDYI.

The lined sea anemone, Edwardsiella lineata, parasitizes the ctenophore Mnemiopsis leidyi, which is one of the most destructive marine invasive species in the world. Mnemiopsis leidyi is known to tolerate a wide range of environmental conditions. However, the environmental tolerances of its most prominent parasite have never been characterized. Here we determined the effects of temperature (18, 22, 26, and 30 C) and salinity (6, 15, 24, and 33 ppt) on the survival and development of E. lineata from a vermiform parasite to a free-living polyp. At higher temperatures and lower salinities, E. lineata experienced significantly higher mortality, and it failed to develop into an adult polyp at the highest temperature (30 C) and lowest salinities we tested (6 ppt or 15 ppt). While such temperature and salinity restrictions would not currently prevent E. lineata from infecting M. leidyi in many of the European waters where it has become a destructive invasive species, these environmental limitations may be reducing overlap between host and parasite within the host's native range, a situation that could be exacerbated by climate change.

Rampant asexual reproduction and limited dispersal in a mangrove population of the coral Porites divaricata.

Corals are critical to marine biodiversity. Reproduction and dispersal are key to their resilience, but rarely quantified in nature. Exploiting a unique system-a fully censused, longitudinally characterized, semi-isolated population inhabiting mangroves-we used 2bRAD sequencing to demonstrate that rampant asexual reproduction most likely via parthenogenesis and limited dispersal enable the persistence of a natural population of thin-finger coral (Porites divaricata). Unlike previous studies on coral dispersal, knowledge of colony age and location enabled us to identify plausible parent-offspring relationships within multiple clonal lineages and develop tightly constrained estimates of larval dispersal; the best-fitting model indicates dispersal is largely limited to a few metres from parent colonies. Our results explain why this species is adept at colonizing mangroves but suggest limited genetic diversity in mangrove populations and limited connectivity between mangroves and nearby reefs. As P. divaricata is gonochoristic, and parthenogenesis would be restricted to females (whereas fragmentation, which is presumably common in reef and seagrass habitats, is not), mangrove populations likely exhibit skewed sex ratios. These findings suggest that coral reproductive diversity can lead to distinctly different demographic outcomes in different habitats. Thus, coral conservation will require the protection of the entire coral habitat mosaic, and not just reefs.

Distinct Phenotypes Associated with Mangrove and Lagoon Habitats in Two Widespread Caribbean Corals, Porites astreoides and Porites divaricata.

AbstractAs coral reefs experience dramatic declines in coral cover throughout the tropics, there is an urgent need to understand the role that non-reef habitats, such as mangroves, play in the ecological niche of corals. Mangrove habitats present a challenge to reef-dwelling corals because they can differ dramatically from adjacent reef habitats with respect to key environmental parameters, such as light. Because variation in light within reef habitats is known to drive intraspecific differences in coral phenotype, we hypothesized that coral species that can exploit both reef and mangrove habitats will exhibit predictable differences in phenotypes between habitats. To investigate how intraspecific variation, driven by either local adaptation or phenotypic plasticity, might enable particular coral species to exploit these two qualitatively different habitat types, we compared the phenotypes of two widespread Caribbean corals, Porites divaricata and Porites astreoides, in mangrove versus lagoon habitats on Turneffe Atoll, Belize. We document significant differences in colony size, color, structural complexity, and corallite morphology between habitats. In every instance, the phenotypic differences between mangrove prop root and lagoon corals exhibited consistent trends in both P. divaricata and P. astreoides. We believe this study is the first to document intraspecific phenotypic diversity in corals occupying mangrove prop root versus lagoonal patch reef habitats. A difference in the capacity to adopt an alternative phenotype that is well suited to the mangrove habitat may explain why some reef coral species can exploit mangroves, while others cannot.

The impact of autotrophic versus heterotrophic nutritional pathways on colony health and wound recovery in corals.

For animals that harbor photosynthetic symbionts within their tissues, such as corals, the different relative contributions of autotrophy versus heterotrophy to organismal energetic requirements have direct impacts on fitness. This is especially true for facultatively symbiotic corals, where the balance between host-caught and symbiont-produced energy can be altered substantially to meet the variable demands of a shifting environment. In this study, we utilized a temperate coral-algal system (the northern star coral, Astrangia poculata, and its photosynthetic endosymbiont, Symbiodinium psygmophilum) to explore the impacts of nutritional sourcing on the host's health and ability to regenerate experimentally excised polyps. For fed and starved colonies, wound healing and total colony tissue cover were differentially impacted by heterotrophy versus autotrophy. There was an additive impact of positive nutritional and symbiotic states on a coral's ability to initiate healing, but a greater influence of symbiont state on the recovery of lost tissue at the lesion site and complete polyp regeneration. On the other hand, regardless of symbiont state, fed corals maintained a higher overall colony tissue cover, which also enabled more active host behavior (polyp extension) and endosymbiont behavior (photosynthetic ability of Symbiondinium). Overall, we determined that the impact of nutritional state and symbiotic state varied between biological functions, suggesting a diversity in energetic sourcing for each of these processes.

Intraspecific variation in oxidative stress tolerance in a model cnidarian: Differences in peroxide sensitivity between and within populations of Nematostella vectensis.

Nematostella vectensis is a member of the phylum Cnidaria, a lineage that includes anemones, corals, hydras, and jellyfishes. This estuarine anemone is an excellent model system for investigating the evolution of stress tolerance because it is easy to collect in its natural habitat and to culture in the laboratory, and it has a sequenced genome. Additionally, there is evidence of local adaptation to environmental stress in different N. vectensis populations, and abundant protein-coding polymorphisms have been identified, including polymorphisms in proteins that are implicated in stress responses. N. vectensis can tolerate a wide range of environmental parameters, and has recently been shown to have substantial intraspecific variation in temperature preference. We investigated whether different clonal lines of anemones also exhibit differential tolerance to oxidative stress. N. vectensis populations are continually exposed to reactive oxygen species (ROS) generated during cellular metabolism and by other environmental factors. Fifteen clonal lines of N. vectensis collected from four different estuaries were exposed to hydrogen peroxide. Pronounced differences in survival and regeneration were apparent between clonal lines collected from Meadowlands, NJ, Baruch, SC, and Kingsport, NS, as well as among 12 clonal lines collected from a single Cape Cod marsh. To our knowledge, this is the first example of intraspecific variability in oxidative stress resistance in cnidarians or in any marine animal. As oxidative stress often accompanies heat stress in marine organisms, resistance to oxidative stress could strongly influence survival in warming oceans. For example, while elevated temperatures trigger bleaching in corals, oxidative stress is thought to be the proximal trigger of bleaching at the cellular level.

Transcription factor NF-κB is modulated by symbiotic status in a sea anemone model of cnidarian bleaching.

Transcription factor NF-κB plays a central role in immunity from fruit flies to humans, and NF-κB activity is altered in many human diseases. To investigate a role for NF-κB in immunity and disease on a broader evolutionary scale we have characterized NF-κB in a sea anemone (Exaiptasia pallida; called Aiptasia herein) model for cnidarian symbiosis and dysbiosis (i.e., "bleaching"). We show that the DNA-binding site specificity of Aiptasia NF-κB is similar to NF-κB proteins from a broad expanse of organisms. Analyses of NF-κB and IκB kinase proteins from Aiptasia suggest that non-canonical NF-κB processing is an evolutionarily ancient pathway, which can be reconstituted in human cells. In Aiptasia, NF-κB protein levels, DNA-binding activity, and tissue expression increase when loss of the algal symbiont Symbiodinium is induced by heat or chemical treatment. Kinetic analysis of NF-κB levels following loss of symbiosis show that NF-κB levels increase only after Symbiodinium is cleared. Moreover, introduction of Symbiodinium into naïve Aiptasia larvae results in a decrease in NF-κB expression. Our results suggest that Symbiodinium suppresses NF-κB in order to enable establishment of symbiosis in Aiptasia. These results are the first to demonstrate a link between changes in the conserved immune regulatory protein NF-κB and cnidarian symbiotic status.

Sex-specific and developmental expression of Dmrt genes in the starlet sea anemone, Nematostella vectensis.

BACKGROUND: The molecular mechanisms underlying sex determination and differentiation in animals are incredibly diverse. The Dmrt (doublesex and mab-3 related transcription factor) gene family is an evolutionary ancient group of transcription factors dating to the ancestor of metazoans that are, in part, involved in sex determination and differentiation in numerous bilaterian animals and thus represents a potentially conserved mechanism for differentiating males and females dating to the protostome-deuterostome ancestor. Recently, the diversity of this gene family throughout animals has been described, but the expression and potential function for Dmrt genes is not well understood outside the bilaterians. RESULTS: Here, we report sex- and developmental-specific expression of all 11 Dmrts in the starlet sea anemone Nematostella vectensis. Nine out of the eleven Dmrts showed significant differences in developmental expression, with the highest expression typically in the adult stage and, in some cases, with little or no expression measured during embryogenesis. When expression was compared in females and males, seven of the eleven Dmrt genes had significant differences in expression with higher expression in males than in females for six of the genes. Lastly, expressions of two Dmrt genes with differential expression in each sex are located in the mesenteries and into the pharynx in polyps. CONCLUSIONS: Our results show that the phylogenetic diversity of Dmrt genes in N. vectensis is matched by an equally diverse pattern of expression during development and in each sex. This dynamic expression suggests multiple functions for Dmrt genes likely present in early diverging metazoans. Detailed functional analyses of individual genes will inform hypotheses regarding the antiquity of function for these transcription factors.

Methods for analyzing the evolutionary relationship of NF-κB proteins using free, web-driven bioinformatics and phylogenetic tools.

Phylogenetic analysis enables one to reconstruct the functional evolution of proteins. Current understanding of NF-κB signaling derives primarily from studies of a relatively small number of laboratory models-mainly vertebrates and insects-that represent a tiny fraction of animal evolution. As such, NF-κB has been the subject of limited phylogenetic analysis. The recent discovery of NF-κB proteins in "basal" marine animals (e.g., sponges, sea anemones, corals) and NF-κB-like proteins in non-metazoan lineages extends the origin of NF-κB signaling by several hundred million years and provides the opportunity to investigate the early evolution of this pathway using phylogenetic approaches. Here, we describe a combination of bioinformatic and phylogenetic analyses based on menu-driven, open-source computer programs that are readily accessible to molecular biologists without formal training in phylogenetic methods. These phylogenetically based comparisons of NF-κB proteins are powerful in that they reveal deep conservation and repeated instances of parallel evolution in the sequence and structure of NF-κB in distant animal groups, which suggest that important functional constraints limit the evolution of this protein.

Production of a reference transcriptome and transcriptomic database (EdwardsiellaBase) for the lined sea anemone, Edwardsiella lineata, a parasitic cnidarian.

BACKGROUND: The lined sea anemone Edwardsiella lineata is an informative model system for evolutionary-developmental studies of parasitism. In this species, it is possible to compare alternate developmental pathways leading from a larva to either a free-living polyp or a vermiform parasite that inhabits the mesoglea of a ctenophore host. Additionally, E. lineata is confamilial with the model cnidarian Nematostella vectensis, providing an opportunity for comparative genomic, molecular and organismal studies. DESCRIPTION: We generated a reference transcriptome for E. lineata via high-throughput sequencing of RNA isolated from five developmental stages (parasite; parasite-to-larva transition; larva; larva-to-adult transition; adult). The transcriptome comprises 90,440 contigs assembled from >15 billion nucleotides of DNA sequence. Using a molecular clock approach, we estimated the divergence between E. lineata and N. vectensis at 215-364 million years ago. Based on gene ontology and metabolic pathway analyses and gene family surveys (bHLH-PAS, deiodinases, Fox genes, LIM homeodomains, minicollagens, nuclear receptors, Sox genes, and Wnts), the transcriptome of E. lineata is comparable in depth and completeness to N. vectensis. Analyses of protein motifs and revealed extensive conservation between the proteins of these two edwardsiid anemones, although we show the NF-κB protein of E. lineata reflects the ancestral structure, while the NF-κB protein of N. vectensis has undergone a split that separates the DNA-binding domain from the inhibitory domain. All contigs have been deposited in a public database (EdwardsiellaBase), where they may be searched according to contig ID, gene ontology, protein family motif (Pfam), enzyme commission number, and BLAST. The alignment of the raw reads to the contigs can also be visualized via JBrowse. CONCLUSIONS: The transcriptomic data and database described here provide a platform for studying the evolutionary developmental genomics of a derived parasitic life cycle. In addition, these data from E. lineata will aid in the interpretation of evolutionary novelties in gene sequence or structure that have been reported for the model cnidarian N. vectensis (e.g., the split NF-κB locus). Finally, we include custom computational tools to facilitate the annotation of a transcriptome based on high-throughput sequencing data obtained from a "non-model system."

Isolation of DNA, RNA and protein from the starlet sea anemone Nematostella vectensis.

Among marine invertebrates, the starlet sea anemone Nematostella vectensis has emerged as an important laboratory model system. One advantage of working with this species relative to many other marine invertebrates is the ease of isolating relatively pure DNA, RNA and protein. Nematostella can be raised at high densities, under clean culture conditions, and it lacks integumentary or skeletal structures that can impede the recovery of DNA, RNA or protein. Here we describe methods used in our lab to isolate DNA, RNA and protein from Nematostella embryos, larvae and adults. The methods described here are less expensive than commercial kits and are more easily scalable to larger tissue amounts. Preparation of DNA can be completed in ∼7 h, RNA preparation in ∼1.5 h and protein preparation in ∼1 h.

Characterizing the spatiotemporal expression of RNAs and proteins in the starlet sea anemone, Nematostella vectensis.

In an effort to reconstruct the early evolution of animal genes and proteins, there is an increasing focus on basal animal lineages such as sponges, cnidarians, ctenophores and placozoans. Among the basal animals, the starlet sea anemone Nematostella vectensis (phylum Cnidaria) has emerged as a leading laboratory model organism partly because it is well suited to experimental techniques for monitoring and manipulating gene expression. Here we describe protocols adapted for use in Nematostella to characterize the expression of RNAs by in situ hybridization using either chromogenic or fluorescence immunohistochemistry (∼1 week), as well as to characterize protein expression by whole-mount immunofluorescence (∼3 d). We also provide a protocol for labeling cnidocytes (∼3 h), the phylum-specific sensory-effector cell type that performs a variety of functions in cnidarians, including the delivery of their venomous sting.

Collecting, rearing, spawning and inducing regeneration of the starlet sea anemone, Nematostella vectensis.

Over the past 20 years, the starlet sea anemone, Nematostella vectensis, a small estuarine animal, has emerged as a powerful model system for field and laboratory studies of development, evolution, genomics, molecular biology and toxicology. Here we describe how to collect Nematostella, culture it through its entire sexual life cycle and induce regeneration for the production of clonal stocks. In less than 1 h at a suitable field site, a researcher on foot can collect hundreds of individual anemones. In a few months, it is possible to establish a laboratory colony that will be reliable in generating hundreds or thousands of fertilized eggs on a roughly weekly schedule. By inducing regeneration roughly every 2 weeks, in less than 6 months, one can establish a clonal stock consisting of hundreds of genetically identical anemones. These results can be achieved very inexpensively and without specialized equipment.

A report from the second Nematostella vectensis research conference.

This report summarizes information discussed at the second Nematostella vectensis research conference, which took place on August 27, 2012 in Boston, MA, USA. The startlet sea anemone Nematostella is emerging as one of leading model organisms among cnidarians, in part because of the extensive genome and transcriptome resources that are becoming available for Nematostella, which were the focus of several presentations. In addition, research was presented on the use of Nematostella in developmental, regeneration, signal transduction, host-symbiont, and gene-environment interaction studies.

NF-κB is required for cnidocyte development in the sea anemone Nematostella vectensis.

The sea anemone Nematostella vectensis (Nv) is a leading model organism for the phylum Cnidaria, which includes anemones, corals, jellyfishes and hydras. A defining trait across this phylum is the cnidocyte, an ectodermal cell type with a variety of functions including defense, prey capture and environmental sensing. Herein, we show that the Nv-NF-κB transcription factor and its inhibitor Nv-IκB are expressed in a subset of cnidocytes in the body column of juvenile and adult anemones. The size and distribution of the Nv-NF-κB-positive cnidocytes suggest that they are in a subtype known as basitrichous haplonema cnidocytes. Nv-NF-κB is primarily cytoplasmic in cnidocytes in juvenile and adult animals, but is nuclear when first detected in the 30-h post-fertilization embryo. Morpholino-mediated knockdown of Nv-NF-κB expression results in greatly reduced cnidocyte formation in the 5 day-old animal. Taken together, these results indicate that NF-κB plays a key role in the development of the phylum-specific cnidocyte cell type in Nematostella, likely by nuclear Nv-NF-κB-dependent activation of genes required for cnidocyte development.

Production of a reference transcriptome and transcriptomic database (PocilloporaBase) for the cauliflower coral, Pocillopora damicornis.

BACKGROUND: Motivated by the precarious state of the world's coral reefs, there is currently a keen interest in coral transcriptomics. By identifying changes in coral gene expression that are triggered by particular environmental stressors, we can begin to characterize coral stress responses at the molecular level, which should lead to the development of more powerful diagnostic tools for evaluating the health of corals in the field. Furthermore, the identification of genetic variants that are more or less resilient in the face of particular stressors will help us to develop more reliable prognoses for particular coral populations. Toward this end, we performed deep mRNA sequencing of the cauliflower coral, Pocillopora damicornis, a geographically widespread Indo-Pacific species that exhibits a great diversity of colony forms and is able to thrive in habitats subject to a wide range of human impacts. Importantly, P. damicornis is particularly amenable to laboratory culture. We collected specimens from three geographically isolated Hawaiian populations subjected to qualitatively different levels of human impact. We isolated RNA from colony fragments ("nubbins") exposed to four environmental stressors (heat, desiccation, peroxide, and hypo-saline conditions) or control conditions. The RNA was pooled and sequenced using the 454 platform. DESCRIPTION: Both the raw reads (n=1, 116, 551) and the assembled contigs (n=70, 786; mean length=836 nucleotides) were deposited in a new publicly available relational database called PocilloporaBase http://www.PocilloporaBase.org. Using BLASTX, 47.2% of the contigs were found to match a sequence in the NCBI database at an E-value threshold of ≤.001; 93.6% of those contigs with matches in the NCBI database appear to be of metazoan origin and 2.3% bacterial origin, while most of the remaining 4.1% match to other eukaryotes, including algae and amoebae. CONCLUSIONS: P. damicornis now joins the handful of coral species for which extensive transcriptomic data are publicly available. Through PocilloporaBase http://www.PocilloporaBase.org, one can obtain assembled contigs and raw reads and query the data according to a wide assortment of attributes including taxonomic origin, PFAM motif, KEGG pathway, and GO annotation.

Two polymorphic residues account for the differences in DNA binding and transcriptional activation by NF-κB proteins encoded by naturally occurring alleles in Nematostella vectensis.

The NF-κB family of transcription factors is activated in response to many environmental and biological stresses, and plays a key role in innate immunity across a broad evolutionary expanse of animals. A simple NF-κB pathway is present in the sea anemone Nematostella vectensis, an important model organism in the phylum Cnidaria. Nematostella has previously been shown to have two naturally occurring NF-κB alleles (Nv-NF-κB-C and Nv-NF-κB-S) that encode proteins with different DNA-binding and transactivation abilities. We show here that polymorphic residues 67 (Cys vs. Ser) and 269 (Ala vs. Glu) play complementary roles in determining the DNA-binding activity of the NF-κB proteins encoded by these two alleles and that residue 67 is primarily responsible for the difference in their transactivation ability. Phylogenetic analysis indicates that Nv-NF-κB-S is the derived allele, consistent with its restricted geographic distribution. These results define polymorphic residues that are important for the DNA-binding and transactivating activities of two naturally occurring variants of Nv-NF-κB. The implications for the appearance of the two Nv-NF-κB alleles in natural populations of sea anemones are discussed.

Evolutionary expression of glucose-dependent-insulinotropic polypeptide (GIP).

Glucose-dependent insulinotropic polypeptide (GIP) is a mammalian incretin hormone released into the circulation following nutrient ingestion. We examined the functional evolution of GIP and its relationship with insulin to delineate their respective roles in promoting nutrient efficiency. Expression patterns were examined in the sea lamprey (Petromyzon marinus), a basal vertebrate lacking a distinct pancreas, and in the zebrafish, Xenopus laevis, chicken, and mouse, organisms possessing extraintestinal pancreata. Although sea lamprey genomic analysis predicted a potential GIP-like gene, transcripts were not detected, and insulin expression was confined to the caudal pancreatic bud. GIP was detected in both the intestine and pancreas of the zebrafish and X. laevis. In contrast, GIP and insulin expression were limited to the intestine and pancreas, respectively, in chicken and mouse. Phylogenetic analysis of the glucagon-like ligands suggested proglucagon as the common ancestor, supporting the theory that GIP arose as a gene duplication of proglucagon. Insulin-secreting cells in the sea lamprey intestine may have obviated the need for an enteroinsular axis, and zebrafish may represent an evolutionary transition where GIP does not yet function as an incretin hormone. These observations are consistent with the hypothesis that GIP and insulin influence survival advantage by enhancing the efficiency of nutrient absorption and energy storage.

Evolution of function of a fused metazoan tRNA synthetase.

The origin and evolution of multidomain proteins are driven by diverse processes including fusion/fission, domain shuffling, and alternative splicing. The 20 aminoacyl-tRNA synthetases (AARS) constitute an ancient conserved family of multidomain proteins. The glutamyl-prolyl tRNA synthetase (EPRS) of bilaterian animals is unique among AARSs, containing two functional enzymes catalyzing ligation of glutamate and proline to their cognate transfer RNAs (tRNAs). The ERS and PRS catalytic domains in multiple bilaterian taxa are linked by variable number of helix-turn-helix domains referred to as WHEP-TRS domains. In addition to its canonical aminoacylation activities, human EPRS exhibits a noncanonical function as an inflammation-responsive regulator of translation. Recently, we have shown that the WHEP domains direct this auxiliary function of human EPRS by interacting with an mRNA stem-loop element (interferon-gamma-activated inhibitor of translation [GAIT] element). Here, we show that EPRS is present in the cnidarian Nematostella vectensis, which pushes the origin of the fused protein back to the cnidarian-bilaterian ancestor, 50-75 My before the origin of the Bilateria. Remarkably, the Nematostella EPRS mRNA is alternatively spliced to yield three isoforms with variable number and sequence of WHEP domains and with distinct RNA-binding activities. Whereas one isoform containing a single WHEP domain binds tRNA, a second binds both tRNA and GAIT element RNA. However, the third isoform contains two WHEP domains and like the human ortholog binds specifically to GAIT element RNA. These results suggest that alternative splicing of WHEP domains in the EPRS gene of the cnidarian-bilaterian ancestor gave rise to a novel molecular function of EPRS conserved during metazoan evolution.

Characterization of the core elements of the NF-κB signaling pathway of the sea anemone Nematostella vectensis.

The sea anemone Nematostella vectensis is the leading developmental and genomic model for the phylum Cnidaria, which includes anemones, hydras, jellyfish, and corals. In insects and vertebrates, the NF-κB pathway is required for cellular and organismal responses to various stresses, including pathogens and chemicals, as well as for several developmental processes. Herein, we have characterized proteins that comprise the core NF-κB pathway in Nematostella, including homologs of NF-κB, IκB, Bcl-3, and IκB kinase (IKK). We show that N. vectensis NF-κB (Nv-NF-κB) can bind to κB sites and activate transcription of reporter genes containing multimeric κB sites or the Nv-IκB promoter. Both Nv-IκB and Nv-Bcl-3 interact with Nv-NF-κB and block its ability to activate reporter gene expression. Nv-IKK is most similar to human IKKε/TBK kinases and, in vitro, can phosphorylate Ser47 of Nv-IκB. Nv-NF-κB is expressed in a subset of ectodermal cells in juvenile and adult Nematostella anemones. A bioinformatic analysis suggests that homologs of many mammalian NF-κB target genes are targets for Nv-NF-κB, including genes involved in apoptosis and responses to organic compounds and endogenous stimuli. These results indicate that NF-κB pathway proteins in Nematostella are similar to their vertebrate homologs, and these results also provide a framework for understanding the evolutionary origins of NF-κB signaling.

A conserved cluster of three PRD-class homeobox genes (homeobrain, rx and orthopedia) in the Cnidaria and Protostomia.

BACKGROUND: Homeobox genes are a superclass of transcription factors with diverse developmental regulatory functions, which are found in plants, fungi and animals. In animals, several Antennapedia (ANTP)-class homeobox genes reside in extremely ancient gene clusters (for example, the Hox, ParaHox, and NKL clusters) and the evolution of these clusters has been implicated in the morphological diversification of animal bodyplans. By contrast, similarly ancient gene clusters have not been reported among the other classes of homeobox genes (that is, the LIM, POU, PRD and SIX classes). RESULTS: Using a combination of in silico queries and phylogenetic analyses, we found that a cluster of three PRD-class homeobox genes (Homeobrain (hbn), Rax (rx) and Orthopedia (otp)) is present in cnidarians, insects and mollusks (a partial cluster comprising hbn and rx is present in the placozoan Trichoplax adhaerens). We failed to identify this 'HRO' cluster in deuterostomes; in fact, the Homeobrain gene appears to be missing from the chordate genomes we examined, although it is present in hemichordates and echinoderms. To illuminate the ancestral organization and function of this ancient cluster, we mapped the constituent genes against the assembled genome of a model cnidarian, the sea anemone Nematostella vectensis, and characterized their spatiotemporal expression using in situ hybridization. In N. vectensis, these genes reside in a span of 33 kb with the same gene order as previously reported in insects. Comparisons of genomic sequences and expressed sequence tags revealed the presence of alternative transcripts of Nv-otp and two highly unusual protein-coding polymorphisms in the terminal helix of the Nv-rx homeodomain. A population genetic survey revealed the Rx polymorphisms to be widespread in natural populations. During larval development, all three genes are expressed in the ectoderm, in non-overlapping territories along the oral-aboral axis, with distinct temporal expression. CONCLUSION: We report the first evidence for a PRD-class homeobox cluster that appears to have been conserved since the time of the cnidarian-bilaterian ancestor, and possibly even earlier, given the presence of a partial cluster in the placozoan Trichoplax. Very similar clusters comprising these three genes exist in Nematostella and diverse protostomes. Interestingly, in chordates, one member of the ancestral cluster (homeobrain) has apparently been lost, and there is no linkage between rx and orthopedia in any of the vertebrates. In Nematostella, the spatial expression of these three genes along the body column is not colinear with their physical order in the cluster but the temporal expression is, therefore, using the terminology that has been applied to the Hox cluster genes, the HRO cluster would appear to exhibit temporal but not spatial colinearity. It remains to be seen whether the mechanisms responsible for the evolutionary conservation of the HRO cluster are the same mechanisms responsible for cohesion of the Hox cluster and other ANTP-class homeobox clusters that have been widely conserved throughout animal evolution.