Skip Navigation


Integrative and Comparative Biology Advance Access originally published online on August 2, 2007
Integrative and Comparative Biology 2007 47(6):854-864; doi:10.1093/icb/icm077
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary data
Right arrow All Versions of this Article:
47/6/854    most recent
icm077v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Matus, D. Q.
Right arrow Articles by Martindale, M. Q.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2007. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology. All rights reserved. For permissions please email: journals.permissions@oxfordjournals.org.

The Hox gene complement of a pelagic chaetognath, Flaccisagitta enflata

David Q. Matus1,*, Kenneth M. Halanych{dagger} and Mark Q. Martindale*
*Kewalo Marine Laboratory, University of Hawaii, 41 Ahui Street, Honolulu, HI 96813, USA; {dagger}Biological Sciences, Auburn University, 101 Rouse, Auburn, AL 36830, USA

Correspondence: 1E-mail: matus{at}hawaii.edu


    Synopsis
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
Chaetognaths are transparent marine animals that are ubiquitous and abundant members of oceanic zooplanktonic communities. Their phylogenetic position within the Metazoa, however, has remained obscure since their discovery. Morphology and embryology have traditionally allied chaetognaths with deuterostomes, but molecular evidence suggests otherwise. Two recent multigene expressed sequence tag (EST) molecular phylogenomic studies suggest that chaetognaths are either sister to the Lophotrochozoa (Matus et al. 2006) or to all protostomes (Marlétaz et al. 2006). We have isolated eight Hox genes, one Parahox gene, and Mox, a related homeodomain gene, from the pelagic chaetognath, Flaccisagitta enflata. Although chaetognath central class Hox genes lack the Lox5 or "spiralian" parapeptide, a diagnostic amino-acid motif that has been utilized previously to assign lophotrochozoan affinity, they do possess a central class Hox gene that has a partial "Ubd-A peptide" found in both ecdysozoan and lophotrochozoan Ubx/Abd-A/Lox2/Lox4 genes. Additionally, we report the presence of two distinct chaetognath posterior Hox genes that possess both ecdysozoan and lophotrochozoan signature amino-acid motifs. The phylogenetic position of chaetognaths, as well as the evolution of the Hox cluster, is discussed in light of these data.


    Introduction
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
Chaetognaths are active predators possessing a torpedo-like, coelomate body with lateral and caudal fins (Fig. 1A). They are hermaphroditic with direct development, lacking both a trochophore-like (i.e., lophotrochozoan) and dipleurula-like (i.e., deuterostome) larval stage. Several features (Fig. 1B) have allied them with deuterostomes (Willmer 1990Go), including a tripartite adult body plan, postanal tail, and developmental features such as radial cleavage, enterocoely, and posterior fate of the blastopore (deuterostomy). Molecular data (Telford and Holland 1993Go; Wada and Satoh 1994Go; Halanych 1996Go; Erber et al. 1998Go; Mallatt and Winchell 2002Go), "total evidence" (Zrzavy et al. 1998Go; Giribet et al. 2000Go; Peterson and Eernisse 2001Go) and mitochondrial genomic studies (Helfenbein et al. 2004Go; Papillon et al. 2004Go) have failed to support the inclusion of the chaetognaths within Deuterostomia. These, and other studies have been unable to discriminate whether chaetognaths are sister to protostomes (Giribet et al. 2000Go; Helfenbein et al. 2004Go) (Fig. 1C), within Ecdysozoa (Halanych 1996Go; Peterson and Eernisse 2001Go) (Fig. 1D) or within Lophotrochozoa (Erber et al. 1998Go; Haase et al. 2001Go; Shimotori and Goto 2001Go) (Fig. 1E). The most recent attempts to place the Chaetognatha have utilized expressed sequence tag (EST) data, both from the benthic chaetognath Spadella cephaloptera (Marlétaz et al. 2006Go) and from the pelagic chaetognath, Flaccisagitta enflata (Matus et al. 2006Go). While taxonomic sampling and gene number varied between these two studies, both definitively showed that chaetognaths are protostomes, although their suggested affinity within the protostomes differed. Based on 79 ribosomal genes, S. cephaloptera was shown to be sister to all protostomes (Marlétaz et al. 2006Go) (Fig. 1B), while analyses based on 72 genes by Phillippe et al. (2005Go), as well as increased taxonomic sampling using the single nuclear gene, tropomyosin, suggested that chaetognaths were sister to lophotrochozoans (Matus et al. 2006Go) (Fig. 1E).


