- Open Access
Surprising complexity of the ancestral apoptosis network
© Zmasek et al; licensee BioMed Central Ltd. 2007
- Received: 20 July 2007
- Accepted: 24 October 2007
- Published: 24 October 2007
Apoptosis, one of the main types of programmed cell death, is regulated and performed by a complex protein network. Studies in model organisms, mostly in the nematode Caenorhabditis elegans, identified a relatively simple apoptotic network consisting of only a few proteins. However, analysis of several recently sequenced invertebrate genomes, ranging from the cnidarian sea anemone Nematostella vectensis, representing one of the morphologically simplest metazoans, to the deuterostomes sea urchin and amphioxus, contradicts the current paradigm of a simple ancestral network that expanded in vertebrates.
Here we show that the apoptosome-forming CED-4/Apaf-1 protein, present in single copy in vertebrate, nematode, and insect genomes, had multiple paralogs in the cnidarian-bilaterian ancestor. Different members of this ancestral Apaf-1 family led to the extant proteins in nematodes/insects and in deuterostomes, explaining significant functional differences between proteins that until now were believed to be orthologous. Similarly, the evolution of the Bcl-2 and caspase protein families appears surprisingly complex and apparently included significant gene loss in nematodes and insects and expansions in deuterostomes.
The emerging picture of the evolution of the apoptosis network is one of a succession of lineage-specific expansions and losses, which combined with the limited number of 'apoptotic' protein families, resulted in apparent similarities between networks in different organisms that mask an underlying complex evolutionary history. Similar results are beginning to surface for other regulatory networks, contradicting the intuitive notion that regulatory networks evolved in a linear way, from simple to complex.
- Additional Data File
- Death Effector Domain
- Card Domain
- Apoptosis Network
- Apoptotic Network
Compared to C. elegans, the vertebrate apoptosis network is extensive, both in the number and in the size of the protein families involved. While C. elegans has one homolog of each (CED-3, CED-4, and CED-9), human has 12 CED-3 (caspase) homologs and 13 CED-9 homologs (Bcl-2-like proteins containing multiple BH motifs) as well as a number of highly divergent proteins that play an analogous role to the EGL-1 protein (BH3 motif only) (three additional caspase related genes, for which confirmation for a role in apoptosis is absent, have been found in C. elegans) [6–8]. All mammals, as well as birds, amphibians, and, to a lesser degree, fish, show somewhat similar expansions of these families . The CED-4/Apaf-1 family is an exception, being the only protein from the core of the apoptosis network that was not duplicated in any of the genomes studied until recently. Therefore, it was logical to expect that the role of this protein is indeed central and unique and that all homologs studied to date represent one-to-one orthologs that have evolved by speciation events only. Such one-to-one orthologs usually tend to display a high level of functional similarity and could be effectively used as functional models of each other . In this context, it was somewhat puzzling that an increasing body of experimental evidence suggested fundamental functional differences between C. elegans CED-4 and Drosophila Dark and their homologs in other species. In vertebrates, cytochrome c binds to Apaf-1 to trigger assembly of the apoptosome , which in turn leads to caspase activation. In contrast, no cytochrome c binding has been recognized for C. elegans CED-4 and remains controversial for Drosophila Dark [5, 11].
With the recent completion of three marine invertebrate genomes, namely two from Deuterostomia (the sea urchin Strongylocentrotus purpuratus and the amphioxus Branchiostoma floridae; unpublished; see Materials and methods) and one from Cnidaria (the sea anemone Nematostella vectensis), we are now able to obtain a more complete picture of how the complex vertebrate apoptosis network might have evolved and how representative the simple networks seen in insects and nematodes are of the systems present in other invertebrate animals [12–15].
Core apoptosis domains in several completed animal genomes
Death domain (DD)
Death effector domain (DED)
H. sapiens (human)
M. musculus (mouse)
C. familiaris (dog)
G. gallus (chicken)
X. tropicalis (western clawed frog)
B. rerio (zebrafish)
F. rubripes (Japanese pufferfish)
T. nigroviridis (green pufferfish)
B. floridae (amphioxus)
C. intestinalis (sea squirt)
S. purpuratus (purple sea urchin)
D. melanogaster (fruit fly)
N. vectensis (starlet sea anemone)
The presence of numerous CED-4/Apaf-1 homologs in the common ancestor of Bilateria and Cnidaria suggests that initially there might have been several mechanisms to activate the intrinsic apoptosis pathways and/or several downstream pathways activated by similar signals and that the mechanism of human Apaf-1 and its vertebrate orthologs presents only one of several possibilities. This also explains why the biochemical/structural mechanism of C. elegans CED-4 and Drosophila Dark can be significantly different from human Apaf-1 .
The functional variations among different branches of the Apaf-1 family are illustrated by their different domain organizations. Human Apaf-1 and its Nematostella, amphioxus, and sea urchin homologs exhibit the same or similar domain organization (CARD [two for Nematostella]-NB-ARC-WD40 repeats). Nematode and most, but not all, insect sequences seem to lack WD40 repeats , suggesting that the loss of the receptor domain of CED-4 is a (relatively) recent event, specific to nematode/insect Apaf-1 homologs. The expanded repertoire of CED-4/Apaf-1 homologs in sea urchin, amphioxus, and Nematostella contains proteins with novel domain combinations. This includes replacement of the single CARD domain at the amino terminus with pairs of CARD domains (Nematostella and amphioxus), death domains (amphioxus and, as previously described in , sea urchin), death effector domains (Nematostella), and TIR domains (amphioxus), all of which function as protein-protein interaction facilitators . At the carboxyl terminus, the WD40 repeats are occasionally missing, replaced by TPR repeats , or supplemented by double death domain repeats. Therefore, it seems that functional differences among CED-4/Apaf-1 homologs could include both the sensing mechanism (carboxy-terminal receptor domains) and the downstream recruitment function (amino-terminal protein-protein interaction domains). While we can only speculate on how such a rich set of domain combinations (as seen in amphioxus) came to be, a correlation between domain versatility and abundance has been observed . Interestingly, the TIR-NB-ARC domain architecture, present in one of the amphioxus proteins, resembles plant disease-resistant (R) genes involved in a process called hypersensitive response , which bears some similarity to apoptosis in animals , suggesting possibly even more distant evolutionary connections.
The final step in apoptosis is proteolysis of a variety of target proteins in the cell by 'effector' caspases, which are activated in a proteolytic cascade by several 'apical' ('initiator') caspases . Both types are clearly present in all animals (Additional data file 1). Yet, again, Nematostella, amphioxus, and sea urchin have representatives in more subtypes (defined by human caspases) than nematodes and insects.
It has been proposed that the invention of apoptosis was an essential requirement for the evolution of multicellular animals , and indeed it has been demonstrated that the apoptotic pathways involving members of the Bcl-2 family are present in the most basal metazoan phylum, the sponges (Porifera) [32, 33]. Our results suggest that the bilaterian-cnidarian ancestor living 650-1,000 million years ago already had an apoptotic regulatory network composed of Apaf-1, Bcl-2 and caspase family members. Surprisingly, this ancient apoptosis network appears to have been more complex than previously thought and the simple networks seen in present day insects and nematodes are the result of significant gene losses. Furthermore, a central protein in the classical apoptosis model, the apoptosome forming Apaf-1 , which exists as a single homolog in all genomes studied so far, has multiple homologs in several morphologically simple invertebrates and many extant Apaf-1 homologs may not be orthologous. This suggests that multiple mechanisms triggering apoptosis, as well as multiple downstream pathways implementing it, may have existed in early organisms. Many gene copy number differences are found that can be explained only by lineage-specific duplications and gene losses. Apparently, different organisms evolved unique apoptosis networks, which interestingly involved essentially the same gene families, hence sometimes providing an appearance of similarity between independently evolved networks. Interestingly, apoptosis regulators are not the only protein families involved in development and disease exhibiting surprising, almost vertebrate-like complexity in Cnidaria, and thus, presumably, the common cnidarian-bilaterian ancestor [34, 35]. Analyses of Nematostella Wnt genes revealed unforeseen ancestral diversity: Nematostella and bilaterians share at least eleven of the twelve known Wnt subfamilies, while five subfamilies appear to be lost in nematodes/insects . Similarly, proteins with innate immunity domains have been found to be expanded in Cnidaria . These results show that biological systems may not (always) evolve linearly from simple to complex. This urges caution in interpreting results from studies of C. elegans and D. melanogaster and indeed any model organisms for understanding apoptosis (or other regulatory pathways) in human. A more prudent approach might be to carefully select specific model systems for each protein family studied in such a way as to minimize the difference between the model and human. Such a selection process ideally should include phylogenetic analysis, thus reinforcing the view that "Nothing in biology makes sense except in the light of evolution." - Theodosius Dobzhansky (1900-1975).
Phylogenetic inference combined with domain composition analysis of Apaf-1, Bcl-2, and caspase proteins - central players in the apoptosis network - reveal a yet unpredicted ancestral complexity within each family. In particular, the relative simplicity of these regulatory networks observed in ecdysozoan species is not the result of a gradual increase in network complexity correlating with morphological complexity, but apparently the result of widespread gene losses. Our results emphasize the importance of explicit phylogenetic analysis covering a sufficiently large sample of species space, not only in the detection of orthologous sequences, but also in model organism selection and in the study of network evolution.
Sequence database searches
N. vectensis and B. floridae 1.0 genome assemblies and protein sets were downloaded from the Joint Genome Institute . The Strongylocentrotus purpuratus assembly Spur_v2.0 and GLEAN3 gene models were obtained from Baylor College of Medicine HGSC . The other genome sequences and corresponding protein sets were downloaded from Ensembl 38 or SWISS-PROT [40, 41]. Several rounds of PSI-TBLASTN searches were performed against each genome by using as seeds human NB-ARC, caspase, CARD, death, and death effector domains as well as Bcl-2 sequences from a variety of genomes . The hits were then mapped to the corresponding genome protein set to acquire the full-length protein sequences (for sea urchin and Nematostella, some of the gene models were in addition predicted by genscan) . All identified genes were checked by reciprocal BLAST analysis, Pfam 21.0 protein searches , Conserved Domain Search (CD-Search), and Reverse PSI-BLAST (RPS-BLAST) .
Multiple sequence alignments and phylogeny reconstructions
To ensure alignment of homologous domains, sequences were trimmed to one Pfam 21.0 model (NB-ARC, Bcl-2, Peptidase_C14 for the caspase domain) . Multiple sequence alignments were produced by PROBCONS 1.11 , MAFFT 5.861 (localpair, maxiterate 1000) , T-COFFEE 4.93 , and hmmalign from HMMER 2.3.2 [49, 50]. Multiple sequence alignment columns with a gap in more than 50% of sequences were deleted. MrBayes 3.1.2 was used with 10,000,000 generations, a sample frequency of 1,000, a mixture of amino-acid models with fixed rate matrices and equal rates, and 25% burn-in . For maximum likelihood approaches, PhyML 2.4.4 was used with the VT (variable time) model and four relative rate substitution categories [52, 53]. Pairwise distances (for the Neighbor Joining and Fitch-Margoliash methods from PHYLIP 3.66 [54–56], and FastME 1.1 ) were calculated by TREE-PUZZLE 5.2 using the VT model . Tree and domain composition diagrams were drawn using ATV 4a1 . All conclusions presented in this work are robust relative to the alignment methods, the alignment processing, the phylogeny reconstruction methods, and the parameters used. All sequence, alignment, and phylogeny files are available upon request.
Domain composition analysis
The following additional data files are available with the online version of this paper. Additional data file 1 is a figure illustrating the evolutionary history of caspase protein family members. Additional data file 2 is a table listing the CED-4/Apaf-1 protein family members used in this study. Additional data file 3 is a table listing the multi-motif Bcl-2 protein family members used in this study. Additional data file 4 is a table listing the caspase protein family members used in this study.
We thank Drs John C Reed, Guy S Salvesen, and Cheryl Bender for discussions and comments on the manuscript. This research was supported by NIH grants AI056324 and GM076221. N. vectensis, B. floridae, and Xenopus tropicalis genome data were produced by the US Department of Energy Joint Genome Institute . S. purpuratus genome data were produced by the Sea Urchin Genome Project at Baylor College of Medicine.
- Meier P, Finch A, Evan G: Apoptosis in development. Nature. 2000, 407: 796-801. 10.1038/35037734.PubMedView ArticleGoogle Scholar
- Opferman JT, Korsmeyer SJ: Apoptosis in the development and maintenance of the immune system. Nature Immunol. 2003, 4: 410-415. 10.1038/ni0503-410.View ArticleGoogle Scholar
- Yuan J, Horvitz HR: A first insight into the molecular mechanisms of apoptosis. Cell. 2004, 116: S53-S56. 10.1016/S0092-8674(04)00028-5.PubMedView ArticleGoogle Scholar
- Koonin EV, Aravind L: Origin and evolution of eukaryotic apoptosis: the bacterial connection. Cell Death Differ. 2002, 9: 394-404. 10.1038/sj.cdd.4400991.PubMedView ArticleGoogle Scholar
- Manoharan A, Kiefer T, Leist S, Schrader K, Urban C, Walter D, Maurer U, Borner C: Identification of a genuine mammalian homolog of nematodal CED-4: is the hunt over or do we need better guns?. Cell Death Differ. 2006, 13: 1310-1317. 10.1038/sj.cdd.4401948.PubMedView ArticleGoogle Scholar
- Adrain C, Brumatti G, Martin SJ: Apoptosomes: protease activation platforms to die from. Trends Biochem Sci. 2006, 31: 243-247. 10.1016/j.tibs.2006.03.004.PubMedView ArticleGoogle Scholar
- Shaham S: Identification of multiple Caenorhabditis elegans caspases and their potential roles in proteolytic cascades. J Biol Chem. 1998, 273: 35109-35117. 10.1074/jbc.273.52.35109.PubMedView ArticleGoogle Scholar
- Abraham MC, Shaham S: Death without caspases, caspases without death. Trends Cell Biol. 2004, 14: 184-193. 10.1016/j.tcb.2004.03.002.PubMedView ArticleGoogle Scholar
- Reed JC: Mechanisms of apoptosis. Am J Pathol. 2000, 157: 1415-1430.PubMedPubMed CentralView ArticleGoogle Scholar
- Eisen JA: Phylogenomics: improving functional predictions for uncharacterized genes by evolutionary analysis. Genome Res. 1998, 8: 163-167.PubMedView ArticleGoogle Scholar
- Kornbluth S, White K: Apoptosis in Drosophila: neither fish nor fowl (nor man, nor worm). J Cell Sci. 2005, 118: 1779-1787. 10.1242/jcs.02377.PubMedView ArticleGoogle Scholar
- Sea Urchin Genome Sequencing Consortium, Sodergren E, Weinstock GM, Davidson EH, Cameron RA, Gibbs RA, Angerer RC, Angerer LM, Arnone MI, Burgess DR, et al: The genome of the sea urchin Strongylocentrotus purpuratus. Science. 2006, 314: 941-952. 10.1126/science.1133609.PubMed CentralView ArticleGoogle Scholar
- Halanych KM: The new view of animal phylogeny. Annu Rev Ecol Evol Systematics. 2004, 35: 229-256. 10.1146/annurev.ecolsys.35.112202.130124.View ArticleGoogle Scholar
- Darling JA, Reitzel AR, Burton PM, Mazza ME, Ryan JF, Sullivan JC, Finnerty JR: Rising starlet: the starlet sea anemone, Nematostella vectensis. BioEssays. 2005, 27: 211-221. 10.1002/bies.20181.PubMedView ArticleGoogle Scholar
- Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, et al: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317: 86-94. 10.1126/science.1139158.PubMedView ArticleGoogle Scholar
- Aravind L, Dixit VM, Koonin EV: Apoptotic molecular machinery: vastly increased complexity in vertebrates revealed by genome comparisons. Science. 2001, 291: 1279-1284. 10.1126/science.291.5507.1279.PubMedView ArticleGoogle Scholar
- Dunn SR, Phillips WS, Spatafora JW, Green DR, Weis VM: Highly conserved caspase and Bcl-2 homologues from the sea anemone Aiptasia pallida: lower metazoans as models for the study of apoptosis evolution. J Mol Evol. 2006, 63: 95-107. 10.1007/s00239-005-0236-7.PubMedView ArticleGoogle Scholar
- Robertson AJ, Croce J, Carbonneau S, Voronina E, Miranda E, McClay DR, Coffman JA: The genomic underpinnings of apoptosis in Strongylocentrotus purpuratus. Dev Biol. 2006, 300: 321-334. 10.1016/j.ydbio.2006.08.053.PubMedView ArticleGoogle Scholar
- van der Biezen EA, Jones JDG: The NB-ARC domain: a novel signalling motif shared by plant resistance gene products and regulators of cell death in animals. Curr Biol. 1998, 8: R226-R228. 10.1016/S0960-9822(98)70145-9.PubMedView ArticleGoogle Scholar
- Inohara N, Chamaillard M, McDonald C, Nunez G: NOD-LRR PROTEINS: Role in host-microbial interactions and inflammatory disease. Annu Rev Biochem. 2005, 74: 355-383. 10.1146/annurev.biochem.74.082803.133347.PubMedView ArticleGoogle Scholar
- Neuwald AF, Aravind L, Spouge JL, Koonin EV: AAA+: a class of chaperone-like ATPases associated with the assembly, operation, and disassembly of protein complexes. Genome Res. 1999, 9: 27-43.PubMedGoogle Scholar
- Kufer TA, Fritz JH, Philpott DJ: NACHT-LRR proteins (NLRs) in bacterial infection and immunity. Trends Microbiol. 2005, 13: 381-388. 10.1016/j.tim.2005.06.004.PubMedView ArticleGoogle Scholar
- Smith TF, Gaitatzes C, Saxena K, Neer EJ: The WD repeat: a common architecture for diverse functions. Trends Biochem Sci. 1999, 24: 181-185. 10.1016/S0968-0004(99)01384-5.PubMedView ArticleGoogle Scholar
- Park HH, Lo Y-C, Lin S-C, Wang L, Yang JK, Wu H: The death domain superfamily in intracellular signaling of apoptosis and inflammation. Annu Rev Immunol. 2007, 25: 561-586. 10.1146/annurev.immunol.25.022106.141656.PubMedPubMed CentralView ArticleGoogle Scholar
- D'Andrea LD, Regan L: TPR proteins: the versatile helix. Trends Biochem Sci. 2003, 28: 655-662. 10.1016/j.tibs.2003.10.007.PubMedView ArticleGoogle Scholar
- Vogel C, Teichmann SA, Pereira-Leal J: The relationship between domain duplication and recombination. J Mol Biol. 2005, 346: 355-365. 10.1016/j.jmb.2004.11.050.PubMedView ArticleGoogle Scholar
- Dangl JL, Jones JDG: Plant pathogens and integrated defence responses to infection. Nature. 2001, 411: 826-833. 10.1038/35081161.PubMedView ArticleGoogle Scholar
- Lacomme C, Santa Cruz S: Bax-induced cell death in tobacco is similar to the hypersensitive response. Proc Natl Acad Sci USA. 1999, 96: 7956-7961. 10.1073/pnas.96.14.7956.PubMedPubMed CentralView ArticleGoogle Scholar
- Fuentes-Prior P, Salvesen GS: The protein structures that shape caspase activity, specificity, activation and inhibition. Biochem J. 2004, 384: 201-232. 10.1042/BJ20041142.PubMedPubMed CentralView ArticleGoogle Scholar
- Krylov DM, Wolf YI, Rogozin IB, Koonin EV: Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution. Genome Res. 2003, 13: 2229-2235. 10.1101/gr.1589103.PubMedPubMed CentralView ArticleGoogle Scholar
- Cikala M, Wilm B, Hobmayer E, Bottger A, David CN: Identification of caspases and apoptosis in the simple metazoan Hydra. Curr Biol. 1999, 9: 959-962. 10.1016/S0960-9822(99)80423-0.PubMedView ArticleGoogle Scholar
- Wiens M, Krasko A, Müller CI, Müller WEG: Molecular evolution of apoptotic pathways: cloning of key domains from sponges (Bcl-2 homology domains and death domains) and their phylogenetic relationships. J Mol Evol. 2000, 50: 520-531.PubMedGoogle Scholar
- Wiens M, Müller WEG: Cell death in Porifera: molecular players in the game of apoptotic cell death in living fossils. Canadian J Zool. 2006, 84: 307-321. 10.1139/z05-165.View ArticleGoogle Scholar
- Technau U, Rudd S, Maxwell P, Gordon PMK, Saina M, Grasso LC, Hayward DC, Sensen CW, Saint R, Holstein TW, et al: 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.PubMedView ArticleGoogle Scholar
- Kortschak RD, Samuel G, Saint R, Miller DJ: EST analysis of the Cnidarian Acropora millepora reveals extensive gene loss and rapid sequence divergence in the model invertebrates. Current Biol. 2003, 13: 2190-2195. 10.1016/j.cub.2003.11.030.View ArticleGoogle Scholar
- Kusserow A, Pang K, Sturm C, Hrouda M, Lentfer J, Schmidt HA, Technau U, von Haeseler A, Hobmayer B, Martindale MQ, et al: Unexpected complexity of the Wnt gene family in a sea anemone. Nature. 2005, 433: 156-160. 10.1038/nature03158.PubMedView ArticleGoogle Scholar
- Miller DJ, Hemmrich G, Ball EE, Hayward DC, Khalturin K, Funayama N, Agata K, Bosch TCG: The innate immune repertoire in Cnidaria - ancestral complexity and stochastic gene loss. Genome Biol. 2007, 8: R59-10.1186/gb-2007-8-4-r59.PubMedPubMed CentralView ArticleGoogle Scholar
- US Department of Energy Joint Genome Institute. [http://www.jgi.doe.gov/]
- Sea Urchin Genome Project. [http://www.hgsc.bcm.tmc.edu/projects/seaurchin/]
- Ensembl. [http://www.ensembl.org/]
- SWISS-PROT. [http://ca.expasy.org/sprot/]
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Burge CB, Karlin S: Finding the genes in genomic DNA. Curr Opin Struct Biol. 1998, 8: 346-354. 10.1016/S0959-440X(98)80069-9.PubMedView ArticleGoogle Scholar
- Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer ELL, et al: The Pfam protein families database. Nucleic Acids Res. 2004, 32: D138-141. 10.1093/nar/gkh121.PubMedPubMed CentralView ArticleGoogle Scholar
- Marchler-Bauer A, Bryant SH: CD-Search: protein domain annotations on the fly. Nucleic Acids Res. 2004, 32: W327-331. 10.1093/nar/gkh454.PubMedPubMed CentralView ArticleGoogle Scholar
- Do CB, Mahabhashyam MSP, Brudno M, Batzoglou S: ProbCons: Probabilistic consistency-based multiple sequence alignment. Genome Res. 2005, 15: 330-340. 10.1101/gr.2821705.PubMedPubMed CentralView ArticleGoogle Scholar
- Katoh K, Kuma K-i, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.PubMedPubMed CentralView ArticleGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-217. 10.1006/jmbi.2000.4042.PubMedView ArticleGoogle Scholar
- Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.PubMedView ArticleGoogle Scholar
- Nuin P, Wang Z, Tillier E: The accuracy of several multiple sequence alignment programs for proteins. BMC Bioinformatics. 2006, 7: 471-10.1186/1471-2105-7-471.PubMedPubMed CentralView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticleGoogle Scholar
- Muller T, Vingron M: Modeling amino acid replacement. J Computational Biol. 2000, 7: 761-776. 10.1089/10665270050514918.View ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Fitch WM, Margoliash E: Construction of phylogenetic trees. Science. 1967, 155: 279-284. 10.1126/science.155.3760.279.PubMedView ArticleGoogle Scholar
- Felsenstein J: PHYLIP - phylogeny inference package. Cladistics. 1989, 5: 164-166.Google Scholar
- Desper R, Gascuel O: Fast and accurate phylogeny reconstruction algorithms based on the minimum-evolution principle. J Computational Biol. 2002, 9: 687-705. 10.1089/106652702761034136.View ArticleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.PubMedView ArticleGoogle Scholar
- Zmasek CM, Eddy SR: ATV: display and manipulation of annotated phylogenetic trees. Bioinformatics. 2001, 17: 383-384. 10.1093/bioinformatics/17.4.383.PubMedView ArticleGoogle Scholar
- Jaroszewski L, Rychlewski L, Li Z, Li W, Godzik A: FFAS03: a server for profile-profile sequence alignments. Nucleic Acids Res. 2005, 33: W284-288. 10.1093/nar/gki418.PubMedPubMed CentralView ArticleGoogle Scholar
- Zdobnov EM, Apweiler R: InterProScan - an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.PubMedView ArticleGoogle 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.