Skip to main content

Novel genes dramatically alter regulatory network topology in amphioxus



Regulation in protein networks often utilizes specialized domains that 'join' (or 'connect') the network through specific protein-protein interactions. The innate immune system, which provides a first and, in many species, the only line of defense against microbial and viral pathogens, is regulated in this way. Amphioxus (Branchiostoma floridae), whose genome was recently sequenced, occupies a unique position in the evolution of innate immunity, having diverged within the chordate lineage prior to the emergence of the adaptive immune system in vertebrates.


The repertoire of several families of innate immunity proteins is expanded in amphioxus compared to both vertebrates and protostome invertebrates. Part of this expansion consists of genes encoding proteins with unusual domain architectures, which often contain both upstream receptor and downstream activator domains, suggesting a potential role for direct connections (shortcuts) that bypass usual signal transduction pathways.


Domain rearrangements can potentially alter the topology of protein-protein interaction (and regulatory) networks. The extent of such arrangements in the innate immune network of amphioxus suggests that domain shuffling, which is an important mechanism in the evolution of multidomain proteins, has also shaped the development of immune systems.


Protein networks are often 'joined' (or 'connected') by specialized protein-protein interaction domains that specifically recognize their targets and thus connect upstream and downstream elements of the network. The group of proteins involved in apoptosis, members of which incorporate the death domain (DD), death effector domain (DED), and caspase recruitment domain (CARD) [1], and the group of proteins involved in innate immunity, members of which incorporate the Toll/interleukin-1 receptor (TIR) domains [2, 3], represent excellent examples of such networks. Genomes of extensively studied organisms, such as Caenorhabditis elegans, Drosophila melanogaster, and human, display strong conservation of many elements of these two networks. In genome evolution, domain recombination events, such as fusion and fission, can create proteins with novel domain combinations that may lead to new functions, including providing new connections inside an existing network or between different networks [4, 5]. Traditionally, it was generally accepted that 'simpler' organisms have less complex networks and that 'more advanced' organisms add new elements to the canonical 'cores' of these networks. However, analyses of recently sequenced genomes, including sea urchin, amphioxus, and sea anemone, challenge this notion [68]. For instance, we have shown that the evolution of the apoptotic regulatory network consists of a succession of lineage-specific expansions and losses, which, combined with the limited number of 'apoptotic' protein families, has resulted in apparent similarities between networks in different organisms that mask an underlying complex evolutionary history [9]. Here, we focus our analysis on the innate immune system and discuss the potential effects of domain rearrangements on network topology.

The innate immune system mediates the primary line of defense against bacterial and viral infection and has distinctive roles in inflammatory diseases as well as in cancer [1012]. In evolutionary terms, innate immunity is very ancient, and several of its mediators can be traced to the basal metazoans (that is, Porifera [13] and Cnidaria [14]). Defense systems that share similarity to animals' innate immunity have also been described in plants, although the exact relationships between these two systems are not clear [15, 16]. The evolutionary history of innate immunity and its relationship to adaptive immune systems is of profound significance to our understanding of immune competence, interrelationships of immune mediators, and immune regulatory networks [17, 18]. The recent sequencing of the amphioxus and sea urchin genomes, which occupy critical positions in the evolution of the deuterostomes (Figure 1), provides a basis for approaching this broad question.

Figure 1
figure 1

Evolutionary relationships of select metazoans. Taxa are arranged in descending order of phylogenetic emergence relative to vertebrates. The protostomes/deuterostomes split is indicated by a red circle. The blue shading is used to distinguish deuterostomes from all other animals. One branch of the deuterostomes includes the chordates (shown against a light blue background) and the other includes the echinoderms (shown against a deep blue background). Times of phylogenetic divergence are not to scale, and the tree branches are intended only to depict general relationships. The phylogenetic relationships between chordates described here are based on the current view that the cephalochordate is the most basal group in the extant chordate lineage [1921].

Sea urchin, an echinoderm, is a representative of one of the two main branches of the deuterostome phylogeny [6]. Amphioxus, a cephalochordate, coming from one of the most basal groups in the extant chordate lineage [1921], represents the other (Figure 1). A large expansion in several multigene families encoding pathogen recognition molecules relative to both vertebrates, such as mammals, and invertebrates, such as C. elegans and D. melanogaster, was reported in sea urchin [22, 23]. Using different bioinformatics resources and tools as well as directed analysis of specific gene transcripts, we studied the innate immune genes in the recently completed amphioxus genome. We found a similar expansion in the numbers of innate receptors; however, unlike sea urchin, much of this expansion in amphioxus consists of genes with novel domain combinations. It is rather unexpected that such radical changes can occur in a relatively conserved network. At this point, amphioxus seems to be unique in the scale of its novel domain rearrangements, although the phenomenon of domain shuffling is likely to be a common mechanism of genome evolution. The extent of such changes in amphioxus highlights the importance of this mechanism in the evolutionary development of the innate immune system.


Large multigene families encoding innate receptors

