Wikisource:WikiProject Open Access/Programmatic import from PubMed Central/Multigene Phylogeny of Choanozoa and the Origin of Animals
Animals are evolutionarily related to fungi and to the predominantly unicellular protozoan phylum Choanozoa, together known as opisthokonts. To establish the sequence of events when animals evolved from unicellular ancestors, and understand those key evolutionary transitions, we need to establish which choanozoans are most closely related to animals and also the evolutionary position of each choanozoan group within the opisthokont phylogenetic tree. Here we focus on Ministeria vibrans, a minute bacteria-eating cell with slender radiating tentacles. Single-gene trees suggested that it is either the closest unicellular relative of animals or else sister to choanoflagellates, traditionally considered likely animal ancestors. Sequencing thousands of Ministeria protein genes now reveals about 14 with domains of key significance for animal cell biology, including several previously unknown from deeply diverging Choanozoa, e.g. domains involved in hedgehog, Notch and tyrosine kinase signaling or cell adhesion (cadherin). Phylogenetic trees using 78 proteins show that Ministeria is not sister to animals or choanoflagellates (themselves sisters to animals), but to Capsaspora, another protozoan with thread-like (filose) tentacles. The Ministeria/Capsaspora clade (new class Filasterea) is sister to animals and choanoflagellates, these three groups forming a novel clade (filozoa) whose ancestor presumably evolved filose tentacles well before they aggregated as a periciliary collar in the choanoflagellate/sponge common ancestor. Our trees show ichthyosporean choanozoans as sisters to filozoa; a fusion between ubiquitin and ribosomal small subunit S30 protein genes unifies all holozoa (filozoa plus Ichthyosporea), being absent in earlier branching eukaryotes. Thus, several successive evolutionary innovations occurred among their unicellular closest relatives prior to the origin of the multicellular body-plan of animals.
Interpretation of molecular phylogeny, morphological and biochemical features suggest that the vast majority of the eukaryote diversity may belong to only six supergroups for review see ref . One of these supergroups, the opisthokonts, is composed of animals, fungi and several smaller groups of unicellular eukaryotes belonging to the phylum Choanozoa –. Recent molecular phylogenetic evidence indicates that animals and fungi evolved independently from different unicellular protozoan choanozoan ancestors –. Choanozoa include both naked phagotrophic protozoa, characterized by long threadlike (filose) cellular projections , –, often involved in feeding, and others with rigid cell walls, which are saprotrophs or parasites. Of these, nucleariid filose amoebae  are probably the closest relatives of fungi , . However, the closest relatives of animals were not previously firmly established. Traditionally, one choanozoan lineage, the choanoflagellates, was regarded as most closely related to animals –,  because some of their tentacles (with a rigid internal skeleton of bundled actin as in animal intestinal microvilli)  are aggregated as a collar surrounding the cilium ( = flagellum) in the choanocytes of sponges . A recent multigene analysis showed that choanoflagellates are more closely related to animals than are ichthyosporean Choanozoa . Single gene phylogenies of a broader diversity of Choanozoa raised the possibility that the recently discovered minute marine protozoan Ministeria vibrans might be closer still to animals , . However those phylogenies were not robust, so it was unclear whether Ministeria was sister to animals, or was instead a highly modified choanoflagellate or else (as we show here for the first time) an even deeper branching choanozoan that might give clues as to the nature of still earlier events in the evolution of the unicellular ancestors of animals. To clarify this important evolutionary question we constructed a cDNA library for M. vibrans from cultures with both aggregated and dispersed cells to favour the highest possible gene diversity and sequenced 4,700 randomly chosen clones. We have used 78 of the most conserved genes to calculate a more extensive and robust multigene phylogenetic tree than hitherto for opisthokonts (Fig. 1; Supporting Table S1). In addition, we have searched our cDNA sequences for evidence of domains and genes for key animal properties that may have originated in the unicellular ancestors of animals rather than during the origin of multicellularity itself, notably for signaling pathways and cell-adhesion both essential in multicellular animal development.
Results and Discussion
Phylogeny of Choanozoa places Ministeria as sister to Capsaspora
Maximum likelihood (ML) and Bayesian inferences (BI) of a multigene alignment composed of 30 taxa and 17,482 amino acid characters were congruent in showing Amoebozoa (the closest outgroup to opisthokonts) , fungi and animals as three distinct clades with maximal statistical support (Fig. 1). Choanozoa in all trees are divided into three distinct lineages with highly supported branching order: 1) choanoflagellates are sisters to animals, 2) the second clade, here designated Filasterea, comprises Ministeria and Capsaspora, and is robustly sister to choanoflagellates plus animals; 3) Ichthyosporea, comprising Sphaeroforma and Amoebidium, which are all walled parasites (unlike the phagotrophic choanoflagellates and Filasterea), are robustly sisters to Filasterea plus choanoflagellates and animals (for taxon and genes used in multigene data, see Table S1). Compared with recent evolutionary studies , , our results greatly improve the resolution and support for the choanozoan branches and show the novel clade Filasterea for the first time. The whole tree is well resolved, with every node having maximal Bayesian support, even among the previously hard to resolve lower fungi where the recent separation of Allomycetes  from other Chytridiomycetes is clearly evident. A fusion between ubiquitin and ribosomal small subunit S30 protein once thought animal-specific  is here identified in all included choanozoan lineages, thus unifying the holozoa (filozoa plus Ichthyosporea) . Our present analysis leaves undetermined the position of one key free-living choanozoan lineage (Corallochytrium, whose cell wall hints at a relationship to Ichthyosoporea, though some trees put it nearer choanoflagellates). No molecular data are available at all for two other putatively choanozoan groups: Fonticulida , social amoebae with flat cristae and filopodia that we accordingly place in Discicristoidia, and Aphelidida , algal parasites with flat cristae that we place in the parasitic Ichthyosporea. However, sequences might reveal either as a distinct lineage of key significance for eukaryote megaphylogeny.
Precursors of animal cell adhesion and signaling pathway components among Choanozoa
Polypeptides involved in cell signaling pathways and cell adhesion are essential components in embryogenesis and development of animal body plans –. The majority of the signaling pathways are present across the animal kingdom, even among the primitive sponges that have differentiated cells, epithelia and connective tissue but lack a nervous system or distinct organs –. Domains related to genes in the animal hedgehog, tyrosine kinase and Notch signaling pathways and many adhesion components have also been identified in choanoflagellates, suggesting that some protein domains involved in signaling and cell adhesion originated in the unicellular ancestors of animals , . Such signaling and cell adhesion genes are not always composed of domains unique for animals, but may include different combinations of animal novelties and ancestral domains ; domains located on the same gene among Eumetazoa can be divided on separate genes in Choanozoa and sponges , . It has been suggested that domain shuffling of signaling genes was important in the transition from unicellular eukaryotes to differentiated multicellular animals , . Thus, in order to identify precursors of genes and domains involved morphogenesis of animals, it is necessary to apply several approaches that can detect homologies even if sequences lack animal-specific domains –, . By using a combination of sequence similarity searches (i.e. BLASTx against NCBInr, NCBIest and Gene ontology databases) and structure analyses against Pfam, InterProScan and Conserved Domain Database, we have annotated several sequences from Ministeria and other Choanozoa as signaling and cell adhesion components (Table 1). As in previous studies ,  many of our gene sequences are incomplete, so in some cases we can only be sure that domains homologous with those in related animal genes are present in Ministeria (for details, see Supporting Table S2). Nonetheless, there is a remarkably high diversity of domains involved in animal early embryonic development – undoubtedly many more precursors of key animal genes could be found in Ministeria if a full genome were sequenced. Annotation of M. vibrans cDNA sequences reveals sequences homologous to domains in Notch receptor and ligand (e.g. Notch 1 and 3), hedgehog (Hint domain), and animal tyrosine kinase receptor (e.g. Ros1 protooncogene). Thus, the Ministeria data contain domains from the precursors of the hedgehog, tyrosine kinase and Notch pathways, congruent with sequences identified from choanoflagellates . In addition, many components involved in cell adhesion such as crumbs, cadherin, focal adhesion kinase and integrin beta have been identified (Table 1 and Table S2). All these have been identified as animal-specific in recent analyses , –, . Among other available choanozoan cDNA library sequences we identified similar components from choanoflagellates and Capsaspora see also refs.  and . Only a few homologues of the signaling and adhesion gene domains could be identified in Ichthyosporea. If a similar paucity is maintained when full genomes are available, this might suggest that most such genes took their modern recognizable form in the last common ancestor of filozoa, after it separated from Ichthyosporea. However, the possibility exists that as Ichthyosporea are parasites they lost some domains secondarily; similar data on the deepest branching Choanozoa, the nucleariid amoeba, are needed to test this. Clearly, some essential prerequisites for development of multicellular organisms and cell differentiation originated earlier than previously thought, distinctly before the origin of the animal kingdom and choanoflagellates.
- "Phylogenetic distribution of annotated animal-like signaling and cell adhesion components within Choanozoa.(10.1371/journal.pone.0002098.t001)"
Ichthyosporea Filasterea Choanoflagellata GO annotation Amoebium paraciticum Sphaeroforma arctica Capsaspora owczarzaki Ministeria vibrans Monosiga brevicollis Monosiga ovata Hedgehog signaling Hedgehog • • • Hyperplastic discs homolog • Receptor tyrosine kinase signaling Ros/insulin family • Fibroblast growth factor • • Ephrin type-A receptor • Ephrin type-B receptor • Non-receptor tyrosine kinase signaling Abl tyrosine kinase • Shark • Notch signaling Notch • • • Delta • ADAM 10 • Cell contact and adhesion proteins Integrin-beta • Crumbs • Caveolin-1 • Focal adhesion kinase 1 • • • • Cadherin • ADAMTS • Selectin • Tetraspanin • Sushi • SHC • • ECM molecules and receptors 40S ribosomal protein SA •
As Filasterea are all unicellular, these components must originally have functioned in single cells and were only later co-opted for new functional pathways in multicellular animals. The identified hedgehog relatives in choanoflagellates , , Ministeria and sponges  have different domain structures and may illustrate such modification , . In Eumetazoa, the hedgehog protein has two main domains, the amino-terminal, signaling domain of hedgehog and the Hint domain, possessing autocatalytic cleavage activity, releasing the signaling domain. Recently, it was found that the hedgehog domains exist in separate genes in the sponge Amphimedon queenslandica and the choanoflagellate M. brevicollis,  (Fig. 2). The Ministeria polypeptide has a structure only containing the Hint domain and the conserved cleavage motive. Thus, the typical hedgehog domain combination of higher animal genes seems to be absent in Choanozoa and sponges. This also implies that hedgehog domains predate animals and choanoflagellates, and that several successive innovations occurred in this gene family before and after the actual origin of these lineages , , , . The fusion of the two domains in the eumetazoan ancestor changed the transmembrane signaling protein into the diffusible hedgehog ligand of higher animals, enabling signaling over a distance of up to 30 cells , a prerequisite for regulation of tissues and organs.
Intriguingly, we could also identify domains partially related to Notch in Ichthyosporea and Filasterea, but the typical N-terminal domains characterized from animal homologues are missing from the sequenced transcripts; even sponges seem to lack many of these domains. If the other Notch domains are present in the genomes of deeply diverging Choanozoa, they may be localized on separate genes similar to what has recently been suggested from genome analysis of choanoflagellates . Our discovery of such a variety of domains that are possible precursors of key animal functions highlights the need for targeted and comparative studies of their functions in Choanozoa and sponges, as recently initiated for a choanoflagellate tyrosine signaling kinase .
Implications on the morphological evolution of early diverging opisthokonts
Fig. 3 summaries how several key innovations can be mapped onto the robust phylogeny for the major groups of Choanozoa provided here. Other molecular innovations will only become apparent with full genome sequences of each choanozoan lineage, but can eventually be mapped onto this phylogeny. As Fig. 3 indicates, the origin of slender filose projections, which apparently occurred during the origin of opisthokonts when their ancestor diverged from Amoebozoa (which typically have broad pseudopodia) was probably a key enabling innovation for opisthokont evolution. There is a distinct difference between the non-tapering tentacles of filozoa with their rigid core of bundled actin and the broader and more flexible tapering branched ‘filopodia’ of nucleariids . During the origin of fungi, which probably originally had branched rhizoids as in chytrids , the cell wall was probably laid down around such filopodia, which sometimes branch in nucleariids, unlike in Filasterea, thereby creating the ancestor of fungal rhizoids and hyphae. By contrast, in Ichthyosporea and Corallochytrium the ancestral filose projections of opisthokonts must have been lost when their cell walls evolved. However, filose tentacles were retained in filozoa for phagotrophic feeding, with a key role, as a filter-feeding collar of tentacles throughout the origin of sponges from choanoflagellate-like ancestors .
The grouping of Ministeria with Capsaspora was unexpected from previous sequence data , . As it is fully congruent with their morphological similarities, since Capsaspora also has long untapered filose tentacles with a microfilamentous skeleton , unlike Ichthyosporea, we establish a new class and separate families for them:
Taxonomy: class Filasterea Cavalier-Smith. Diagnosis: naked, unicellular, uninucleate aerobic protozoa, with Golgi dictyosome, flat mitochondrial cristae and very slender long non-tapering, projecting tentacles with an internal skeleton of microfilaments, but which are not organized into a periciliary collar as in choanoflagellates. Trophic phase without wall. Sole order Ministeriida Cavalier-Smith 1997. Family Ministeriidae Cavalier-Smith. Diagnosis: free-living marine filastereans with symmetric radiating tentacles and sometimes vestigial cilium; cysts unknown. Family Capsasporidae Cavalier-Smith. Diagnosis: animal symbionts with lateral tentacles and specialized feeding peduncle; walled resting cyst.
Materials and Methods
cDNA library construction and contig assembly
M. vibrans (ATCC 50519) cells were grown in ATCC medium 1525 at 17°C. Cells were scraped and harvested by centrifugation and flash frozen in liquid nitrogen. Frozen samples were shipped on dry ice to Agencourt Bioscience Corporation (Beverly MA, USA). The mRNA was isolated and treated with a oligo(dT) -primed procedure for first strand cDNA synthesis. Double strand DNA fragments were directionally ligated into Agencourt pAGEN-1 vector. About 4700 randomly picked clones were sequenced from the 5′-end with the average read length above 600 bp. Sequence contigs were constructed by a Phred/phrap pipeline at the Bioportal service at University of Oslo ([http:www.bioportal.uio.no http:www.bioportal.uio.no]).
Single- and multigene phylogenetic analyses
Resulting contigs and singletons were screened for sequences usable for phylogeny by BLASTx searches against NCBI nr databases and independent BLASTx searches against manually curated single gene alignments used previously . Sequences from Ministeria vibrans and publicly available sequences from the nr- and dbEST databases (http://www.ncbi.nlm.nih.gov/blast/blast_databases.shtml) were retrieved from BLAST results and aligned to a previously published alignment  using the Mafft program . For each single gene alignment orthologous gene copies were selected on the basis of phylogenetic relationship in maximum likelihood trees and bootstrap values (i.e. with 100 pseudoreplicates and one tree search for each) inferred with PhyML . Taxa with several almost identical sequences, only the sequence displaying the shortest branch length on the tree was used in subsequent analyses. Among all constructed single gene alignments, we selected 78 genes that had the highest possible coverage of Ministeria sequences and that contained at least one sequences from one of choanozoan (i.e. Amoebidium parasiticum, Sphaeroforma arctica, Capsaspora owczarzaki, Ministeria vibrans, Proterospongia sp. and Monosiga brevicollis) lineages. Taxon-sampling were chosen to reflect the phylogenetic range within the opisthokonts (Amoebozoa used as outgroup). All ambiguously aligned sites were deleted before single- and multigene phylogenetic reconstructions. The final concatenated multi-gene matrix included 30 taxa and 17.482 amino acid characters. Each of the included taxa has maximum 70% missing characters (for further information about the taxa and included genes, see supporting Table 1).
Maximum likelihood phylogeny of the concatenated data was inferred with RAxML MPI version 2.2.3 . The rtREV+F evolutionary model was preferred by the ProtTest program . Among site-rate variation were accounted for by using the PROTMIX approximation (initial CAT rate categories set to 25 and final optimization with 4 gamma shape categories) and PROTGAMMA (gamma shaped distribution of site-rates with 4 rate categories) in two separate tree searches. Tree searches were done with 100 randomly generated starting trees. The branching order was identical in both trees inferred with PROTMIX and PROTGAMMA. Furthermore, because of the highly congruent branch pattern from these two runs, the bootstrap analysis was done with the most efficient PROTCAT method (as recommended recently ) and rtREV+F on 100 pseudoreplicates and one random starting tree for each replicate.
Bayesian inference was done with PhyloBayes version 2.3 . The CAT evolutionary model was used in together with a gamma distributed across-site variation (4 discrete rate categories). The evolution of the log-likelihood as a function of time was used to estimate if the 2 parallel chains had reached the stationary-state. This was then used to set the burn-in and compare the frequency of the bipartitions between several independent runs. The largest discrepancy (maxdiff) between the bipartitions was less than 0.1 and therefore we considered the Markov chain Monte Carlo chains to have converged. The tree and posterior probability values presented in Fig. 1 are a consensus of the cold chains from the 2 independent runs. All phylogenetic analyses were performed on the freely available Bioportal at Universtity of Oslo (http://www.bioportal.uio.no).
Identification of precursors of animal signaling and cell adhesion domains among Choanozoa
All contigs and singletons generated from the M. vibrans cDNA library and other available Choanozoa sequences were annotated in order to identify homologs of previously defined signaling and cell adhesion protein components that are involved in morphogenesis of animals and only identified in the animal kingdom and the choanoflagellates , –, . Thus, although components involved in for instance tyrosine kinase signaling pathways have been identified outside the Animalia, we only searched for types of components that have been suggested as animal-specific , –, . First, the sequences were screened against the NCBInr protein databases using BLASTx. All query sequences that retrieved relevant genes were further annotated by searches against the Gene Ontology (GO) database (http://amigo.geneontology.org/cgi-bin/amigo/go.cgi). A similar procedure was recently applied for sponge cDNA sequences . For each search only the GO descriptions with highest e-values were considered and assembled in a local database (Table 1 and Table S2). These sequences were subsequently searched against several motif and domain databases (Fig. 2) including InterProScan (http://www.ebi.ac.uk/InterProScan/) and Conserved Domain database (http://www.ncbi.nlm.nih.gov/sites/entrezdbcdd). Furthermore, in order to uncover protein domains in other genes than those identified in Table 1, single domains from human Hedgehog, Notch and Delta genes were used as query in BLAST searches (i.e. BLASTp and tBLASTn) against the NCBI databases (i.e. NCBInr and NCBIest) and local databases comprising the Ministeria vibrans cDNA and sequences downloaded from the Protist EST Project (PEP; http://www.bch.umontreal.ca/pepdb/pep_main.html).
(0.96 MB PDF) (data file pone.0002098.s001.pdf)
(0.08 MB PDF) (data file pone.0002098.s002.pdf) We thank Ema E. Chao and Hans Erik Karlsen for assistance with cell culturing and infrastructure, and the Bioportal at University of Oslo for computational resources. The Protist EST Project generated sequence data from some of the species included in the phylogenetic analyses.Competing Interests: The authors have declared that no competing interests exist.Funding: TC-S thanks NERC and the Canadian Institute for Advanced Research Evolutionary Biology Program for fellowship support and NERC for research grants, and KSJ the Norwegian Research Council for a research project and scholarships to KS, MAM and ME.
- Keeling, PJ; Burger, G; Durnford, DG; Lang, BF; Lee, RW; Pearlman, RE; Roger, AJ & Gray, MW (2005), “The tree of eukaryotes.”, Trends Ecol Evol 20: 670–676, 16701456, pmid:16701456
- Cavalier-Smith, T; Rayner, ADM; Brasier, CM & Moore, D (1987), “The origin of Fungi and pseudofungi.”, Evolutionary biology of the Fungi. (Cambridge University Press): 339–353
- Ruiz-Trillo, I; Burger, G; Holland, PW; King, N & Lang, BF (2007), “The origins of multicellularity: a multi-taxon genome initiative.”, Trends Genet 23: 113–118, 17275133, pmid:17275133
- Cavalier-Smith, T; Coombs, GxGH; Vickerman, K; Sleigh, MA & Warren, A (1998), “Neomonada and the origin of animals and fungi.”, Evolutionary relationships among Protozoa. (London: Chapman and Hall): 375–407
- Tong, SM (1997), “Heterotrophic flagellates and other protists from Southampton Water, UK.”, Ophelia 47: 71–131
- Cavalier-Smith, T (2003), “Protist phylogeny and the high-level classification of Protozoa.”, Eur J Protistol 39: 338–348
- Dykova, I; Veverkova, M; Fiala, I; Machackova, B & Peckova, H (2003), “Nuclearia pattersoni sp. n. (Filosea), a new species of amphizoic amoeba isolated from gills of roach (Rutilus rutilus), and its rickettsial endosymbiont.”, Folia Parasitol (Praha) 50: 161–170, 14535340, pmid:14535340
- Ruiz-Trillo, I; Lane, CE; Archibald, JM & Roger, AJ (2006), “Insights into the evolutionary origin and genome architecture of the unicellular opisthokonts Capsaspora owczarzaki and Sphaeroforma arctica.”, J Eukaryot Microbiol 53: 379–384, 16968456, pmid:16968456
- Medina, M; Collins, AG; Taylor, JW; Valentine, JW & Lipps, JH (2003), “Phylogeny of Opistokonta and the evolution of multicellularity and complexity in Fungi and Metazoa.”, Int J of Astrobiol 2(3): 203
- Cavalier-Smith, T & Chao, EE (2003), “Phylogeny of Choanozoa, Apusozoa, and other protozoa and early eukaryote megaevolution.”, J Mol Evol 56: 540–563, 12698292, pmid:12698292
- Jimenez-Guri, E; Philippe, H; Okamura, B & Holland, PW (2007), “Buddenbrockia is a cnidarian worm.”, Science 317: 116–118, 17615357, pmid:17615357
- Karpov, SA & Leadbeater, BS (1998), “Cytoskeleton structure and composition in choanoflagellates.”, J Euk Microbiol 45: 361–367
- James-Clark, H (1868), “On the Spongiae Ciliatae as Infusoria Flagellata: or observations on the structure, animality, and relationship of Leucosolenia botryoides, Bowerbank.”, Memoirs Boston Soc Nat Hist 1: 305–340 plus plates 9, 10
- Ruiz-Trillo, I; Roger, AG; Burger, G; Gray, MW & Lang, BF (2008), “A Phylogenomic Investigation into the Origin of Metazoa.”, Mol Biol Evol 25: 664–672, 18184723 doi:10.1093/molbev/msn006, pmid:18184723
- Steenkamp, ET; Wright, J & Baldauf, SL (2006), “The protistan origins of animals and fungi.”, Mol Biol Evol 23: 93–106, 16151185, pmid:16151185
- Richards, TA & Cavalier-Smith, T (2005), “Myosin domain evolution and the primary divergence of eukaryotes.”, Nature 436: 1113–1118, 16121172, pmid:16121172
- King, K; Westbrook, MJ; Young, SL; Kuo, A & Abedin, M (2008), “The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans.”, Nature 451: 783–788, 18273011, pmid:18273011
- Cavalier-Smith, T (1998), “A revised six-kingdom system of life.”, Biol Rev Camb Philos Soc 73: 203–266, 9809012, pmid:9809012
- Baker, RT; Williamson, NA & Wettenhall, RE (1996), “The yeast homolog of mammalian ribosomal protein S30 is expressed from a duplicated gene without a ubiquitin-like protein fusion sequence. Evolutionary implications.”, J Biol Chem 271: 13549–13555, 8662789, pmid:8662789
- Lang, BF; O'Kelly, C; Nerad, T; Gray, MW & Burger, G (2002), “The closest unicellular relatives of animals.”, Curr Biol 12: 1773–1778, 12401173, pmid:12401173
- Gromov, BV (2000), “Algal parasites of the genera Aphelidium, Amoeboaphelidium and Pseudoaphelidium, from the Cienkovski's “Monadea” group of as representatives of a new class.”, Zool Zh 79: 517–525 (in Russian)
- Pires-daSilva, A & Sommer, RJ (2003), “The evolution of signalling pathways in animal development.”, Nat Rev Genet 4: 39–49, 12509752, pmid:12509752
- Tepass, U; Godt, D & Winklbauer, R (2002), “Cell sorting in animal development: signalling and adhesive mechanisms in the formation of tissue boundaries.”, Curr Opin Genet Dev 12: 572–582, 12200163, pmid:12200163
- Nichols, SA; Dirks, W; Pearse, JS & King, N (2006), “Early evolution of animal cell signaling and adhesion genes.”, Proc Natl Acad Sci U S A 103: 12451–12456, 16891419, pmid:16891419
- Putnam, NH; Srivastava, M; Hellsten, U; Dirks, B & Chapman, J (2007), “Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization.”, Science 317: 86–94, 17615350, pmid:17615350
- Adamska, M; Matus, DQ; Adamski, M; Green, K & Rokhsar, DS (2007), “The evolutionary origin of hedgehog proteins.”, Curr Biol 17: R836–R837, 17925209, pmid:17925209
- King, N; Hittinger, CT & Carroll, SB (2003), “Evolution of key cell signaling and adhesion protein families predates animal origins.”, Science 301: 361–363, 12869759, pmid:12869759
- Snell, EA; Brooke, NM; Taylor, WR; Casane, D & Philippe, H (2006), “An unusual choanoflagellate protein released by Hedgehog autocatalytic processing.”, Proc Biol Sci 273: 401–407, 16615205, pmid:16615205
- Bayly, RD; Ngo, M; Aglyamova, GV & Agarwala, S (2007), “Regulation of ventral midbrain patterning by Hedgehog signaling.”, Development 134(11): 2115–2124, 17507412, pmid:17507412
- Segawa, Y; Suga, H; Iwabe, N; Oneyama, C & Akagi, T (2006), “Functional development of Src tyrosine kinases during evolution from a unicellular ancestor to multicellular animals.”, Proc Natl Acad Sci U S A 103: 12021–12026, 16873552, pmid:16873552
- Cavalier-Smith, T & McLaughlin, DJ (2001), “What are Fungi?.”, The Mycota VII Part A: Systematics and Evolution. (Berlin, Heidelberg: Springer-Verlag): 3–37
- Hertel, LA; Bayne, CJ & Loker, ES (2002), “The symbiont Capsaspora owczarzaki, nov. gen. nov. sp., isolated from three strains of the pulmonate snail Biomphalaria glabrata is related to members of the Mesomycetozoea.”, Int J Parasitol 32: 1183–1191, 12117501, pmid:12117501
- Rodriguez-Ezpeleta, N; Brinkmann, H; Burey, SC; Roure, B & Burger, G (2005), “Monophyly of primary photosynthetic eukaryotes: green plants, red algae, and glaucophytes.”, Curr Biol 15: 1325–1330, 16051178, pmid:16051178
- Katoh, K; Kuma, K; Toh, H & Miyata, T (2005), “MAFFT version 5: improvement in accuracy of multiple sequence alignment.”, Nucleic Acids Res 33: 511–518, 15661851, pmid:15661851
- Guindon, S & Gascuel, O (2003), “A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.”, Syst Biol 52: 696–704, 14530136, pmid:14530136
- Stamatakis, A (2006), “RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.”, Bioinformatics 22: 2688–2690, 16928733, pmid:16928733
- Abascal, F; Zardoya, R & Posada, D (2005), “ProtTest: Selection of best-fit models of protein evolution.”, Bioinformatics 21: 2104–2105, 15647292, pmid:15647292
- Lartillot, N & Philippe, H (2004), “A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process.”, Mol Biol Evol 21: 1095–1109, 15014145, pmid:15014145