The homeodomain complement of the ctenophore Mnemiopsis leidyi suggests that Ctenophora and Porifera diverged prior to the ParaHoxozoa
© The article is a work of the United States Government. 2010
Received: 4 June 2010
Accepted: 4 October 2010
Published: 4 October 2010
The much-debated phylogenetic relationships of the five early branching metazoan lineages (Bilateria, Cnidaria, Ctenophora, Placozoa and Porifera) are of fundamental importance in piecing together events that occurred early in animal evolution. Comparisons of gene content between organismal lineages have been identified as a potentially useful methodology for phylogenetic reconstruction. However, these comparisons require complete genomes that, until now, did not exist for the ctenophore lineage. The homeobox superfamily of genes is particularly suited for these kinds of gene content comparisons, since it is large, diverse, and features a highly conserved domain.
We have used a next-generation sequencing approach to generate a high-quality rough draft of the genome of the ctenophore Mnemiopsis leidyi and subsequently identified a set of 76 homeobox-containing genes from this draft. We phylogenetically categorized this set into established gene families and classes and then compared this set to the homeodomain repertoire of species from the other four early branching metazoan lineages. We have identified several important classes and subclasses of homeodomains that appear to be absent from Mnemiopsis and from the poriferan Amphimedon queenslandica. We have also determined that, based on lineage-specific paralog retention and average branch lengths, it is unlikely that these missing classes and subclasses are due to extensive gene loss or unusually high rates of evolution in Mnemiopsis.
This paper provides a first glimpse of the first sequenced ctenophore genome. We have characterized the full complement of Mnemiopsis homeodomains from this species and have compared them to species from other early branching lineages. Our results suggest that Porifera and Ctenophora were the first two extant lineages to diverge from the rest of animals. Based on this analysis, we also propose a new name - ParaHoxozoa - for the remaining group that includes Placozoa, Cnidaria and Bilateria.
Basic Local Alignment Search Tool
expressed sequence tag
polymerase chain reaction
phylogenetic inferences using ML
randomized accelerated ML
Ctenophores are a phylum of marine metazoans with uncertain phylogenetic affinity. Their signature morphological features include a set of eight ciliated comb rows that are used for swimming; these are controlled by an aborally located statocyst called the apical sense organ. Most ctenophores have a pair of feeding tentacles that contain specialized adhesive cells called colloblasts. Ctenophores are extremely fragile and difficult to culture and, as such, we know very little about their biology relative to other metazoans .
The unique ctenophore body plan has made it difficult to untangle its phylogenetic position in relation to other animal phyla. The earliest comparative classifications by Cuvier allied ctenophores with cnidarians and echinoderms in the Radiata . Later, Leuckart grouped ctenophores with sponges and cnidarians in the Coelenterata . Associations proposed between ctenophores and other taxa include groupings with Platyhelminthes, trochozoans, bilaterians and subsets of sponges and cnidarians (see  for a review). Many of the early molecular studies using 18 s ribosomal RNA (rRNA) sequences placed ctenophores sister to Placozoa, Cnidaria and Bilateria (for example, Wainright et al., 1993; Smothers et al. 1994; Bridge et al. 1995; Collins 1998; Kim et al. 1999 [5–9]).
The 18 S rRNA placement of ctenophores has recently been challenged using data generated in expressed sequence tag (EST)-based phylogenomic studies. One study re-allies ctenophores with the cnidarians , while another has ctenophores branching at the base of the animal tree [11, 12]. Yet another study combined morphological, structural and sequence data, leading to the placement of ctenophores in a clade with all other non-bilaterians sister to the Bilateria . These conflicting results could be due to the use of different methods, the inherent incomplete nature of transcript sequencing (in the case of the EST-based studies) or for other reasons. This series of studies have left most investigators waiting for further evidence to tilt the consensus convincingly in one direction or another.
Phylogenomic approaches hold the promise of reconstructing the true tree of life (reviewed in ). Thus far, most phylogenomic efforts that have included data from all four of the early branching phyla have been restricted to the aforementioned (and conflicting) EST-based analyses (for example [10–12, 15]). Methods based on whole-genome content rather than large concatenated data matrices can provide an independent assessment of current phylogenetic hypotheses and, due to the rarity of events measured, may arguably be more appropriate in this context [14, 16]. One such rare genomic change that has previously been employed involves the presence or absence of gene duplications of homeobox genes [17–19].
Homeobox genes encode transcription factor proteins characterized by the presence of a helix-loop-helix DNA-binding domain called the homeodomain . Homeobox genes were present in the last common ancestor of plants, animals and fungi and underwent extensive independent diversification in each of these lineages [21, 22]. In animals, the homeobox superfamily has been separated into 11 classes and more than 125 gene families [21, 23, 24].
Examination of the homeobox complement of species from early-branching metazoan phyla (such as Cnidaria [25, 26], Placozoa [27, 28] and Porifera [29, 30]) has been an especially fertile area of research, one that has been fuelled by the recent availability of full genomic sequence data from several non-bilaterian genomes. The last remaining non-bilaterian phylum lacking a species with a sequenced genome (and, therefore, a completely examined homeobox repertoire) was Ctenophora.
We have used a next-generation sequencing approach to sequence and assemble the ~150 MB of the lobate ctenophore, Mnemiopsis leidyi. Here, we present the first whole-genome investigation of the ctenophore homeobox superfamily. Our results expand on previous studies that have utilized degenerate polymerase chain reaction (PCR) approaches to identify ctenophore homeoboxes from the PRD, ANTP and SINE classes [22, 31–36]. In addition to these classes, we show that the POU, LIM and TALE classes were also present prior to the divergence of the ctenophores from the rest of Metazoa.
This is also the first study to compare the complete homeobox catalogue of species from all of the non-bilaterian phyla, along with that of the two major bilaterian lineages (Protostomia and Deuterostomia) where complete genomic sequence data is available. As such, this work provides a major missing piece of evidence that is critical to understanding the makeup of the homeodomain superfamily in early metazoan history. With these data in hand, we evaluate the congruency of the homeodomain data with the recently proposed phylogenetic relationships of the early branching phyla.
Overview of homeobox genes in Mnemiopsis
Mnemiopsis homeobox genes.
Most of these class-level assignments are confirmed by the presence of secondary domains, sequence signatures, and/or class-specific introns (Table 1). To all of these classes (with the exception of the 18 homeodomains that remained unclassified), we added corresponding homeodomain sequence data from the demosponge Amphimedon queenslandica , the placozoan Trichoplax adhaerens , the cnidarian Nematostella vectensis and the choanoflagellate Monosiga brevis : we then performed class-specific phylogenetic analyses. We named Mnemiopsis homeodomains that showed a strong affiliation for a particular family accordingly; otherwise, the name of the class is used in conjunction with a preliminary number that was originally assigned to the homeodomain.
ANTP class NKL subclass
The following groupings are consistent in both trees and have support values over 50 in our best Bayesian tree: (1) MlANTP65 with the Tlx family and (2) MlANTP22 with the Human Ventx gene. MlANPT25, MlANTP35, MlANTP67 and MlANTP78 group with the Dlx family consistently in both trees but, in the full tree, the Dlx clade also includes MlANTP66. Similarly, MlANTP19, MlANTP68, MlANTP71 and MlANTP72 form clades positioned sister to the Barx, Bsx, Dbx and Hlx familes. However, the relationships between Mnemiopsis homeodomains is inconsistent between these trees. The other NKL homeodomains identified are MlANTP21, MlANTP23, MlANTP25, MlANTP35, MlANTP37, MlANTP47, MlANTP48, MlANTP51, MlANTP63, MlANTP66, MlANTP67 and MlANTP78.
Consistent with our analysis, a previous study classified MlANTP65 as a Tlx-like homeodomain . The same study also associated MlANTP66 with the Dlx family, MlANTP67 with the Barh family and MlANTP68 with the Bsx family, an observation that was not consistently reproduced in our trees.
Evidence in the form of diagnostic residues can provide additional support to classifications . The following homeodomains all contain the diagnostic residues associated with the NKL subclass ([AKST][DENPS][LAST][Q][V] at positions 41-45): MlANTP19, MlANTP37, MlANTP48, MlANTP68, MlANTP71 and MlANTP72. The only other Mnemiopsis homeodomain with the NKL signature is the SINE class homeodomain MlSIX36. (The position of MlSIX36 on the tree in Figure 1 and its upstream SIX domain led to its SINE class designation.) In addition to the NKL signature, the MlANTP37 homeodomain also contains the HOXL signature ([KT][IV]WFQNRR[AMV]K[DEHKLMQWY][KR][KR] at positions 46-58) and the MlANTP48 homeodomains contains the HOXL2 signature (LE[AGKNR]E at positions 16-19) (Table 1).
ANTP class HOXL-related
Four paralogous Mnemiopsis ANTP homeodomains (MlANTP03a, MlANTP03b, MlANTP03c and MlANTP03d) group with the engrailed family in our superfamily tree (Figure 1) and with the Evx family in the ANTP tree (Figure 2). Despite the engrailed family being assigned to the NKL subclass in HomeoDB , engrailed has been historically allied with the extended Hox subclass based on synteny [41, 42] and phylogeny . Evx is also considered a member of the extended Hox subclass. While it is difficult to pin down the exact relationship of the MlANTP03 homeodomains, it does appear that they are the most likely descendants of the homeodomain that gave rise to the HOXL genes in the lineage leading to Placozoa, Cnidaria and Bilateria. Consistent with this classification, the MlANTP03a and MlANTP03b genes both contain the HOXL2 diagnostic residue signature. There are no clear ParaHox or Hox genes in Mnemiopsis.
Of the seven Mnemiopsis PRD class homeodomains, six have a Q at position 50. The exception (MlPRD43) is missing sequence information at that position. (Note: just prior to the submission of this manuscript, a new assembly has revealed the likely 3' end of this homeodomain that includes a Q at position 50). We did not find any Mnemiopsis genes with an S at position 50. The only other Mnemiopsis genes with a K at position 50 are the 18 SINE class genes that, like the K50 PRD class genes, also characteristically have a K residue at position 50.
Consistent with the absence of lysine or serine residues at position 50 in Mnemiopsis and Amphimedon PRD homeodomains, we see no grouping of Mnemiopsis or Amphimedon homeodomains with S50 or K50 clades, with the following exceptions: (1) MlPRD16 groups with the Nematostella S50 homeodomain NvPRD074, albeit with virtually no support (ML bootstrap = 2, Bayesian posterior probability distribution = 2), within a larger clade of Q50 homedomains; and (2) the Amphimedon homeodomain AqQ50a groups with the highly divergent HsDUXBI (ML bootstrap = 27, Bayesian posterior probability distribution = 95), also within a larger clade of Q50 homedomains (Figure 1). The overwhelming evidence suggests that Mnemiopsis and Amphimedon are devoid of S50 and K50 PRD class homeodomains. Conversely, Nematostella and Trichoplax both have clear K50 and S50 homeodomains. The phylogenetic distribution of Q50, S50 and K50 PRD homeodomains in our study is consistent with the hypothesis that Q50 homeodomains were the founders of the PRD class .
Five of the seven Mnemiopsis PRD class homeodomains contain the diagnostic residues (L[EINQRV][^DGHMPTVWY][^CDGKMNPQR][FL][^CFILPTWY][AEFHKQRV][ADEGKNSTW][CHKMPQR][FHY]P at positions 16-26) associated with paired homeodomains in bilaterians: MlPRD10a, MlPRD10b, MlPRD16, MlPRD44, MlPRD61 (Table 1). No other Mnemiopsis homeodomains display this pattern. MlPRD10b, MlPRD16 and MlPRD44 had been identified as Paired class genes in a previous study and were named Prd3, Prd1 and Prd2 respectively . MlPRD44 also contains the HOXL2 diagnostic residues (Table 1).
MlPRD16 and MlPRD61 have clear octapeptide sequences upstream of the homeodomain (SSISSLLS and HSIDDILG, respectively), a hallmark characteristic of a subset of the PRD class homeodomains. MlPRD43 and MlPRD50 have less-conserved but possible octapeptides as well (QRILGILS and YNIEGLLG, respectively). There are no paired domains associated with any Mnemiopsis homeodomains, but there are two independent paired domain sequences that appear to be direct orthologs of the two identified in the ctenophore Coeloplana willeyi .
Like most PRD class homedomains , all but one of the Mnemiopsis PRD homeodomains have an intron that occurs in the vicinity of the 46th and 47th codons. The one exception, MlPRD10a, has a single intron that interrupts the 37th codon. This might be the result of a retrotransposition event involving a transcript from its paralog (MlPRD10b) followed by an intron gain event. There are additional introns in the N-termini of the homeodomains of MlPRD10b, MlPRD50, and MlPRD61.
MlPbx, MlMeis and MlPknox, like other TALE class homeodomains, have a three amino acid insertion in the loop between the first and second alpha-helices (Table 1). MlPbx, MlMeis and MlPknox consistently group with the Pbx, Meis and Pknox families, respectively, in both trees with moderate support (Figures 1 and 4B). In all three cases, the phylogenetic assignment of these homeodomains is reinforced by the identification of several conserved motifs outside of the homeodomain, as well as by conserved intron positions (Table 1).
Like other Pbx genes (and unlike other TALE genes), MlPbx has a glycine residue at position 50 of the homeodomain. In addition, a Basic Local Alignment Search Tool (BLAST) search to the contig containing MlPBX shows significant similarity to the PBC domain  located ~1.5 KB upstream of the homeodomain, as assessed by BLAST [percent identity (ID)= 25/67, expectation (E)-value = 2 × 10-7). Like the cnidarian and human PBX genes, MlPbx has an intron that interrupts the second codon and one that falls between the 47th and 48th codon of the 63-codon TALE homeobox.
Meis homeodomain proteins have several conserved motifs in addition to the homeodomain . A GENSCAN prediction containing the MlMeis homeodomain shows similarity to the upstream MEIS A domain (ID = 19/69, E-value = 0.005), as well as weaker similarity to the MEIS D domain downstream of the homeodomain (ID = 16/48, E-value = 0.014). Similar to bilaterians and cnidarians, MlMeis has two introns. One falls between the 25th and 26th codons, while another interrupts the 51st codon.
The GENSCAN-predicted peptide that contains the MlPknox homeodomain also includes the abbreviated MEIS A domain that is characteristic of the Pknox family, as well as the MEIS B motif (ID = 33/139, E-value = 8 × 10-5). MlPknox, like the human PKNOX1 and PKNOX2 genes, has an intron that separates the 25th and 26th codons and one that interrupts the 51st codon of the homeobox.
We were unable to identify an Irx homeodomain in Mnemiopsis, despite there being Irx family members from Amphimedon, Trichoplax and Nematostella. Also absent was the Tgif homeodomain, found only in cnidarians and bilaterians.
MlIsl, MlLhx1.5, MlLhx3.4 and MlLmx make up the four LIM class homeodomains of Mnemiopsis (Figure 4C). We assigned these four Mnemiopsis homeodomains to the Isl, Lhx1/5, Lhx3/4 and Lmx families, respectively, based on the consistency between tree runs (Figures 1 and 4C) and moderate support in the full homeodomain tree (Figure 1). BLAST searches of the genomic scaffolds containing Mnemiopsis LIM homeodomains reveal LIM-type zinc finger domains immediately upstream of these four homeodomains. Additional BLAST searches also reveal traces of LIM domains independent of homeodomains in the Mnemiopsis genome (data not shown), suggesting the existence of LIM domain transcription regulator genes.
Eighteen SINE class homeodomains representing seven distinct SINE lineages were recovered from the Mnemiopsis genome (Figure 4D, Table 1). Of these, all but one have the characteristic lysine at position 50 (as described in ). The exception is MlSIX41, for which we are missing the sequence information from the C-terminus of the homeodomain (including position 50). Additionally, 17 of the 18 SINE class homeodomains have a SIX domain upstream of the homeodomain. The exception, MlSIX32a, is situated on the N-terminal end of a small scaffold in our current assembly, so its absence may be due to the resolution of our assembly.
The SINE class is monophyletic in our superfamily tree except for a clade of five Mnemiopsis homeodomains (MlSIX59a, MlSIX59b, MlSIX59c, MlSIX59 d and MlSIX59e), which group with Zhx/Homez (Figure 1). This exception is perhaps not completely unexpected given that, like the Zhx/Homez genes, MlSIX59 homeodomains are quite divergent; they are five of only six homeodomains in our entire Mnemiopsis set that do not include a tryptophan at position 48, which is characteristic of the typical homeodomain. The other homeodomain, MlSIX45, is also a member of the SINE class.
The Mnemiopsis SINE class homeodomains do not clearly separate into the three families recognized in bilaterians. Only two of the 18 maintain the four family-defining diagnostic residues (positions 3-6) in the homeodomain . MlSIX41 and MlSIX27 have the SIX1/2 family 'ETSY' pattern in positions 3-6 of the homeodomain. However, neither MlSIX41 nor MlSIX27 group convincingly with the Six1/2 group. The Mnmemiopsis SIX class is the result of extensive ctenophore-specific diversification. A more in-depth phylogenetic analysis that includes SIX domains may provide additional insight into these relationships.
Unclassified Mnemiopsis homeodomains
Two clades consisting of 18 Mnemiopsis homeodomains appear as separate offshoots in our superfamily tree (Figure 1, Table 1). None of these 18 homeodomains have introns, or any of the known class signatures, that would hint that they might belong to an existing class. MlHD60 and MlHD79 have insertions that interrupt the homeodomain but these insertions are unlike the known insertions seen in the TALE, HNF and PROS classes. The MlHD60 insertion consists of two amino acids that occur in the third alpha-helix. The other insertion occurs in the loop region between the first and second alpha-helices but, unlike the TALE insertions, it consists only of a single amino acid. The average branch length of the homeodomains in these clades is 5% shorter than for the other Mnemiopsis homeodomains, confirming that these unclassified Mnemiopsis homeodomains do not simply comprise a clade of unusually long branches.
There are no Mnemiopsis homeodomains that grouped with HNF, CUT, PROS, or CERS classes in our analyses. Consistent with this result, no Mnemiopsis homeodomains exhibit insertions between the second and third helices, like those seen in the bilaterian HNF and PROS class homeoboxes. Besides the five apparent SINE class homeodomains, no other Mnemiopsis homeodomains group with zinc finger (ZF) homeodomains.
Evolutionary dynamics of the Mnemiopsis homeodomain superfamily
Paralog count and estimated branch lengths of seven species.
No. of HDs in this test
No. of HDs in species-specific clades
No. of species-specific clades
Average branch length
For each species, we recorded the number of species-specific clades that included more than one homeodomain, as well as the total number of homeodomains in those species-specific clades (Table 2). These numbers give us an approximation of the number of lineage-specific homeodomain duplications that have been preserved in a specific lineage since it split from its closest relative in the analysis . A species that has recently undergone extensive genome reduction would be expected to harbour less species-specific clades than a genome that has experienced a recent genomic expansion. Our data shows that very few paralogous homeodomains exist in the Amphimedon (7) and Trichoplax (0) genomes, whereas the human genome has a remarkably high level of paralogs (197). Mnemiopsis (45) and Nematostella (66) are both very close to the mean (Table 2).
Branch lengths provide a means of measuring the level of divergence for a particular homeodomain. Longer branches correspond to higher levels of divergence. We rooted the same neighbor-joining tree described above at its midpoint and determined the average branch lengths for each species' set of homeodomains (Table 2). In our tree, the Mnemiopsis branches tend to be longer than for all the other species except for C. elegans, which is known to have very long branches . The Mnemiopsis average branch length is slightly closer to the mean than it is to the C. elegans average, suggesting that the Mnemiopsis homeodomains are moderately divergent. The trees used in this analysis are included as supplemental material (Additional Files 2 and 3) and branch lengths can be visualized directly using a tree-viewing program such as Figtree .
The apparent absence of Hox/ParaHox, S50 and K50 PRD and HNF class homeodomains in Amphimedon and Mnemiopiopsis supports an organismal clade that unites cnidarians, bilaterians and placozoans. Although previous analyses have strongly suggested a monophyletic relationship of these three phyla (particularly Collins ), we are not aware of any studies to date that have given this subkingdom a name. We propose the name ParaHoxozoa for the clade containing Bilateria, Placozoa and Cnidaria.
Studies have shown that the ancestor of the ParaHoxozoa had at least one Hox/ParaHox gene [26, 54] and that these genes are missing from both ctenophores and sponges ( and this study). While cnidarians and bilaterians have homeodomains that are clearly descended from both ancestral Hox and ParaHox genes, Trichoplax, has only a single Hox/ParaHox gene, Trox2 . This homeodomain consistently forms a moderately supported clade with the ParaHox GSX family ([28, 56] and this study) and may be a true ParaHox gene. However, it has been postulated that the Trox2 gene may be a direct descendant of an ancestral 'ProtoHox' gene rather than a proper ParaHox gene [28, 54, 56]. The name ParaHoxozoa was chosen since, based strictly on the trees themselves, a ParaHox clearly unites this group.
Evidence distinguishing Mnemiopsis from the ParaHoxozoa
The presence of ample paralogs in the Mnemiopsis genome suggests that it has not undergone extensive genome reduction. Currently, we cannot rule out the possibility that recent duplication events have masked more ancient gene losses or that the paralog count in the homeobox superfamily is not typical of most superfamilies. Regardless, the level of paralog retention does give some initial insight into the evolutionary dynamics of the Mnemiopsis genome, while also lowering the probability that the Mnemiopsis lineage at one time had (and subsequently lost) these missing homeodomains. The availability of other ctenophore genomes for future studies and analysis of other gene superfamilies will help to better-resolve the evolutionary dynamics of the Mnemiopsis genome.
The unusually high level of conservation of the homeobox genes and their vital roles in early development further lowers the possibility that the loss of entire classes and subclasses of these genes could be tolerated. There are numerous examples of homeobox families missing from a wide array of phylogenetically disparate lineages. For example, Ciona intestinalis and multiple species of parasitic platyhelminthes appear to have lost several Hox genes [57, 58]. In addition, the genomes of Drosophila and human are each missing several homeobox families . However, it is important to note that no species examined thus far has been shown to be missing entire classes or subclasses of homeobox genes.
Long branches do have the ability to distort phylogenies and the presence of longer branches in our Mnemiopsis data set was initially a concern. However, if ctenophores are among the earliest metazoan branches and substantial expansion of the homeodomain superfamily occurred along the lineage leading to Mnemiopsis, these long branches are not entirely unexpected. Conversely, the somewhat shorter branches in Amphimedon homeodomains, despite its early-branching phylogenetic position, are likely due to it having far fewer paralogs.
Another possible concern is that our paralog results are an artifact caused by the attraction of long branches and that these long branches represent class members that appear to be missing but are perhaps unrecognizable in Mnemiopsis. Two pieces of evidence undermine this assertion. First, C. elegans, which has 9% longer branches than Mnemiopsis and 67% more homeodomains, would be expected to have more species-specific homeodomain clades, but in fact has 23% fewer (Table 2). Second, when we remove the C. elegans homeodomains from the dataset (which vastly reduces the number of long branches in our tree) and redo the analysis, the number of Mnemiopsis -specific clades does not change (Additional File 3).
Incongruence with the 'Coelenterate' hypothesis
A scenario grouping cnidarians and ctenophores into a 'coelenterate' clade consistent with the phylogenomic study by Philippe et al.  seems unlikely based on the number of clades that include a Nematostella homeodomain to the exclusion of a Mnemiopsis representative. For example, Nematostella has representatives of the ANTP homeodomain families Mox, Emx, Gbx, Ro, Mnx, Vax, Not, Nk1, NK2 and NK3, as well as two ParaHox and three Hox-related homeoboxes (Figure 2). If Ctenophora and Cnidaria were, in fact, sister taxa, an extraordinary number of gene losses would have been required in the Mnemiopsis lineage, given the observation that Mnemiopsis is lacking these families. This amount of loss seems very unlikely given the pattern of paralog retention in Mnemiopsis. A more recent study by this group also failed to recover the 'coelenterate' clade .
Mnemiopsis extended Hox?
The origin of the four paralogous Mnemiopsis ANTP homeodomains (MlANTP03a, MlANTP03b, MlANTP03c and MlANTP03d) that group with the engrailed family in the superfamily tree (Figure 1) and with the Evx family in the ANTP tree is difficult to interpret. The presence of the HOXL signature in two of these four homeodomains, combined with their tendency to form a clade with hox-related genes, suggests that these four homeodomains may have descended from an ancestral homeodomain; this ancestral homeodomain may then, in turn, have led to the formation of the extended Hox (Evx, Gbx, Meox, Mnx, Rough), Hox and ParaHox classes through a series of duplications in the ParaHoxozoa stem.
Functional analysis of the MlANTP03 genes may give further insight into the role of this ancestral gene, which would be particularly interesting given the roles that extended Hox genes play in critical biological processes such as neurogenesis, myogenesis, axial patterning, segmentation, gastrulation and photoreception.
Identifying the basal branch
Despite strong evidence uniting the ParaHoxozoa to the exclusion of Porifera and Ctenophora, it does not appear that homeodomain sequences alone (at least based on the current set of available sequences) will be powerful enough to incontrovertibly determine the basal metazoan branch. There are several unambiguous clades containing Mnemiopsis homeodomains grouping with parahoxozoan clades to the exclusion of Amphimedon representatives (for example, Lmx, Hlx and Pbx). The many occurrences of this pattern insinuate that sponges branched off at the base of the Metazoa. This pattern is consistent with many of the 18 S phylogenies (see Table 1 in  and references therein), as well as a recent EST-based phylogeny . However, there is substantial evidence suggesting the alternative relationship.
A number of parahoxozoan clades include Amphimedon homeodomains but lack a corresponding Mnemiopsis homeodomain (for example, NK2 and Irx). The scarcity of paralogous homeodomains could be a sign that the Amphimedon genome is reduced. This, in turn, increases the chance that multiple homeodomain families were lost in Amphimedon. In addition, the lack of a homeobox associated with a paired domain in Mnemiopsis opens the possibility that the fusion of the homeobox and paired domain postdates the divergence of ctenophores from the rest of Metazoa. This would be consistent with other phylogenomic studies [11, 12].
It is possible that similar analyses of additional gene families with the same set of taxa might be able to resolve the basal branch. Alternatively, it could be that sequence data from additional sponge and ctenophore genomes will be required to satisfactorily settle the debate. The current study does make a strong case that this uncertainty might eventually be resolved using a gene family-type approach.
Implications for the evolution of homeobox clusters
Extensive genomic clustering of ANTP homeoboxes in multiple metazoan genomes suggests that the ANTP class homeobox genes were formed by a series of tandem duplications (reviewed in ). A recent study showed that six of the eight Amphimedon NKL homeoboxes are clustered  and are likely descendants of the ancestral ANTP megacluster. One of the four Mnemiopsis homeobox clusters - specifically, the cluster containing MlANTP19 and MlANTP47 - is potentially a remnant of this ANTP megacluster as well (Figure 5). In our phylogenetic analyses, MlANTP47 consistently groups with the Amphimedon NK5.7b homeodomain with moderate support (ML bootstrap = 24, Bayesian posterior probability distribution = 95; Figure 2). The family-level identity of MlANTP19 is uncertain but it does not appear to be paralogous with MlANTP47.
It is very difficult to draw any conclusions as to the implications of the other three linked homeoboxes, particularly given the difficulty of assigning family-level orthology to the Mnemiopsis homeodomains. The process for detecting ancestral/functional linkages involves identifying orthologs in multiple evolutionarily disparate genomes. Therefore, it may be possible that these linkages are significant and, perhaps, representative of some ancestral cluster, but the current phylogenetic resolution of the Mnemiopsis homeoboxes may make this relationship difficult to detect.
Distribution of classes and subclasses of homeobox genes among early branching taxa.
The expansion of the homeobox superfamily has played a major role in the evolution of animal phyla . An understanding of this expansion, along with an accurate animal phylogeny, is critical to understanding metazoan evolution. With this new dataset and phylogeny, as well as the help of functional genomic techniques, we can start piecing together the evolutionary steps that led to such astounding evolutionary feats such as the development of nervous systems, muscular systems, and complex symmetry.
Sequencing and assembly
Genomic DNA was isolated from the larvae of two separate self-fertilizing hermaphroditic individuals. A library from one source was sequenced to 10× coverage using 454 sequencing: 8.1 million reads totaling 2.7 gigabases were assembled into 29,877 contigs (contig-N50 = 11KB) using the Phusion assembler . Subsequently, we constructed and sequenced a paired-end DNA library with insert sizes of ~4kb from the second genomic source using Illumina sequencing. These 2.8 million paired end reads were used to compile the contigs into 10,106 scaffolds (scaffold-N50 = 123 KB) bringing the physical coverage to ~50×.
Retrieval of Mnemiopsis homeodomains
A TBLASTN search of the Mnemiopsis assembly was conducted using a set of bilaterian homeodomains downloaded from the Homeodomain Resource . This set was compared to and supplemented by previous Mnemiopsis homeodomains that were generated from degenerate PCR and RACE .
Gene isolation via RACE PCR
MlPRD86, MlPbx, MlSIX59b, MlSIX41, MlSIX13f, MlPRD10b, MlANTP65, MlANTP67, MlPRD76b, MlPRD60, MlPRD77b and MlPRD86 homeoboxes had been isolated but not published prior to the sequencing of the genome and are included in this study. These were isolated as previously described .
The Mnemiopsis homeodomains were aligned by eye to the dataset used by Holland et al. 2007, which consisted of all human homeodomains and a representative set of protostome set consisting mostly of homeodomains from Drosophila melanogaster . We supplemented these sequences with eight Branchiostoma floridae homeodomains that are known to be missing from humans. Insertions in the loop-region of the Mnemiopsis homeodomains were removed as done in  and other studies. In determining which amino acids to remove from atypical homeodomains, we realized that there was inconsistency in the Holland set as to which three amino acids were removed between the Drosophila and the human sequences. As with the human sequences, we removed the 23rd, 24th and 25th amino acids from the Mnemiopsis atypical TALE sequences and adjusted the Drosophila sequences so they conformed to this rule. Alignment is available as supplementary material (Additional File 1).
The superclass alignment was divided into six separate alignments, (1) ANTP, (2) PRD, (3) SIX, (4) LIM, (5) POU and (6) TALE based on the best superclass tree and the class membership of the bilaterian homeodomains determined from HomeoDB . The three amino acid insertions were reinserted into the TALE alignment. To each of these datasets we added homeodomains from other non-bilaterian species based on published classifications. We added 32 homeodomains from the demosponge Amphimedon queenslandica , 37 homeodomains from the placozoan Trichoplax adhaerens , 127 homedomains from the cnidarian Nematostella vectensis  and two homeodomains from the choanoflagellate Monosiga brevis . These alignments are available as supplemental material (Additional File 1).
The Perl script proteinModelSelection.pl (available from the RaxML  web site) was used to determine the best scoring amino acid substitution model for our supertree alignment (RTREV +GAMMA). All subsequent analyses used this model.
Three independent runs of RaxML version 7.0.4  were conducted. Two runs used random starting trees with the following command line (raxmlHPC-MPI -m PROTGAMMARTREV -s ALN.phy -#10 -n NAME -k). One run used a neighbor joining starting tree that was generated with default parameters in Phylip version 3.6a3 . The command line for this run was (raxmlHPC -m PROTGAMMARTREV -s ALN.phy -t NJ.tre -n NAME -k).
Three independent runs of PhyML version 3.0  were conducted. Two runs used random starting with the following command line (phyml -i ALN.phy -d aa -m RtRev -a e -q --rand_start -s SPR --r_seed 'cat FILE_W_RANDSEED'). One run used a neighbor joining starting tree with the following command line (beorun phyml -i ALN.phy -d aa -m RtRev -a e -q -s SPR).
Two independent runs with the MPI version of Mr. Bayes version 3.1.2  were conducted with the following execution block (prset aamodelpr = fixed(RTREV); lset rates = gamma; mcmcp mcmcdiagn = no nruns = 1 ngen = 5000000 printfreq = 5000 samplefreq = 500 nchains = 5 savebrlens = yes; mcmc; sumt filename = tale_w_insert.nex nRuns = 1 Relburnin = YES BurninFrac =.25 Contype = Allcompat;). Log likelihood values were plotted and their progression was visually examined over time. All runs were found to be asymptotic before the.25 burnin fraction.
Likelihood values for all runs (3 RaxML, 3 PhyML, 2 Bayes and 1 neighbor joining) were generated using PhyMl version 3.0  with the following command line (phyml -i ALN.phy -c 4 -m RtREV -a e -o lr -f d -u TREE.tre -d aa -b 0 -s NNI). The tree with the highest likelihood value was chosen for all downstream analysis. The best trees were as follows: Supertree = RaxML-randomstart, ANTP = RaxML-randomstart, PRD = RaxML-randomstart, SIX = PhyML-randomstart, LIM = RaxML-randomstart, POU = RaxML-NJstart, and TALE = RaxML-randomstart. The best trees are available as supplemental material (Additional File 4).
Support was assessed by 100 replicates of the bootstrap using the method that provided the highest likelihood value. Bootstrap values greater than 50 were applied to the best tree. To these three trees, Bayesian posterior probabilities were also added from the Bayesian tree with the highest likelihood score. The best trees with support values are available as supplemental material (Additional File 4).
Paralog retention analysis
From our initial superfamily alignment, we removed all sequences that were not Mnemiopsis, Drosophila, or human. To this, we added all Trichoplax, Nematostella and Amphimedon sequences that were used in our class-level phylogenies. Finally we added all the C. elegans homeodomains that were used in . We ran a neighbor-joining analysis using default parameters in Phylip . This tree was then parsed with a Perl script (count_species_specific_clades.pl) that identified occurrences of homogeneous clades of taxa with identical two-letter prefixes. This Perl script, tree and alignment are included as supplemental material (Additional File 2). A version of this analysis was run without the C. elegans data (see Additional File 3).
Branch length analysis
The tree used in the paralog analysis was subsequently used to estimate average branch lengths (Table 2). This tree was opened in FigTree v.1.2.3 , rooted at the midpoint, and saved. This rooted tree was then opened in TreeStat v.1.2  and the Root-Tip Lengths were calculated. A Perl script (calculate_average_root_tip.pl) was used to parse the output of TreeStat and calculate the average branch lengths. This Perl script, tree, and alignment are included as supplemental material (Additional File 2). In order to compare lengths of subsets of Mnemiopsis data, we ran a separate instance of this analysis, substituting the two-letter 'Ml' prefix of the subset to be analysed with 'Zz'. In order to assess the feasibility of using this technique to assess the relative branch lengths from trees based on 60 amino acid matrices, we performed a simulation study, and the details of this simulation study are available as supplemental material (Additional File 5).
This work was supported by the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health. Additional support was provided by the NSF Graduate Research Fellowship Program to KP and grants from NASA and NSF to MQM. We would like to thank Alice Young, Brian Schmidt, Natalie Gurson, Richelle Legaspi and Betsy Novotny, who were largely responsible for the sequencing work performed at NISC. JR would like to thank Arjun Prasad, Adam Reitzel, Bhavesh Borate and Tyra Wolfsberg for insightful conversations regarding this work. We are also thankful to the three anonymous reviewers whose thoughtful suggestions, corrections and insightful comments greatly improved this finished manuscript.
- Pang K, Martindale MQ: Ctenophores. Curr Biol 2008, 18: R1119–1120.PubMedView Article
- Cuvier G: Le Règne Animal Distribué Selon son Organisation, pour Servir de Base à l'Histoire Naturelle des Animaux et d'Introduction à l'Anatomie Comparée. Paris: Deterville; 1817.
- Leuckart R: Ueber die Morphologie und die Verwandtschaftsverhältnisse der wirbellosen thiere. Ein Beitrag zur Charakteristik und Classification der thierischen Formen. Braunschweig,: F. Vieweg und Sohn;; 1848.
- Harbison GR: On the classification and evolution of the Ctenophora. In The Origins and Relationships of Lower Iinvertebrates. Volume 28. Edited by: Conway Morris S, George JD, Gibson R, Platt HM. London, UK: The Systematics Assocation; 1985:78–100.
- Bridge D, Cunningham CW, DeSalle R, Buss LW: Class-level relationships in the phylum Cnidaria: molecular and morphological evidence. Mol Biol Evol 1995, 12: 679–689.PubMed
- Collins AG: Evaluating multiple alternative hypotheses for the origin of Bilateria: an analysis of 18 S rRNA molecular evidence. Proc Natl Acad Sci USA 1998, 95: 15458–15463.PubMedView Article
- Kim J, Kim W, Cunningham CW: A new perspective on lower metazoan relationships from 18 S rDNA sequences . Molec Biol Evolution 1999, 16: 423–427.
- Smothers JF, von Dohlen CD, Smith LH, Spall RD: Molecular evidence that the myxozoan protists are metazoans. Science 1994, 265: 1719–1721.PubMedView Article
- Wainright PO, Hinkle G, Sogin ML, Stickel SK: Monophyletic origins of the metazoa: an evolutionary link with fungi. Science 1993, 260: 340–342.PubMedView Article
- Philippe H, Derelle R, Lopez P, Pick K, Borchiellini C, Boury-Esnault N, Vacelet J, Renard E, Houliston E, Queinnec E, et al.: Phylogenomics revives traditional views on deep animal relationships. Curr Biol 2009, 19: 706–712.PubMedView Article
- Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, et al.: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature 2008, 452: 745–749.PubMedView Article
- Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Baguna J, Bailly X, Jondelius U, et al.: Assessing the root of bilaterian animals with scalable phylogenomic methods. Proc Biol Sci 2009, 276: 4261–4270.PubMedView Article
- Schierwater B, Kolokotronis SO, Eitel M, Desalle R: The Diploblast-Bilateria Sister hypothesis: Parallel revolution of a nervous systems may have been a simple step. Commun Integr Biol 2009, 2: 403–405.PubMedView Article
- Delsuc F, Brinkmann H, Philippe H: Phylogenomics and the reconstruction of the tree of life. Nat Rev Genet 2005, 6: 361–375.PubMedView Article
- Pick KS, Philippe H, Schreiber F, Erpenbeck D, Jackson DJ, Wrede P, Wiens M, Alie A, Morgenstern B, Manuel M, Worheide G: Improved phylogenomic taxon sampling noticeably affects non-bilaterian relationships. Mol Biol Evol 2010, 27: (9): 1983–1987.PubMedView Article
- Rokas A, Holland PW: Rare genomic changes as a tool for phylogenetics. Trends Ecol Evol 2000, 15: 454–459.PubMedView Article
- Balavoine G: Are Platyhelminthes Coelomates without a Coelom? An argument based on the evolution of hox genes. Amer Zool 1998, 38: 843–858.
- Cook CE, Smith ML, Telford MJ, Bastianello A, Akam M: Hox genes and the phylogeny of the arthropods. Curr Biol 2001, 11: 759–763.PubMedView Article
- 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–776.PubMedView Article
- McGinnis W, Levine MS, Hafen E, Kuroiwa A, Gehring WJ: A conserved DNA sequence in homoeotic genes of the Drosophila Antennapedia and bithorax complexes. Nature 1984, 308: 428–433.PubMedView Article
- Bürglin TR: Homeodomain Proteins. In Encyclopedia or Molecular Cell Biology and Molecular Medicine. 2nd edition. Edited by: Meyers RA. Weinheim: Wiley-VCH Verlag; 2005.
- Derelle R, Lopez P, Le Guyader H, Manuel M: Homeodomain proteins belong to the ancestral molecular toolkit of eukaryotes. Evol Dev 2007, 9: 212–219.PubMedView Article
- Holland PW, Booth HA, Bruford EA: Classification and nomenclature of all human homeobox genes. BMC Biol 2007, 5: 47.PubMedView Article
- Zhong YF, Butts T, Holland PW: HomeoDB: a database of homeobox gene diversity. Evol Dev 2008, 10: 516–518.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, Kuo A, Mitros T, Salamov A, Carpenter ML, et al.: The Trichoplax genome and the nature of placozoans. Nature 2008, 454: 955–960.PubMedView Article
- Schierwater B, Kamm K, Srivastava M, Rokhsar D, Rosengarten RD, Dellaporta SL: The early ANTP gene repertoire: insights from the placozoan genome. PLoS One 2008, 3: e2457.PubMedView Article
- Degnan BM, Vervoort M, Larroux C, Richards GS: Early evolution of metazoan transcription factors. Curr Opin Genet Dev 2009, 19: 591–599.PubMedView Article
- Larroux C, Luke GN, Koopman P, Rokhsar DS, Shimeld SM, Degnan BM: Genesis and expansion of metazoan transcription factor gene classes. Mol Biol Evol 2008, 25: 980–996.PubMedView Article
- Bebenek IG, Gates RD, Morris J, Hartenstein V, Jacobs DK: Sine oculis in basal Metazoa. Devel Genes Evolution 2004, 214: 342–351.
- Finnerty JR: Homeobox genes in the Ctenophora: identification of paired-type and Hox homologues in the atentaculate ctenophore, Beroe ovata . Molec Marine Biol Biotech 1996, 5: 249–258.
- Hoshiyama D, Iwabe N, Miyata T: Evolution of the gene families forming the Pax/Six regulatory network: isolation of genes from primitive animals and molecular phylogenetic analyses. FEBS Lett 2007, 581: 1639–1643.PubMedView Article
- Lee SE, Gates RD, Jacobs DK: Gene fishing: the use of a simple protocol to isolate multiple homeodomain classes from diverse invertebrate taxa. J Mol Evol 2003, 56: 509–516.PubMed
- Martinelli C, Spring J: T-box and homeobox genes from the ctenophore Pleurobrachia pileus: comparison of Brachyury, Tbx2/3 and Tlx in basal metazoans and bilaterians. FEBS Lett 2005, 579: 5024–5028.PubMedView Article
- Pang K, Martindale MQ: Developmental expression of homeobox genes in the ctenophore Mnemiopsis leidyi . Dev Genes Evol 2008, 218: 307–319.PubMedView Article
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22: 2688–2690.PubMedView Article
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003, 52: 696–704.PubMedView Article
- King N, Westbrook MJ, Young SL, Kuo A, Abedin M, Chapman J, Fairclough S, Hellsten U, Isogai Y, Letunic I, et al.: The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans. Nature 2008, 451: 783–788.PubMedView Article
- Fonseca NA, Vieira CP, Holland PW, Vieira J: Protein evolution of ANTP and PRD homeobox genes. BMC Evol Biol 2008, 8: 200.PubMedView Article
- Castro LF, Holland PW: Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes. Evol Dev 2003, 5: 459–465.PubMedView Article
- Pollard SL, Holland PW: Evidence for 14 homeobox gene clusters in human genome ancestry. Curr Biol 2000, 10: 1059–1062.PubMedView Article
- Banerjee-Basu S, Baxevanis AD: Molecular evolution of the homeodomain family of transcription factors. Nucleic Acids Res 2001, 29: 3258–3269.PubMedView Article
- Galliot B, de Vargas C, Miller D: Evolution of homeobox genes: Q50 Paired-like genes founded the Paired class. Dev Genes Evol 1999, 209: 186–197.PubMedView Article
- Larroux C, Fahey B, Liubicich D, Hinman VF, Gauthier M, Gongora M, Green K, Worheide G, Leys SP, Degnan BM: Developmental expression of transcription factor genes in a demosponge: insights into the origin of metazoan multicellularity. Evol Dev 2006, 8: 150–173.PubMedView Article
- Groger H, Callaerts P, Gehring WJ, Schmid V: Characterization and expression analysis of an ancestor-type Pax gene in the hydrozoan jellyfish Podocoryne carnea . Mech Dev 2000, 94: 157–169.PubMedView Article
- Ryan AK, Rosenfeld MG: POU domain family values: flexibility, partnerships and developmental codes. Genes Dev 1997, 11: 1207–1225.PubMedView Article
- Bürglin TR: Analysis of TALE superclass homeobox genes (MEIS, PBC, KNOX, Iroquois, TGIF) reveals a novel domain conserved between plants and animals. Nucleic Acids Res 1997, 25: 4173–4180.PubMedView Article
- Mukherjee K, Bürglin TR: Comprehensive analysis of animal TALE homeobox genes: new conserved motifs and cases of accelerated evolution. J Mol Evol 2007, 65: 137–153.PubMedView Article
- Seo HC, Curtiss J, Mlodzik M, Fjose A: Six class homeobox genes in drosophila belong to three distinct families and are involved in head development. Mech Dev 1999, 83: 127–139.PubMedView Article
- Fitch WM: Distinguishing homologous from analogous proteins. Syst Zool 1970, 19: 99–113.PubMedView Article
- Fitch DHA: Introduction to nematode evolution and ecology (August 31 2005). [http://www.wormbook.org/] The C. elegans Reseach Community, WormBook;
- Rambaut A, ed: FigTree. In Book FigTree. Edinburgh: Institute of Evolutionary Biology, University of Edinburgh; 2006.
- Jakob W, Sagasser S, Dellaporta S, Holland P, Kuhn K, Schierwater B: The Trox-2 Hox/ParaHox gene of Trichoplax (Placozoa) marks an epithelial boundary. Dev Genes Evol 2004, 214: 170–175.PubMedView Article
- Larroux C, Fahey B, Degnan SM, Adamski M, Rokhsar DS, Degnan BM: The NK homeobox gene cluster predates the origin of Hox genes. Curr Biol 2007, 17: 706–710.PubMedView Article
- Monteiro AS, Schierwater B, Dellaporta SL, Holland PW: A low diversity of ANTP class homeobox genes in Placozoa. Evol Dev 2006, 8: 174–182.PubMedView Article
- Wada S, Tokuoka M, Shoguchi E, Kobayashi K, Di Gregorio A, Spagnuolo A, Branno M, Kohara Y, Rokhsar D, Levine M, et al.: A genomewide survey of developmentally relevant genes in Ciona intestinalis . II. Genes for homeobox transcription factors. Dev Genes Evol 2003, 213: 222–234.PubMedView Article
- Koziol U, Lalanne AI, Castillo E: Hox genes in the parasitic platyhelminthes Mesocestoides corti , Echinococcus multilocularis and Schistosoma mansoni : evidence for a reduced Hox complement. Biochem Genet 2009, 47: 100–116.PubMedView Article
- Ryan JF, Baxevanis AD: Hox, Wnt and the evolution of the primary body axis: insights from the early-divergent phyla. Biol Direct 2007, 2: 37.PubMedView Article
- Garcia-Fernandez J: The genesis and evolution of homeobox gene clusters. Nat Rev Genet 2005, 6: 881–892.PubMedView Article
- Valentine JW: Why no new phyla after the cambrian? Genome and ecospace hypotheses revisited. PALAIOS 1995, 10: 190–194.View Article
- Mullikin JC, Ning Z: The phusion assembler. Genome Res 2003, 13: 81–90.PubMedView Article
- Moreland RT, Ryan JF, Pan C, Baxevanis AD: The homeodomain resource: a comprehensive collection of sequence, structure, interaction, genomic and functional information on the homeodomain protein family. Database (Oxford) 2009, 2009: bap004.
- Felsenstein J: PHYLIP (Phylogeny Inference Package). Version 3.63. 2004. [Distributed by the author: Department of Genetics, University of Washington, Seattle, USA]
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19: 1572–1574.PubMedView Article
- Nam J, Nei M: Evolutionary change of the numbers of homeobox genes in bilateral animals. Mol Biol Evol 2005, 22: 2386–2394.PubMedView Article
- Rambaut A, Drummond AJ, eds: [http://tree.bio.ed.ac.uk/software/treestat/] TreeStat. v1.1:Tree Statstic Calculaton Tool.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.