Innate immune responses depend on several families of pattern-recognition receptors that recognize pathogen-associated molecular patterns and cellular danger signals, which originate from invading pathogens or are released by dying or injured cells. Two families of pattern-recognition receptors, the transmembrane Toll-like receptors (TLRs) [2426] and the intracellular NOD-like receptors (NLRs) [2729], are of particular interest because of their role in a number of diseases. Major differences in the numbers of the above pattern-recognition receptors, as well as in other receptors, such as scavenger receptor cysteine-rich (SRCR) proteins [30], have been reported in sea urchin relative to both vertebrates and other invertebrates [22, 23]. A similar expansion in these families is seen in the amphioxus genome (Table 1; Additional data file 3). The several-fold increases in the number of genes in these families in both sea urchin and amphioxus over other known invertebrates and vertebrates suggest that there is considerably more specificity in innate recognition in the former two species. It appears as if expansion of innate receptors is a shared characteristic of representatives of both arms of deuterostome evolution (Figure 1). From the standpoint of mammalian immunity, the findings in amphioxus are most interesting as the phenomena along the chordate arm of evolution has been lost in higher vertebrates; relatively few members of these families of innate receptors are found in vertebrate genomes.

Table 1 Expansion of protein families with innate immunity domains in amphioxus

The domain content of innate receptors in amphioxus is unique

TLRs consist of multiple leucine-rich repeats (LRRs) at the amino terminus and a TIR domain at the carboxyl terminus that recruits TIR domain-containing adaptors for downstream signaling [2, 31] (Figure 2a); examples (in human) are myeloid differentiation factor 88 (MyD88), TIR domain-containing adaptor protein (TIRAP), TIR domain-containing adaptor inducing interferon-β (TRIF), TRIF-related adaptor molecule (TRAM), and sterile α and HEAT-Armadillo motifs containing protein (SARM). Approximately eight domain combinations containing the TIR domain occur in mammals, five in Drosophila, and three in C. elegans (Figure 3; Additional data file 4). TIR domain combinations seen in Drosophila and C. elegans are also found in human. In contrast, 20 (out of a total of 28) domain combinations containing a TIR domain in amphioxus are specific to this organism. The difference with sea urchin is of particular note, since only about six TIR domain combinations exist in sea urchin, although the number of proteins containing TIR domains in sea urchin is even larger than in amphioxus (Table 1).

Figure 2
figure 2

The diversification of the innate immune arsenal in amphioxus. (a) A simplified model of extracellular and intracellular innate immune signaling in human. TLR signaling involves recruitment of a number of TIR domain-containing adaptors, including myeloid differentiation factor 88 (MyD88), TIR domain-containing adaptor protein (TIRAP), TIR domain-containing adaptor inducing interferon-β (TRIF), TRIF-related adaptor molecule (TRAM), and sterile α and HEAT-Armadillo motifs containing protein (SARM), which in turn activates transcription factors such as nuclear factor-κB (NF-κB) and interferon regulatory factors (IRFs) that ultimately lead to tumor necrosis factor (TNF) and type I interferon (IFN) production. NLR signaling can also stimulate inflammatory responses via the NF-κB pathway. Also, NLRs can form the inflammasome with apoptosis-associated speck-like protein (ASC) and procaspase-1, leading to the generation of the active form of interleukin (IL)-1β and IL-18. (b) The diversity of the innate immune system in amphioxus. Novel domain architectures as well as significant expansion in receptor number are evident. Selected 'direct connection' gene models are shown against a pink background. The cellular localization of amphioxus TLR proteins is still unclear; some of them could be localized in endosome in a manner equivalent to that seen in mammals. Domains: BIR, baculovirus inhibitor of apoptosis repeat domain [1]; CARD, caspase recruitment domain [1]; CASPASE, caspase [1]; DD, death domain [1]; DED, death effector domain [1]; IPAF, ICE (IL-1β converting enzyme) protease activating factor; LRR, leucine-rich repeat [24]; NACHT, NAIP, CIITA, HET-E, and TP1 [28]; NALP, NACHT, LRR, and PYRIN-domain-containing protein; NB-ARC, nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4 [42]; PYRIN, amino-terminal domain of protein pyrin [1]; TIR, Toll/interleukin-1 receptor [3, 26]; TNFR, tumor necrosis factor receptor [59]; WD40, Trp-Asp 40 [60].

Figure 3
figure 3

Different domain combinations in innate immunity receptor families. Numbers of different domains that combine with an individual TIR or NACHT domain in each designated genome are displayed. 'Average of all domains' (purple bars) means the average of domain combinations over all domains found in a genome. A detailed list of partner domains that combine with TIR or NACHT in each genome is given in Additional data file 4. The absolute numbers differ slightly when different Ensembl protein datasets or thresholds are used, but the relative fluctuations between different genomes are the same.

NLRs contain a nucleotide binding NACHT (domain present in neuronal apoptosis inhibitory protein (NAIP), CIITA, HET-E, and TP1) domain and are members of a distinct subfamily of the AAA+ (ATPase associated with diverse cellular activities) family [32]. In vertebrates, NLRs possess one of several types of linker domains (CARD, PYRIN/PAAD [amino-terminal domain of protein pyrin/pyrin, AIM (absent-in-melanoma), ASC (apoptosis-associated speck-like protein), and DD-like], or BIR (baculovirus inhibitor of apoptosis repeat)) at the amino terminus and multiple LRRs at the carboxyl terminus that effect pathogen recognition [3, 28] (Figure 2a). Upon activation, NLRs are believed to assemble into complexes (inflammasomes) and recruit and activate additional proteins, such as caspase-1 and caspase-5 [33]. In amphioxus, approximately 21 different domain combinations involve NACHT domains, whereas approximately 5 are predicted in mammals (Figure 3; Additional data file 4). The NACHT domain is absent in Drosophila and C. elegans. Finally, it is noteworthy that in amphioxus SRCR-containing proteins, the SRCR domain - another domain related to the innate immune system [30] - is also combined with a greater diversity of other domains than in comparable proteins of sea urchins and other animals (Additional data file 3), similar to observations noted about TIR and NACHT domains.

