- Research article
- Open Access
Origin and diversification of the basic helix-loop-helix gene family in metazoans: insights from comparative genomics
- Elena Simionato†1,
- Valérie Ledent†2,
- Gemma Richards3,
- Morgane Thomas-Chollier2, 4,
- Pierre Kerner1,
- David Coornaert2,
- Bernard M Degnan3 and
- Michel Vervoort1, 5Email author
© Simionato et al; licensee BioMed Central Ltd. 2007
- Received: 21 November 2006
- Accepted: 02 March 2007
- Published: 02 March 2007
Molecular and genetic analyses conducted in model organisms such as Drosophila and vertebrates, have provided a wealth of information about how networks of transcription factors control the proper development of these species. Much less is known, however, about the evolutionary origin of these elaborated networks and their large-scale evolution. Here we report the first evolutionary analysis of a whole superfamily of transcription factors, the basic helix-loop-helix (bHLH) proteins, at the scale of the whole metazoan kingdom.
We identified in silico the putative full complement of bHLH genes in the sequenced genomes of 12 different species representative of the main metazoan lineages, including three non-bilaterian metazoans, the cnidarians Nematostella vectensis and Hydra magnipapillata and the demosponge Amphimedon queenslandica. We have performed extensive phylogenetic analyses of the 695 identified bHLHs, which has allowed us to allocate most of these bHLHs to defined evolutionary conserved groups of orthology.
Three main features in the history of the bHLH gene superfamily can be inferred from these analyses: (i) an initial diversification of the bHLHs has occurred in the pre-Cambrian, prior to metazoan cladogenesis; (ii) a second expansion of the bHLH superfamily occurred early in metazoan evolution before bilaterians and cnidarians diverged; and (iii) the bHLH complement during the evolution of the bilaterians has been remarkably stable. We suggest that these features may be extended to other developmental gene families and reflect a general trend in the evolution of the developmental gene repertoires of metazoans.
- bHLH Protein
- bHLH Gene
- bHLH Domain
- Nematostella Vectensis
- bHLH Family
The basic helix-loop-helix (bHLH) protein superfamily constitutes an ancient class of eukaryotic transcription factors that are found in fungi, plants and metazoans [1, 2]. The bHLH transcription factors are named after their highly conserved domain (about 60 amino acids long) that consists of a DNA-binding basic region (b) followed by two α-helices separated by a variable loop region (HLH) . Interaction between the helix regions of two different proteins leads to the formation of homodimeric or heterodimeric complexes, and the basic region of each partner recognizes and binds to a core hexanucleotide DNA sequence. Many bHLH proteins also include additional domains that are important for their activity as transcriptional regulators, such as 'leucine zipper', 'PAS' or 'orange' domains, which are mainly involved in protein-protein interactions [3–5]. In unicellular eukaryotes, such as Saccharomyces cerevisiae, bHLH proteins mainly regulate metabolic pathways [1, 6]. In contrast, in metazoans and plants, the bHLH proteins are mainly involved in controlling developmental processes, in regulating the cell cycle, and in sensing environmental signals [1, 2, 7, 8].
Because of the important functions that they display in various organisms, bHLH proteins have been the subject of a number of studies aimed at the identification of their full complement encoded by completely sequenced genomes. The putative full set of genes encoding bHLH proteins ('bHLH genes') has been reported for Saccharomyces cerevisiae (8 bHLH genes), Drosophila melanogaster (58), Caenorhabditis elegans (39), Homo sapiens (125), Ciona intestinalis (46), Arabidopsis thaliana (118 to 147), and Oryza sativa (131 to 167) [4, 6, 9–16]. In most of these studies, phylogenetic analyses of the amino acid sequences of the bHLHs were used to define orthologous families (that is groups of genes that derive from a common ancestor). In addition, these phylogenetic studies have enabled the definition of higher-order groups (named A, B, C, D, E and F) within the bHLH superfamily, which comprise evolutionarily related families of orthologs that share structural and biochemical properties [1, 2, 4, 17, 18]. In brief, groups A and B include bHLH proteins that bind core DNA sequences referred to as E boxes (CANNTG); respectively CACCTG or CAGCTG (group A) and CACGTG or CATGTTG (group B). Group C corresponds to the families of bHLH proteins known as 'bHLH-PAS', as they contain a 'PAS' domain in addition to the bHLH. They bind to ACGTG or GCGTG core sequences. Group D corresponds to HLH proteins that lack a basic domain and are hence unable to bind DNA. These proteins act as antagonists of group A bHLH proteins. Group E includes proteins related to the Drosophila Hairy and Enhancer of split bHLH (HER) proteins. These proteins bind preferentially to sequences referred to as N boxes (CACGCG or CACGAG). Most group E proteins also contain two characteristic domains in addition to the bHLH, the 'orange' domain and a WRPW peptide in their carboxyterminal part. Group F corresponds to the COE proteins, which lack a basic domain and are characterized by the presence of an additional domain, the COE domain, involved in both dimerization and DNA binding. Yeast and plant bHLHs are all included in group B .
Derivation of the putative full set of bHLH genes of 12 metazoan species
We first developed a list of bHLH genes encoded by the three genomes (Drosophila melanogaster, Caenorhabditis elegans, and Homo sapiens) we had previously studied [4, 11], using the current version of their genome assembly (Figure 1, additional file 1). While we found the same 39 bHLHs as in our previous study in Caenorhabditis elegans, we identified, in Drosophila melanogaster, one additional bHLH as compared to our published analysis, raising the total number of Drosophila bHLHs to 59. In Homo sapiens, our previous study was made on a draft assembly of the genome and a few bHLH genes that we identified cannot be found in the current assembly or correspond to pseudogenes. We found a total number of 118 bHLH genes in the current assembly of the genome of Homo sapiens.
We then used the sequences of all the bHLHs from Homo sapiens and Drosophila melanogaster to identify, through similarity searches using BLAST algorithm, the bHLHs encoded in the genomes of 9 additional species, the demosponge Amphimedon queenslandica (formerly named Reniera sp.), the cnidarians Hydra magnipapillata (a hydrozoan) and Nematostella vectensis (an anthozoan), the annelid Capitella sp. I (formerly known as Capitella capitata), the mollusk Lottia gigantea, the arthropods Daphnia pulex (a crustacean) and Tribolium castaneum (an insect), the echinoderm Strongylocentrotus purpuratus, and the chordate Branchiostoma floridae (a cephalochordate). These species were chosen because, together with the four metazoan species in which the putative full set of bHLHs has been previously determined (Drosophila melanogaster, Caenorhabditis elegans, Ciona intestinalis, and Homo sapiens; [4, 11, 12]), they provide a significant coverage of the main metazoan evolutionary lineages. The phylogenetic relationships between these species and the total number of bHLHs that have been found for each species are shown in Figure 1. Additional informations about the studied species and the status of the corresponding genome projects can be found in additional file 2. We identified a total number of 695 bHLH genes in the 12 aforementioned species. All the identified sequences can be found in additional file 3. In most cases, we were able to retrieve the complete bHLH domain. Given our extensive searches, we are confident that we obtained a significant coverage of the bHLH genes present in the genomes of all the studied species. We have used this very rich sampling of metazoan bHLH genes to better understand the evolution of this superfamily through phylogenetic analyses.
Phylogenetic analyses of the bHLH genes
Carrying out evolutionary analyses of multigene families requires that orthologs, which have evolved by vertical descent from a common ancestor, are distinguished from paralogs, which arise by duplication and domain shuffling within a genome . We therefore constructed phylogenetic trees to define metazoan families of orthologs. Given the large number of sequences, we independently analysed the sequences from the different species. We made 9 different multiple alignments, each comprising all the identified bHLHs from Homo sapiens and Drosophila melanogaster plus those from one of the nine other species. We also made two additional multiple alignments, one with the bHLHs from Homo sapiens and Drosophila melanogaster plus those from the two cnidarian species, Nematostella vectensis and Hydra magnipapillata and the other with the bHLHs from Homo sapiens and Drosophila melanogaster plus those from the two arthropod species, Tribolium castaneum and Daphnia pulex. These alignments allowed the study of specific relationships of bHLHs within cnidarians and arthropods. The multiple alignments were then used to construct phylogenetic trees. We used four different methods of phylogenetic reconstruction: distance (neighbour-joining; NJ), maximum parsimony (MP), maximum likelihood (ML), and Bayesian inference (BI). In general, the phylogenetic trees obtained by the different methods were congruent and displayed very similar topologies.
These phylogenetic trees were then used to define families of metazoan bHLH orthologs, using the same criterion as in our previous studies. That is, we defined bHLH families of orthologs as monophyletic groups which include sequences from different species and whose monophyly is consistent across the different phylogenetic methods and supported by bootstrap values and posterior marginal probabilities superior to 50% [4, 11]. This criterion was relaxed for the Mesp, Myc, and Hairy/E(spl) families as discussed previously [4, 11]. Genes that cannot be confidently assigned to any families (because for example their inclusion in a given family is weakly supported or not found in all the trees constructed by the different phylogenetic methods) are categorized as 'orphan' genes. The results of our analysis are consistent with those we previously obtained by using only the sequences from Drosophila melanogaster, Caenorhabditis elegans, and Homo sapiens . The additional sequences included in the present study allowed us to define 2 additional metazoan families of orthologs, Delilah and MyoRa and b previously grouped into the MyoR family (additional file 1). We moreover grouped together two previously defined bilaterian families of orthologs, Hairy and E(spl), into a single one as there is no clear support for the monophyly of either of these families in our current analysis (a detailed phylogenetic study of these families will be published elsewhere). These analyses identified 45 metazoan bHLH families of orthologs (see additional file 1 for the full list and for the number of members in each species; the data that have been used to construct additional file 1 and Figure 1 can be found in additional files 4 to 15, where each file corresponds to one species and contains a table with all the bHLHs from this species, the family to which each bHLH belongs, and the statistical support for its inclusion in the corresponding family in the trees constructed using the four different phylogenetic methods).
In the case of the bilaterian species, at the exclusion of Caenorhabditis elegans (in which the bHLH complement is much derived; ), we were able to confidently allocate most of the bHLHs (more than 90%) to defined metazoan families of orthologs (additional file 1). In addition, we found, in each species, at least one member for most families (39 to 44 of the 45 families), confirming the very good coverage of the bHLH genes present in the genome of these species (Figure 1, additional file 1). This also confirms one of the main conclusions of our previous work, the fact that most bHLH families were already present in Urbilateria, the last common ancestor of all bilaterians (see Discussion) .
Nematostella-specific duplications of bHLH genes
bHLHs from non bilaterians do not cluster with those of fungi
A recent report has shown that cnidarians may have retained some ancestral eukaryotic genes that have been lost in bilaterians . We were interested to investigate whether some of the 'orphan' cnidarian (and by extension Amphimedon queenslandica) bHLHs may be such genes. We therefore retrieved a large set of bHLHs from fungi, which are among the closest relatives of eukaryotes (reviewed in ) and for which many completely sequenced genomes are available. We identified 98 bHLHs from 13 different species as summarized in the additional file 18. We first performed phylogenetic analyses on a multiple alignment of these sequences with those of Homo sapiens and Drosophila melanogaster. We observed that the fungi bHLHs form several fungi-specific families of orthologs that do not include bilaterian genes (not shown). We nevertheless noted that two of these families have a tendency to cluster with the MITF and SREBP families but with very low statistical support (not shown). We then made multiple alignments of the fungi bHLHs with those of Amphimedon queenslandica or those of Nematostella vectensis and Hydra magnipapillata; all analyses also included Drosophila bHLHs (used here as representative of bilaterians). We did not find any monophyletic groups that included both fungi and cnidarian and/or sponge bHLHs (not shown). We therefore conclude that the 'orphan' bHLHs from the cnidarians and from Amphimedon queenslandica are probably not ancestral eukaryotic bHLHs that were lost in bilaterians.
We report in this article the in silico identification of 695 putative bHLHs encoded by 12 sequenced metazoan genomes. Our choice of species represents a significant sampling of metazoan diversity, as we include representatives of the main metazoan evolutionary lineages (Figure 1). Our dataset therefore allows us to study the evolution of the bHLH superfamily at the scale of the whole metazoan kingdom. This is, to our knowledge, the first time that an entire gene superfamily has been studied at such a large scale. We must however caution that our analysis has probably not been carried out on a fully comprehensive dataset of bHLHs sequences. In most cases, our search for bHLH genes has been with unannotated genome assemblies or even on unassembled whole genome shotgun reads. It is therefore possible that we may have missed some bHLHs and/or that we may have included bHLH domains from some pseudogenes. Nevertheless, these data are sufficient for the purpose of this study, i.e. to obtain a qualitatively accurate assessment of the metazoan bHLH complement and of its evolution within this lineage.
Ancestrality and stability of the bHLH complement in bilaterians
In contrast, the total number of bHLHs representing each family of orthologs is much more variable between different species (ranging from 50 to 118). For example, Branchiostoma floridae has 28 bHLHs more than Strongylocentrotus purpuratus (78 versus 50) but each species has approximately the same number of families (44 versus 43). This difference in the total number of bHLHs is due to the specific expansion of a few families in Branchiostoma floridae: Hairy/E(spl) and Hey (17 genes in Branchiostoma floridae versus 5 in Strongylocentrotus purpuratus); MyoD (4 versus 2); MyoRa (4 versus 1); and PTFb (3 versus 1). Another significant example is that of Drosophila melanogaster, which has more bHLHs than Daphnia pulex (59 versus 57) but in which fewer families are represented (39 versus 41). The larger number of bHLHs in Drosophila melanogaster as compared to Daphnia pulex is mainly due the expansion of one family, Hairy/E(spl) (11 or 12 versus 5 or 6). We can thus conclude that the number of bHLHs found in a given species is a poor indicator of the diversity of its bHLH complement, as some families of orthologs have undergone species- or clade-specific expansions. This holds true for Homo sapiens, although in this case most of the families have undergone varying degrees of expansion; 36 of the 44 families (about 80%) have more than one member in Homo sapiens while in other bilaterians the mean value is about 18% (range 7% to 26%). The bHLH complement in other vertebrates [; unpublished observations] is quite similar to that of Homo sapiens, indicating that an expansion of the bHLH complement occurred at some early stage of vertebrate evolution, probably related to large-scale duplication events, possibly of the whole genome (reviewed in ).
Cnidarian bHLHs: a complex ancestral complement, lineage-specific duplications and losses
We have found that 29 to 33 of the 44 conserved bilaterian bHLH families of orthologs have cnidarian members (Figure 1). Following the same reasoning as we used to infer the minimal complement of bHLHs in Urbilateria, we can therefore conclude that the last common ancestor of cnidarians and bilaterians (Ureumetazoa) possessed at least 29 to 33 bHLHs (Figure 5). The estimate given here for the bHLH complement of Ureumetazoa almost certainly represents an underestimate because (i) it is likely that additional bHLHs will be found in the future (we may have missed some genes due to the fact that our searches were made within unfinished cnidarian genome assemblies); (ii) some cnidarian sequences that were included in this analysis could not be placed unambiguously into specific metazoan families of orthologs (and may represent divergent members of additional ancestral families); (iii) our analyses were done on the two cnidarian genomes available and we cannot exclude the possibility that some ancestral bHLHs may have been lost in both species. Ureumetazoa therefore possessed a high number and diversity of bHLHs, similar to that found in bilaterians. Although initially surprising, this is fully consistent with recent studies performed on other gene families, which indicate a complex eumetazoan ancestral complement, for example, of homeobox genes (at least 56; ) and Wnt genes (at least 11; ). Taken together, we infer that Ureumetazoa possessed an elaborate complement of developmental genes and as such may have been morphogenetically complex.
Besides these speculations on Ureumetazoa, what appears clear is that cnidarians, and in particular anthozoans such as Nematostella vectensis, are neither 'simple', nor 'primitive' organisms, as is still often assumed. Nematostella vectensis possesses 68 bHLHs, i.e significantly more than most invertebrate bilaterians (Figure 1). Several recent studies have shown that Nematostella vectensis also possesses 130 to 139 homeobox genes including probably 7 Hox genes [21, 22, 28], 14 Sox (Sry-related HMG-box) genes , 15 Fox (Forkhead domain) genes , 12 WNT genes , and 6 TGF-β genes and numerous antagonists of these [30–32]. For several gene families, the number of genes found in Nematostella vectensis is in fact higher than that found in invertebrate models such as Drosophila melanogaster and Caenorhbaditis elegans [21, 22, 27, 33, 34]. Interestingly, many of these genes are developmentally expressed in Nematostella vectensis, and these patterns show some similarities to those found in bilaterians (e.g. [27–32]). So far, very few bHLH encoding genes have been studied at the expression level in cnidarians [35, 36]. It would thus be of particular interest to study whether the bHLH genes we have identified in Nematostella vectensis and Hydra magnipapillata are expressed during development and how these expression patterns may be related to those of their bilaterian counterparts.
It is also becoming clear, from this study and others, that in both the anthozoan Nematostella vectensis and the hydrozoan Hydra magnipapillata, several lineage specific genomic modifications have occurred. As already mentioned, the high number of bHLHs encoded by the genome of Nematostella vectensis is due to two factors, the inheritance of numerous (about 30) bHLH genes from Ureumetazoa ('ancestral complexity'), but also the expansion of some of the bHLH families through duplications ('lineage-specific duplications'). These two factors have contributed to the high level of complexity of the Nematostella vectensis genome, as has been previously suggested based on the analysis of a large set of expressed sequence tags (ESTs)  and is now clearly confirmed for two large multigenic superfamilies, the bHLH (this study) and the homeobox [21, 22, 34]. For both superfamilies, the lineage-specific duplicates are often physically linked, suggesting that they have been produced by tandem duplication, an apparently frequent event in the evolutionary lineage leading to Nematostella vectensis. It would now be interesting to see whether the duplications found in Nematostella vectensis are shared by other cnidarians. This is clearly not the case in Hydra magnipapillata, where only two metazoan families of orthologs (ASCa and Myc) contain more than one member (additional file 1). Although these families do also contain more than one member in Nematostella vectensis, the bHLHs from the ASCa and Myc families of the two cnidarians do not cluster in phylogenetic trees, suggesting independent duplications. Hydra magnipapillata has far fewer bHLHs than Nematostella vectensis (33 versus 68), again a similar situation to that reported for the homeobox genes (53 versus 139) . From our analysis, we can conclude that Hydra magnipapillata has lost at least 11 ancestral types of bHLHs (Mesp, Twist, Hand, SCL, Mnt, SREBP, Mlx, TF4, Clock, Bmal, and Hey; additional file 1). By comparison, only a single metazoan family of orthologs (Ahr) has a member in Hydra magnipapillata but not in Nematostella vectensis, therefore indicating that this family has been lost in Nematostella vectensis (additional file 1). Loss of conserved developmental genes may be a general trend in Hydra magnipapillata, which could reflect a general simplification of the genome in this lineage.
The Amphimedon queenslandica bHLH complement may reveal an intermediate step in the diversification of bHLHs in metazoans
The number of bHLHs and the number of represented families in Amphimedon queenslandica are markedly less than those found in cnidarians and bilaterians (Figure 1). This observation made on bHLHs can be extended to most developmental gene superfamilies: the genome of Amphimedon queenslandica encodes a considerable diversity of transcription factors and cell-cell molecules, but displays fewer members for most orthologs than cnidarians and bilaterians [; B.M.D. et al., unpublished observations]. Prior to making any evolutionary interpretations from these observations, in particular in terms of reconstruction of ancestral properties, we must caution that (i) we only have data on one sponge species and we can therefore not rule out the possibility that Amphimedon queenslandica may have secondarily lost many genes (e.g. as we observed in cnidarians for Hydra magnipapillata) and (ii) the sponges represent a very complex and diversified phyla that is probably paraphyletic [38–40]. Therefore Amphimedon queenslandica (a demosponge from the Haplosclerida group) cannot be considered as representative for all sponges. In particular, calcisponges may be the sister-group of eumetazoans (Figure 5) [38–40], (i.e. more akin to eumetazoans than to the other sponges), and might significantly differ, in terms of genomic content, from demosponges such as Amphimedon queenslandica. Demosponges themselves constitute a complex group that is also probably not monophyletic. In particular the homoscleromorphs do not clearly cluster with the other demosponges in a recent analysis , and sponges from this group possess eumetazoan-like ultrastructural features, including an epithelium characterized by closely apposed cells with an underlying basement membrane and regularly distributed cell-cell junctions . Interestingly, a report of a set of ESTs from a homoscleromorph species, Oscarella carmela, has been recently published . We have identified 4 bHLH genes in these ESTs, with two corresponding to the Myc and Hey families also found in Amphimedon queenslandica and two to the Mad and Emc families for which no member can be found in Amphimedon queenslandica (unpublished observations). This highlights the importance for broad sampling of sponges in future studies.
Keeping this in mind, we can nevertheless infer that the last common ancestor of Amphimedon queenslandica and eumetazoans, and thus probably that of all metazoans (Urmetazoa) possessed a minimal number of 10 to 14 bHLHs (Figure 5). Indeed, we have found that 10 bHLH orthologs do contain clear Amphimedon queenslandica members (one member in each family; additional file 1), indicating that these families are ancestral to metazoans. In addition, 5 Amphimedon queenslandica bHLHs can be related to several metazoan families of orthologs but cannot be confidently allocated to any of these families: one to the ARNT and Bmal families; one to three different families (Hif, Sim, and Trh); one to many group A families, which constitute the so-called Atonal and Twist superfamilies; and two to the ASCa and ASCb families (Figure 3; additional file 1). These data suggest that at least 4 additional bHLH families are ancestral to metazoans, raising the probable minimal number of bHLHs of Urmetazoa to 14 (Figure 5).
Noteworthy, there is a single bHLH in Amphimedon queenslandica (Amq1; additional file 15, Figure 3) that may represent the prototype of the numerous atonal and twist-related bHLHs found in eumetazoans, which are tissue-specific transcription factors involved in key developmental processes, such as neurogenesis and myogenesis (e.g. [4, 44] and references therein). Deciphering the expression of the Amphimedon queenslandica relative would be therefore of particular interest. More generally, the reduced number of bHLHs found in Amphimedon queenslandica and the fact that gene expression analyses can be performed during the development of this species , offers an altogether unique opportunity to study the different functions assumed by this entire gene superfamily in a non-eumetazoan species. Most of the bHLH genes that we have been identified from the genome traces have now been cloned from developmental cDNAs, indicating that they are expressed during embryogenesis, and the analysis of their expression patterns is currently being performed [G.R and B.M.D, unpublished data].
A general model for the diversification of the bHLH complement in metazoans
The model is summarized in Figure 5. Our phylogenetic analysis shows that there are no bHLH families of orthologs which are shared between fungi and metazoans, indicating an independent diversification of the bHLH complement in these two clades. The minimal number of bHLH genes in the last common ancestor of fungi and metazoans (Uropisthokonta) is therefore one. This single bHLH was most likely from higher-order B, as all of the fungi bHLHs are from this group (not shown). From our analysis of the bHLHs from Amphimedon queenslandica, we have inferred that the last common ancestor to all modern metazoans (i.e. Urmetazoa) possessed at least 10–14 different bHLHs, including members from 5 of the 6 higher-order groups of bHLHs. We can therefore infer that a first important expansion (both in terms of number and diversity) of the bHLH complement occurred after the divergence of fungal and metazoan lineages, but before the divergence of demosponges from the other metazoans. Although it is tempting to speculate that this expansion of the bHLH complement (and probably of many other gene families) has occurred in the early evolution of the metazoans and is related to the acquisition of multicellularity, we cannot rule out that it might have occurred even earlier. Indeed, fungi do not appear to be the sister-group to the Metazoa. There are other opisthokont eukaryotes more closely related, in particular the choanoflagellates . Choanoflagellates are eukaryotes that share some similarities with one of the most prominent cell-type of sponges (the choanocyte) and are usually considered as unicellular, although some of them display some sort of multicellularity/colonialism . It has been shown that choanoflagellates express relatives of a number of cell and adhesion protein families that have not previously been isolated from nonmetazoans, including cadherins, C-type lectins, several tyrosine kinases, and tyrosine kinase pathway components . It is therefore possible that at least part of the diversification of the bHLH complement might have occurred in a common ancestor of the choanoflagellates and the metazoans. The completion of the genome projects conducted on two choanoflagellate species, Monosiga brevicollis and Monosiga ovata, will help to address this issue.
From our analysis of the bHLHs from the cnidarians Nematostella vectensis and Hydra magnipapillata, we have inferred that the ureumetazoan genome contained at least 29–33 different bHLHs, with representatives from each of the 6 higher-order groups of bHLHs. A second phase of expansion of the bHLH complement has therefore occurred after the demosponge lineage split from other animals but before the divergence of cnidarians and bilaterians. As demosponges probably do not represent the sister-group of the cnidarian + bilaterian clade, it would be worthwhile to study the bHLH complement in organisms with an intermediate position between demosponges and eumetazoans, such as calcareous sponges and ctenophores [38–40]. This second expansion particularly concerned the group A bHLH families, which contain tissue-specific transcription factors that control the determination, the specification and the differentiation of many cell types in bilaterians. We therefore suggest that this second expansion of the bHLH complement can be correlated with the increase in the diversity of cell types that occurred before eumetazoan cladogenesis.
Finally, a third phase of expansion, which almost exclusively concerns the group A bHLHs, may have occurred after the divergence of cnidarians and other eumetazoans. Indeed, there are several bHLH families of orthologs that have bilaterian but not cnidarian members, and several related cnidarian bHLHs that cannot be allocated to any of the defined families (Figure 2). Although we cannot rule out the possibility that the cnidarian bHLHs may be derived members of the aforementioned bHLH families of orthologs, the existence of these cnidarian 'orphan' bHLHs may indicate an independent expansion of the bHLH complement in bilaterian and cnidarian lineages.
In this study, we identified the putative full set of bHLHs encoded by the newly sequenced genomes of 12 different species representative of the main metazoan evolutionary lineages, including three non bilaterian species, two cnidarians and a demosponge. Phylogenetic analysis of the sequence of the 695 identified bHLHs allowed us to conclude that (i) a first diversification of the bHLH superfamily has occurred in the pre-Cambrian and prior to metazoan cladogenesis; (ii) a second expansion of the bHLHs has occurred early in metazoan evolution before bilaterians and cnidarians diverged; and (iii) the bHLH complement during the evolution of the bilaterians has been remarkably stable. We suggest that these features may be extended to other developmental gene families and reflect a general trend in the evolution of the developmental gene repertoires of metazoans.
Retrieval of the bHLH sequences
As a starting point, we used the list of bHLH proteins we previously identified in Homo sapiens, Drosophila melanogaster, and Caenorhabditis elegans . In order to confirm this list, we performed similarity searches for each sequence using TBLASTN and BLASTP algorithms  on the current assembly and the predicted proteins of the three aforementioned genomes maintained by the National Center for Biotechnology Information (NCBI) . The whole set of Homo sapiens and Drosophila melanogaster bHLHs were then used for all the similarity searches using BLAST algorithm in the other species. The similarity searches using BLAST algorithm were performed at low stringency in order to obtain all possible bHLHs sequences, including divergent members relative to those of Homo sapiens and Drosophila melanogaster. All the retrieved sequences were used as queries to undertake similarity searches using BLAST algorithm against the NCBI protein database to ascertain that they corresponded to genuine bHLH domains.
In the case of Strongylocentrotus purpuratus and Tribolium castaneum, we used the BLAST server dedicated to these two species (on the web site of the Human Genome Sequencing Center, Baylor College of Medicine) . We performed similarity searches using both TBLASTN (on the WGS assembly and on the individual WGS reads) and BLASTP (on the protein predictions) algorithms to retrieve the maximum number of bHLHs encoded by these two genomes. For all the other studied metazoan species, no genome assembly was available. We therefore downloaded the publicly available shotgun traces through the Trace Archive v3.0 at the NCBI , made banks of the traces for each species, and performed similarity searches using TBLASTN algorithm against these banks, all of them being hosted by the Belgian EMBnet Node . We typically retrieved 6–7 overlapping traces that corresponded to the same aminoacid sequence. The nucleotide sequences of the traces were retrieved using the Ensembl Trace server  or from the banks hosted by the Belgian EMBnet node. Contigs were made for each set of overlapping traces by using CAP3  at the web site of the Pôle de Bioinformatique de Lyon (PBIL) . Aminoacid sequences were subsequently predicted using both Geneid  (through the web server at the Genome BioInformatics Research Lab, ), Genscan  (through the web server at the Institut Pasteur, ) and TBLASTN against the NCBI nr database. In some cases, only a part of the bHLH domain was retrieved. In these cases, the genomic sequences were extended through repeated steps of discontiguous Mega BLAST (at the NCBI) by using as query the left and right part of the previously assembled contig and then assembling new contigs with all the retrieved trace sequences (by using CAP3). The obtained contigs were analysed with Geneid, Genscan, and TBLASTN. This process was pursued until the whole bHLH domain was obtained or until it became impossible to further extend the contig. This workflow was successful in most cases. In the case of Amphimedon queenslandica, we systematically made larger gene assemblies using the same type of workflow. Preliminary phylogenetic analyses were then performed on the whole dataset for each species, and if some of the bHLH families of orthologs were found to lack a member in a given species, additional similarity searches using TBLASTN algorithm were performed specifically using members of the missing families as queries. This allowed us, in a very few cases, to retrieve additional bHLH domains missed in the first BLAST screen. Finally, each of the retrieved bHLH sequences were used to make similarity searches using BLASTP algorithm against the NCBI nr database (in order to see whether the corresponding gene may already have been cloned) and TBLASTN algorithm against the NCBI Expressed Sequence Tag (EST) database (in order to see whether there were ESTs corresponding to the identified bHLHs).
In the case of Nematostella vectensis, we made our initial similarity search using TBLASTN algorithm on the WGS reads. However, in a later phase of our work, genome assemblies became available on the web site of the US Department of Energy Joint Genome Institute  and on the StellaBase web site [60, 61]. We therefore made additional similarity searches using BLAST algorithm on these genome assemblies and we assessed the presence of our identified bHLHs in the genome assemblies and retrieved the corresponding transcript models using TBLASTN algorithm. In order to study the putative physical linkages between some of the cnidarian bHLH genes, we used the genome browser available on the web site of the DOE Joint genome institute . We also downloaded from the same source the nucleotide sequence of the corresponding scaffolds and used Geneid and Genscan [55, 57] to predict the size and position of the putative genes contained in these scaffolds. In the case of the searches for the fungi bHLHs, we used the bHLHs from Saccharomyces cerevisiae that we had previously identified  and a selection of human bHLHs as queries for similarity searches using TBLASTN and BLASTP algorithms on the fungi WGS assemblies and protein predictions, respectively, that are available at the NCBI.
All the nucleotide sequences (traces, contigs, ESTs) we have isolated are available upon request.
Multiple alignments were performed with Clustal W  using the ClustalW web server at the Bioinformatics Center of the Kyoto University  or using ClustalW at the Belgian EMBnet Node and they were subsequently manually improved. Handling of the multiple alignments was done using SEAVIEW  or GeneDoc . Unweighted maximum-parsimony (MP) and neighbour-joining (NJ) reconstructions were performed with the PAUP 4.0 program . NJ analyses were done using the BioNJ algorithm  and 10,000 bootstrap replicates. MP analyses were performed with the following settings: heuristic search of over 250 bootstrap replicates; MAXTREES set at 3000, and other parameters set at default values. Maximum likelihood (ML) analyses were performed with PHYML . PHYML analyses were performed using the Jones-Taylor-Thornton (JTT) amino-acid substitution model , the frequencies of amino acids being estimated from the data set, and rate heterogeneity across sites being modelled by two rate categories (one constant and eight γ-rates). The amino acid substitution model was chosen using ModelGenerator . Statistical support for the different internal branches was assessed by bootstrap resampling (150 bootstrap replicates), as implemented in PHYML . Bayesian inference was performed using the Markov chain Monte Carlo method as implemented in the MRBAYES (version 3) package [71, 72]. We used the JTT substitution frequency matrix  with among-sites rate variation modelled by means of a discrete γ distribution with four equally probable categories. Two independent Markov chains were run, each containing from 1,500,000 to 3,000,000 Monte Carlo steps (depending on the number of steps required to get chain convergence). One out of every 250 trees was saved. The trees obtained in the two runs were meshed and the first 25% of the trees were discarded as 'burnin'. Majority consensus of the obtained trees was computed by means of the PAUP 4.0 program. Marginal probabilities at each internal branche were taken as a measure of statistical support.
All the alignments and the trees are available upon request.
We are extremely grateful to the Department of Energy (DoE) Joint Genome Institute, the Baylor College of Medicine (BCM-HGSC), the J. Craig Venter Institute, and the National Human Genome Research Institute for sequencing the genomes of the different species used in this study and for making these sequences publicly available. We are also much grateful to the scientists who set up and led these projects, in particular D. Rokhsar, K. Remington, S. Richards, M. Martindale, and J. Finnerty. Our work has been supported by the CNRS and the Ministère Français de la Recherche through the ACI "Jeunes chercheurs et jeunes chercheuses" (to M.V.), and the Australian Research Council to B.M.D. V.L and D.C. were financially supported by the Belgian Science Policy and M.T.-C. by the Vrije Universiteit Brussel (Geconcerteerde Onderzoeksactie 29). E.S. holds a fellowship from the Ministère Français de la Recherche and P.K. a "Bourse pour Docteur-Ingénieur" from the CNRS. We thank Dr M. Adamski for providing assembly tools and expertise for the Amphimedon queenslandica genomic and EST sequences.
- Massari ME, Murre C: Helix-Loop-Helix proteins: regulators of transcription in eucaryotic organisms. Mol Cell Biol. 2000, 20: 429-440. 10.1128/MCB.20.2.429-440.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones S: An overview of the basic helix-loop-helix proteins. Genome Biol. 2004, 5: 226-10.1186/gb-2004-5-6-226.PubMed CentralView ArticlePubMedGoogle Scholar
- Pontig CP, Aravind L: PAS: a multifunctional domain family comes to light. Curr Biol. 1997, 7: R674-677.View ArticleGoogle Scholar
- Ledent V, Vervoort M: The basic Helix-loop-Helix protein family: comparative genomics and phylogenetic analysis. Genome Res. 2001, 11: 754-770. 10.1101/gr.177001.PubMed CentralView ArticlePubMedGoogle Scholar
- Amoutzias GD, Robertson DL, Oliver SG, Bornberg-Bauer E: Convergent evolution of gene networks by single-gene duplications in higher eukaryotes. EMBO Rep. 2004, 5: 274-279. 10.1038/sj.embor.7400096.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson KA, Lopez JM: Saccharomyces cerevisiae basic Helix-Loop-Helix proteins regulate diverse biological processes. Nucleic Acids Res. 2000, 28: 1499-1505. 10.1093/nar/28.7.1499.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu YZ, Hogenesch JB, Bradfield CA: The PAS superfamily: sensors of environmental and developmental signals. Annu Rev Pharmacol Toxicol. 2000, 40: 519-561. 10.1146/annurev.pharmtox.40.1.519.View ArticlePubMedGoogle Scholar
- Luscher B: Function and regulation of the transcription factors of the Myc/Max/Mad network. Gene. 2001, 277: 1-14. 10.1016/S0378-1119(01)00697-7.View ArticlePubMedGoogle Scholar
- Moore AW, Barbel S, Jan LY, Jan YN: A genomewide survey of basic helix-loop-helix factors in Drosophila. Proc Natl Acad Sci. 2000, 97: 10436-10441. 10.1073/pnas.170301897.PubMed CentralView ArticlePubMedGoogle Scholar
- Peyrefitte S, Kahn D, Haenlin M: New members of the Drosophila Myc transcription factor subfamily revealed by a genome-wide examination for basic Helix-Loop-Helix genes. Mech Dev. 2001, 104: 99-104. 10.1016/S0925-4773(01)00360-4.View ArticlePubMedGoogle Scholar
- Ledent V, Paquet O, Vervoort M: Phylogenetic analysis of the human basic helix-Loop-Helix proteins. Genome Biol. 2002, 3: 0030.1-0030.18. 10.1186/gb-2002-3-6-research0030.View ArticleGoogle Scholar
- Satou Y, Imai KS, Levine M, Kohara Y, Rokhsar D, Satoh N: A genomewide survey of developmentally relevant genes in Ciona intestinalis. I. Genes for bHLH transcription factors. Dev Genes Evol. 2003, 213: 213-221. 10.1007/s00427-003-0319-7.View ArticlePubMedGoogle Scholar
- Toledo-Ortiz G, Huq E, Quail PH: The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003, 5: 1749-1770. 10.1105/tpc.013839.View ArticleGoogle Scholar
- Heim MA, Jakoby M, Werber M, Martin C, Weisshaar B, Bailey PC: The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity. Mol Biol Evol. 2003, 20: 735-747. 10.1093/molbev/msg088.View ArticlePubMedGoogle Scholar
- Buck MJ, Atchley WR: Phylogenetic analysis of plant basic helix-loop-helix proteins. J Mol Evol. 2003, 56: 742-750. 10.1007/s00239-002-2449-3.View ArticlePubMedGoogle Scholar
- Li X, Duan X, Jiang H, Sun Y, Tang Y, Yuan Z, Guo J, Liang W, Chen L, Yin J, Ma H, Wang J, Zhang D: Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis. Plant Physiol. 2006, 141: 1167-1184. 10.1104/pp.106.080580.PubMed CentralView ArticlePubMedGoogle Scholar
- Atchley WR, Fitch WM: A natural classification of the basic helix-loop-helix class of transcription factors. Proc Natl Acad Sci USA. 1997, 94: 5172-5176. 10.1073/pnas.94.10.5172.PubMed CentralView ArticlePubMedGoogle Scholar
- Dubois L, Vincent A: The COE – Collier/Olf1/EBF – transcription factors: structural conservation and diversity of developmental functions. Mech Dev. 2001, 108: 3-12. 10.1016/S0925-4773(01)00486-5.View ArticlePubMedGoogle Scholar
- Fitch WM: Distinguishing homologous from analogous proteins. Syst Zool. 1970, 19: 99-113. 10.2307/2412448.View ArticlePubMedGoogle Scholar
- Technau U, Rudd S, Maxwell P, Gordon PM, Saina M, Grasso LC, Hayward DC, Sensen CW, Saint R, Holstein TW, Ball EE, Miller DJ: Maintenance of ancestral complexity and non-metazoan genes in two basal cnidarians. Trends Genet. 2005, 21: 633-639. 10.1016/j.tig.2005.09.007.View ArticlePubMedGoogle Scholar
- Chourrout D, Delsuc F, Chourrout P, Edvardsen RB, Rentzsch F, Renfer E, Jensen MF, Zhu B, de Jong P, Steele RE, Technau U: Minimal ProtoHox cluster inferred from bilaterian and cnidarian Hox complements. Nature. 2006, 442: 684-687. 10.1038/nature04863.View ArticlePubMedGoogle Scholar
- Ryan JF, Burton PM, Mazza ME, Kwong GK, Mullikin JC, Finnerty JR: The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: Evidence from the starlet sea anemone, Nematostella vectensis. Genome Biol. 2006, 7: R64-10.1186/gb-2006-7-7-r64.PubMed CentralView ArticlePubMedGoogle Scholar
- Simpson AG, Roger AJ: The real 'kingdoms' of eukaryotes. Curr Biol. 2004, 14: R693-696. 10.1016/j.cub.2004.08.038.View ArticlePubMedGoogle Scholar
- Peterson KJ, Lyons JB, Nowak KS, Takacs CM, Wargo MJ, McPeek MA: Estimating metazoan divergence times with a molecular clock. Proc Natl Acad Sci USA. 2004, 101: 6536-6541. 10.1073/pnas.0401670101.PubMed CentralView ArticlePubMedGoogle Scholar
- Douzery EJ, Snell EA, Bapteste E, Delsuc F, Philippe H: The timing of eukaryotic evolution: does a relaxed molecular clock reconcile proteins and fossils?. Proc Natl Acad Sci USA. 2004, 101: 15386-15391. 10.1073/pnas.0403984101.PubMed CentralView ArticlePubMedGoogle Scholar
- Panopoulou G, Poustka AJ: Timing and mechanism of ancient vertebrate genome duplications – the adventure of a hypothesis. Trends Genet. 2005, 21: 559-567. 10.1016/j.tig.2005.08.004.View ArticlePubMedGoogle Scholar
- Kusserow A, Pang K, Sturm C, Hrouda M, Lentfer J, Schmidt HA, Technau U, von Haeseler A, Hobmayer B, Martindale MQ, Holstein TW: Unexpected complexity of the Wnt gene family in a sea anemone. Nature. 2005, 433: 156-160. 10.1038/nature03158.View ArticlePubMedGoogle Scholar
- Kamm K, Schierwater B, Jakob W, Dellaporta SL, Miller DJ: Axial patterning and diversification in the cnidaria predate the Hox system. Curr Biol. 2006, 16: 920-926. 10.1016/j.cub.2006.03.036.View ArticlePubMedGoogle Scholar
- Magie CR, Pang K, Martindale MQ: Genomic inventory and expression of Sox and Fox genes in the cnidarian Nematostella vectensis. Dev Genes Evol. 2005, 215: 618-630. 10.1007/s00427-005-0022-y.View ArticlePubMedGoogle Scholar
- Matus DQ, Thomsen GH, Martindale MQ: Dorso/ventral genes are asymmetrically expressed and involved in germ-layer demarcation during cnidarian gastrulation. Curr Biol. 2006, 16: 499-505. 10.1016/j.cub.2006.01.052.View ArticlePubMedGoogle Scholar
- Matus DQ, Pang K, Marlow H, Dunn CW, Thomsen GH, Martindale MQ: Molecular evidence for deep evolutionary roots of bilaterality in animal development. Proc Natl Acad Sci USA. 2006, 103: 11195-11200. 10.1073/pnas.0601257103.PubMed CentralView ArticlePubMedGoogle Scholar
- Rentzsch F, Anton R, Saina M, Hammerschmidt M, Holstein TW, Technau U: Asymmetric expression of the BMP antagonists chordin and gremlin in the sea anemone Nematostella vectensis: Implications for the evolution of axial patterning. Dev Biol. 2006, 296: 375-387. 10.1016/j.ydbio.2006.06.003.View ArticlePubMedGoogle Scholar
- Ball EE, Hayward DC, Saint R, Miller DJ: A simple plan – cnidarians and the origins of developmental mechanisms. Nat Rev Genet. 2004, 5: 567-577. 10.1038/nrg1402.View ArticlePubMedGoogle Scholar
- Kamm K, Schierwater B: Ancient complexity of the non-Hox ANTP gene complement in the anthozoan Nematostella vectensis. Implications for the evolution of the ANTP superclass. J Exp Zoolog B Mol Dev Evol. 2006, 306: 589-596. 10.1002/jez.b.21123.View ArticleGoogle Scholar
- Martindale MQ, Pang K, Finnerty JR: Investigating the origins of triploblasty: 'mesodermal' gene expression in a diploblastic animal, the sea anemone Nematostella vectensis (phylum, Cnidaria; class, Anthozoa). Development. 2004, 131: 2463-2474. 10.1242/dev.01119.View ArticlePubMedGoogle Scholar
- Pang K, Matus DQ, Martindale MQ: The ancestral role of COE genes may have been in chemoreception: evidence from the development of the sea anemone, Nematostella vectensis (Phylum Cnidaria; Class Anthozoa). Dev Genes Evol. 2004, 214: 134-138. 10.1007/s00427-004-0383-7.View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1525-142X.2006.00086.x.View ArticlePubMedGoogle Scholar
- Cavalier-Smith T, Allsopp MT, Chao EE, Boury-Esnault N, Vacelet J: Sponge phylogeny, animal monophyly, and the origin of the nervous system: 18S rRNA evidence. Can J Zool. 1996, 74: 2031-2045.View ArticleGoogle Scholar
- Zrzavy J, Mihulka S, Kepka P, Bezdek A: Phylogeny of the Metazoa based on morphological and 18S ribosomal DNA evidence. Cladistics. 1998, 14: 249-285.View ArticleGoogle Scholar
- Borchiellini C, Manuel M, Alivon E, Boury-Esnault N, Vacelet J, Le Parco Y: Sponge paraphyly and the origin of Metazoa. J Evol Biol. 2001, 14: 171-179. 10.1046/j.1420-9101.2001.00244.x.View ArticleGoogle Scholar
- Borchiellini C, Chombard C, Manuel M, Alivon E, Vacelet J, Boury-Esnault N: Molecular phylogeny of Demospongiae: implications for classification and scenarios of character evolution. Mol Phylogenet Evol. 2004, 32: 823-837. 10.1016/j.ympev.2004.02.021.View ArticlePubMedGoogle Scholar
- Boury-Esnault N, Ereskovsky A, Bezac C, Tokina D: Larval development in the Homoscleromorpha (Porifera, Demospongiae). Invertebr Biol. 2003, 122: 187-202.View ArticleGoogle Scholar
- Nichols SA, Dirks W, Pearse JS, King N: Early evolution of animal cell signaling and adhesion genes. Proc Natl Acad Sci USA. 2006, 103: 12451-12456. 10.1073/pnas.0604065103.PubMed CentralView ArticlePubMedGoogle Scholar
- Hassan BA, Bellen HJ: Doing the MATH: is the mouse a good model for fly development?. Genes Dev. 2000, 14: 1852-65.PubMedGoogle Scholar
- King N: Choanoflagellates. Curr Biol. 2005, 15: R113-4. 10.1016/j.cub.2005.02.004.View ArticlePubMedGoogle Scholar
- King N: The unicellular ancestry of animal development. Dev Cell. 2004, 7: 313-325. 10.1016/j.devcel.2004.08.010.View ArticlePubMedGoogle Scholar
- King N, Hittinger CT, Carroll SB: Evolution of key cell signaling and adhesion protein families predates animal origins. Science. 2003, 301: 361-363. 10.1126/science.1083853.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSIBLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- The National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/]
- The Human Genome Sequencing Center, Baylor College of Medicine. [http://www.hgsc.bcm.tmc.edu/]
- The Belgian EMBnet Node. [http://www.be.embnet.org/]
- The Ensembl Trace server. [http://trace.ensembl.org/]
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle Scholar
- The Pôle de Bioinformatique de Lyon CAP3 server. [http://pbil.univ-lyon1.fr/cap3.php]
- Parra G, Blanco E, Guigo R: GeneID in Drosophila. Genome Res. 2000, 10: 511-515. 10.1101/gr.10.4.511.PubMed CentralView ArticlePubMedGoogle Scholar
- The Genome BioInformatics Research Lab. [http://www1.imim.es/geneid.html]
- Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.View ArticlePubMedGoogle Scholar
- The Institut Pasteur. [http://bioweb.pasteur.fr/seqanal/interfaces/genscan]
- The DOE Joint genome institute Nematostella vectensis home page. [http://genome.jgi-psf.org/Nemve1/Nemve1.home.html]
- Sullivan JC, Ryan JF, Watson JA, Webb J, Mullikin JC, Rokhsar D, Finnerty JR: StellaBase: the Nematostella vectensis Genomics Database. Nucleic Acids Res. 2006, D495-499. 10.1093/nar/gkj020. 34 DatabaseGoogle Scholar
- The StellaBase web site. [http://evodevo.bu.edu/stellabase/]
- Thompson JD, Higgins JD, Gibson TJ: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- The Bioinformatics Center of the Kyoto University. [http://align.genome.jp/]
- Galtier N, Gouy M, Gautier C: SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996, 12: 543-548.PubMedGoogle Scholar
- Nicholas KB, Nicholas HB, Deerfield DW: GeneDoc: Analysis and Visualization of Genetic Variation. EMBnet.NEWS. 1997, 4: 14-Google Scholar
- Swofford DL: PAUP: Phylogenetic Analysis Using Parsimony (and Other Methods), Version 4. 1998, Sunderland, MA: SinauerGoogle Scholar
- Gascuel O: BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997, 14: 685-695.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.PubMedGoogle Scholar
- Keane TM, Creevey CJ, Pentony MM, Naughton TJ, Mclnerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006, 6: 29-10.1186/1471-2148-6-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Adoutte A, Balavoine G, Lartillot N, Lespinet O, Prud'homme B, de Rosa R: The new animal phylogeny: reliability and implications. Proc Natl Acad Sci USA. 2000, 97: 4453-4456. 10.1073/pnas.97.9.4453.PubMed CentralView ArticlePubMedGoogle Scholar
- Delsuc F, Brinkmann H, Chourrout D, Philippe H: Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature. 2006, 439: 965-968. 10.1038/nature04336.View ArticlePubMedGoogle Scholar
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.