BACKGROUND
genes that belong to the homeobox family are characterized by the ability to code for a protein that contains a recognizable, although very variable, 'homeodomain', usually  <dig> amino acids in length  <cit> .
many of these genes are transcription factors that play important roles in the embryonic development of bilaterian and non-bilaterian animals.
changes in homeobox gene content and deployment during evolution may have contributed to the evolution of body plan differences in animals  <cit> .
therefore, comparison of homeobox gene sets from different animals may shed light on the evolutionary events that gave rise to animal body plan diversity.
nevertheless, gene orthology is not always easy to establish when comparing divergent animals  <cit> .
for this purpose, phylogenetic analysis, conservation of synteny, paralogy within the human genome, insertions within the homeodomain, key amino acid residues, and several motifs outside of the homeodomain can all be used.
these features can also be used to classify homeobox genes into classes, subclasses and families.
in the latest revision, holland et al.
<cit>  classified all  <dig> human homeobox genes into  <dig> classes  and  <dig> gene families.
the antp class is further divided into hoxl and nkl subclasses.
it should be noted that the only protein region that can be aligned in all  <dig> gene families or  <dig> genes is the homeodomain  <cit> .
here, we report amino acid patterns typical of bilaterian hoxl, nkl, and prd genes that can be used to quickly and efficiently retrieve amino acid sequences belonging to these classes and subclasses, among hundreds of other homeodomain sequences.
retrieving a given class or subclass of sequences from many animal genomes may be thus an easier task than previously thought.
however, we show that these typical amino acid patterns should be cautiously used when sequences come from non-bilaterian animals.
since phylogenetic analysis was one of the primary sources of evidence used to establish the different classes and sub-classes, it is likely that most of homeobox gene classes and subclasses represent monophyletic lineages.
hence, it is conceivable that amino acid changes important for protein function in a given lineage may be revealed as fixed differences between classes and subclasses.
previous attempts have been made to classify genes within hox families  <cit>  but not at the level of whole classes or subclasses.
here, we show that, because of their chemical properties, amino acid usage is different in antp and prd classes at five positions.
furthermore, at nine positions, amino acid usage is different between hoxl and nkl subclasses.
our findings support the notion that many chemically important changes happened early in the evolution of homeobox genes, and that these changes can be used as additional evidence to establish gene orthology.
these results can be helpful for experimental studies aimed onto investigating the biochemical functions of key homeodomain residues in different gene classes as well.
methods
data
we have used the hand-curated data set of holland et al.
<cit> , the prosite homeodomain data set  <cit>  and the ncbi database  <cit> .
identification of amino acid patterns typical of hoxl, nkl and prd genes
in order to find amino acid patterns that distinguish genes from the antp subclasses, as well as genes from the prd class, from genes belonging to other classes, a fast word discovery program  was first used to find statistically interesting words.
the minimum word length used was two  and the minimum number of sequences where the word should occur  was set to 10% of the number of sequences in the data set being considered.
the words found were then filtered to discard those words that occurred more than three times in more than one dataset.
the software bioredx  was then used to find the largest and more general amino acid patterns that distinguish two sets of sequences.
this program accepts as input two sequence files and a word seed.
based on the seed given, tries to find patterns that occur in one file  and not in the other .
this software was used to find patterns that occurred less than four times in the negative file and considered patterns up to  <dig> amino acid residues larger than the initial word seed.
in the reported amino acid patterns, amino acids listed within brackets are those allowed at a given position.
for compactness of representation it is also possible to negate the class.
the negation is denoted by "^".
in this case, the amino acid residues listed are the ones that are not allowed at that particular position.
the approach here used to find amino acid patterns does not use a set of aligned sequences.
therefore, in order to make sure that the derived amino acid patterns are found in the same region of the homeodomain sequence, all sequences belonging to a given class and sub-classes were aligned using the pam  <dig> scoring matrix  <cit> .
identification of key amino acid residues
there are many schemes for comparing and grouping amino acids.
nevertheless, none of them can possibly capture the vast number of contexts in which amino acids are found within proteins.
the scheme here used is that of livingstone and barton  <cit> , based on the amino acid properties size, polarity, hydrophobicity, charge, aliphaticity and aromaticity.
for each position and property we calculated the probability of gaining or losing a constraint using the distribution of the amino acid property's observed at that position in all sequences considered.
only changes in the properties values with a probability of occurrence lower than 5% were considered.
phylogenetic analyses
a neighbour-joining tree, using pair-wise deletion, as implemented in the mega software  <cit>  was constructed in order to classify the set of  <dig> homeodomain sequences from a variety of animal species , into the  <dig> classes and two sub-classes scheme proposed by holland et al.
<cit> .
RESULTS
in order to establish a method for fast retrieval of sequences belonging to a given class or subclass, amino acid patterns characteristic of different homeobox classes and subclasses were sought.
as a starting point we used the hand-curated human prd, nkl and hoxl data sets compiled by holland et al.
<cit> .
the homeodomains of human en <dig>  en <dig>  dlx <dig>  dlx <dig>  dlx <dig>  dlx <dig>  dlx <dig>  dlx <dig>  noto, and hopx are difficult to classify and were not used in these initial analyses.
therefore, the prd, nkl and hoxl data sets contain  <dig>   <dig> and  <dig> sequences, respectively.
within each of the three data sets, the amino acid sequences show considerable variation, suggesting that the human sequences include the majority of amino acid variability allowed at a given position along the sequence.
this assumption was later tested .
1) amino acid patterns characteristic of the different human gene classes and sub-classes
many characteristic amino acid patterns were detected using the approach described in material and methods.
here, we describe the seven patterns that are most pertinent, taking into account their coverage across a group of genes and absence in the other groups.
the pattern  wfqnrr k    , named hoxl  <dig> pattern, is present in all human hoxl homeodomains, and in none of the human nkl and prd genes.
the pattern le e , named hoxl  <dig> pattern, is present in all human hoxl homeodomains except mnx <dig>  only two human prd genes  encode this amino acid pattern.
none of the human nkl genes encode this amino acid pattern.
the pattern      , named nkl pattern, is present in  <dig> out of  <dig> human nkl sequences, in one human hoxl sequence , and in none of the human prd sequences.
the human nkl sequences that do not present this pattern are: hhex, nanog, and ventx.
the pattern l         p , named prd pattern, is present in all human prd homeodomains and never present in human hoxl and nkl sequences.
the patterns    , named antp  <dig> pattern, and  w qn r k, named antp  <dig> pattern,  are present in all human hoxl and nkl genes and never present in human prd homeodomains.
the pattern w      , named antp-prd pattern, is found in all human prd, nkl and hoxl homeodomains.
2) generality of the amino acid patterns found
in order to test the generality of the amino acid patterns derived in the previous section, we used the  <dig> homeodomain sequences classified by holland et al.
<cit> , which include the  <dig> human sequences used above.
these sequences include genes classified as hoxl, nkl and prd, plus other homeobox gene classes .
it should be noted, that in these analyses, the difficult to classify genes that were excluded above were now included, and these were classified as tentatively suggested by holland et al.
<cit> .
using amino acid sequences, a neighbour-joining tree  was built with the  <dig> sequences, plus  <dig> homeodomain sequences from a variety of animal species.
the latter  <dig> sequences were obtained from the file ps <dig> available at prosite  <cit>  after removing all non-animal sequences from the file.
sequences from clusters that are supported by a bootstrap value of 80% or higher, and that include at least one sequence that has been classified by holland et al.
<cit> , as belonging to a given class or subclass, were classified as belonging to that class or subclass.
using this phylogenetic argument  <dig> sequences could be classified as hoxl,  <dig> as nkl,  <dig> as prd,  <dig> as lim,  <dig> as pou,  <dig> as hnf,  <dig> as sine,  <dig> as tale,  <dig> as cut,  <dig> as pros,  <dig> as zf, and  <dig> as cers.
nevertheless,  <dig> sequences remained unclassified.
uncl.
– unclassified sequences.
the total number of sequences in each class is shown in parentheses.
the hoxl  <dig> pattern is not as specific as the hoxl  <dig> pattern, since it is observed in about 46% of the genes belonging to the pou class and in a few nkl  and prd  genes.
all  <dig> prd genes showing the hoxl  <dig> pattern belong to the pax4/ <dig> gene family, thus it is likely that this is a case of convergent evolution.
the human pou homeodomains showing hoxl  <dig> pattern are pou1f <dig>  pou3f <dig>  pou3f <dig>  pou3f <dig>  and pou3f <dig>  the latter four genes are closely related, although the pou1f <dig> gene is distantly related to those genes  <cit> .
nevertheless, this is also likely a case of convergent evolution, since the alternative hypothesis  implies many independent losses.
therefore we argue that the necessity of using the amino acid combinations implied by hoxl  <dig> pattern is a derived feature that appeared early in the hoxl lineage.
only three sequences  classified as hoxl using a phylogenetic argument, do not show the hoxl  <dig> pattern, namely: drosophila melanogaster exex , human mnx <dig> and mouse mnx <dig> .
mnx <dig> genes are peculiar in other ways.
for instance, they show both the nkl and hoxl  <dig> pattern .
the nkl pattern is highly specific but not widely representative of all nkl genes.
the pattern is found in 57% of the sequences classified as belonging to the nkl class , but in only ten other classified genes .
in addition,  <dig> of the  <dig> phylogenetically 'unclassified' genes show the nkl pattern, and of these  <dig> have already been classified as belonging to the nkl class by prosite .
the three hoxl genes with the nkl pattern are the three mnx <dig> genes .
mnx <dig> genes also show the hoxl  <dig> pattern that is highly specific for hoxl genes.
thus, it is unlikely that mnx <dig> genes have been misclassified.
therefore, the presence of the nkl pattern in mnx <dig> sequences may be the result of convergent evolution.
the seven sine sequences where the nkl pattern is found are six <dig> from oryzias latipes, gallus gallus, mus musculus and homo sapiens and six <dig> from gallus gallus, mus musculus and homo sapiens.
all six human genes cluster together with a high bootstrap value .
nevertheless, the nkl pattern is only observed in the six <dig> and six <dig> genes, two closely related genes.
it is thus, likely a case of convergent evolution.
the requirement to use one of the amino acid combinations implied by the nkl pattern is a derived feature that appeared likely early in the nkl lineage.
most of the classified sequences where the nkl pattern is not found have been annotated by holland et al.
<cit> , and prosite as belonging to the nanog, noto, ventx, en, dlx, and barx, families.
the en, dlx and noto genes are difficult to classify.
thus, information on these genes was not used to derive the nkl amino acid pattern.
the prd pattern is highly specific, being found only in sequences classified as prd, and is also highly representative being found in 86% of the sequences classified as prd.
the prd genes that do not show this pattern are d. melanogaster odsh, otp, phdp, ptx <dig>  ip <dig>  hopx genes from danio rerio, homo sapiens, bos taurus, rattus norvegicus, mus musculus, sus scrofa, and gallus gallus , human pax <dig>  pax <dig> and pax <dig> , hydra vulgaris dmbx, otp genes from heliocidaris erythrogramma, heliocidaris tuberculata, lytechinus variegates, paracentrotus lividus, and saccoglossus kowalevskii , otx from strongylocentrotus purpuratus , alx from strongylocentrotus purpuratus , anf genes  from gallus gallus, and xenopus laevis , and the prop gene from canis familiaris .
the prd sequences that do not show this pattern are not related in any particular way.
therefore, the absence of this pattern in these sequences is likely the result of several independent losses.
therefore, it seems likely that the requirement to use one of the amino acid combinations implied by the prd pattern is a derived feature that appeared early in the evolution of prd genes.
since the prd-specific pattern is highly specific, it is likely that the  <dig> unclassified sequences that show this pattern are also prd sequences.
according to prosite classification these genes are gsc from danio rerio, xenopus laevis , gallus gallus, mus musculus, saguinus labiatus, gorilla gorilla, pongo pygmaeus, pan paniscus, and pan troglodytes, alx <dig> from mus musculus and bos taurus, and unc- <dig> and ceh- <dig> genes from caenorhabditis elegans .
according to prosite,  <dig> of these genes belong to the prd class; the status of the ceh- <dig> gene is unknown.
the antp  <dig> pattern is found in most hoxl  and nkl  gene sequences and is almost absent in prd  gene sequences.
nevertheless, such a pattern is also found in sequences from genes belonging to all other classes except the pros class .
it should, however, be noted that the sample size of pros genes is very small .
the antp sequences that do not show this pattern are drosophila melanogaster eve, lbe, and lbl, homo sapiens noto, heterodontus francisci evx <dig> , gallus gallus hme <dig> , fugu rubripes hxdbb, and caenorhabditis elegans vab <dig> and hm <dig> .
the hxb <dig> gene from gallus gallus  is incomplete, and thus it is not possible to determine whether it shows this pattern.
the hopx sequence  from danio rerio is the only prd sequence that shows antp  <dig> pattern.
although the antp  <dig> pattern is short , the broad distribution indicates that, very likely, it is not the result of convergent evolution.
it is more likely that the ability to use the amino acid combinations implied by this pattern is an ancestral feature of homeobox containing genes.
for some likely functional reason, antp genes have retained such a pattern, in contrast with prd genes  where only  <dig> % of all sequences show this pattern.
the antp  <dig> pattern is found in most hoxl  but in only 68% of the nkl sequences.
this pattern is not found in sequences from other classes.
it is thus, highly specific, although less representative across the nkl subclass of antp class genes.
the hoxl genes that do not show this pattern are  hmd <dig> from gallus gallus, hxaba from danio rerio, hxb <dig> from salmo salar and hxbb from fugu rubripes.
these sequences are not related in any particular way.
it is not possible to determine whether the hxb <dig> gene from gallus gallus  shows this pattern, since this is a partial sequence.
nkl sequences not showing the antp  <dig> pattern belong to the dlx , en  and noto gene families that have been previously found to be difficult to classify.
hence, information on these genes was not used to derive the antp  <dig> pattern.
once again, genes from these families stand out as oddities .
the antp-prd pattern identifies most hoxl , nkl  and prd  sequences used.
nevertheless, it is also present in all sine sequences and in 79% of the lim sequences.
the  <dig> antp plus prd sequences that do not show the expected pattern are the hopx genes from homo sapiens, mus musculus, rattus norvegicus, sus scrofa, bos taurus, gallus gallus, and danio rerio, the hm <dig> gene from caenorhabditis elegans, and the artemia sanfranciscana hmen gene.
it should be noted that hopx are difficult to classify, that in the latest revision  <cit>  were classified as prd genes, and these were not used to derive the amino acid patterns being tested.
the pax <dig>  pax <dig> and pax <dig> sequences from homo sapiens, as well as the hxb <dig> gene from gallus gallus, the hxdbb gene from fugu rubripes, and the dlx <dig> and dlx <dig> genes from eleutherodactylus coqui, are partial, thus it is not possible to determine whether they show the antp-prd amino acid pattern.
3) further characterization of nkl genes
given the failure to identify about 43% of the classified nkl sequences using the nkl pattern, we performed additional analyses to see whether the nkl pattern could be refined to accommodate the members of the nanog, noto, ventx, en, dlx, and barx, families.
the results are shown in table  <dig>  when the nkl and refined patterns are considered, 98% of all nkl sequences could be classified as such.
it should be noted that these refined nkl patterns are almost completely absent in the hoxl, prd, lim, pou, hnf, sine, tale, cut, pros, zf and cers classes.
the  <dig> hoxl sequences that contain the en pattern but not the nkl pattern all belong to the hox <dig> gene family.
this happens because the en pattern allows for an asparagine at position  <dig>  five out of the six genes that show the nanog pattern but not the nkl pattern belong to the hoxl gsx gene family.
these sequences are found because the pattern allows for a lysine at position  <dig>  the  <dig> hoxl sequences that show a hit when using the ventx pattern but not the nkl pattern all belong to the gbx family.
this happens because a valine is now allowed at position  <dig>  the three cut genes that show the en and ventx patterns are the homo sapiens onecut <dig> gene and the orthologues in mus musculus and rattus norvegicus.
these sequences show a hit because both patterns allow for an isoleucine at position  <dig>
uncl.
– unclassified sequences.
the total number of sequences in each class is shown in parentheses.
dlx-barx –     ; en –     ;
noto –     ; nanog –     ;
ventx –
4) non-bilaterian homeobox sequences
only a few non-bilaterian homeobox sequences  are contained in the data sets used.
therefore, we collected from the ncbi database a set of  <dig> non-redundant non-bilaterian homeobox sequences that encompass the regions where the amino acid patterns here reported are located, and that showed the antp-prd pattern derived above .
most bilaterian antp and prd sequences show this antp-prd pattern ; furthermore, this pattern is only found in antp, prd, lim and sine sequences .
therefore, by imposing the presence of the antp-prd pattern we hoped to enrich the data set for non-bilaterian antp and prd sequences.
table  <dig> summarizes the results.
as expected, of the retrieved gene sequences only one seems to belong to gene classes other than antp, prd, lim and sine.
antp1; antp2;
antp1; antp2;
antp1; antp2;
antp1; antp2;
antp1; antp2;
antp1; antp2;
only those non-bilaterian sequences that encompass the regions where the amino acid patterns here described are located, and that showed the antp-prd pattern were used.
in total  <dig> non-redundant sequences were used .
cases where non-bilaterian genes are clearly misidentified when using amino acid patterns are shown underlined .
uncl.
– unclassified sequences.
* this signature can also be considered characteristic of nkl genes since only 68% of bilaterian nkl sequences show the antp <dig> pattern
most non-bilaterian genes showing the full nkl or prd signature have already been classified as such.
furthermore, the antp demox class found in demospongiae and apparently absent in all other animals  <cit> , is characterized by the presence of the antp-prd, and antp  <dig> and  <dig> patterns, as expected for antp genes that do not belong to the hoxl or nkl lineage .
although there are no true hox genes in sponges  <cit> ,  <dig> sequences show the full hoxl signature.
interestingly, all such sequences are non-annotated, thus they must be hard to classify gene sequences.
one ephydatia fluviatilis gene has been labeled as the msx  gene.
our pattern analyses show that the corresponding sequence shows a pattern exclusively found in the mnx <dig>  gene family, thus, also suggesting that non-bilaterian sequences with full signatures may be wrongly identified.
about 53% of all non-bilaterian sequences used did not show an easily recognizable signature.
we conclude that the reported amino acid patterns must be cautiously used when applied to sequences from non-bilaterian animals.
5) identification of single key amino acid changes in bilaterian prd, nkl and hoxl genes
key amino acid changes that occurred during the evolution of homeobox containing genes may be revealed as changes affecting a single amino acid position.
thus, we compared, at each amino acid position, the large hoxl , nkl  and prd  data sets described above, for the following chemical properties: size, charge, polarity, hydrophobicity, aromaticity, and aliphapaty.
for completeness, the  <dig>   <dig> and  <dig> sequences that could not be classified using a phylogenetic argument but that show amino acid patterns typical of hoxl, nkl and prd genes were also used.
chemical properties observed in more than 99% of the sequences belonging to one class, and observed in less than three-quarters of sequences belonging to another class are shown in tables  <dig>   <dig>  and  <dig>  rules are described relative to the inferred common ancestor situation of the two categories being compared.
for instance, when comparing prd and antp class genes the rule "not negatively charged at position 27" means that it was inferred  that the common ancestor to these two gene classes could use a negatively charged amino acid at position 27; in this example, antp genes subsequently almost completely lost the ability to use a negatively charged amino acid, but 94% of prd genes retained the use of such an amino acid at this position .
an almost complete change in amino acid usage regarding polarity and hydrophobicity is also observed for amino acid position  <dig>  furthermore, about 30% of the prd genes show a charged amino acid at position  <dig>  hoxl and nkl genes never use charged amino acids at this position, and this is also the case for lim genes, here used as an outgroup to antp and prd genes.
therefore, the possibility of using a charged amino acid at position  <dig> seems to be a derived feature that appeared in the prd lineage.
this amino acid position has been previously used to sub-classify prd genes into three categories  <cit> .
a charged amino acid is found in sequences from caenorhabditis elegans, strongylocentrotus purpuratus, and hydra vulgaris, among others, thus this switch is an event that happened early in the evolution of prd genes.
in brackets is indicated the total number of sequences analysed.
phy – set of sequences classified using a phylogenetic approach
unc – set of sequences not classified using a phylogenetic approach.
these sequences, nevertheless, show amino acid patterns typical of prd, hoxl and nkl genes.
ne – negatively charged amino acids; p – positively charged amino acids; a -aromatic amino acids; lp – less-polar amino acids; np – non-polar amino acids; lh – less hydrophobic amino acids; vh – very hydrophobic amino acids; s- small amino acids; t – tinny amino acids; al – alyphatic amino acids.
in brackets is indicated the total number of sequences analysed.
legend as in table 4
in brackets is indicated the total number of sequences analysed.
legend as in table 4
genes belonging to the hoxl lineage show derived constraints at amino acid positions  <dig>   <dig>   <dig>   <dig>   <dig> and  <dig>  relative to the inferred ancestral state for antp class genes .
only three sequences classified as hoxl do not follow the general pattern for these genes, namely, the hma <dig> gene  from helobdella triserialis that shows an aromatic amino acid at position  <dig>  and the hxb <dig> and pal <dig> genes  from sus scrofa and caenorhabditis elegans, respectively, that show small amino acids at position  <dig>  no derived single amino acid constraints were found for the nkl lineage.
discussion
the amino acid patterns and single key amino acid changes here identified shed light on some of the major likely functional changes that occurred during the evolution of the hoxl, nkl and prd genes.
our results show that important amino acid changes happened very early in the evolution of these genes, and thus it is possible to identify an archetype for bilaterian prd, antp, hoxl and nkl genes.
as with every generalization, however, some genes do not fit the archetype.
experimental studies are now needed in order to understand why the archetypes possess such chemical properties.
we propose the following archetypes for homeodomain amino acid sequences showing the antp-prd pattern .
prd genes show the l         p  prd pattern, use a negatively charged amino acid at position  <dig>  and a less polar or non-polar amino acid  at position  <dig>  about 80% of all animal bilaterian prd genes follow this description.
antp  genes show the     antp  <dig> pattern.
in contrast to prd genes, antp genes do not use a negatively charged amino acid at position  <dig> .
furthermore, at position  <dig>   <dig> % of all bilaterian antp genes use a polar amino acid .
97% of all antp bilaterian genes follow this description .
hoxl genes are characterized by the presence of the hoxl  <dig> pattern.
96% of all bilaterian hoxl genes show the  wfqnrr k     hoxl  <dig> pattern.
nkl genes show the nkl pattern  or a nkl derived pattern.
about 98% of the bilaterian nkl genes follow this pattern.
if the nkl pattern is generalized to include the derived nkl patterns then many hits are observed in other homeodomain classes .
this is why several nkl-related patterns are reported.
the above definitions suggest that the region in between amino acid positions 16– <dig> and 44– <dig> are the most important for function specificity of genes belonging to different classes and sub-classes.
the former region corresponds to the end of helix  <dig>  the inter helix 1– <dig> region and to the first positions of helix  <dig>  <cit> .
the region in between amino acid positions 44– <dig> correspond mainly to helix  <dig>  this helix, also called recognition helix is essential for successful and specific dna binding  <cit> .
the situation observed for nkl genes suggest that specific amino acid patterns may exist for groups of genes within major classes, as previously suggested  <cit> .
it also indicates that not all motifs implied by an amino acid pattern occur in a given homeodomain class.
on the other hand, the situation observed for non-bilaterian vs bilaterian genes suggests that we may have failed to identify all relevant chemical changes.
for instance, it is conceivable that important amino acid changes are observed when the protostome and the deuterostome lineages are compared; such distinctions were not addressed in this work because we grouped genes from a given homeodomain class together, irrespective of the organism.
it is tempting to use these features to classify homeobox containing genes.
nevertheless, given the very old age of the gene families considered, the possibility of convergent evolution must be considered when analyzing a given amino acid sequence.
therefore, for the purpose of gene classification, the features here described should not be decisive but rather used as additional evidence.
it should be noted, that the reconstruction of gene genealogies is a hard problem, namely when proteins belonging to the same family share some of the same protein interaction partners, thus facing a similar selective environment .
any additional piece of evidence that may shed light on the correct classification of genes should therefore be used.
CONCLUSIONS
in this work we report a method for the fast retrieval of bilaterian antp and prd sequences.
given the availability of a sufficiently large curated data set this method can be applied to any group of proteins.
furthermore, we report some of the main amino acid changes that occurred early in the hoxl, nkl and prd lineages.
these features can be used for the classification of gene sequences, although, as shown, convergent evolution must be considered as an explanation for the presence of a given pattern in a sequence.
the possibility that the region of the protein that allows distinguishing the different classes, also allows the distinction of different families within classes, has important practical and evolutionary consequences and must be explored in more detail.
authors' contributions
naf, jv and cpv conceived the design of the study.
naf developed the software needed for the large-scale analyses, while jv and cpv collected, aligned the sequence data, and performed the phylogenetic analyses.
jv drafted the manuscript.
all authors participated in the results discussion and helped writing the final version of the manuscript.
all authors read and approved the final manuscript.
appendix
gi numbers of the  <dig> non-bilaterian non-redundant sequences used for the analysis presented in table  <dig>  all sequences show the antp-prd amino acid pattern.
sequences showing patterns antp-prd, antp <dig>  antp <dig>  hoxl <dig> and hoxl2
nematostella vectensis: 32263856; 32263873; 74039490; 74039492; 82621559; 82621587; 82621611; 82621663; 83356315; 156363224; 156363226; 156387518; 156397205
sequences showing patterns antp-prd, antp <dig>  antp <dig>  hoxl <dig> and nkl
ephydatia fluviatilis: 438584
sequences showing patterns antp-prd, antp <dig>  antp <dig> and nkl
acropora millepora: 117581722; 117581724; 117581726; 117581728; anemonia erythraea: 158936936; ephydatia fluviatilis: 438585; 510584; halichondria bowerbanki: 33641771; nematostella vectensis: 32816231; 78190373; 82621509; 82621525; 82621533; 82621555; 82621557; 82621585; 82621591; 82621609; 82621625; 82621657; 82621665; 82621667; 110339021; 110339029; 110339121; 156353243; 156363798; 156367335; 156372762; 156397943; 168693291;suberites domuncula: 34786940; 49659003; sycon raphanus: 11066243
sequences showing patterns antp-prd, antp <dig> and nkl
anthopleura japonica: 144369330; 144369363; hydra vulgaris: 2331219; 7635735; nematostella vectensis: 82570553; 82621543; 82621677; 110339077; 110339133; 156376845; 156402818; 156406963; 156407174
sequences showing patterns antp-prd, antp <dig>  antp <dig> and en
nematostella vectensis: 82621563; 82621631; 82621655; 156366927; 156402173
sequences showing patterns antp-prd, antp <dig>  antp <dig> and noto
nematostella vectensis: 82621579; 82621647; 82621651; 110339063; 156388033; 156400908; 156407394
sequences showing patterns antp-prd, antp <dig>  antp <dig> and ventx
nematostella vectensis: 110339027
sequences showing patterns antp-prd, antp <dig> and antp2
baikalospongia intermédia: 62238212; ephydatia fluviatilis: 1438870; ephydatia muelleri: 3184520; hydra vulgaris: 7635740; nematostella vectensis: 82621517; 82621615; 82621619; 82621623; 82621635; 82621653; 156375827; 156375891; 156376334; 156401296; 156402692; podocoryne carnea: 28188799; potamolepis sp.
: 62238208; suberites domuncula: 49659001
sequences showing patterns antp-prd and prd
anthopleura japónica: 144369334; cladonema californicum: 9652040; hydra vulgaris: 3021450; 3021452; 6503072; nematostella vectensis: 38569883; 78190377; 82395396; 82395402; 82395404; 82395406; 82395408; 82395412; 82395414; 82570541; 82570557; 82570559; 82570563; 110339183; 110339213; 110339215; 110339225; 113120207; 156377162; 156615306;tripedalia cystophora: 33391193
sequences showing other pattern combinations
acropora formosa: 228960; acropora millepora: 7335704; 7595811; 7595813; 13506878; 117581730; 117581732;anthopleura japónica: 144369323; 144369326; 144369366; aurelia aurita: 50841484; cassiopea xamachana: 4894653; 4894655; 4894659; cladonema radiatum: 47155918; 47155920; 47155922; eleutheria dichotoma: 1147626; 91982983; 91982989; hydra littoralis: 2102728; hydra magnipapillata: 630481; 630482; 89242120; 144369369; 144369375; hydra viridis: 7120; 7124; 7130; 83763566; hydra vulgaris: 4433647; 4838455; 7635731; 7635733; 7635737; 7635742; 9945022; hydractinia symbiolongicarpus: 2980868; 83272159; nematostella vectensis: 5081328; 32816237; 74039494; 78190375; 82395398; 82395410; 82570519; 82570521; 82570527; 82570529; 82570537; 82570539; 82570545; 82570549; 82570555; 82570561; 82621507; 82621511; 82621513; 82621515; 82621519; 82621523; 82621527; 82621531; 82621535; 82621539; 82621541; 82621547; 82621549; 82621551; 82621565; 82621567; 82621571; 82621573; 82621575; 82621589; 82621595; 82621601; 82621603; 82621605; 82621613; 82621621; 82621627; 82621629; 82621633; 82621637; 82621643; 82621645; 82621671; 82621675; 83356317; 99030986; 110338989; 110339023; 110339061; 110339099; 110339101; 110339115; 110339171; 110339191; 110339211; 110339247; 113120203; 156356964; 156358580; 156361301; 156361303; 156364678; 156368221; 156371265; 156371439; 156372678; 156374167; 156377158; 156377160; 156377164; 156389076; 156390886; 156393340; 156395806; 156396978; 156398319; 156407912; podocoryne carnea: 976094; 6118056; 6465862; 7649930; 9964019; 47155914; 62002543; sarsia sp.
: 9988771; scolionema suvaense: 158936914; suberites domuncula: 49659005; trichoplax adhaerens: 38565482

 genes that belong to the homeobox family are characterized by the ability to code for a protein that contains a recognizable, although very variable,'homeodomain', usually amino acids in length. many of these genes are transcription factors that play important roles in the embryonic development of bilaterian and non-bilaterian animals. changes in homeobox gene content and deployment during evolution may have contributed to the evolution of body plan differences in animals. therefore, comparison of homeobox gene sets from different animals may shed light on the evolutionary events that gave rise to animal body plan diversity. nevertheless, gene orthology is not always easy to establish when comparing divergent animals. for this purpose, phylogenetic analysis, conservation of synteny, paralogy within the human genome, insertions within the homeodomain, key amino acid residues, and several motifs outside of the homeodomain can all be used. these features can also be used to classify homeobox genes into classes, subclasses and families. in the latest revision, holland et al. classified all human homeobox genes into classes and gene families. the antp class is further divided into hoxl and nkl subclasses. it should be noted that the only protein region that can be aligned in all gene families or genes is the homeodomain. here, we report amino acid patterns typical of bilaterian hoxl, nkl, and prd genes that can be used to quickly and efficiently retrieve amino acid sequences belonging to these classes and subclasses, among hundreds of other homeodomain sequences. retrieving a given class or subclass of sequences from many animal genomes may be thus an easier task than previously thought. however, we show that these typical amino acid patterns should be cautiously used when sequences come from non-bilaterian animals. since phylogenetic analysis was one of the primary sources of evidence used to establish the different classes and sub-classes, it is likely that most of homeobox gene classes and subclasses represent monophyletic lineages. hence, it is conceivable that amino acid changes important for protein function in a given lineage may be revealed as fixed differences between classes and subclasses. previous attempts have been made to classify genes within hox families but not at the level of whole classes or subclasses. here, we show that, because of their chemical properties, amino acid usage is different in antp and prd classes at five positions. furthermore, at nine positions, amino acid usage is different between hoxl and nkl subclasses. our findings support the notion that many chemically important changes happened early in the evolution of homeobox genes, and that these changes can be used as additional evidence to establish gene orthology. these results can be helpful for experimental studies aimed onto investigating the biochemical functions of key homeodomain residues in different gene classes as well. methods data we have used the hand-curated data set of holland et al., the prosite homeodomain data set and the ncbi database. identification of amino acid patterns typical of hoxl, nkl and prd genes in order to find amino acid patterns that distinguish genes from the antp subclasses, as well as genes from the prd class, from genes belonging to other classes, a fast word discovery program was first used to find statistically interesting words. the minimum word length used was two and the minimum number of sequences where the word should occur was set to 10% of the number of sequences in the data set being considered. the words found were then filtered to discard those words that occurred more than three times in more than one dataset. the software bioredx was then used to find the largest and more general amino acid patterns that distinguish two sets of sequences. this program accepts as input two sequence files and a word seed. based on the seed given, tries to find patterns that occur in one file and not in the other. this software was used to find patterns that occurred less than four times in the negative file and considered patterns up to amino acid residues larger than the initial word seed. in the reported amino acid patterns, amino acids listed within brackets are those allowed at a given position. for compactness of representation it is also possible to negate the class. the negation is denoted by "^". in this case, the amino acid residues listed are the ones that are not allowed at that particular position. the approach here used to find amino acid patterns does not use a set of aligned sequences. therefore, in order to make sure that the derived amino acid patterns are found in the same region of the homeodomain sequence, all sequences belonging to a given class and sub-classes were aligned using the pam scoring matrix. identification of key amino acid residues there are many schemes for comparing and grouping amino acids. nevertheless, none of them can possibly capture the vast number of contexts in which amino acids are found within proteins. the scheme here used is that of livingstone and barton, based on the amino acid properties size, polarity, hydrophobicity, charge, aliphaticity and aromaticity. for each position and property we calculated the probability of gaining or losing a constraint using the distribution of the amino acid property's observed at that position in all sequences considered. only changes in the properties values with a probability of occurrence lower than 5% were considered. phylogenetic analyses a neighbour-joining tree, using pair-wise deletion, as implemented in the mega software was constructed in order to classify the set of homeodomain sequences from a variety of animal species, into the classes and two sub-classes scheme proposed by holland et al..  in order to establish a method for fast retrieval of sequences belonging to a given class or subclass, amino acid patterns characteristic of different homeobox classes and subclasses were sought. as a starting point we used the hand-curated human prd, nkl and hoxl data sets compiled by holland et al.. the homeodomains of human en en dlx dlx dlx dlx dlx dlx noto, and hopx are difficult to classify and were not used in these initial analyses. therefore, the prd, nkl and hoxl data sets contain and sequences, respectively. within each of the three data sets, the amino acid sequences show considerable variation, suggesting that the human sequences include the majority of amino acid variability allowed at a given position along the sequence. this assumption was later tested. 1) amino acid patterns characteristic of the different human gene classes and sub-classes many characteristic amino acid patterns were detected using the approach described in material and methods. here, we describe the seven patterns that are most pertinent, taking into account their coverage across a group of genes and absence in the other groups. the pattern wfqnrr k, named hoxl pattern, is present in all human hoxl homeodomains, and in none of the human nkl and prd genes. the pattern le e, named hoxl pattern, is present in all human hoxl homeodomains except mnx only two human prd genes encode this amino acid pattern. none of the human nkl genes encode this amino acid pattern. the pattern, named nkl pattern, is present in out of human nkl sequences, in one human hoxl sequence, and in none of the human prd sequences. the human nkl sequences that do not present this pattern are: hhex, nanog, and ventx. the pattern l p, named prd pattern, is present in all human prd homeodomains and never present in human hoxl and nkl sequences. the patterns, named antp pattern, and w qn r k, named antp pattern, are present in all human hoxl and nkl genes and never present in human prd homeodomains. the pattern w, named antp-prd pattern, is found in all human prd, nkl and hoxl homeodomains. 2) generality of the amino acid patterns found in order to test the generality of the amino acid patterns derived in the previous section, we used the homeodomain sequences classified by holland et al., which include the human sequences used above. these sequences include genes classified as hoxl, nkl and prd, plus other homeobox gene classes. it should be noted, that in these analyses, the difficult to classify genes that were excluded above were now included, and these were classified as tentatively suggested by holland et al.. using amino acid sequences, a neighbour-joining tree was built with the sequences, plus homeodomain sequences from a variety of animal species. the latter sequences were obtained from the file ps available at prosite after removing all non-animal sequences from the file. sequences from clusters that are supported by a bootstrap value of 80% or higher, and that include at least one sequence that has been classified by holland et al., as belonging to a given class or subclass, were classified as belonging to that class or subclass. using this phylogenetic argument sequences could be classified as hoxl, as nkl, as prd, as lim, as pou, as hnf, as sine, as tale, as cut, as pros, as zf, and as cers. nevertheless, sequences remained unclassified. uncl. – unclassified sequences. the total number of sequences in each class is shown in parentheses. the hoxl pattern is not as specific as the hoxl pattern, since it is observed in about 46% of the genes belonging to the pou class and in a few nkl and prd genes. all prd genes showing the hoxl pattern belong to the pax4/ gene family, thus it is likely that this is a case of convergent evolution. the human pou homeodomains showing hoxl pattern are pou1f pou3f pou3f pou3f and pou3f the latter four genes are closely related, although the pou1f gene is distantly related to those genes. nevertheless, this is also likely a case of convergent evolution, since the alternative hypothesis implies many independent losses. therefore we argue that the necessity of using the amino acid combinations implied by hoxl pattern is a derived feature that appeared early in the hoxl lineage. only three sequences classified as hoxl using a phylogenetic argument, do not show the hoxl pattern, namely: drosophila melanogaster exex, human mnx and mouse mnx. mnx genes are peculiar in other ways. for instance, they show both the nkl and hoxl pattern. the nkl pattern is highly specific but not widely representative of all nkl genes. the pattern is found in 57% of the sequences classified as belonging to the nkl class, but in only ten other classified genes. in addition, of the phylogenetically'unclassified' genes show the nkl pattern, and of these have already been classified as belonging to the nkl class by prosite. the three hoxl genes with the nkl pattern are the three mnx genes. mnx genes also show the hoxl pattern that is highly specific for hoxl genes. thus, it is unlikely that mnx genes have been misclassified. therefore, the presence of the nkl pattern in mnx sequences may be the result of convergent evolution. the seven sine sequences where the nkl pattern is found are six from oryzias latipes, gallus gallus, mus musculus and homo sapiens and six from gallus gallus, mus musculus and homo sapiens. all six human genes cluster together with a high bootstrap value. nevertheless, the nkl pattern is only observed in the six and six genes, two closely related genes. it is thus, likely a case of convergent evolution. the requirement to use one of the amino acid combinations implied by the nkl pattern is a derived feature that appeared likely early in the nkl lineage. most of the classified sequences where the nkl pattern is not found have been annotated by holland et al., and prosite as belonging to the nanog, noto, ventx, en, dlx, and barx, families. the en, dlx and noto genes are difficult to classify. thus, information on these genes was not used to derive the nkl amino acid pattern. the prd pattern is highly specific, being found only in sequences classified as prd, and is also highly representative being found in 86% of the sequences classified as prd. the prd genes that do not show this pattern are d. melanogaster odsh, otp, phdp, ptx ip hopx genes from danio rerio, homo sapiens, bos taurus, rattus norvegicus, mus musculus, sus scrofa, and gallus gallus, human pax pax and pax, hydra vulgaris dmbx, otp genes from heliocidaris erythrogramma, heliocidaris tuberculata, lytechinus variegates, paracentrotus lividus, and saccoglossus kowalevskii, otx from strongylocentrotus purpuratus, alx from strongylocentrotus purpuratus, anf genes from gallus gallus, and xenopus laevis, and the prop gene from canis familiaris. the prd sequences that do not show this pattern are not related in any particular way. therefore, the absence of this pattern in these sequences is likely the result of several independent losses. therefore, it seems likely that the requirement to use one of the amino acid combinations implied by the prd pattern is a derived feature that appeared early in the evolution of prd genes. since the prd-specific pattern is highly specific, it is likely that the unclassified sequences that show this pattern are also prd sequences. according to prosite classification these genes are gsc from danio rerio, xenopus laevis, gallus gallus, mus musculus, saguinus labiatus, gorilla gorilla, pongo pygmaeus, pan paniscus, and pan troglodytes, alx from mus musculus and bos taurus, and unc- and ceh- genes from caenorhabditis elegans. according to prosite, of these genes belong to the prd class; the status of the ceh- gene is unknown. the antp pattern is found in most hoxl and nkl gene sequences and is almost absent in prd gene sequences. nevertheless, such a pattern is also found in sequences from genes belonging to all other classes except the pros class. it should, however, be noted that the sample size of pros genes is very small. the antp sequences that do not show this pattern are drosophila melanogaster eve, lbe, and lbl, homo sapiens noto, heterodontus francisci evx, gallus gallus hme, fugu rubripes hxdbb, and caenorhabditis elegans vab and hm. the hxb gene from gallus gallus is incomplete, and thus it is not possible to determine whether it shows this pattern. the hopx sequence from danio rerio is the only prd sequence that shows antp pattern. although the antp pattern is short, the broad distribution indicates that, very likely, it is not the result of convergent evolution. it is more likely that the ability to use the amino acid combinations implied by this pattern is an ancestral feature of homeobox containing genes. for some likely functional reason, antp genes have retained such a pattern, in contrast with prd genes where only % of all sequences show this pattern. the antp pattern is found in most hoxl but in only 68% of the nkl sequences. this pattern is not found in sequences from other classes. it is thus, highly specific, although less representative across the nkl subclass of antp class genes. the hoxl genes that do not show this pattern are hmd from gallus gallus, hxaba from danio rerio, hxb from salmo salar and hxbb from fugu rubripes. these sequences are not related in any particular way. it is not possible to determine whether the hxb gene from gallus gallus shows this pattern, since this is a partial sequence. nkl sequences not showing the antp pattern belong to the dlx, en and noto gene families that have been previously found to be difficult to classify. hence, information on these genes was not used to derive the antp pattern. once again, genes from these families stand out as oddities. the antp-prd pattern identifies most hoxl, nkl and prd sequences used. nevertheless, it is also present in all sine sequences and in 79% of the lim sequences. the antp plus prd sequences that do not show the expected pattern are the hopx genes from homo sapiens, mus musculus, rattus norvegicus, sus scrofa, bos taurus, gallus gallus, and danio rerio, the hm gene from caenorhabditis elegans, and the artemia sanfranciscana hmen gene. it should be noted that hopx are difficult to classify, that in the latest revision were classified as prd genes, and these were not used to derive the amino acid patterns being tested. the pax pax and pax sequences from homo sapiens, as well as the hxb gene from gallus gallus, the hxdbb gene from fugu rubripes, and the dlx and dlx genes from eleutherodactylus coqui, are partial, thus it is not possible to determine whether they show the antp-prd amino acid pattern. 3) further characterization of nkl genes given the failure to identify about 43% of the classified nkl sequences using the nkl pattern, we performed additional analyses to see whether the nkl pattern could be refined to accommodate the members of the nanog, noto, ventx, en, dlx, and barx, families. the results are shown in table when the nkl and refined patterns are considered, 98% of all nkl sequences could be classified as such. it should be noted that these refined nkl patterns are almost completely absent in the hoxl, prd, lim, pou, hnf, sine, tale, cut, pros, zf and cers classes. the hoxl sequences that contain the en pattern but not the nkl pattern all belong to the hox gene family. this happens because the en pattern allows for an asparagine at position five out of the six genes that show the nanog pattern but not the nkl pattern belong to the hoxl gsx gene family. these sequences are found because the pattern allows for a lysine at position the hoxl sequences that show a hit when using the ventx pattern but not the nkl pattern all belong to the gbx family. this happens because a valine is now allowed at position the three cut genes that show the en and ventx patterns are the homo sapiens onecut gene and the orthologues in mus musculus and rattus norvegicus. these sequences show a hit because both patterns allow for an isoleucine at position uncl. – unclassified sequences. the total number of sequences in each class is shown in parentheses. dlx-barx –; en –; noto –; nanog –; ventx – 4) non-bilaterian homeobox sequences only a few non-bilaterian homeobox sequences are contained in the data sets used. therefore, we collected from the ncbi database a set of non-redundant non-bilaterian homeobox sequences that encompass the regions where the amino acid patterns here reported are located, and that showed the antp-prd pattern derived above. most bilaterian antp and prd sequences show this antp-prd pattern; furthermore, this pattern is only found in antp, prd, lim and sine sequences. therefore, by imposing the presence of the antp-prd pattern we hoped to enrich the data set for non-bilaterian antp and prd sequences. table summarizes the results. as expected, of the retrieved gene sequences only one seems to belong to gene classes other than antp, prd, lim and sine. antp1; antp2; antp1; antp2; antp1; antp2; antp1; antp2; antp1; antp2; antp1; antp2; only those non-bilaterian sequences that encompass the regions where the amino acid patterns here described are located, and that showed the antp-prd pattern were used. in total non-redundant sequences were used. cases where non-bilaterian genes are clearly misidentified when using amino acid patterns are shown underlined. uncl. – unclassified sequences. * this signature can also be considered characteristic of nkl genes since only 68% of bilaterian nkl sequences show the antp pattern most non-bilaterian genes showing the full nkl or prd signature have already been classified as such. furthermore, the antp demox class found in demospongiae and apparently absent in all other animals, is characterized by the presence of the antp-prd, and antp and patterns, as expected for antp genes that do not belong to the hoxl or nkl lineage. although there are no true hox genes in sponges, sequences show the full hoxl signature. interestingly, all such sequences are non-annotated, thus they must be hard to classify gene sequences. one ephydatia fluviatilis gene has been labeled as the msx gene. our pattern analyses show that the corresponding sequence shows a pattern exclusively found in the mnx gene family, thus, also suggesting that non-bilaterian sequences with full signatures may be wrongly identified. about 53% of all non-bilaterian sequences used did not show an easily recognizable signature. we conclude that the reported amino acid patterns must be cautiously used when applied to sequences from non-bilaterian animals. 5) identification of single key amino acid changes in bilaterian prd, nkl and hoxl genes key amino acid changes that occurred during the evolution of homeobox containing genes may be revealed as changes affecting a single amino acid position. thus, we compared, at each amino acid position, the large hoxl, nkl and prd data sets described above, for the following chemical properties: size, charge, polarity, hydrophobicity, aromaticity, and aliphapaty. for completeness, the and sequences that could not be classified using a phylogenetic argument but that show amino acid patterns typical of hoxl, nkl and prd genes were also used. chemical properties observed in more than 99% of the sequences belonging to one class, and observed in less than three-quarters of sequences belonging to another class are shown in tables and rules are described relative to the inferred common ancestor situation of the two categories being compared. for instance, when comparing prd and antp class genes the rule "not negatively charged at position 27" means that it was inferred that the common ancestor to these two gene classes could use a negatively charged amino acid at position 27; in this example, antp genes subsequently almost completely lost the ability to use a negatively charged amino acid, but 94% of prd genes retained the use of such an amino acid at this position. an almost complete change in amino acid usage regarding polarity and hydrophobicity is also observed for amino acid position furthermore, about 30% of the prd genes show a charged amino acid at position hoxl and nkl genes never use charged amino acids at this position, and this is also the case for lim genes, here used as an outgroup to antp and prd genes. therefore, the possibility of using a charged amino acid at position seems to be a derived feature that appeared in the prd lineage. this amino acid position has been previously used to sub-classify prd genes into three categories. a charged amino acid is found in sequences from caenorhabditis elegans, strongylocentrotus purpuratus, and hydra vulgaris, among others, thus this switch is an event that happened early in the evolution of prd genes. in brackets is indicated the total number of sequences analysed. phy – set of sequences classified using a phylogenetic approach unc – set of sequences not classified using a phylogenetic approach. these sequences, nevertheless, show amino acid patterns typical of prd, hoxl and nkl genes. ne – negatively charged amino acids; p – positively charged amino acids; a -aromatic amino acids; lp – less-polar amino acids; np – non-polar amino acids; lh – less hydrophobic amino acids; vh – very hydrophobic amino acids; s- small amino acids; t – tinny amino acids; al – alyphatic amino acids. in brackets is indicated the total number of sequences analysed. legend as in table 4 in brackets is indicated the total number of sequences analysed. legend as in table 4 genes belonging to the hoxl lineage show derived constraints at amino acid positions and relative to the inferred ancestral state for antp class genes. only three sequences classified as hoxl do not follow the general pattern for these genes, namely, the hma gene from helobdella triserialis that shows an aromatic amino acid at position and the hxb and pal genes from sus scrofa and caenorhabditis elegans, respectively, that show small amino acids at position no derived single amino acid constraints were found for the nkl lineage. discussion the amino acid patterns and single key amino acid changes here identified shed light on some of the major likely functional changes that occurred during the evolution of the hoxl, nkl and prd genes. our results show that important amino acid changes happened very early in the evolution of these genes, and thus it is possible to identify an archetype for bilaterian prd, antp, hoxl and nkl genes. as with every generalization, however, some genes do not fit the archetype. experimental studies are now needed in order to understand why the archetypes possess such chemical properties. we propose the following archetypes for homeodomain amino acid sequences showing the antp-prd pattern. prd genes show the l p prd pattern, use a negatively charged amino acid at position and a less polar or non-polar amino acid at position about 80% of all animal bilaterian prd genes follow this description. antp genes show the antp pattern. in contrast to prd genes, antp genes do not use a negatively charged amino acid at position. furthermore, at position % of all bilaterian antp genes use a polar amino acid. 97% of all antp bilaterian genes follow this description. hoxl genes are characterized by the presence of the hoxl pattern. 96% of all bilaterian hoxl genes show the wfqnrr k hoxl pattern. nkl genes show the nkl pattern or a nkl derived pattern. about 98% of the bilaterian nkl genes follow this pattern. if the nkl pattern is generalized to include the derived nkl patterns then many hits are observed in other homeodomain classes. this is why several nkl-related patterns are reported. the above definitions suggest that the region in between amino acid positions 16– and 44– are the most important for function specificity of genes belonging to different classes and sub-classes. the former region corresponds to the end of helix the inter helix 1– region and to the first positions of helix. the region in between amino acid positions 44– correspond mainly to helix this helix, also called recognition helix is essential for successful and specific dna binding. the situation observed for nkl genes suggest that specific amino acid patterns may exist for groups of genes within major classes, as previously suggested. it also indicates that not all motifs implied by an amino acid pattern occur in a given homeodomain class. on the other hand, the situation observed for non-bilaterian vs bilaterian genes suggests that we may have failed to identify all relevant chemical changes. for instance, it is conceivable that important amino acid changes are observed when the protostome and the deuterostome lineages are compared; such distinctions were not addressed in this work because we grouped genes from a given homeodomain class together, irrespective of the organism. it is tempting to use these features to classify homeobox containing genes. nevertheless, given the very old age of the gene families considered, the possibility of convergent evolution must be considered when analyzing a given amino acid sequence. therefore, for the purpose of gene classification, the features here described should not be decisive but rather used as additional evidence. it should be noted, that the reconstruction of gene genealogies is a hard problem, namely when proteins belonging to the same family share some of the same protein interaction partners, thus facing a similar selective environment. any additional piece of evidence that may shed light on the correct classification of genes should therefore be used.  in this work we report a method for the fast retrieval of bilaterian antp and prd sequences. given the availability of a sufficiently large curated data set this method can be applied to any group of proteins. furthermore, we report some of the main amino acid changes that occurred early in the hoxl, nkl and prd lineages. these features can be used for the classification of gene sequences, although, as shown, convergent evolution must be considered as an explanation for the presence of a given pattern in a sequence. the possibility that the region of the protein that allows distinguishing the different classes, also allows the distinction of different families within classes, has important practical and evolutionary consequences and must be explored in more detail. authors' contributions naf, jv and cpv conceived the design of the study. naf developed the software needed for the large-scale analyses, while jv and cpv collected, aligned the sequence data, and performed the phylogenetic analyses. jv drafted the manuscript. all authors participated in the results discussion and helped writing the final version of the manuscript. all authors read and approved the final manuscript. appendix gi numbers of the non-bilaterian non-redundant sequences used for the analysis presented in table all sequences show the antp-prd amino acid pattern. sequences showing patterns antp-prd, antp antp hoxl and hoxl2 nematostella vectensis: 32263856; 32263873; 74039490; 74039492; 82621559; 82621587; 82621611; 82621663; 83356315; 156363224; 156363226; 156387518; 156397205 sequences showing patterns antp-prd, antp antp hoxl and nkl ephydatia fluviatilis: 438584 sequences showing patterns antp-prd, antp antp and nkl acropora millepora: 117581722; 117581724; 117581726; 117581728; anemonia erythraea: 158936936; ephydatia fluviatilis: 438585; 510584; halichondria bowerbanki: 33641771; nematostella vectensis: 32816231; 78190373; 82621509; 82621525; 82621533; 82621555; 82621557; 82621585; 82621591; 82621609; 82621625; 82621657; 82621665; 82621667; 110339021; 110339029; 110339121; 156353243; 156363798; 156367335; 156372762; 156397943; 168693291;suberites domuncula: 34786940; 49659003; sycon raphanus: 11066243 sequences showing patterns antp-prd, antp and nkl anthopleura japonica: 144369330; 144369363; hydra vulgaris: 2331219; 7635735; nematostella vectensis: 82570553; 82621543; 82621677; 110339077; 110339133; 156376845; 156402818; 156406963; 156407174 sequences showing patterns antp-prd, antp antp and en nematostella vectensis: 82621563; 82621631; 82621655; 156366927; 156402173 sequences showing patterns antp-prd, antp antp and noto nematostella vectensis: 82621579; 82621647; 82621651; 110339063; 156388033; 156400908; 156407394 sequences showing patterns antp-prd, antp antp and ventx nematostella vectensis: 110339027 sequences showing patterns antp-prd, antp and antp2 baikalospongia intermédia: 62238212; ephydatia fluviatilis: 1438870; ephydatia muelleri: 3184520; hydra vulgaris: 7635740; nematostella vectensis: 82621517; 82621615; 82621619; 82621623; 82621635; 82621653; 156375827; 156375891; 156376334; 156401296; 156402692; podocoryne carnea: 28188799; potamolepis sp.: 62238208; suberites domuncula: 49659001 sequences showing patterns antp-prd and prd anthopleura japónica: 144369334; cladonema californicum: 9652040; hydra vulgaris: 3021450; 3021452; 6503072; nematostella vectensis: 38569883; 78190377; 82395396; 82395402; 82395404; 82395406; 82395408; 82395412; 82395414; 82570541; 82570557; 82570559; 82570563; 110339183; 110339213; 110339215; 110339225; 113120207; 156377162; 156615306;tripedalia cystophora: 33391193 sequences showing other pattern combinations acropora formosa: 228960; acropora millepora: 7335704; 7595811; 7595813; 13506878; 117581730; 117581732;anthopleura japónica: 144369323; 144369326; 144369366; aurelia aurita: 50841484; cassiopea xamachana: 4894653; 4894655; 4894659; cladonema radiatum: 47155918; 47155920; 47155922; eleutheria dichotoma: 1147626; 91982983; 91982989; hydra littoralis: 2102728; hydra magnipapillata: 630481; 630482; 89242120; 144369369; 144369375; hydra viridis: 7120; 7124; 7130; 83763566; hydra vulgaris: 4433647; 4838455; 7635731; 7635733; 7635737; 7635742; 9945022; hydractinia symbiolongicarpus: 2980868; 83272159; nematostella vectensis: 5081328; 32816237; 74039494; 78190375; 82395398; 82395410; 82570519; 82570521; 82570527; 82570529; 82570537; 82570539; 82570545; 82570549; 82570555; 82570561; 82621507; 82621511; 82621513; 82621515; 82621519; 82621523; 82621527; 82621531; 82621535; 82621539; 82621541; 82621547; 82621549; 82621551; 82621565; 82621567; 82621571; 82621573; 82621575; 82621589; 82621595; 82621601; 82621603; 82621605; 82621613; 82621621; 82621627; 82621629; 82621633; 82621637; 82621643; 82621645; 82621671; 82621675; 83356317; 99030986; 110338989; 110339023; 110339061; 110339099; 110339101; 110339115; 110339171; 110339191; 110339211; 110339247; 113120203; 156356964; 156358580; 156361301; 156361303; 156364678; 156368221; 156371265; 156371439; 156372678; 156374167; 156377158; 156377160; 156377164; 156389076; 156390886; 156393340; 156395806; 156396978; 156398319; 156407912; podocoryne carnea: 976094; 6118056; 6465862; 7649930; 9964019; 47155914; 62002543; sarsia sp.: 9988771; scolionema suvaense: 158936914; suberites domuncula: 49659005; trichoplax adhaerens: 38565482 
our results help focus experimental studies onto investigating the biochemical functions of key homeodomain residues in different gene classes. although homeobox genes have been the subject of many studies, little is known about the main amino acid changes that occurred early in the evolution of genes belonging to different classes. the reported protein features will facilitate the robust classification of homeobox genes from newly sequenced bilaterian genomes. furthermore, we look for diagnostic amino acid residues that can be used to distinguish hoxl, nkl and prd genes. in this study, we report a method for the fast and efficient retrieval of sequences belonging to the antp and prd classes.
 although homeobox genes have been the subject of many studies, little is known about the main amino acid changes that occurred early in the evolution of genes belonging to different classes.  in this study, we report a method for the fast and efficient retrieval of sequences belonging to the antp and prd classes. furthermore, we look for diagnostic amino acid residues that can be used to distinguish hoxl, nkl and prd genes.  the reported protein features will facilitate the robust classification of homeobox genes from newly sequenced bilaterian genomes. nevertheless, in non-bilaterian genomes our findings must be cautiously applied. in principle, as long as a good manually curated data set is available the approach here described can be applied to non-bilaterian organisms as well. our results help focus experimental studies onto investigating the biochemical functions of key homeodomain residues in different gene classes. 