Unique domain combinations imply unique topology of innate receptors

Activation of downstream host-defense mechanisms occurs via specialized signal transduction pathways that are mediated by a number of specific protein domains [3, 34]. Domain shuffling can create multidomain proteins with new domain architectures and functions, including proteins serving as novel connectors in regulatory pathways [5]. Organisms differ not only in the sizes of protein families, but also in their domain architectures - the combination of different domains in multidomain proteins. To study such differences, we have previously developed the Comparative Analysis of Protein Domain Organization (CADO) software package [35], which provides a tool that can visualize and analyze domain combinations of proteins in a given genome. CADO defines protein organization as a graph in which protein domains are represented as nodes, and domain combinations, defined as instances of two domains found in one protein, are represented as edges (lines). Using CADO, domain graphs of two (or more) genomes can be compared, identifying similarities and differences both in individual domain combinations and in general topology of the domain graph [35, 36].

CADO-based analysis was applied in order to determine if the expansion of the innate immunity receptor families also resulted in changes to the overall topology of the innate immune network in terms of unique domain combinations. Based on the comparison of amphioxus, human, and sea urchin genomes, the TIR domain combination repertoire of sea urchin is very close to that seen in human (Figure 4a), although the copy number of TIR-containing sequences between human and sea urchin differs approximately 10-fold (Table 1). Almost all the TIR domain combinations present in human and sea urchin can also be identified in amphioxus, which are shown by gray lines in Figure 4b,c; however, amphioxus has many more unique TIR domain combinations. Most of the domain combinations seen in amphioxus are specific to this organism (red lines in Figure 4b,c).

Figure 4
figure 4

Difference between protein domain networks involving the TIR domain in amphioxus, human, and sea urchin. (a) A comparison by CADO of the domain network anchored by the TIR domain in human and sea urchin. (b) CADO picture anchored by the TIR domain between human and amphioxus. (c) CADO picture anchored by the TIR domain between amphioxus and sea urchin. A line connecting two domains indicates a predicted single protein domain combination. Common domain combinations between the selected genomes are shown in gray; amphioxus-specific combinations are shown in red; human-specific combinations are shown in blue; and sea urchin-specific combinations are shown in green. Please note that to simplify the graphical representation, Pfam clans are adopted for some Pfam domains. The CADO picture may differ slightly when different thresholds are used, for instance, the Ig-TIR domain combination can be found in sea urchin when using SMART domain definitions.

Similar observations have been made for NLRs. In this case, most of the differences reside in the amino-terminal domain. Instead of a vertebrate-specific PYRIN/PAAD domain, amphioxus can have CARD, DD, or DED as connector domains (Figure 2b). The DD-NACHT and DED-NACHT direct domain combinations seen in NLRs have not been seen in vertebrates but are found in sea urchin [6, 23] and Nematostella vectensis [7, 37]. Because the amino-terminal prodomain in amphioxus caspases can be any of the DD, DED, or CARD types, these hybrid intracellular pathogen recognition receptors may directly trigger the apoptosis response (Figure 2b), rather than function through an ASC-like 'hub'.

Other types of hybrid genes, including those encoding tumor necrosis factor receptor (TNFR)-caspase, LRRs-caspase, TIR-NACHT, TIR-[NB-ARC]-WD40s (NB-ARC is nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4; WD40 is Trp-Asp 40), TIR-sterile alpha motif (SAM), TIR-Laminin and so on, which potentially could mediate immune-related functions, have also been identified in the amphioxus genome.

The unique predicted hybrid genes are expressed

Despite the presence of unusually complex patterns of repetitive DNA, the current assembly of the amphioxus genome is generally highly reliable [19]; notwithstanding this high level of confidence in the hybrid gene predictions, it is essential to note that cDNA transcripts of many of the predicted hybrid proteins have been recovered. The TNFR-caspase domain protein (Joint Genome Institute (JGI) model: Brafl1_82667) represents one of the shortcut pathways of particular interest (Figure 2b; Additional data file 6 part a). This predicted transmembrane protein contains an extracellular TNFR domain and an intracellular caspase domain and presumably provides a shortcut between inflammatory-type signals and cell death. cDNA analyses not only validate this domain architecture but also have identified other related gene sequences, including more than one type of both TNFR and caspase domains. These transcripts are the products of three genetic regions on scaffolds: _41, _114, and _457. Other examples include cDNAs encoding: the death-caspase domain combination predicted in model Brafl1_105741 (fgenesh2_pg.scaffold_505000014); the death-NACHT domain combination in model Brafl1_82459 (fgenesh2_pg.scaffold_111000114) and Brafl1_89453 (fgenesh2_pg.scaffold_187000018); the DED-NACHT combination in model Brafl1_98233 (fgenesh2_pg.scaffold_317000043); and the TIR-SAM combination in model Brafl1_131196 (estExt_fgenesh2_pg.C_5050026), which are described in Additional data file 5.

The recovery of transcripts corresponding to the 'direct connector' genes is, in itself, important as many of these genes most likely exhibit developmental stage-specific expression, may be expressed in relatively low abundance, and/or are transcribed in cells that are present in relatively low numbers or are undergoing apoptosis. Efforts to locate the expression of hybrid genes are currently underway.