Figure 1
View larger version (45K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1 The chaetognath body plan and possible phylogenetic affinities. (A) Micrograph of an 18 mm-long Flaccisagitta enflata collected from Kaneohe Bay, Oahu, (scale bar = 1 mm) and diagram of Sagitta sp. showing the deuterostome-like tripartite body, postanal tail, and other major anatomical features. (BE) Phylogenetic scenarios showing the possible affinities of the Chaetognatha. (B) Chaetognaths have been traditionally allied with the deuterostomes (Willmer 1990Go). (C) Chaetognaths as the sister group to Protostomia (Giribet et al. 2000Go; Helfenbein et al. 2004Go; Marlétaz et al. 2006Go). (D) Chaetognaths within or sister to the Ecdysozoa (Halanych 1996Go; Peterson and Eernisse 2001Go). (E) Chaetognaths within or sister to the Lophotrochozoa (Erber et al. 1998Go; Haase et al. 2001Go; Shimotori and Goto 2001Go; Matus et al. 2006Go).

 
In both of these recent EST analyses (Marlétaz et al. 2006Go; Matus et al. 2006Go) taxonomic sampling, particularly within the protostomes is sparse. As with many recent EST phylogenies (Philippe et al. 2005Go; Bourlat et al. 2006Go), Lophotrochozoa is only represented by Mollusca, Annelida, and Platyhelminthes, while Ecdysozoa is represented by Nematoda, Tardigrada, and Arthropoda. This represents only six of the ~26 phyla that make up these clades. Increased taxonomic sampling from within these groups might result in better resolution of the chaetognath phylogenetic position.

Another set of genes that have been used to gain insight into phylogenetic relationships is Hox genes. Hox genes are developmental regulatory genes involved in body-plan formation that are believed to have diversified before the last common bilaterian ancestor (Brooke et al. 1998Go; de Rosa et al. 1999Go; Kourakis and Martindale 2000Go; Balavoine et al. 2002Go). Their highly conserved nature has made them useful for revealing relationships of several enigmatic metazoan groups, e.g., brachiopods and priapulids (de Rosa et al. 1999Go) and dicyemids (Kobayashi et al. 1999Go). Orthological assignments of certain Hox genes can be determined by phylogenetic analysis and/or the presence of key diagnostic amino-acid residues (Sharkey et al. 1997Go; Anderson et al. 1998Go; de Rosa et al. 1999Go; Kobayashi et al. 1999Go; Kourakis and Martindale 2000Go; Telford 2000Go; Balavoine et al. 2002Go), which results in assignment to specific paralog groups (PG). While studying chaetognath development, we used degenerate PCR to sample the genome of the pelagic chaetognath, F. enflata for Hox genes to elucidate phylogenetic placement of arrow worms within the Metazoa. Chaetognaths possess definitive members from all four primordial bilaterian Hox paralogy classes (anterior, PG 3, central, and posterior), consistent with an ancestral bilaterian cluster of 8–10 genes (Fig. 2). Although genomic linkage data are not yet available, there is no evidence of cluster duplication in F. enflata.


Figure 2
View larger version (50K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2 The Metazoan Hox cluster, including data from the chaetognaths. (A) PCR survey for chaetognath Hox genes resulted in the identification of nine unique gene fragments. The tree on the left represents bilaterian phylogeny (Aguinaldo et al. 1997Go; de Rosa et al. 1999Go; Peterson and Eernisse 2001Go; Balavoine et al. 2002Go). B = Bilaterian common ancestor; E = Ecdysozoan stem; L = Lophptrochozoan stem; P = Protostome common ancestor; D = Deuterostome common ancestor. Horizontal black lines represent genomic linkage (when available). Uncertain orthology relationships are noted by question marks, dashed boxes indicate short PCR fragments, solid boxes indicate complete, or nearly complete, homeodomains. Hox data for two chaetognaths, Spadella cephaloptera (Papillon et al. 2003Go, 2005Go) and Flaccisagitta enflata (this study) suggest that chaetognaths possess at least ten Hox genes.

 

    Methods
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
Adult chaetognaths (F. enflata,) were obtained from Kaneohe Bay, Kaneohe, Hawaii, and kindly identified by Erik V. Thuesen (The Evergreen State College, Washington, USA). Hox and Parahox genes were initially isolated by PCR from genomic DNA using degenerate primers (see, Supplementary Materials). Flanking sequence was obtained with gene-specific primers by nested 5' and 3' RACE PCR using embryonic and adult cDNA as a template. Sequences have been reported to GenBank under accession numbers EF490810–EF490820. Putative orthological assignments were made via BLASTX searches of the GenBank. An amino-acid alignment of the 60 amino-acid homeodomain (HD) and the 12 amino acids directly flanking the 3' end of the homeodomain was made using representative sequences from acoel (Symsagittifera roscoffensis), ecdysozoan (Tribolium castaneum: beetle), lophotrochozoan (Nereis virens and Capitella capitata sp. I: annelids), and deuterostome (Branchiostoma floridae: cephalochordate) representatives (see, Supplementary Materials). A Bayesian phylogenetic analysis was conducted utilizing the MPI version of MrBayes 3.1.2 (Huelsenbeck and Ronquist 2001Go) on the University of Hawaii Dell computing cluster, using a mixed amino-acid model with gamma, which selected RtRev with a 100% posterior probability with 10,000,000 generations sampled every 100 generations with four chains over four independent runs. A summary tree was produced from the last 3,96,000 trees representing 39,600,000 stationary generations. Additionally, Neighbor Joining (NJ) (using mean AA distances) was conducted using phylogenetic analysis using parsimony (*and other methods; PAUP*) v4.0b10 (Swofford 2000Go). ProtTest (Abascal et al. 2005Go) also selected the rtrev + G model, which was used for maximum likelihood (ML) analyses conducted using RAXML v2.2.1 (Stamatakis 2006Go). An initial search was conducted in RAXML v2.2.1 (Stamatakis 2006Go) using 500 searches using the rtrev + G model. A consensus of this search is reported. ML bootstrap analysis was also conducted with RAXML v2.2.1 with 500 iterations and 15 independent searches, bootstrapped with SeqBoot (Felsenstein 1989Go), for each iteration. Additional details of phylogenetic analyses, GenBank accession numbers, alignments, and primer sequences are available in Supplementary Materials as well as being available upon request.


    Results and discussion
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
Previous work on the benthic chaetognath, S. cephaloptera, identified a chaetognath Hox cluster of at least six genes (Papillon et al. 2003Go), although for many of the genes only partial HDs were recovered. Due to the strong sequence conservation of the HD throughout metazoan evolution, however, obtaining flanking-sequence information is necessary to make more definitive orthological and paralogy assignments. Our initial PCR survey on F. enflata genomic DNA, utilizing primers specific to Hox-class HD genes, resulted in the identification of nine unique Hox gene fragments as well as the identification of the posterior ParaHox gene, caudal, and a HD gene related to the Hox cluster, mox (Minguillon and Garcia-Fernandez 2003Go). Using RACE PCR from adult and embryonic cDNA templates we isolated additional sequence information from eight of the nine Hox genes isolated in the initial PCR surveys from F. enflata genomic DNA. Based upon both predicted amino-acid sequences (Fig. 3) and a phylogenetic analysis of the HD and 3' flanking sequence (Fig. 4), F. enflata possesses orthologs to anterior (Fen Hox1), PG3 (Fen Hox3), four central class genes (Fen Hox4, Hox6, Hox7, and Hox8), two posterior class genes (Fen PostA and Fen PostB), and one gene (Fen MedPost) that appears to be the ortholog of the benthic chaetognath, S. cephaloptera MedPost gene. Failure to recover PG2 (Hox2 and proboscipedia) genes from the initial genomic DNA survey might not be surprising due to the presence of an intron in the HD of PG2 genes in many taxa (Kmita-Cunisse et al. 1998Go; Callaerts et al. 2002Go). Additionally, while flanking sequence for the PG3 gene, FenHox3, was unable to be recovered by RACE PCR, it shares 100% amino-acid identity and 85% nucleotide identity (12 changes, all third codon) with the ortholog isolated from the benthic chaetognath S. cephaloptera (Papillon et al. 2003Go) as well as identity with other bilaterian PG3 genes suggesting that it is correctly classified as a PG3 gene (SF1).


Figure 3
View larger version (66K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3 The chaetognath Hox gene fragments share amino-acid motifs found in ecdysozoan and lophotrochozoan Hox genes. Inferred amino-acid sequences of homeodomains (in bold) and N-terminal and C-terminal flanking regions for central class Paralog Group (PG) 4 (A), PG 6/7 (B), PG 8 (C) and posterior class (D) Hox genes. Chaetognath sequences are in red, ecdysozoan sequences are in black, lophotrochozoan genes are in dark blue, acoel flatworm sequences are in grey, and deuterostome sequences are in green. (A) The chaetognath Hox4 gene possesses the 3' LPTNK motif found in all PG4/Hox4 genes. (B) Chaetognath central class Hox genes (Fen Hox6 Fen Hox7 and SceMed4) (in red) lack the lophotrochozoan "lox 5 parapeptide" (de Rosa et al. 1999Go; Balavoine et al. 2002Go) (shaded in pink) found in all Lox5 genes. (C) Chaetognath MedPost genes lack the "UbdA parapeptide" (shaded in orange) (de Rosa et al. 1999Go; Balavoine et al. 2002Go), while Fen Hox8 shares four of seven amino acids in common with the "UbdA parapeptide." Both representative lophotrochozoan PG8 genes (Lox2 and Lox4) and ecdysozoan PG8 genes (Ubx and AbdA) are shown in the alignment. All chaetognath central class Hox genes share conserved central-class amino-acid motifs within their homeodomains (shaded in yellow), including the two chaetognath MedPost genes. (D) The chaetognath MedPost genes also share diagnostic amino-acid residues characteristic of posterior class Hox genes (shaded in green). Chaetognaths also possess two definitive posterior Hox genes (Fen PostA and PostB) that also share posterior-class diagnostic amino-acid residues (shaded in green). Both share a threonine at position 24 found in lophotrochozoan Post-1 and an acoel posterior gene (shaded in light blue). The chaetognath PostA gene shares an arginine at position 36 found in all ecdysozoan AbdB genes, an acoel flatworm posterior gene, and two hemichordate posterior genes (shaded in dark blue).

 

Figure 4
View larger version (37K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 4 Phylogenetic tree of metazoan Hox genes. Numbers above the branches represent a Bayesian analysis based on 3,96,000 trees from 39,600,000 stationary generations. Posterior probabilities are shown above the branches. An initial ML analysis was conducted using 500 separate searches to find the tree with the best likelihood score. A consensus was calculated from these searches. Support greater than 98% is indicated by a black circle at each node. Bootstrap support values (500 iterations) for the ML analysis are shown below the branches.

 
The phylogenetic analyses suggest orthology for all of the chaetognath Hox genes isolated in this study, with 100% posterior probability (pp) for Fen Mox, Fen Hox1, Fen MedPost, Fen Cdx, and the two posterior genes, Fen PostA and Fen PostB with known bilaterian PG and 95% pp for Fen Hox6 with a polychaete Lox5 gene (Fig. 4). While the remaining chaetognath Hox genes are not supported by >95% pp, their suggested orthology seems likely, as their position is supported by other methods of phylogenetic analysis, including ML and NJ. The Fen Hox4 clusters with other PG4 genes with 88% pp. The Fen Hox7 clusters with other protostome PG7 genes (Tribolium antp and Nereis Hox7 with 81% pp) and together is sister to deuterostome Branchiostoma Hox6 and Hox7 genes (83% pp). The Fen Hox8 is sister group to protostome and deuterostome PG8 genes (Ubx, Abd-A, Lox2, Lox4, and Branchiostoma Hox8; 79% pp). Finally, the two chaetognath MedPost genes, from F. enflata and S. cephaloptera cluster together with 100% pp, forming the sister group to posterior Hox and Parahox genes (81% pp) (Fig. 4).

PG4 genes (Dfd and Hox4) are characterized by the presence of a diagnostic amino-acid motif found immediately 3' to the HD (LPNTK). The chaetognath gene (Fen Hox4) possesses this motif, making this a definitive PG4 gene (Fig. 3A). Since this amino-acid motif is present in all metazoan PG4 genes surveyed to date (protostomes and deuterostomes), this motif has not been useful in assessing phylogenetic relatedness relative to other lineages. As it has already been established that the chaetognaths are protostomes (Helfenbein et al. 2004Go; Papillon et al. 2004Go; Marlétaz et al. 2006Go; Matus et al. 2006Go), of particular interest are those Hox genes that have been used to assign taxa to either Lophotrochozoa or Ecdysozoa, namely Lox5, Ubx, AbdA, Lox2, and Lox4 (de Rosa et al. 1999Go; Telford 2000Go; Balavoine et al. 2002Go).

The five central class Hox genes recovered from F. enflata, as well as those isolated from the benthic S. cephaloptera (Papillon et al. 2003Go), lack the "Lox5" or "spiralian" (de Rosa et al. 1999Go) peptide motif found in the 3' flanking region of all lophotrochozoan Lox5 genes, a central class gene considered to be orthologous to either ecdysozoan Fushi-tarazu (PG6) genes (Telford 2000Go) or to Antennapedia (PG7) genes (de Rosa et al. 1999Go; Balavoine et al. 2002Go) (Fig. 3B).

Another diagnostic amino-acid motif, the Ubd-A parapeptide, is found 3' to the HD in ecdysozoan Ubx and AbdA genes as well as in lophotrochozoan lox2 and lox4 genes (PG8) (de Rosa et al. 1999Go; Balavoine et al. 2002Go). One chaetognath central class Hox gene, Fen Hox 8, possesses four of the seven amino acids (QxIxEMx) found in the Ubd-A parapeptide (Figs. 3C and 5), which suggests that it may be an ortholog of protostome ubx/abdA/lox2/lox4 genes.


Figure 5
View larger version (25K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 5 Alignment of the consensus sequence of the protostome UbdA parapeptide. The consensus sequence from ecdysozoan Ubx and Abd-A genes, lophotrochozoan Lox2 and Lox4 genes, and the chaetognath Fen Hox8 sequence is shown. Fen Hox8 shares four of seven amino acids with other protostome UbdA parapeptides (shaded in grey).

 
Previous work on the benthic chaetognath, S. cephaloptera, identified at least six chaetognath Hox genes, including a PG3 gene, four central class genes of unknown orthology/paralogy, and a gene which the authors believed to possess both central and posterior class diagnostic residues (Spadella MedPost), but no definitive posterior class Hox genes (Papillon et al. 2003Go). When the amino-acid sequences of these six S. cephaloptera Hox gene fragments are compared to those recovered from F. enflata in this study (SF1 and Fig. 3) it appears that SceMed1 is an ortholog of FenHox7, SceMed2 is an ortholog of FenHox6, and SceMed3 is an ortholog of FenHox8. An ortholog of SceMed4, which has been suggested to be related to PG5 (sex combs reduced) genes (Papillon et al. 2005Go) does not appear to have been identified in F. enflata. This is consistent with phylogenetic analyses, which failed to identify a PG5 gene as well (Fig. 4).

Based upon predicted amino-acid sequence (Fig. 3D) and phylogenetic analyses (Fig. 4), the S. cephaloptera MedPost gene appears to be orthologous to the F. enflata MedPost gene and may be a chaetognath-specific Hox gene. Both MedPost genes share the central class diagnostic amino acids: QT (positions 6–7 in the HD), LTR(R/K)RRI (positions 26–32) and a glutamic acid (E) at position 59 (Papillon et al. 2003Go) (Fig. 3C and D). 5' of the HD, Fen MedPost has a definitive hexapeptide (SIYPWM), an amino-acid motif important functionally for interaction with extradenticle/pbx/homothorax proteins (Berthelsen et al. 1999Go; Ryoo et al. 1999Go), but missing from all posterior class genes. However, both chaetognath MedPost genes also possess diagnostic posterior class residues (K3, A14, R18, Y20, V21, and Q36), which are thought to be exclusive to posterior class genes with the exception of the Q36 in the priapulid Ubx gene (Fig. 3C and D).

Bayesian, ML, and NJ phylogenetic analyses place these chaetognath MedPost genes as sister to posterior class and cdx genes (Fig. 4), although with differing statistical support. While Bayesian posterior probability for this grouping is 81%, below the accepted 95% confidence level, it is still five times more likely than any other grouping for the chaeotognath MedPost genes. ML bootstrap support for this node was only 30% (see, Supplementary Materials). However, in the initial RAXML v 2.1.2 (Stamatakis 2006Go) searches, 99% of the 500 replicates performed identified a ML tree that placed MedPost genes as the sister group to posterior class Hox genes (Fig. 4 and Supplementary Materials). When those amino acids that are diagnostic of central class Hox genes were removed, both ML and Bayesian analyses place the MedPost genes within the posterior class (see, Supplementary Material). When posterior class diagnostic amino acids were removed, the MedPost genes remained as the sister group to the posterior genes (see, Supplementary Materials). Together, the results of the phylogenetic analyses suggest that chaetognath MedPost genes are definitive central class genes, but are consistently placed as the sister group to the posterior class Hox genes.

The MedPost gene has been suggested to retain features of an ancestral central/posterior class gene (Papillon et al. 2003Go). This was based largely on the absence of definitive posterior class genes from chaetognaths. We have isolated two posterior class genes using degenerate primers specific to posterior Hox genes. The two chaetognath posterior Hox genes (Fen PostA and Fen PostB) cluster with posterior genes with 100% pp in Bayesian phylogenetic analyses (Fig. 4). Additionally, they possess amino acids diagnostic of posterior class Hox genes, including K3 in PostA, and Y20 in PostA and PostB. Interestingly, both chaetognath posterior genes share a threonine at position 24, which is found in lophotrochozoan posterior1 and in the acoel, Symsagitiffera posterior Hox gene (Fig. 3D). The chaetognath PostA gene shares an arginine at position 36 that has been suggested to be ecdysozoan-specific (AbdominalB genes) (de Rosa et al. 1999Go; Balavoine et al. 2002Go; Cook et al. 2004Go), but further taxonomic sampling has shown it to be found in an acoel (Symsagitiffera Post) and two hemichordate posterior genes (Ptychodera Hox9/10 and Hox11/13a) but is absent from lophotrochozoan posterior class genes (Fig. 3D). Thus, the two chaetognath posterior Hox genes share plesiomorphic characters with lophotrochozoans, ecdysozoans, and deuterostomes making phylogenetic interpretations difficult. The presence of definitive posterior class genes in chaetognaths, however, suggests that the MedPost genes are more likely to be a chaetognath-specific innovation rather than evolutionary antecedents of central and posterior class genes, as previously suggested by others (Papillon et al. 2003Go). Additionally, the identification of posterior class Hox genes in both cnidarians (Ryan et al. 2007Go) and acoels (Cook et al. 2004Go), organisms believed to have diverged prior to chaetognaths (Collins et al. 2005Go; Jondelius et al. 2002Go; Ruiz-Trillo et al. 1999Go, 2002Go) during metazoan evolution, further bolsters the conclusion that chaetognath MedPost genes are not directly descended from the Hox posterior-central class gene ancestor.

The absence of a definitive Lox5 motif and the presence of a partial UbdA parapeptide in a chaetognath have implications for both the evolution of the bilaterian Hox cluster and the phylogenetic position of the chaetognaths and can be interpreted in two ways. First, the chaetognaths may have diverged prior to the lophotrochozoan/ecdysozoan split (Fig. 6A). Alternatively, because the plesiomorphic condition of both the Lox5 and the UbdA peptides is unclear, we cannot rule out that both the lack of diagnostic residues (Lox5) and divergence (UbdA) could be due to secondary changes unique to the chaetognath lineage (Fig. 6B). In such a case, Chaetognatha may be sister to Lophotrochozoa but the phylogenetic signal within their Hox genes has been erased. Additionally, due to the nature of PCR-biased sampling, one cannot rule out that a second gene containing a UbdA parapeptide exists in chaetognaths. However, the presence of one gene that only possesses a partial UbdA parapeptide suggests that the protostome PG8 genes, Ubx and AbdA (ecdysozoan) and Lox2 and Lox4 (lophotrochozoan) duplicated after their divergence from the chaetognaths (Fig. 6A).


Figure 6
View larger version (16K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 6 Hypothetical scenarios for the evolution of the Hox cluster based upon alternative positions of the Chaetognatha. The presence of anterior (PG1 and PG2) and posterior class Hox genes in the Cnidaria (Ryan et al. 2007Go) suggests that the cnidarian–bilaterian ancestor possessed at least two Hox genes (1). The earliest branching lineage to have definitive central class genes are the acoel flatworms, which have a Hox cluster of at least three genes (2). The protostome–deuterostome ancestral Hox cluster underwent diversification of central class genes, including a gene possessing the Hox4 LPTNK motif (absent from acoels) (3). (A) Chaetognaths are the sister group to the rest of the protostomes (Marlétaz et al. 2006Go). The chaetognath–protostome ancestral Hox cluster would have consisted of a minimum of eight genes, and included a central class gene possessing a UbdA parapeptide (4). Following a split from the Chaetognatha, the genes possessing a UbdA parapeptide duplicated in their respective lineages to yield Ubx/AbdA and Lox2/Lox4 (5). A central class gene containing the Lox5/spiralian parapeptide arose in the lophotrochozoan ancestor (6). Chaetognaths evolved a MedPost gene with central and posterior class attributes after splitting from the rest of the protostomes (7). (B) Chaetognaths are the sister group to the lophotrochozoans (Matus et al. 2006Go). The ecdysozoan–lophotrochozoan ancestor has a Hox cluster of at least eight genes including a gene containing the UbdA parapeptide. This duplication of this gene occurs separately in lophotrochozoans into Lox2 and Lox4 after diverging from chaetognaths (8) and ecdysozoans into Ubx and AbdA (9). The Lox5 parapeptide evolves after the split between the Chaetognatha and the rest of the Lophotrochozoa (6).

 
Given that recent phylogenetic EST analyses places chaetognaths as either sister to protostomes (Marlétaz et al. 2006Go) or sister to lophotrochozoans (Matus et al. 2006Go), the question arises as to whether the data from chaetognath Hox genes sheds any light on these alternative hypotheses. Chaetognaths clearly possess a near complete set of Hox genes, with members from all four classes (anterior, PG3, central, and posterior) as well as a potentially unique class, the MedPost gene, which supports previous theories of an ancestral bilaterian Hox cluster of 8–10 genes (de Rosa et al. 1999Go; Kourakis and Martindale 2000Go; Balavoine et al. 2002Go). One explanation for the data presented here is a sister-group relationship between chaetognaths and all other protostomes. Cnidarians have been shown to possess both anterior class (PG1 and PG2) and posterior class genes (Ryan 2007Go: 2025), suggesting that the cnidarian–bilaterian ancestor possessed at least three Hox genes (two anterior and one posterior) (Fig. 6). Acoels are the earliest branching metazoan that possesses definitive central class Hox genes (Cook et al. 2004Go). The protostome–deuterostome ancestor would have then undergone expansion of central and anterior class Hox genes, including the innovation of a definitive Hox4/dfd paralog with the LPTNK motif found 3' of all Hox4 genes (but absent from the acoel central class gene) (Fig. 3A). The chaetognath–protostome ancestor would then have possessed a Hox cluster of at least eight genes including a central class gene possessing an Ubd-A parapeptide. The gene containing this motif would then have duplicated into Ubx/AbdA and Lox2/Lox4 following the split from the chaetognaths. The Lox5 peptide would have evolved in lophotrochozoans following the split from ecdysozoans. Chaetognaths likely evolved a MedPost gene specific to their lineage since a clear ortholog of this gene has not been recovered from any other metazoan (Fig. 6A).

If chaetognaths are actually the sister group to lophotrochozoans as suggested by EST analysis, mitochondrial genomic analyses, and tropomyosin phylogeny (Matus et al. 2006Go), this would suggest the scenario presented in Fig. 6B. A central Hox gene containing the "Ubd-A parapeptide" would have evolved at the base of the protostomes, and genes containing this motif would have diverged subsequently in ecdysozoans (Ubx and AbdA) and the remainder of the lophotrochozoans (Lox2 and Lox4), to the exclusion of chaetognaths. The Lox5 parapeptide would have evolved after the chaetognath–lophotrochozoan split (Fig. 6B). Unfortunately, due to the nature of degenerate PCR, it is possible that we did not isolate all of the Hox genes in the chaetognath genome, which could have profound impacts on the interpretation of these results. Total genome sequencing of a chaetognath will greatly aid in the characterization of the chaetognath Hox cluster and perhaps shed further light on the evolution of Hox genes in the protostomes.

Placement as either sister to protostomes or sister to Lophotrochozoa raises interesting issues about several key morphological and embryological innovations. The unique set of chaetognath morphological and developmental characters, such as radial cleavage, enterocoely, fate of the blastopore, brachyury staining (Takada et al. 2002Go), direct development, a tripartite body plan, and protostome-like central nervous system (Hyman 1959Go; Willmer 1990Go), may provide insight into the protostomian or lophotrochozoan ancestor. If these features are bilaterian plesiomorphies, the conflict between chaetognaths possessing "deuterostome-like" morphology and embryology can be reconciled with modern molecular phylogenetic approaches. The observation that chaetognaths may possess a single complete complement of Hox genes and that adults and embryos are readily available provides an opportunity to better understand morphological, genomic, and developmental evolution at a key node in the Metazoa.


    Supplementary data
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
Supplementary data are available at ICB online.


    Acknowledgments
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
We wish to thank Allison Sweeney and Sonke Johnson for the opportunity to participate in this symposium. We would also like to thank Deirdre Killebrew and the University of Hawaii Dell Cluster for computational resources, and Casey Dunn and Joe Ryan for help with phylogenetic analyses. This work was supported by NSF support to M.Q.M. (IBN-0236565) and K.H.M. (IBN-0333843). This work is AU Marine Biology Program contribution 26.

Conflict of interest: None declared.


    Footnotes
 
From the symposium "Integrative Biology of Pelagic Invertebrates" presented at the annual meeting of the Society for Integrative and Comparative Biology, January 3–7, 2007, at Phoenix, Arizona.


    References
 Top
 Synopsis
 Introduction
 Methods
 Results and discussion
 Supplementary data
 Acknowledgments
 References
 
Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics (2005) 21::2104–5.[Abstract/Free Full Text]

Aguinaldo AA, Turbeville JM, Linford LS, Rivera MC, Garey JR, Raff RA, Lake JA. Evidence for a clade of nematodes, arthropods and other moulting animals. Nature (1997) 387::489–93.[CrossRef][Medline]

Anderson CL, Canning EU, Okamura B. A triploblast origin for Myxozoa? Nature (1998) 392::346–7.[CrossRef][Medline]

Balavoine G, de Rosa R, Adoutte A. Hox clusters and bilaterian phylogeny. Mol Phylogenet Evol (2002) 24::366–73.[CrossRef][Web of Science][Medline]

Berthelsen J, Kilstrup-Nielsen C, Blasi F, Mavilio F, Zappavigna V. The subcellular localization of PBX1 and EXD proteins depends on nuclear import and export signals and is modulated by association with PREP1 and HTH. Genes Dev (1999) 13::946–53.[Abstract/Free Full Text]

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–8.[CrossRef][Medline]

Brooke NM, Garcia-Fernandez J, Holland PWH. The ParaHox gene cluster is an evolutionary sister of the Hox gene cluster. Nature (1998) 392::920–2.[CrossRef][Medline]

Callaerts P, Lee PN, Hartmann B, Farfan C, Choy DW, Ikeo K, Fischbach KF, Gehring WJ, de Couet HG. HOX genes in the sepiolid squid Euprymna scolopes: implications for the evolution of complex body plans. Proc Natl Acad Sci USA (2002) 99::2088–93.[Abstract/Free Full Text]

Collins AG, Cartwright P, McFadden CS, Scheirwater B. Phylogenetic context and Basal Metazoan Model systems. Integr Comp Biol (2005) 45::585–94.[Abstract/Free Full Text]

Cook CE, Jimenez E, Akam M, Salo E. The Hox gene complement of acoel flatworms, a basal bilaterian clade. Evol Dev (2004) 6::154–63.[CrossRef][Web of Science][Medline]

de Rosa R, Grenier JK, Andreeva T, Cook CE, Adoutte A, Akam M, Carroll SB, Balavoine G. Hox genes in brachiopods and priapulids and protostome evolution. Nature (1999) 399::772–6.[CrossRef][Medline]

Erber A, Riemer D, Bovenschulte M, Weber K. Molecular phylogeny of metazoan intermediate filament proteins. J Mol Evol (1998) 47::751–62.[CrossRef][Web of Science][Medline]

Felsenstein J. PHYLIP-phylogeny inference package. Cladistics (1989) 5::164–6.

Giribet G, Distel DL, Polz M, Sterrer W, Wheeler WC. Triploblastic relationships with emphasis on the acoelomates and the position of gnathostomulidam cycliophora, plathelminthes, and chaetognatha: a combined approach of 18S rDNA sequences and morphology. Syst Biol (2000) 49::539–62.[Abstract/Free Full Text]

Haase A, Stern M, Wächtler K, Bicker G. A tissue-specific marker of Ecdysozoa. Dev Genes Evol (2001) 211::428–33.[CrossRef][Web of Science][Medline]

Halanych KM. Testing hypotheses of chaetognath origins: long branches revealed by 18S ribosomal DNA. Syst Biol (1996) 45::223–46.[Abstract/Free Full Text]

Helfenbein KG, Fourcade HM, Vanjani RG, Boore JL. The mitochondrial genome of Paraspadella gotoi is highly reduced and reveals that chaetognaths are a sister group to protostomes. Proc Natl Acad Sci USA (2004) 101::10639–43.[Abstract/Free Full Text]

Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics (2001) 17::754–5.[Abstract/Free Full Text]

Hyman LH. The enterocoelus coelomates - phylum chaetognatha. The invertebrates: smaller coelomate groups. (1959) 5:. New York: McGraw-Hill. 1–66.

Jondelius U, Ruiz-Trillo I, Baguña J, Riutort M. The nemertodermatids are basal bilaterians and not members of the pllatyhelminthes. Zool Scripta (2002) 31::201–15.[CrossRef][Web of Science]

Kmita-Cunisse M, Loosli F, Bierne J, Gehring WJ. Homeobox genes in the ribbonworm Lineus sanguineus: evolutionary implications. Proc Natl Acad Sci USA (1998) 95::3030–5.[Abstract/Free Full Text]

Kobayashi M, Furuya H, Holland PWH. Dicyemids are higher animals. Nature (1999) 401::762.[Medline]

Kourakis MJ, Martindale MQ. Combined-method phylogenetic analysis of Hox and Parahox genes of the metazoa. J Exp Zool (Mol Dev Evol) (2000) 288::175–91.[CrossRef]

Mallatt J, Winchell CJ. Testing the new animal phylogeny: first use of combined large-subunit and small-subunit rRNA gene sequences to classify the protostomes. Mol Biol Evol (2002) 19::289–301.[Abstract/Free Full Text]

Marlétaz F, Martin E, Perez Y, Papillon D, Caubit X, Lowe CJ, Freeman B, Fasano L, Dossat C, Wincker P, Weissenbach J, Le Parco Y. Chaetognath phylogenomics: a protostome with deuterostome-like development. Curr Biol (2006) 16::R577–8.[CrossRef][Web of Science][Medline]

Matus DQ, Copley RR, Dunn CW, Hejnol A, Eccleston H, Halanych KM, Martindale MQ, Telford MJ. Broad taxon and gene sampling indicate that chaetognaths are protostomes. Curr Biol (2006) 16::R575–6.[CrossRef][Web of Science][Medline]

Minguillon C, Garcia-Fernandez J. Genesis and evolution of the Evx and Mox genes and the extended Hox and ParaHox gene clusters. Genome Biol (2003) 4::R12.[CrossRef][Medline]

Papillon D, Perez Y, Caubit X, Le Parco Y. Identification of chaetognaths as protostomes is supported by the analysis of their mitochondrial genome. Mol Biol Evol (2004) 21::2122–9.[Abstract/Free Full Text]

Papillon D, Perez Y, Caubit X, Le Parco Y. Hox gene survey in the chaetognath Spadella cephaloptera: evolutionary implications. Dev Genes Evol (2003) 213::142–8.[Web of Science][Medline]

Papillon D, Perez Y, Fasano L, Le Parco Y, Caubit X. Restricted expression of a median Hox gene in the central nervous system of chaetognaths. Dev Genes Evol (2005) 215::369–73.[CrossRef][Web of Science][Medline]

Peterson KJ, Eernisse DJ. Animal phylogeny and the ancestry of bilaterians: inferences from morphology and 18S rDNA gene sequences. Evol Dev (2001) 3::170–205.[CrossRef][Web of Science][Medline]

Philippe H, Lartillot N, Brinkmann H. Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia. Mol Biol Evol (2005) 22::1246–53.[Abstract/Free Full Text]

Ruiz-Trillo I, Riutort M, Littlewood DTJ, Herniou EA, Baguñà J. Acoel flatworms: earliest extant bilaterian metazoans, not members of Platyhelminthes. Science (1999) 283::1919–23.[Abstract/Free Full Text]

Ruiz-Trillo I, Paps J, Loukota M, Ribera C, Jondelius U, Baguñà J, Riutort M. A phylogenetic analysis of Myosin heavy chain type II sequences corroborates that Acoela and Nemertodermatida are basal bilaterians. Proc Natl Acad Sci USA (2002) 99::11246–51.[Abstract/Free Full Text]

Ryan JF, Mazza ME, Pang K, Matus DQ, Baxevanis AD, Martindale MQ, Finnerty JR. Pre-bilaterian origins of the Hox cluster and the Hox code: evidence from the sea anemone, Nematostella vectensis. PLoS ONE (2007) 2::e153.[CrossRef]

Ryoo HD, Marty T, Casares F, Affolter M, Mann RS. Regulation of Hox target genes by a DNA bound homothorax/hox/extradenticle complex. Development (1999) 126::5137–48.[Abstract]

Sharkey M, Graba Y, Scott MP. Hox genes in evolution: protein surfaces and paralog groups. Trends Genet (1997) 13::145–51.[CrossRef][Web of Science][Medline]

Shimotori T, Goto T. Developmental fates of the first four blastomeres of the chaetognath Paraspadella gotoi: relationship to protostomes. Dev Growth Differ (2001) 43::371–82.[CrossRef][Web of Science][Medline]

Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics (2006) 22::2688–90.[Abstract/Free Full Text]

Swofford DL. PAUP* Phylogenetic analysis using parsimony *and other methods. (2000) Sunderland, MA: Sinauer Associates.

Takada N, Goto T, Satoh N. Expression pattern of the Brachyury gene in the arrow worm Paraspadella gotoi (Chaetognatha). Genesis (2002) 32::240–5.[CrossRef][Web of Science][Medline]

Telford MJ. Turning "Hox signatures" into synapomorphies. Evol Dev (2000) 2::360–4.[CrossRef][Web of Science][Medline]

Telford MJ, Holland PH. The phylogenetic affinities of the chaetognaths: a molecular analysis. Mol Biol Evol (1993) 10::660–76.[Abstract]

Wada H, Satoh N. Details of the evolutionary history from invertebrates to vertebrates, as deduced from the sequences of 18S rDNA. Proc Natl Acad Sci USA (1994) 91::1801–4.[Abstract/Free Full Text]

Willmer P. Invertebrate relationships: patterns in animal evolution. (1990) Cambridge: Cambridge University Press.

Zrzavy J, Mihulka S, Kepka P, Bezdek A, Tietz D. Phylogeny of Metazoa Based on Morphological and 18S Ribosomal DNA Evidence. Cladistics (1998) 14::249–85.[Web of Science]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary data
Right arrow All Versions of this Article:
47/6/854    most recent
icm077v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Matus, D. Q.
Right arrow Articles by Martindale, M. Q.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?