- Open Access
The innate immune repertoire in Cnidaria - ancestral complexity and stochastic gene loss
Genome Biologyvolume 8, Article number: R59 (2007)
Characterization of the innate immune repertoire of extant cnidarians is of both fundamental and applied interest - it not only provides insights into the basic immunological 'tool kit' of the common ancestor of all animals, but is also likely to be important in understanding the global decline of coral reefs that is presently occurring. Recently, whole genome sequences became available for two cnidarians, Hydra magnipapillata and Nematostella vectensis, and large expressed sequence tag (EST) datasets are available for these and for the coral Acropora millepora.
To better understand the basis of innate immunity in cnidarians, we scanned the available EST and genomic resources for some of the key components of the vertebrate innate immune repertoire, focusing on the Toll/Toll-like receptor (TLR) and complement pathways. A canonical Toll/TLR pathway is present in representatives of the basal cnidarian class Anthozoa, but neither a classic Toll/TLR receptor nor a conventional nuclear factor (NF)-κB could be identified in the anthozoan Hydra. Moreover, the detection of complement C3 and several membrane attack complex/perforin domain (MAC/PF) proteins suggests that a prototypic complement effector pathway may exist in anthozoans, but not in hydrozoans. Together with data for several other gene families, this implies that Hydra may have undergone substantial secondary gene loss during evolution. Such losses are not confined to Hydra, however, and at least one MAC/PF gene appears to have been lost from Nematostella.
Consideration of these patterns of gene distribution underscores the likely significance of gene loss during animal evolution whilst indicating ancient origins for many components of the vertebrate innate immune system.
The innate immune system is the first line of defense against pathogens, and in non-chordates is assumed to be the sole means by which any non-self cells are detected and either killed or contained . Innate immunity in vertebrates is essentially a two-tier system consisting on one hand of phagocyte activation by the interaction of specialized surface receptors with pathogens or pathogen-derived components, and on the other of the direct opsonization and lysis of pathogens via the complement cascade. Whilst the vertebrate innate immune system has been the subject of intense investigation and is relatively well understood, studies of invertebrate immunity, which have focused primarily on the arthropods Drosophila and various horseshoe crab species [2–4], have revealed some striking similarities. For example, in both Drosophila and vertebrates, the Toll/Toll-like receptor (TLR) mediates the activation of appropriate response genes to microbial challenge [5, 6].
Toll and the TLRs are transmembrane proteins with a characteristic domain structure consisting of an extracellular amino-terminal domain containing leucine-rich repeats (LRRs) responsible for pattern recognition and an intracellular Toll interleukin receptor (TIR) domain that mediates signal transmission. Although the Toll and TLR families of arthropods and mammals are thought to have independently diversified [7, 8], all Tolls and TLRs signal via a common pathway that is conserved between Drosophila and mammals. The ultimate step in this pathway is translocation of nuclear factor (NF)-κB or its fly counterpart (the Dif/Rel heterodimer) into the nucleus, where it stimulates transcription of appropriate response genes. The immune repertoire of the horseshoe crab Carcinoscorpius includes a complex complement pathway that has both opsonic and lytic effector functions . Horseshoe crab complement C3 is functionally homologous with mammalian C3, mediating phagocytosis of bacteria (by hemocytes) in a strikingly similar manner.
Whilst these specific studies imply that at least some innate immune mechanisms have been conserved, broader comparative studies highlight the extent of gene loss and divergence in various metazoan lineages. For example, although Carcinoscorpius clearly uses a vertebrate-like complement system, none of the central components of the cascade (C2, C3, C4, C5) are encoded by the genomes of the ecdysozoans Drosophila, Caenorhabditis or Anopheles. Moreover, the sole Toll/TLR in Caenorhabditis elegans and C. brigssae is not known to function in the context of immunity, nor does that reported in the horseshoe crab Tachypleus tridentatus . There are also important differences between the Toll/TLR systems of Drosophila and mammals. For example, some mammalian TLRs themselves act as pattern recognition receptors (PRRs) upon microbial challenge, whereas in fly this is not the case . Moreover, whereas most of the ten or so vertebrate TLRs function primarily in immunity, only one of the nine fly (and ten mosquito) Tolls functions in this context. The others play a role in development , most famously in controlling differentiation in the dorsal/ventral axis.
The significance of gene loss in animal evolution has recently been brought into focus by preliminary expressed sequence tag (EST) and genomic analyses of some 'basal' animals (Figure 1), particularly the anthozoan cnidarians Acropora millepora and Nematostella vectensis [12, 13] and the planarian Dugesia japonica . Paradoxically, the genomes of these morphologically simple animals contain many genes previously thought to have evolved much later in the context of vertebrate complexity, and most of the complexity of signaling pathways and transcription factors associated with higher animals is represented in the anthozoan datasets [13, 15–17]. In contrast to Drosophila and Caenorhabditis, which have undergone substantial gene loss, for at least some groups of genes Acropora and Nematostella appear to have preserved much of the genetic complexity of the common metazoan ancestor. For example, whereas fly and worm have each lost approximately half of the ancestral Wnt complement, all but one of the 12 known Wnt subfamilies is represented in Nematostella . The emerging cnidarian EST and genomic datasets are, therefore, potentially highly informative with respect to the ancestral immunological repertoire. In addition to this basic evolutionary significance, understanding the bases of cnidarian immunity is of major applied significance in light of the dramatic decline in coral health that has occurred on a global scale over the past 20 years. Increasing human activity in coastal zones throughout the world has led to declines in water quality with increased sediment, nutrient and heavy metal concentrations. These have all had detrimental effects on corals with an associated increase in the prevalence of disease, and perhaps led to some new diseases, although this is less certain, as coral diseases are very poorly understood. Most are named for their symptoms, for example, 'white band disease', 'black band disease', and 'rapid wasting syndrome' and causative agents are frequently unknown. In the face of this uncertainty it is important that coral defense mechanisms should be better understood.
Although cnidarians have no specialized immune cells, at least some display highly specific allorecognition characteristics. Allorecognition, xenorecognition, and killing mechanisms have been demonstrated in several hydrozoans [18–20] and anthozoans [21–23]. Allorecognition is thought to protect colonial cnidarians from fusion with genetically different individuals and to prevent germ line parasitism. The effector mechanisms range from contact avoidance involving chemical sensing, to barrier formation, or usage of nematocysts. For example, the sea anemone Anthopleura xanthogrammmica will 'tolerate' adjacent clonal individuals, but will attempt to 'reject' heterogenic clones with which it comes into contact [24, 25]. In the anthozoans, Stylophora pistillata and Montipora verrucosa branches within one colony will readily fuse while branches of genetically different individuals never undergo fusion [22, 26, 27]. Fusion of two conspecific individuals is occasionally referred to as 'natural transplantation'. Observations in sea anemones (Anthopeura elegantissima, Phymactis clematis) and gorgonians (Eunicella stricta) indicate that individual colonies possess unique sets of histocompatibility elements, which are recognized as nonself by all other conspecific colonies [23, 28]. In Hydrozoa, the same phenomenon was reported for Millepora dichotoma  and studied in great detail in the colonial marine hydroid Hydractinia echinata . Hydractinia in fact was among the first invertebrates shown to display a genetically based system of intolerance against allogenic tissue: for more than 50 years it has been known that allorecognition and the inability to fuse stolons of different colonies is under the control of one polymorphic gene [31, 32]. Recent efforts using defined genetic lines of the hydroid Hydractinia symbiolongicarpus have confirmed this and shown that the single chromosomal region contains at least two loci .
The availability of whole genome sequences of two basal cnidarians at respectably high levels of redundancy - the hydrozoan Hydra magnipapillata (>6-fold coverage) and the anthozoan Nematostella vectensis (>7.6-fold coverage ) - together with large-scale EST datasets for these and for the coral Acropora millepora potentially offer new perspectives on the origins of mammalian immune functions. Here we report the results of a screen of the available genomic and EST resources for the cnidarian counterparts of key components of the vertebrate innate immune repertoire.
The toll receptor and other proteins containing the TIR domain
Searching the Hydra predicted protein collection using the available hidden Markov models (HMMs) identified only four TIR domain-containing proteins, two of which are clearly related to MyD88, which functions downstream of TLRs (Table 1). Consistent with their assignment as MyD88 family members, both of these Hydra proteins also contain the characteristic DEATH domain. The two other Hydra TIR proteins are atypical transmembrane proteins in having relatively short extracellular domains that are devoid of the LRR domains that characterize Toll and the TLRs (Figure 2). cDNAs encoding these proteins have previously been isolated by the Bosch laboratory (unpublished data) and their functions are presently under investigation; these proteins are known as HyTRR-1 and HyTRR-2. Surprisingly, extensive searching of the Hydra genome and all available EST/cDNA resources failed to identify any proteins having the canonical Toll/TLR structure, characterized by possession of both LRR and TIR domains.
Whereas only four TIR proteins are present in Hydra, substantially more could be identified amongst the predicted proteins from Nematostella using HMM-based search methods. Five of them were sufficiently complete to be included in the analyses presented here. These include a single MyD88 homolog (NvMyD88) and a protein (NvTLR-1) clearly related to members of the Toll/TLR family (Figure 2). Whereas the mammalian TLRs, and some members of the fly Toll/TLR family, have only a carboxy-terminal cysteine-rich motif flanking the LRRs proximal to the membrane, Nematostella NvTLR-1 is predicted to contain both carboxy- and amino-terminal-flanking cysteine-rich motifs in the extracellular part of the protein (Figure 2). This suggests that fly and anemone Toll more closely reflect the ancestral domain structure than do the mammalian TLRs. Moreover, a phylogenetic analysis (Figure 3) groups the TIR in Nematostella NvTLR-1 with its fly and human counterparts, with strong bootstrap support. Surprisingly, three more of the predicted Nematostella TIR proteins also contain multiple immunoglobulin (Ig) domains (Figure 2), and thus reflect the domain structure of mammalian interleukin 1 receptors (IL-1Rs). NvIL-1R1 and NvIL-1R2 each contain three Ig domains, and NvIL-1R3 contains two predicted Ig domains (Figure 2) but may be incomplete. In the phylogenetic analysis based on TIR domains the Nematostella IL-1R-like proteins form a clade distinct from both the MyD88 and Toll/TLR types (Figure 3), although these cnidarian TIRs appear to be distinct from those in the vertebrate IL-1Rs (data not shown). Several other TIR proteins were detected amongst the sequences of Nematostella (Additional data file 1), but were not subjected to further analysis as the TIR domains were incomplete or the sequences were judged likely to be artifactual. Two complete TIRs were identified by searching the available coral datasets. The trace archive yielded one TIR from Acropora palmata (ApGenomic) and a second was encoded by an A. millepora EST (AmTIR-1). These two coral TIRs are most similar to those in the Nematostella IL-1R-like proteins (Figure 3), but no linked domains have yet been identified in these cases.
The Müller group recently reported the identification of MyD88 in a demosponge, Suberites domuncula . However, whilst phylogenetic analyses clearly grouped the TIR in this sponge sequence with those present in unambiguous MyD88 orthologs (Figure 3), domain searching indicates that the predicted sponge protein may not have a functional DEATH domain.
The Toll/TLR pathway is ancestral but some components are missing or highly divergent in Hydra
Most of the intracellular mediators of Toll/TLR signaling could be identified in Nematostella and Acropora, but some key components appear to have either been lost or diverged beyond recognition in Hydra (Table 1). The absence of a Toll/TLR protein sensu stricto from Hydra is discussed above, but in addition only a single highly derived Rel domain could be found in Hydra whereas unambiguous NF-κB homologs are present in both Nematostella and Acropora (Table 1). In addition to the pathway leading to nuclear localization of NF-κB, Toll/TLR signaling can activate the Jun N-terminal kinase (JNK) and p38 mitogen-activated protein kinase (MAPK) pathways, leading to transcription of a range of target genes via the AP1 (activating protein 1)/ATF (activating transcriptionfactor 1) factors. Toll/TLR signaling via JNK/MAPK requires the participation of the ECSIT (evolutionarily conserved signaling intermediate in Toll pathways) adaptor protein , which also provides a link between the Toll/TLR and TGF-b (transforming growth factor-beta)/BMP (bone morphogenic protein) pathways . The presence of ECSIT as well as the key components of the JNK/MAPK pathway in the cnidarian datasets (Table 1, Figure 4) indicates an early origin for this variant of Toll/TLR signaling. The discovery of a conserved predicted ECSIT coding sequence in the fresh water sponge Ephydatia fluviatilis (N Funayama, personal observation) additionally supports this view.
Wiens et al.  have suggested that a sponge-specific cell surface protein known as SLIP (sponge LPS-interacting protein) functions as a pattern recognition receptor and effectively substitutes for Toll/TLR in antimicrobial defence. The sponge MyD88-related protein nominally functions downstream of SLIP in the proposed pathway . In support of the idea that sponges lack Toll/TLRs, the authors cite a lack of TLRs in a screen of 15,000 ESTs. As sponge 'MyD88' and SLIP are co-immunoprecipitated by the reciprocal antibodies , they clearly can interact in vitro. Some TLR pathway components could also be identified in the available sponge data, but the equivocal status of the sponge MyD88 related protein means that it is unclear at this time whether sponges have a canonical Toll/TLR pathway.
Cnidarian complement C3 and related proteins
The complement component C3 has recently been reported in another anthozoan cnidarian, the octocoral Swiftia , and the corresponding gene has recently been cloned from Acropora (Hayward, unpublished data). The Acropora C3 (C3-Am) gene is first expressed strongly in the endoderm of the planula as it elongates following gastrulation (Figure 5a). The endodermal expression is not uniform, being most intense in a subset of dark staining cells that have not yet been characterized. As the planula elongates expression becomes somewhat weaker, with the strongest expression localized to the aboral endoderm (Figure 5b). Post-settlement (Figure 5c-e) expression is limited to the endoderm and is particularly strong in the endoderm of the polyp as it rises from the calcifying platform at its base (for example, Figure 5d).
C3 has a complex domain structure. Whilst anthozoan C3s resemble their deuterostome counterparts both in domain structure (Figure 5f) and sequence, not only could no corresponding gene be identified in Hydra, but also some of the domains characteristic of C3 (ANATO, C345C; Figure 5f) could not be detected in any Hydra protein. Although lacking a canonical C3, Hydra contains a gene encoding A2M related domains. Interestingly, in situ hybridization in Hydra using a probe covering these typical A2M-related domains (Figure 5f; A2M-comp/A2M-recep) showed expression restricted to the endodermal epithelium (Figure 5g), as was the case with Acropora C3.
MAC/PF domain containing proteins in Cnidaria
Searching for other components of the complement cascade, we identified proteins containing a membrane attack complex/perforin domain (MAC/PF) similar to that present in complement component C6 and related proteins. HMM searching identified just two MAC/PF domain-containing proteins in Hydra (Table 1), whereas four proteins were identified in Nematostella. Two MAC/PF proteins were also identified amongst the Acropora ESTs. Database searches and analyses of predicted domain structures revealed that most of the cnidarian MAC/PF sequences are likely to fall into three groups corresponding to the known proteins types MPEG, TX-60A and apextrin (Table 1, Figure 5h).
TBlastN-based searches of the Nematostella genome identified a gene matching strongly to the human macrophage expressed protein 1 (MPEG1; GenBank:XP_166227) and its abalone homolog abMPEG1 (GenBank:AAR82936) . A clearly related gene in S. domuncula has recently been implicated as an effector in a hypothetical sponge innate immune defence pathway . Recombinant Suberites MPEG has anti-bacterial activity against Gram-negative bacteria, and is up-regulated after lipopolysaccharide (LPS) treatment . The MPEG1 family clearly has an ancient evolutionary history (the sponge and human sequences have 28% identity and 46% similarity)  but only in Suberites has any functional characterization been done. Despite the presence of MPEG1 in the sponge and an anthozoan, no corresponding gene could be identified in Hydra.
The nematocyst venom of at least some anthozoans contains the protein TX-60A , and two of the Nematostella MAC/PF proteins and one of the Acropora ESTs clearly correspond to this protein type (Table 1). TX-60A has an epidermal growth factor (EGF) domain immediately carboxy-terminal of the MAC/PF domain. In Hydra, this domain structure can be found in Hy-MAC, one of the two Hydra MAC/PF proteins (Figure 5h, Table 1). However, it is unclear whether the Hydra and anthozoan sequences are orthologous, as overall sequence identity is low. In situ hybridization analysis shows that expression of Hy-MAC is restricted to gland cells that are interspersed throughout the endoderm of Hydra (Figure 5i). Since endodermal gland cells and nematocysts are terminally differentiated , this pattern of expression is not easy to reconcile with a common function for the venom TX-60A and Hy-MAC.
Apextrin, a gene lost from Nematostella
The third class of cnidarian MAC/PF proteins represented in the Hydra and Acropora ESTs (Figure 5h) contains no identifiable domains other than MAC/PF. These proteins have moderate overall similarity to the echinoderm apextrins [41, 42] and to the apicomplexan protein family to which the Plasmodium membrane attack ookinete protein (MAOP)  belongs. MAOP is responsible for rupture of epithelial cells in the insect host by the ookinete stage of the parasite. Surprisingly, apextrin seems to be a case of gene loss from Nematostella as, despite clearly related genes being present in Hydra and Acropora, extensive searching of both the predicted protein collection and the anemone genome using a variety of tools failed to identify an apextrin-related gene (Table 1).
To explore the significance of this case of apparent gene loss, the expression pattern of an apextrin gene was examined by in situ hybridization in Acropora (Figure 5k-o). Apextrin-Am expression first appears in scattered ectodermal cells at the oral (blastopore) end of the embryo as it begins to elongate following blastopore closure (Figure 5k). As the embryo continues to elongate expression increases in intensity and the zone of expression spreads toward the aboral end, initially still in scattered cells (Figure 5l), but as elongation continues apparently in all ectodermal cells (Figure 5m), as is clear in transverse section (Figure 5n). As the planula settles, expression becomes less obvious at the oral end, eventually becoming limited to a belt marking the transition between the tissue that was formerly aboral (bottom), and the end bearing the oral pore, which subsequently forms the mouth of the polyp (Figure 5o). These expression data suggest that the primary function of apextrin-Am is in the ectoderm leading up to metamorphosis; hence, secondary loss of the corresponding gene from Nematostella may be explicable in terms of the very different modes of development of these two animals (see Discussion).
In adult Hydra, whole mount in situ hybridization showed an expression of the apextrin-like gene in groups of ectodermal cells arising from the interstitial cell lineage (Figure 5j), which may reflect a possible functional shift in an organism where metamorphosis is absent.
In addition to the above, the Nematostella dataset yielded MAC/PF domain-containing proteins having high similarity to the neural cell adhesion molecule spondin-1 (Additional data file 1).
These preliminary analyses of the newly available genomic and EST datasets indicate that a surprising number of key components of the innate immune system, including the Toll/TLR pathway and some complement cascade components, were in place at the base of the Eumetazoa. Also represented in the datasets are proteins strongly matching both the tumor necrosis factor (TNF) and Nod-like receptors and with the same domain structures (data not shown).
The analyses presented here are consistent with the idea of a genetically complex common metazoan ancestor [12, 13]; clearly the Toll/TLR, MyD88 and IL-1R protein families were distinct prior to the divergence of the Cnidaria from the Bilateria, and a Toll/TLR pathway may predate even the Porifera/Eumetazoa split. The discovery of proteins with the same domain structure as the IL-1R in Nematostella indicates that this receptor type predates chordate origins and that its original ligands may not have been interleukins. However, as the TIR domains in the cnidarian IL-1R-like and mammalian IL-1R proteins are divergent, separate evolutionary origins cannot yet be ruled out despite the similarities at the level of domain structure.
It is also likely that a prototypic complement effector pathway involving C3 and multiple MAC/PF proteins was in place in Ureumetazoa. Comparisons of the TIR and MAC/PF complements of Hydra and Nematostella highlight the likely extent of gene loss and sequence divergence in the former - not only does Hydra appear to have lost the Toll-receptor, but NF-κB has also either been lost or has diverged beyond recognition. Moreover, only a few TIR proteins could be identified in Hydra compared to the rich representation of these in the anthozoans. In addition, Hydra appears to have lost a number of MAC/PF proteins and lacks an equivalent of the ancestral complement C3 protein, implying that although anthozoans most likely have a prototype complement effector pathway, this has undergone degeneration in Hydra. Hydra is an efficient producer of a number of potent antimicrobial peptides (Bosch, unpublished data), indicating that at least some of the functions of the complement effector system may have been subsumed by pathways involved in synthesis of these peptides.
The observation that genes of two of the three classes of cnidarian immune-repertoire genes examined here are expressed in the endoderm (Figure 2) was unexpected, but is consistent with the primary sites of expression of such genes in many animals. In Drosophila, the fat body (an endodermal derivative, and the functional homolog of the mammalian liver) is the predominant source of the antimicrobial peptides whose synthesis is under the control of the Toll and Imd pathways  and in Caenorhabditis the gut is the primary source of antimicrobials [45, 46]. Although mammalian skin contains immune sensory cells (Langerhans cells, dendritic cells and so on), these are also endodermal derivatives, hence there is a deep evolutionary link between the immune system and the endoderm.
Other specific cases of gene loss from Hydra have been documented  (Seneca, unpublished data), hence there are precedents for the apparent loss of immune repertoire components reported here. More surprising is the implication that gene losses may also have occurred in Nematostella, as evidenced by the apparent absence of the apextrin gene that is present in both Hydra and Acropora. This loss may be associated with differences between Nematostella and Acropora during the planula to polyp transition. Metamorphosis from the motile planula to the sessile polyp involves dramatic tissue remodeling  and gene expression changes in Acropora (Grasso, unpublished data), and the massive up-regulation of apextrin leading up to and during metamorphosis (Figure 5k-o) is consistent with a role in this process. It is worthy of note that apextrin is expressed in the tissue that is least remodeled at the time of metamorphosis. By contrast, it appears that there are no dramatic events comparable to Acropora metamorphosis occurring during this same period of Nematostella development; rather, Nematostella undergoes a simple and continuous transition from planula to polyp and appears to show some neotenous features, continuing to glide over the bottom aboral end first for some time after metamorphosis , and never forming a pedal disc . Thus, loss of this gene may be related to the different modes of development of these two animals. This comparison between two anthozoans with different modes of development is clearly reminiscent of the original identification of apextrin as a gene specifically expressed in the 'direct' developing sea-urchin Heliocidaris erythrogramma during metamorphosis, but which is not expressed during the development of H. tuberculata, a congeneric 'indirect' developer [41, 42]. Because indirect development is considered to be ancestral in sea urchins, these authors hypothesized that apextrin had been coopted to function in metamorphosis in H. erythrogramma. The specific ectodermal expression of homologous proteins during metamorphosis in a cnidarian and an echinoderm is either a remarkable example of convergence or reflects conservation of function at some level. This latter possibility carries with it the implication that the 'simple' continuous planula to polyp transition seen in Nematostella might represent a derived state, and that the more dramatic metamorphosis seen in Acropora might more closely reflect the ancestral condition.
An important general implication that flows from these data is that gene loss may occur stochastically. If the genes in a pathway function only in that pathway, then one might expect that entire pathways would disappear following loss of one key component. However, the Hydra data appear to contradict this; most of the intracellular intermediates of Toll/TLR signaling are present despite loss of all of the major receptor types (both Toll/TLR and the IL-1R types known from more basal cnidarians are apparently absent from Hydra). Moreover, simple comparisons between related animals (for example, Nematostella versus Acropora) are unlikely to be informative in terms of understanding the origins of genes  - loss of a gene from Nematostella and loss of a specific function in the indirect developing urchin seems more likely than the independent co-option of the same protein to related roles in Acropora and a direct developing urchin. Stochastic loss underlying the distribution patterns of genes across the Metazoa may account for some cases of assumed lateral gene transfer - chance retention of an ancestral gene in one or a few animal lineages might easily be confused with lateral transfer. Reconstructing the genome of the common animal ancestor will not, therefore, be a simple undertaking; it will require whole genome data for a wide range of 'lower' as well as 'higher' animals.
Materials and methods
Genomic and EST sequence data were downloaded from public databases at NCBI (dbEST, Trace archive) or originate from private investigators (D Miller, K Agata) and were stored on a local comparative genomics analysis platform . The raw datasets included 10,171,402 genomic reads and 163,221 ESTs for H. magnipapillata, 5,996,730 genomic reads and 146,976 ESTs for N. vectensis, 14,625 genomic reads and 10,232 ESTs for A. millepora, 11,025 genomic reads for A. palmata and 11,450 genomic reads for Porites lobata. A set of 36,820 ESTs from the fresh water sponge E. fluviatilis is maintained in the lab of K Agata. For the corals A. millepora, A. palmata and P. lobata only pilot genomic sequencing data are available to date. Genomic trace data were maintained in the original raw format, whereas ESTs were clustered and assembled using TGICL tools from TIGR . The ESTscan software  was used to infer unigene and predicted peptide sequences from the assembled ESTs.
Database search, sequence analysis and phylogenetic methods
For database searches a local Blast-platform  and the public Blast platform at NCBI  were used. Domain searches were performed using SMART  and a local install of HMMer . Profile HMMs for the investigated domains were obtained from PFAM  and Superfamily  databases. Seqtools  and BioEdit  were used for general sequence analysis. Protein sequence alignments were created using ClustalW  and the HMMalign script included in the HMMer package. Maximum likelihood phylogenetic analyses were undertaken using MolPhy version 2.3  using the Dayhoff substitution matrix and local rearrangement search mode.
Whole mount in situ hybridization on Hydra polyps was performed as previously described by Grens et al. . Acropora embryos were fixed and processed for in situ hybridization as described in de Jong et al.  following the detailed in situ protocol given in Hayward et al. . For all genes shown, sense controls gave no hybridization signals.
Additional data files
The following additional data are available with the online version of this paper. Additional data file lists accession numbers for additional sequences identified within the database searches that were not further characterized in the present study.
Beutler B: Innate immunity: an overview. Mol Immunol. 2004, 40: 845-859. 10.1016/j.molimm.2003.10.005.
Litman GW, Cannon JP, Dishaw LJ: Reconstructing immune phylogeny: new perspectives. Nat Rev Immunol. 2005, 5: 866-879. 10.1038/nri1712.
Royet J, Reichhart JM, Hoffman JA: Sensing and signaling during infection in Drosophila. Curr Opin Immunol. 2005, 17: 11-17. 10.1016/j.coi.2004.12.002.
Flajnik MF, Du Pasquier L: Evolution of innate and adaptive immunity: Can we draw a line?. Trends Immunol. 2004, 25: 640-644. 10.1016/j.it.2004.10.001.
Lemaitre B, Reichhart JM, Hoffman JA: Drosophila host defense: differential induction of antimicrobial peptide genes after infection by various classes of microorganisms. Proc Natl Acad Sci USA. 1997, 94: 14614-14619. 10.1073/pnas.94.26.14614.
Medzhitov R, Preston-Hurlburt P, Janeway CA: A human homologue of the Drosophila Toll protein signals activation of adaptive immunity. Nature. 1997, 388: 394-397. 10.1038/41131.
Rock FL, Hardiman G, Timans JC, Kastelein RA, Bazan JF: A family of human receptors structurally related to Drosophila Toll. Proc Natl Acad Sci USA. 1998, 95: 588-593. 10.1073/pnas.95.2.588.
Luo C, Zheng L: Independent evolution of Toll and related genes in insects and mammals. Immunogenetics. 2000, 51: 92-98. 10.1007/s002510050017.
Zhu Y, Thangamani S, Ho B, Ding JL: The ancient origin of the complement system. EMBO J. 2005, 24: 382-394. 10.1038/sj.emboj.7600533.
Inamori K, Ariki S, Kawabata S: A toll-like receptor in horseshoe crabs. Immunol Rev. 2004, 198: 106-115. 10.1111/j.0105-2896.2004.0131.x.
Michel T, Reichhart JM, Hoffmann JA, Royet J: Drosophila Toll is activated by Gram-positive bacteria through a circulating peptidoglycan recognition protein. Nature. 2001, 414: 756-759. 10.1038/414756a.
Kortschak RD, Samuel G, Saint R, Miller DJ: EST analysis of the cnidarian, Acropora millepora, reveals extensive gene loss and rapid sequence divergence in the model invertebrates. Curr Biol. 2003, 13: 2190-2195. 10.1016/j.cub.2003.11.030.
Technau U, Rudd S, Maxwell P, Gordon PM, Saina M, Grasso LC, Hayward DC, Sensen CW, Saint R, Holstein TW, et al: Maintenance of ancestral complexity and non-metazoan genes in two basal cnidarians. Trends Genet. 2005, 21: 633-639. 10.1016/j.tig.2005.09.007.
Mineta K, Nakazawa M, Cebria F, Ikeo K, Agata K, Gojobori T: Origin and evolutionary process of the CNS elucidated by comparative genomics analysis of planarian ESTs. Proc Natl Acad Sci USA. 2003, 100: 7666-7671. 10.1073/pnas.1332513100.
Kusserow A, Pang K, Sturm C, Hrouda M, Lentfer J, Schmidt HA, Technau U, von Haeseler A, Hobmayer B, Martindale MQ, Holstein TW: Unexpected complexity of the Wnt gene family in a sea anemone. Nature. 2005, 433: 156-160. 10.1038/nature03158.
Kamm K, Schierwater B: Ancient complexity of the non-Hox ANTP gene complement in the anthozoan Nematostella vectensis. Implications for the evolution of the ANTP superclass. J Exp Zoolog B Mol Dev Evol. 2006, 306: 589-596. 10.1002/jez.b.21123.
Ryan JF, Burton PM, Mazza ME, Kwong GK, Mullikin JC, Finnerty JR: The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: evidence from the starlet sea anemone Nematostella vectensis. Genome Biol. 2006, 7: R64-10.1186/gb-2006-7-7-r64.
Campbell RD, Bibb C: Transplantation in coelenterates. Transplant Proc. 1970, 2: 202-211.
Leddy SV, Green DR: Historecognition in the Cnidaria. Phylogenesis of Immune Functions. Edited by: Warr GW, Cohen N. 1991, Boca Raton: CRC Press, 103-116.
Bosch TCG, David CN: Immunocompetence in Hydra: Epithelial cells recognize self-nonself and react against it. J Exp Zool. 1986, 238: 225-234. 10.1002/jez.1402380212.
Hildemann WH, Bigger CH, Johnston IS: Histoincompatibility reactions and allogeneic polymorphism among invertebrates. Transplant Proc. 1979, 11: 1136-1142.
Hildemann WH, Jokiel PL, Bigger CH, Johnston IS: Allogeneic polymorphism and alloimmune memory in the coral, Montipora verrucosa. Transplantation. 1980, 30: 297-301.
Lubbock R: Clone-specific cellular recognition in a sea anemone. Proc Natl Acad Sci USA. 1980, 77: 6667-6669. 10.1073/pnas.77.11.6667.
Francis L: Clone specific segregation in sea anemone Anthopleura elegantissima. Biol Bull. 1973, 144: 64-72. 10.2307/1540147.
Sebens KP: Agonistic behaviour in the inter-tidal sea anemone Anthopleura xanthogrammica. Biol Bull. 1984, 166: 457-472. 10.2307/1541154.
Müller WE, Müller I, Zahn RK, Maidhof A: Intraspecific recognition system in scleractinian corals: morphological and cytochemical description of the autolysis mechanism. J Histochem Cytochem. 1984, 32: 285-288.
Chadwick-Furman NE, Rinkevich B: A complex allorecognition in a reef building coral: delayed responses, reversals and nontransitive hierarchies. Coral Reefs. 1994, 13: 57-63. 10.1007/BF00426436.
Meinardi E, Florin-Christensen M, Paratcha G, Azcurra JM, Florin-Christensen J: The molecular basis of the self/nonself selectivity of a coelenterate toxin. Biochem Biophys Res Commun. 1995, 216: 348-354. 10.1006/bbrc.1995.2630.
Frank U, Rinkevich B: Alloimmune memory is absent in the Red Sea hydrocoral Millepora dichotoma. J Exp Zool. 2001, 291: 25-29. 10.1002/jez.3.
Frank U, Leitz T, Muller WA: The hydroid Hydractinia: a versatile, informative cnidarian representative. Bioessays. 2001, 23: 963-971. 10.1002/bies.1137.
Hauenschild CV: Genetische und entwicklungsphysiologische Untersuchungen über Intersexualität und Gewebeverträglichkeit bei Hydractinia echinata. Flem Wilhem Roux' Archiv. 1954, 147: 1-41. 10.1007/BF00576821.
Hauenschild CV: Über die Vererbung einer Gewebeverträglichkeitseigenschaft bei dem Hydroidpolypen Hydractinia echinata. Z Naturforsch. 1956, 11b: 132-138.
Cadavid LF, Powell AE, Nicotra ML, Moreno M, Buss LW: An invertebrate histocompatibility complex. Genetics. 2004, 167: 357-365. 10.1534/genetics.167.1.357.
Wiens M, Korzhev M, Krasko A, Thakur NL, Perovic-Ottstadt S, Breter HJ, Ushijima H, Diehl-Seifert B, Muller IM, Muller WE: Innate immune defense of the sponge Suberites domuncula against bacteria involves a MyD88-dependent signaling pathway. Induction of a peforin-like molecule. J Biol Chem. 2005, 280: 27949-27959. 10.1074/jbc.M504049200.
Kopp E, Medzhitov R, Carothers J, Xiao C, Douglas I, Janeway CA, Ghosh S: ECSIT is an evolutionarily conserved intermediate in the Toll/IL-1 signal transduction pathway. Genes Dev. 1999, 13: 2059-2071.
Xiao C, Shim JH, Kluppel M, Zhang SS, Dong C, Flavell RA, Fu XY, Wrana JL, Hogan BL, Ghosh S: Ecsit is required for BMP signaling and mesoderm formation during mouse embryogenesis. Genes Dev. 2003, 17: 2933-2949. 10.1101/gad.1145603.
Dishaw LJ, Smith SL, Bigger CH: Characterisation of a C3-like cDNA in a coral: phylogenetic implications. Immunogenetics. 2005, 57: 535-548. 10.1007/s00251-005-0005-1.
Mah SA, Moy GW, Swanson WJ, Vacquier VD: A perforin-like protein from a marine mollusk. Biochem Biophys Res Comm. 2004, 316: 468-475. 10.1016/j.bbrc.2004.02.073.
Oshiro N, Kobayashi C, Iwanaga S, Nozaki M, Namikoshi M, Spring J, Nagai H: A new membrane-attack complex/perforin (MAC/PF) domain lethal toxin from the nematocyst venom of the Okinawan sea anemone Actineria villosa. Toxicon. 2004, 43: 225-228. 10.1016/j.toxicon.2003.11.017.
Bosch TCG: Symmetry breaking in stem cells of the basal metazoan Hydra. Progress in Molecular and Subcellular Biology: Asymmetric Cell Division. Edited by: Macieira-Coelho A. 2006, Heidelberg: Springer, 61-78.
Haag ES, Raff RA: Isolation and characterization of three mRNAs enriched in embryos of the direct-developing sea urchin Heliocidaris erythrogramma: evolution of larval ectoderm. Dev Genes Evol. 1998, 208: 188-204. 10.1007/s004270050173.
Haag ES, Sly BJ, Andrews ME, Raff RA: Apextrin, a novel extracellular protein associated with larval ectoderm evolution in Heliocidaris erythrogramma. Dev Biol. 1999, 211: 77-87. 10.1006/dbio.1999.9283.
Kadota K, Ishino T, Matsuyama T, Chinzei Y, Yuda M: Essential role of membrane-attack protein in malarial transmission to mosquito host. Proc Natl Acad Sci USA. 2004, 101: 16310-16315. 10.1073/pnas.0406187101.
Mylonakis E, Aballay A: Worms and flies as genetically tractable animal models to study host-pathogen interactions. Infect Immun. 2005, 73: 3833-3841. 10.1128/IAI.73.7.3833-3841.2005.
Mallo GV, Kurz CL, Couillault C, Pujol N, Granjeaud S, Kohara Y, Ewbank JJ: Inducible antibacterial defense system in C. elegans. Curr Biol. 2002, 12: 1209-1214. 10.1016/S0960-9822(02)00928-4.
Millet AC, Ewbank JJ: Immunity in Caenorhabditis elegans. Curr Opin Immunol. 2004, 16: 4-9. 10.1016/j.coi.2003.11.005.
Chourrout D, Delsuc F, Chourrout P, Edvardsen RB, Rentzsch F, Renfer E, Jensen MF, Zhu B, de Jong P, Steele RE, Technau U: Minimal protohox cluster inferred from bilaterian and cnidarian Hox complements. Nature. 2006, 442: 684-687. 10.1038/nature04863.
Ball EE, Hayward DC, Reece-Hoyes JS, Hislop NR, Samuel G, Saint R, Harrison PL, Miller DJ: Coral development: from classical embryology to molecular control. Int J Dev Biol. 2002, 46: 671-678.
Hand C, Uhlinger KR: The culture, sexual and asexual reproduction, and growth of the sea anemone Nematostella vectensis. Biol Bull. 1992, 182: 169-176. 10.2307/1542110.
Collins AG, Cartwright P, McFadden CS, Schierwater B: Phylogenetic context and basal metazoan model systems. Integr Comp Biol. 2005, 45: 585-594. 10.1093/icb/45.4.585.
Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.
Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999, 7: 138-148.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP, Bork P: SMART 4.0: towards genomic data integration. Nucleic Acids Res. 2004, 32: D142-D144. 10.1093/nar/gkh088.
Sonnhammer EL, Eddy SR, Birney E, Bateman A, Durbin R: Pfam - multiple sequence alignments and HMM-profiles of protein domains. Nucleic Acids Res. 1998, 26: 320-322. 10.1093/nar/26.1.320.
Gough J, Karplus K, Hughey R, Chothia C: Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J Mol Biol. 2001, 313: 903-919. 10.1006/jmbi.2001.5080.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Adachi J, Hasegawa M: MOLPHY version 2.3: program for molecular phylogenetics based on maximum likelihood. Comput Sci Monogr. 1996, 28: 1-150.
Grens A, Gee L, Fisher DA, Bode HR: CnNK-2, an NK-2 homeobox gene, has a role in patterning the basal end of the axis in hydra. Dev Biol. 1996, 180: 473-488. 10.1006/dbio.1996.0321.
de Jong DM, Hislop NR, Hayward DC, Reece-Hoyes JS, Pontynen PC, Ball EE, Miller DJ: Components of both major axial patterning systems of the Bilateria are differentially expressed along the primary axis of a 'radiate' animal, the anthozoan cnidarian Acropora millepora. Dev Biol. 2006, 298: 632-643. 10.1016/j.ydbio.2006.07.034.
Hayward DC, Catmull J, Reece-Hoyes JS, Berghammer H, Dodd H, Hann SJ, Miller DJ, Ball EE: Gene structure and larval expression of cnox2Am from the coral Acropora millepora. Dev Genes Evol. 2001, 211: 10-19. 10.1007/s004270000112.
This work was supported by Grants from the Deutsche Forschungsgemeinschaft (DFG/SFB617) to TCGB and from the Australian Research Council (ARC) directly to DJM and EEB (Grants A00105431, DP0209460 and DP0344483) and via both the Centre for the Molecular Genetics of Development and the Centre of Excellence for Coral Reef Studies. The authors thank H Bode, R Steele and D Rokshar for their efforts in heading the Hydra and Nematostella Genome projects.
David J Miller, Georg Hemmrich contributed equally to this work.