The large-scale expansion of several families of innate receptors in amphioxus parallels that seen in sea urchin and is a shared feature of both sides of the deuterostome split. The phenomenon of lineage-specific gene expansion has also been reported for protein families in other genomes [38]. Further sequencing efforts are required to establish if the large numbers of novel domain architectures in innate immune-related genes are specific only to amphioxus, are specific only to deuterostomes, or represent a more general mechanism. We stress that the exact functions of these genes from amphioxus remain unknown and that further experimental work is needed; however, it is reasonable to hypothesize that the wide variety of domain combinations reported here likely expands the functions of the innate immune system in amphioxus. It is tempting to speculate that perhaps functionality of the amphioxus specific genes is provided by other regulatory mechanisms in vertebrates and that better understanding of the functions of novel amphioxus genes may help in discovering these mechanisms.

Many of the domain combinations in amphioxus are present in separate proteins in vertebrates that are interconnected by multistep signaling pathways (examples shown in Figure 2b and Additional data file 6). As such, the amphioxus proteins can be viewed as shortcuts between two endpoints. The presence of such shortcuts would change the topology of the network in a way that can be described as a difference between 'hub-and-spoke' versus 'direct connection' networks [39]. For instance, a TIR-NACHT architecture, present in amphioxus but absent in vertebrates, is a shortcut that directly connects the extra- and intracellular pathogenic pattern-recognition pathways (Figure 2b). In human, these two pathways are likely connected 'indirectly' by transforming growth factor-β activated kinase 1 (TAK1), receptor-interacting protein 2 (RIP2), and/or other molecules, although the detailed relationships of this functional integration are not resolved [3, 34, 40]. Proteins composed of LRRs or TNFR domains that directly connect to the caspase domain could provide direct links between pathogen recognition and apoptosis (Figure 2b; Additional data file 6). All these proteins contain the conserved QACXG (where X is R, Q, or G) pentapeptide active-site motif [41] in their caspase domains and, thus, likely have proteolytic function (Additional data file 7). Amphioxus proteins that combine a TIR domain with an NB-ARC domain [42] and WD40 repeats share features with Apaf-1 (apoptotic protease activating factor 1; a central regulator of apoptosis in animals, which consists of a CARD domain, an NB-ARC domain, and multiple WD40 repeats). The association of these structures with an amino-terminal TIR domain suggests a direct link between the innate immunity and apoptosis networks.

In general, the innate immunity and apoptosis networks, which interact through a complex system of signaling pathways in human and other vertebrates, are closely intertwined in amphioxus through multiple direct connection proteins. It is possible that the close relationship between these two major systems represents an important innovation at the base of the deuterostome lineage that has been preserved throughout the vertebrates, albeit implemented through different mechanisms. It has been shown that the artificial joining of domains in novel combinations [4345] create new signaling pathways. Specifically, the chimeric adaptor proteins, which contain a DED with a phosphotyrosine-binding (PTB) or Src homology 2 (SH2) domain, can redirect tyrosine kinase signaling from survival and cell growth to apoptosis [45]. In another example, it has been shown that caspase can be activated by the chemically inducible dimerization (CID) signal, resulting in apoptosis when its catalytic domain is artificially fused to CID-binding domains [43]. These directed studies lend considerable support for potential functions of the multiple shortcut proteins that have been identified in amphioxus. Furthermore, the results suggest that engineering of constructs corresponding to the amphioxus chimeric molecules represents a viable approach for gaining a better understanding of how these molecules function in innate immunity. The presence of direct connectors has important consequences for the flexibility of the network. In the hub-and-spoke model, the number of possible connections is exponential, even with the linear growth of the number of proteins. A very large number of different 'direct connections' would be required to provide equivalent flexibility.

Although not characterized at the transcription level, some of the 'hub' domains and connections that are present in human can also be found in the cnidarian N. vectensis [14, 46], such as the NACHT domain, the death-TIR connection, the Ig-TIR connection, and so on. Thus, the 'hub-and-spoke' model could be considered ancestral and was reduced in the arthropod and nematode lineage by eliminating some 'destinations' and/or even 'hubs' (for example, C. elegans has only one Toll-like receptor, TOL-l [47], and one SARM-like TIR domain containing adaptor, TIR-1 [48]; the NACHT domain is absent in both C. elegans and Drosophila (Table 1)). Taken together with the observations reported here, expansion appears to have occurred at the base of deuterostomes, and further evolution may well have proceeded independently in the echinoderm and cephalochordate branches. Although proteins with novel domain combinations also have been found in sea urchin [23, 49], the extent of such direct connections appears to be far greater in amphioxus. It is reasonable to assume that some direct connections could have been lost with the emergence of the vertebrate adaptive immune system or effectively replaced by additional 'hub' molecules, such as the ASC in the vertebrate lineage [33]. In light of these changes, the topology of the network would become closer to that of the common ancestor. The coexistence of both shortcut and conventional pathways in an extant species is exceptional and underscores the potential relevance of amphioxus for understanding the selective advantages of such arrangements.


Two aspects of genome architecture and complexity influence innate immunity in amphioxus. First, large-scale gene expansion, a characteristic shared with sea urchin, creates a greater level of potential specificity in several families of innate immune receptors than is found in species with adaptive immune systems and could result in refinement of immune function. Second, novel domain architectures and, in particular, direct connections (shortcuts) in regulatory pathways can introduce a more refined level of functional integration of networks than would likely be achieved by the simple duplication and subsequent divergence of genes encoding immune receptors. A model for expansion and the possibility of topology change of a network is presumed in the analyses of the amphioxus genome presented here. A corollary issue raised by these observations is whether specific features of the amphioxus genome, such as the extraordinary level of site variation and unusually complex patterns of repetitive DNA, factor in such changes. Irrespective of their origins, genes with novel architectures in amphioxus could potentially serve as a pathway-level 'Rosetta stone' for elucidating new regulatory connections in the innate systems of contemporary vertebrates, similar to approaches that are used to elucidate protein and regulatory complexes in prokaryotic genomes [50]. Assuming that such shortcuts impart selective advantage, there is reason to look for signaling alternatives that may emulate the predicted distinct function implicit in these unique hybrid structures.

Materials and methods


The v.1.0 genome assembly and related gene models of amphioxus (Branchiostoma floridae) were obtained from the JGI [51] as were the genome assembly 1.0 and related protein set of the sea anemone (N. vectensis). The genome assembly Spur_v2.0 and the GLEAN3 gene models for the sea urchin (Strongylocentrotus purpuratus) were obtained from the Baylor College of Medicine Human Genome Sequencing Center [52]. The other genome sequences and corresponding protein sets, including human, mouse, dog, chicken, Xenopus, zebrafish, fugu, tetraodon, ciona, nematode (C. elegans), and fruit fly (D. melanogaster) were downloaded from Ensembl [53].

Database search and sequence analysis

Several rounds of PSITBLASTN [54] searches were performed against each genome using known human TIR or NACHT domain amino acid sequences as seeds. Hits were mapped to the corresponding genome protein set in order to obtain the full-length protein sequences (for sea urchin and sea anemone, some of the gene models were in addition predicted by GenScan [55]). All identified genes were checked using: first, reciprocal BLAST analysis; second, Pfam protein searches, performed either locally or at the Pfam website [56], which also address the issue of family specificity, such as distinguishing NACHT domain from NB-ARC domain based on different hidden Markov models; third, NCBI CD-Search [57] and local RPS-BLAST search; and fourth, multiple sequence-alignment and phylogeny analysis.

Domain combination analysis

Different combinations of innate immune domains identified in the aforementioned genomes were compared using the CADO [35] approach.

RT-PCR confirmation of select modular transcripts

JGI-predicted models were used to develop PCR strategies for identifying cDNA transcripts. The predicted transcripts were placed onto the current assembly (v.1.0) using local BLAST (v.2.2.11) to verify genomic organization (for example, exon/intron structure and gene copy number). Primers were designed (from visual alignments or with Primer3 [58]) to span domain combinations and specific exon/intron boundaries. Primer design accommodated variations due to genetic polymorphism and haplotype complexity, a significant confounding aspect of this type of analysis. Total RNA was isolated from 30 animals using RNA-Bee (Tel-Test, Inc., Friendswood, TX, USA), and cDNA synthesis was primed using either poly-A or random hexamer strategies (SuperScriptIII, Invitrogen, Carlsbad, CA, USA). cDNAs were combined and served as templates for PCR amplification. Certain transcripts could be detected only after two rounds of nested PCR. Transcribed sequences with the expected length were sequenced to confirm the predicted gene models. The verified amphioxus gene models in this study have been deposited in the GenBank database under accession numbers [GenBank:EU049583] to [GenBank:EU049596] and [GenBank:EU279424] to [GenBank:EU279425] (Additional data file 5).

Additional data files

The following additional data files are available. Additional data file 1 is a table listing the TIR domain containing sequences in different genomes. Additional data file 2 is a table listing the NACHT domain containing sequences in different genomes. Additional data file 3 is a table listing the SRCR domain combinations in different genomes. Additional data file 4 is a table listing partner domains that combine with individual TIR or NACHT domains in different genomes. Additional data file 5 is a table listing the selected JGI-predicted amphioxus gene models that have been verified by RT-PCR. Additional data file 6 is a figure showing examples of novel domain combinations in amphioxus that represent the shortcuts between two or more proteins present in human. Additional data file 7 is a figure showing alignment of sequences in the vicinity of the catalytic center of the caspase domain from human caspases and amphioxus proteins with TNFR-caspase or LRRs-caspase architectures.



apoptosis-associated speck-like protein


Comparative Analysis of Protein Domain Organization


caspase recruitment domain


chemically inducible dimerization


death domain


death effector domain


Joint Genome Institute


leucine-rich repeat


domain present in NAIP, CIITA, HET-E, and TP1


neuronal apoptosis inhibitory protein


nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4


NOD-like receptor


pyrin, AIM (absent-in-melanoma), ASC, and DD-like


amino-terminal domain of protein pyrin


sterile alpha motif


sterile α and HEAT-Armadillo motifs containing protein


scavenger receptor cysteine-rich


Toll/interleukin-1 receptor


Toll-like receptor


tumor necrosis factor receptor


Trp-Asp 40.


  1. Reed JC, Doctor KS, Godzik A: The domains of apoptosis: a genomics perspective. Sci STKE. 2004, 2004: re9-10.1126/stke.2392004re9.

    PubMed  Google Scholar 

  2. O'Neill LA, Bowie AG: The family of five: TIR-domain-containing adaptors in Toll-like receptor signalling. Nat Rev Immunol. 2007, 7: 353-364. 10.1038/nri2079.

    Article  PubMed  Google Scholar 

  3. Werts C, Girardin SE, Philpott DJ: TIR, CARD and PYRIN: three domains for an antimicrobial triad. Cell Death Differ. 2006, 13: 798-815. 10.1038/sj.cdd.4401890.

    Article  PubMed  CAS  Google Scholar 

  4. Fong JH, Geer LY, Panchenko AR, Bryant SH: Modeling the evolution of protein domain architectures using maximum parsimony. J Mol Biol. 2007, 366: 307-315. 10.1016/j.jmb.2006.11.017.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Patthy L: Modular assembly of genes and the evolution of new functions. Genetica. 2003, 118: 217-231. 10.1023/A:1024182432483.

    Article  PubMed  CAS  Google Scholar 

  6. Sea Urchin Genome Sequencing Consortium, Sodergren E, Weinstock GM, Davidson EH, Cameron RA, Gibbs RA, Angerer RC, Angerer LM, Arnone MI, Burgess DR, Burke RD, Coffman JA, Dean M, Elphick MR, Ettensohn CA, Foltz KR, Hamdoun A, Hynes RO, Klein WH, Marzluff W, McClay DR, Morris RL, Mushegian A, Rast JP, Smith LC, Thorndyke MC, Vacquier VD, Wessel GM, Wray G, Zhang L, Elsik CG, et al: The genome of the sea urchin Strongylocentrotus purpuratus. Science. 2006, 314: 941-952. 10.1126/science.1133609.

    Article  PubMed Central  Google Scholar 

  7. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, Jurka J, Genikhovich G, Grigoriev IV, Lucas SM, Steele RE, Finnerty JR, Technau U, Martindale MQ, Rokhsar DS: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317: 86-94. 10.1126/science.1139158.

    Article  PubMed  CAS  Google Scholar 

  8. Holland LZ, Albalat R, Azumi K, Benito-Gutierrez E, Blow MJ, Bronner-Fraser M, Brunet F, Butts T, Candiani S, Dishaw LJ, Ferrier DE, Garcia-Fernandez J, Gibson-Brown JJ, Gissi C, Godzik A, Hallbook F, Hirose D, Hosomichi K, Ikuta T, Inoko H, Kasahara M, Kasamatsu J, Kawashima T, Kimura A, Kobayashi M, Kozmik Z, Kubokawa K, Laudet V, Litman GW, McHardy AC, et al: The amphioxus genome illuminates vertebrate origins and cephalochordate biology. Genome Res. 2008, 18: 1100-1111. 10.1101/gr.073676.107.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Zmasek CM, Zhang Q, Ye Y, Godzik A: Surprising complexity of the ancestral apoptosis network. Genome Biol. 2007, 8: R226-10.1186/gb-2007-8-10-r226.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Medzhitov R: Recognition of microorganisms and activation of the immune response. Nature. 2007, 449: 819-826. 10.1038/nature06246.

    Article  PubMed  CAS  Google Scholar 

  11. Beutler B: Innate immunity: an overview. Mol Immunol. 2004, 40: 845-859. 10.1016/j.molimm.2003.10.005.

    Article  PubMed  CAS  Google Scholar 

  12. Lawton JA, Ghosh P: Novel therapeutic strategies based on toll-like receptor signaling. Curr Opin Chem Biol. 2003, 7: 446-451. 10.1016/S1367-5931(03)00077-2.

    Article  PubMed  CAS  Google Scholar 

  13. Müller WEG, Müller IM: Origin of the metazoan immune system: identification of the molecules and their functions in sponges. Integr Comp Biol. 2003, 43: 281-292. 10.1093/icb/43.2.281.

    Article  PubMed  Google Scholar 

  14. Miller DJ, Hemmrich G, Ball EE, Hayward DC, Khalturin K, Funayama N, Agata K, Bosch TC: The innate immune repertoire in cnidaria - ancestral complexity and stochastic gene loss. Genome Biol. 2007, 8: R59-10.1186/gb-2007-8-4-r59.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Nürnberger T, Brunner F, Kemmerling B, Piater L: Innate immunity in plants and animals: striking similarities and obvious differences. Immunol Rev. 2004, 198: 249-266. 10.1111/j.0105-2896.2004.0119.x.

    Article  PubMed  Google Scholar 

  16. Ausubel FM: Are innate immune signaling pathways in plants and animals conserved?. Nat Immunol. 2005, 6: 973-979. 10.1038/ni1253.

    Article  PubMed  CAS  Google Scholar 

  17. Hoffmann JA, Kafatos FC, Janeway CA, Ezekowitz RA: Phylogenetic perspectives in innate immunity. Science. 1999, 284: 1313-1318. 10.1126/science.284.5418.1313.

    Article  PubMed  CAS  Google Scholar 

  18. Litman GW, Cannon JP, Dishaw LJ: Reconstructing immune phylogeny: new perspectives. Nat Rev Immunol. 2005, 5: 866-879. 10.1038/nri1712.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutierrez EL, Dubchak I, Garcia-Fernandez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov VV, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin-I T, et al: The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008, 453: 1064-1071. 10.1038/nature06967.

    Article  PubMed  CAS  Google Scholar 

  20. Bourlat SJ, Juliusdottir T, Lowe CJ, Freeman R, Aronowicz J, Kirschner M, Lander ES, Thorndyke M, Nakano H, Kohn AB, Heyland A, Moroz LL, Copley RR, Telford MJ: Deuterostome phylogeny reveals monophyletic chordates and the new phylum Xenoturbellida. Nature. 2006, 444: 85-88. 10.1038/nature05241.

    Article  PubMed  CAS  Google Scholar 

  21. Blair JE, Hedges SB: Molecular phylogeny and divergence times of deuterostome animals. Mol Biol Evol. 2005, 22: 2275-2284. 10.1093/molbev/msi225.

    Article  PubMed  CAS  Google Scholar 

  22. Rast JP, Smith LC, Loza-Coll M, Hibino T, Litman GW: Genomic insights into the immune system of the sea urchin. Science. 2006, 314: 952-956. 10.1126/science.1134301.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Hibino T, Loza-Coll M, Messier C, Majeske AJ, Cohen AH, Terwilliger DP, Buckley KM, Brockton V, Nair SV, Berney K, Fugmann SD, Anderson MK, Pancer Z, Cameron RA, Smith LC, Rast JP: The immune gene repertoire encoded in the purple sea urchin genome. Dev Biol. 2006, 300: 349-365. 10.1016/j.ydbio.2006.08.065.

    Article  PubMed  CAS  Google Scholar 

  24. Medzhitov R: Toll-like receptors and innate immunity. Nat Rev Immunol. 2001, 1: 135-145. 10.1038/35100529.

    Article  PubMed  CAS  Google Scholar 

  25. Beutler B, Jiang Z, Georgel P, Crozat K, Croker B, Rutschmann S, Du X, Hoebe K: Genetic analysis of host resistance: Toll-like receptor signaling and immunity at large. Annu Rev Immunol. 2006, 24: 353-389. 10.1146/annurev.immunol.24.021605.090552.

    Article  PubMed  CAS  Google Scholar 

  26. West AP, Koblansky AA, Ghosh S: Recognition and signaling by toll-like receptors. Annu Rev Cell Dev Biol. 2006, 22: 409-437. 10.1146/annurev.cellbio.21.122303.115827.

    Article  PubMed  CAS  Google Scholar 

  27. Fritz JH, Ferrero RL, Philpott DJ, Girardin SE: Nod-like proteins in immunity, inflammation and disease. Nat Immunol. 2006, 7: 1250-1257. 10.1038/ni1412.

    Article  PubMed  CAS  Google Scholar 

  28. Martinon F, Tschopp J: NLRs join TLRs as innate sensors of pathogens. Trends Immunol. 2005, 26: 447-454. 10.1016/

    Article  PubMed  CAS  Google Scholar 

  29. Ting JP, Kastner DL, Hoffman HM: CATERPILLERs, pyrin and hereditary immunological disorders. Nat Rev Immunol. 2006, 6: 183-195. 10.1038/nri1788.

    Article  PubMed  CAS  Google Scholar 

  30. Sarrias MR, Gronlund J, Padilla O, Madsen J, Holmskov U, Lozano F: The Scavenger receptor cysteine-rich (SRCR) domain: an ancient and highly conserved protein module of the innate immune system. Crit Rev Immunol. 2004, 24: 1-37. 10.1615/CritRevImmunol.v24.i1.10.

    Article  PubMed  CAS  Google Scholar 

  31. McGettrick AF, O'Neill LA: The expanding family of MyD88-like adaptors in Toll-like receptor signal transduction. Mol Immunol. 2004, 41: 577-582. 10.1016/j.molimm.2004.04.006.

    Article  PubMed  CAS  Google Scholar 

  32. Martinon F, Tschopp J: Inflammatory caspases: linking an intracellular innate immune system to autoinflammatory diseases. Cell. 2004, 117: 561-574. 10.1016/j.cell.2004.05.004.

    Article  PubMed  CAS  Google Scholar 

  33. Martinon F, Tschopp J: Inflammatory caspases and inflammasomes: master switches of inflammation. Cell Death Differ. 2007, 14: 10-22. 10.1038/sj.cdd.4402038.

    Article  PubMed  CAS  Google Scholar 

  34. Delbridge LM, O'Riordan MX: Innate recognition of intracellular bacteria. Curr Opin Immunol. 2007, 19: 10-16. 10.1016/j.coi.2006.11.005.

    Article  PubMed  CAS  Google Scholar 

  35. Ye Y, Godzik A: Comparative analysis of protein domain organization. Genome Res. 2004, 14: 343-353. 10.1101/gr.1610504.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Wuchty S: Scale-free behavior in protein domain networks. Mol Biol Evol. 2001, 18: 1694-1702.

    Article  PubMed  CAS  Google Scholar 

  37. Darling JA, Reitzel AR, Burton PM, Mazza ME, Ryan JF, Sullivan JC, Finnerty JR: Rising starlet: the starlet sea anemone, Nematostella vectensis. Bioessays. 2005, 27: 211-221. 10.1002/bies.20181.

    Article  PubMed  CAS  Google Scholar 

  38. Lespinet O, Wolf YI, Koonin EV, Aravind L: The role of lineage-specific gene family expansion in the evolution of eukaryotes. Genome Res. 2002, 12: 1048-1059. 10.1101/gr.174302.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Poole RW, Butler V: Airline deregulation: the unfinished revolution. Regulation. 1999, 22: 44-51.

    Google Scholar 

  40. Kufer TA, Sansonetti PJ: Sensing of bacteria: NOD a lonely job. Curr Opin Microbiol. 2007, 10: 62-69. 10.1016/j.mib.2006.11.003.

    Article  PubMed  CAS  Google Scholar 

  41. Cohen GM: Caspases: the executioners of apoptosis. Biochem J. 1997, 326: 1-16.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. van der Biezen EA, Jones JD: The NB-ARC domain: a novel signalling motif shared by plant resistance gene products and regulators of cell death in animals. Curr Biol. 1998, 8: R226-R227. 10.1016/S0960-9822(98)70145-9.

    Article  PubMed  CAS  Google Scholar 

  43. Fan L, Freeman KW, Khan T, Pham E, Spencer DM: Improved artificial death switches based on caspases and FADD. Hum Gene Ther. 1999, 10: 2273-2285. 10.1089/10430349950016924.

    Article  PubMed  CAS  Google Scholar 

  44. Lim YM, Wong S, Lau G, Witte ON, Colicelli J: BCR/ABL inhibition by an escort/phosphatase fusion protein. Proc Natl Acad Sci USA. 2000, 97: 12233-12238. 10.1073/pnas.210253497.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. Howard PL, Chia MC, Del Rizzo S, Liu FF, Pawson T: Redirecting tyrosine kinase signaling to an apoptotic caspase pathway through chimeric adaptor proteins. Proc Natl Acad Sci USA. 2003, 100: 11267-11272. 10.1073/pnas.1934711100.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  46. 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.

    Article  PubMed  CAS  Google Scholar 

  47. Pujol N, Link EM, Liu LX, Kurz CL, Alloing G, Tan MW, Ray KP, Solari R, Johnson CD, Ewbank JJ: A reverse genetic analysis of components of the Toll signaling pathway in Caenorhabditis elegans. Curr Biol. 2001, 11: 809-821. 10.1016/S0960-9822(01)00241-X.

    Article  PubMed  CAS  Google Scholar 

  48. Couillault C, Pujol N, Reboul J, Sabatier L, Guichou JF, Kohara Y, Ewbank JJ: TLR-independent control of innate immunity in Caenorhabditis elegans by the TIR domain adaptor protein TIR-1, an ortholog of human SARM. Nat Immunol. 2004, 5: 488-494. 10.1038/ni1060.

    Article  PubMed  CAS  Google Scholar 

  49. Robertson AJ, Croce J, Carbonneau S, Voronina E, Miranda E, McClay DR, Coffman JA: The genomic underpinnings of apoptosis in Strongylocentrotus purpuratus. Dev Biol. 2006, 300: 321-334. 10.1016/j.ydbio.2006.08.053.

    Article  PubMed  CAS  Google Scholar 

  50. Marcotte EM, Pellegrini M, Ng HL, Rice DW, Yeates TO, Eisenberg D: Detecting protein function and protein-protein interactions from genome sequences. Science. 1999, 285: 751-753. 10.1126/science.285.5428.751.

    Article  PubMed  CAS  Google Scholar 

  51. US Department of Energy Joint Genome Institute. []

  52. Sea Urchin Genome Project. []

  53. Ensembl. []

  54. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Burge CB, Karlin S: Finding the genes in genomic DNA. Curr Opin Struct Biol. 1998, 8: 346-354. 10.1016/S0959-440X(98)80069-9.

    Article  PubMed  CAS  Google Scholar 

  56. Pfam. []

  57. NCBI CD-Search. []

  58. Primer3. []

  59. Locksley RM, Killeen N, Lenardo MJ: The TNF and TNF receptor superfamilies: integrating mammalian biology. Cell. 2001, 104: 487-501. 10.1016/S0092-8674(01)00237-9.

    Article  PubMed  CAS  Google Scholar 

  60. Smith TF, Gaitatzes C, Saxena K, Neer EJ: The WD repeat: a common architecture for diverse functions. Trends Biochem Sci. 1999, 24: 181-185. 10.1016/S0968-0004(99)01384-5.

    Article  PubMed  CAS  Google Scholar 

  61. Verification of JGI-predicted amphioxus gene models by cDNA sequencing. []

Download references


We thank J Rast for discussions and comments and B Pryor for editorial assistance. This work was supported by grants from the National Institutes of Health (AI056324 to QZ, 23338 to GWL, and GM076221 to CMZ and AG). B. floridae and N. vectensis genome data, including gene models and annotations, were produced by the US Department of Energy Joint Genome Institute and downloaded from their Web site. The authors acknowledge the JGI for their efforts in sequencing, assembling, and annotating the amphioxus genome. S. purpuratus genome data were produced by the Sea Urchin Genome Project at the Baylor College of Medicine.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Gary W Litman or Adam Godzik.

Additional information

Authors' contributions

QZ performed the sequence and domain analyses and prepared the figures. CMZ performed phylogenetic analyses. LJD developed approaches for identifying hybrid transcripts. MGM cloned and sequenced hybrid transcripts. YY contributed to the domain analyses of the predicted proteins. GWL interpreted immunology concepts. AG formulated the problem and planned the work. All authors contributed to the interpretation of the results and to writing of the paper.

Qing Zhang, Christian M Zmasek contributed equally to this work.

Electronic supplementary material

Additional data file 1: TIR domain containing sequences in different genomes. (PDF 67 KB)

Additional data file 2: NACHT domain containing sequences in different genomes. (PDF 70 KB)

Additional data file 3: SRCR domain combinations in different genomes. (PDF 17 KB)


Additional data file 4: Partner domains that combine with individual TIR or NACHT domains in different genomes. (PDF 17 KB)

Additional data file 5: Selected JGI-predicted amphioxus gene models that have been verified by RT-PCR. (PDF 14 KB)


Additional data file 6: Examples of novel domain combinations in amphioxus that represent the shortcuts between two or more proteins present in human. (DOC 196 KB)


Additional data file 7: Alignment of sequences in the vicinity of the catalytic center of the caspase domain from human caspases and amphioxus proteins with TNFR-caspase or LRRs-caspase architectures. (DOC 1 MB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Zhang, Q., Zmasek, C.M., Dishaw, L.J. et al. Novel genes dramatically alter regulatory network topology in amphioxus. Genome Biol 9, R123 (2008).

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: