Systematic searches for molecular synapomorphies in model metazoan genomes give some support for Ecdysozoa after accounting for the idiosyncrasies of Caenorhabditis elegans

The Genome of the Softshell Clam Mya arenaria and the Evolution of Apoptosis

Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

The origin of the animals and a 'Savannah' hypothesis for early bilaterian evolution

The earliest evolution of the animals remains a taxing biological problem, as all extant clades are highly derived and the fossil record is not usually considered to be helpful. The rise of the bilaterian animals recorded in the fossil record, commonly known as the 'Cambrian explosion', is one of the most significant moments in evolutionary history, and was an event that transformed first marine and then terrestrial environments. We review the phylogeny of early animals and other opisthokonts, and the affinities of the earliest large complex fossils, the so-called 'Ediacaran' taxa. We conclude, based on a variety of lines of evidence, that their affinities most likely lie in various stem groups to large metazoan groupings; a new grouping, the Apoikozoa, is erected to encompass Metazoa and Choanoflagellata. The earliest reasonable fossil evidence for total-group bilaterians comes from undisputed complex trace fossils that are younger than about 560 Ma, and these diversify greatly as the Ediacaran-Cambrian boundary is crossed a few million years later. It is generally considered that as the bilaterians diversified after this time, their burrowing behaviour destroyed the cyanobacterial mat-dominated substrates that the enigmatic Ediacaran taxa were associated with, the so-called 'Cambrian substrate revolution', leading to the loss of almost all Ediacara-aspect diversity in the Cambrian. Why, though, did the energetically expensive and functionally complex burrowing mode of life so typical of later bilaterians arise? Here we propose a much more positive relationship between late-Ediacaran ecologies and the rise of the bilaterians, with the largely static Ediacaran taxa acting as points of concentration of organic matter both above and below the sediment surface. The breaking of the uniformity of organic carbon availability would have signalled a decisive shift away from the essentially static and monotonous earlier Ediacaran world into the dynamic and burrowing world of the Cambrian. The Ediacaran biota thus played an enabling role in bilaterian evolution similar to that proposed for the Savannah environment for human evolution and bipedality. Rather than being obliterated by the rise of the bilaterians, the subtle remnants of Ediacara-style taxa within the Cambrian suggest that they remained significant components of Phanerozoic communities, even though at some point their enabling role for bilaterian evolution was presumably taken over by bilaterians or other metazoans. Bilaterian evolution was thus an essentially benthic event that only later impacted the planktonic environment and the style of organic export to the sea floor.

A transcriptome approach to ecdysozoan phylogeny

The monophyly of Ecdysozoa, which comprise molting phyla, has received strong support from several lines of evidence. However, the internal relationships of Ecdysozoa are still contended. We generated expressed sequence tags from a priapulid (penis worm), a kinorhynch (mud dragon), a tardigrade (water bear) and five chelicerate taxa by 454 transcriptome sequencing. A multigene alignment was assembled from 63 taxa, which comprised after matrix optimization 24,249 amino acid positions with high data density (2.6% gaps, 19.1% missing data). Phylogenetic analyses employing various models support the monophyly of Ecdysozoa. A clade combining Priapulida and Kinorhyncha (i.e. Scalidophora) was recovered as the earliest branch among Ecdysozoa. We conclude that Cycloneuralia, a taxon erected to combine Priapulida, Kinorhyncha and Nematoda (and others), are paraphyletic. Rather Arthropoda (including Onychophora) are allied with Nematoda and Tardigrada. Within Arthropoda, we found strong support for most clades, including monophyletic Mandibulata and Pancrustacea. The phylogeny within the Euchelicerata remained largely unresolved. There is conflicting evidence on the position of tardigrades: While Bayesian and maximum likelihood analyses of only slowly evolving genes recovered Tardigrada as a sister group to Arthropoda, analyses of the full data set, and of subsets containing genes evolving at fast and intermediate rates identified a clade of Tardigrada and Nematoda. Notably, the latter topology is also supported by the analyses of indel patterns.

The C. elegans rab family: identification, classification and toolkit construction

Rab monomeric GTPases regulate specific aspects of vesicle transport in eukaryotes including coat recruitment, uncoating, fission, motility, target selection and fusion. Moreover, individual Rab proteins function at specific sites within the cell, for example the ER, golgi and early endosome. Importantly, the localization and function of individual Rab subfamily members are often conserved underscoring the significant contributions that model organisms such as Caenorhabditis elegans can make towards a better understanding of human disease caused by Rab and vesicle trafficking malfunction. With this in mind, a bioinformatics approach was first taken to identify and classify the complete C. elegans Rab family placing individual Rabs into specific subfamilies based on molecular phylogenetics. For genes that were difficult to classify by sequence similarity alone, we did a comparative analysis of intron position among specific subfamilies from yeast to humans. This two-pronged approach allowed the classification of 30 out of 31 C. elegans Rab proteins identified here including Rab31/Rab50, a likely member of the last eukaryotic common ancestor (LECA). Second, a molecular toolset was created to facilitate research on biological processes that involve Rab proteins. Specifically, we used Gateway-compatible C. elegans ORFeome clones as starting material to create 44 full-length, sequence-verified, dominant-negative (DN) and constitutive active (CA) rab open reading frames (ORFs). Development of this toolset provided independent research projects for students enrolled in a research-based molecular techniques course at California State University, East Bay (CSUEB).

Deep genomic-scale analyses of the metazoa reject Coelomata: evidence from single- and multigene families analyzed under a supertree and supermatrix paradigm

Solving the phylogeny of the animals with bilateral symmetry has proven difficult. Morphological studies have suggested a variety of alternative hypotheses, of which, Hyman’s Coelomata hypothesis has become the most established. Studies based on 18S rRNA have failed to endorse Coelomata, supporting instead the rearrangement of the protostomes into two new clades: the Lophotrochozoa (including, e.g., the molluscs and the annelids) and the Ecdysozoa (including the Panarthropoda and most pseudocoelomates, such as the nematodes and priapulids). Support for this new animal phylogeny has been attained from expressed sequence tag studies, although these generally have a limited gene sampling. In contrast, deep genomic-scale analyses have often supported Coelomata. However, these studies are problematic due to their limited taxonomic sampling, which could exacerbate tree reconstruction artifacts.

Here, we address both of these sampling limitations; we study the effect of long-branch attraction (LBA) in deep genomic-scale analyses and provide convincing evidence, using both single- and multigene families, that Coelomata is an artifact. We show that optimal outgroup selection is key in avoiding LBA and identify the use of inadequate outgroups as the reason previous deep genomic-scale analyses found strong support for Coelomata.

From parasite genomes to one healthy world: Are we having fun yet?

In 1990, the Human Genome Sequencing Project was established. This laid the ground work for an explosion of sequence data that has since followed. As a result of this effort, the first complete genome of an animal, Caenorhabditis elegans was published in 1998. The sequence of Drosophila melanogaster was made available in March, 2000 and in the following year, working drafts of the human genome were generated with the completed sequence (92%) being released in 2003. Recent advancements and next-generation technologies have made sequencing common place and have infiltrated every aspect of biological research, including parasitology. To date, sequencing of 32 apicomplexa and 24 nematode genomes are either in progress or near completion, and over 600k nematode EST and 200k apicomplexa EST submissions fill the databases. However, the winds have shifted and efforts are now refocusing on how best to store, mine and apply these data to problem solving. Herein we tend not to summarize existing X-omics datasets or present new technological advances that promise future benefits. Rather, the information to follow condenses up-to-date-applications of existing technologies to problem solving as it relates to parasite research. Advancements in non-parasite systems are also presented with the proviso that applications to parasite research are in the making.

The origin and evolution of arthropods

The past two decades have witnessed profound changes in our understanding of the evolution of arthropods. Many of these insights derive from the adoption of molecular methods by systematists and developmental biologists, prompting a radical reordering of the relationships among extant arthropod classes and their closest non-arthropod relatives, and shedding light on the developmental basis for the origins of key characteristics. A complementary source of data is the discovery of fossils from several spectacular Cambrian faunas. These fossils form well-characterized groupings, making the broad pattern of Cambrian arthropod systematics increasingly consensual.

The evolution of the Ecdysozoa

Ecdysozoa is a clade composed of eight phyla: the arthropods, tardigrades and onychophorans that share segmentation and appendages and the nematodes, nematomorphs, priapulids, kinorhynchs and loriciferans, which are worms with an anterior proboscis or introvert. Ecdysozoa contains the vast majority of animal species and there is a great diversity of body plans among both living and fossil members. The monophyly of the clade has been called into question by some workers based on analyses of whole genome datasets. We review the evidence that now conclusively supports the unique origin of these phyla. Relationships within Ecdysozoa are also controversial and we discuss the molecular and morphological evidence for a number of monophyletic groups within this superphylum.

Homoplasy in genome-wide analysis of rare amino acid replacements: the molecular-evolutionary basis for Vavilov's law of homologous series

BACKGROUND

Rare genomic changes (RGCs) that are thought to comprise derived shared characters of individual clades are becoming an increasingly important class of markers in genome-wide phylogenetic studies. Recently, we proposed a new type of RGCs designated RGC_CAMs (after Conserved Amino acids-Multiple substitutions) that were inferred using genome-wide identification of amino acid replacements that were: i) located in unambiguously aligned regions of orthologous genes, ii) shared by two or more taxa in positions that contain a different, conserved amino acid in a much broader range of taxa, and iii) require two or three nucleotide substitutions. When applied to animal phylogeny, the RGC_CAM approach supported the coelomate clade that unites deuterostomes with arthropods as opposed to the ecdysozoan (molting animals) clade. However, a non-negligible level of homoplasy was detected.

RESULTS

We provide a direct estimate of the level of homoplasy caused by parallel changes and reversals among the RGC_CAMs using 462 alignments of orthologous genes from 19 eukaryotic species. It is shown that the impact of parallel changes and reversals on the results of phylogenetic inference using RGC_CAMs cannot explain the observed support for the Coelomata clade. In contrast, the evidence in support of the Ecdysozoa clade, in large part, can be attributed to parallel changes. It is demonstrated that parallel changes are significantly more common in internal branches of different subtrees that are separated from the respective common ancestor by relatively short times than in terminal branches separated by longer time intervals. A similar but much weaker trend was detected for reversals. The observed evolutionary trend of parallel changes is explained in terms of the covarion model of molecular evolution. As the overlap between the covarion sets in orthologous genes from different lineages decreases with time after divergence, the likelihood of parallel changes decreases as well.

CONCLUSION

The level of homoplasy observed here appears to be low enough to justify the utility of RGC_CAMs and other types of RGCs for resolution of hard problems in phylogeny. Parallel changes, one of the major classes of events leading to homoplasy, occur much more often in relatively recently diverged lineages than in those separated from their last common ancestor by longer time intervals of time. This pattern seems to provide the molecular-evolutionary underpinning of Vavilov's law of homologous series and is readily interpreted within the framework of the covarion model of molecular evolution.

The tardigrade Hypsibius dujardini, a new model for studying the evolution of development

Studying development in diverse taxa can address a central issue in evolutionary biology: how morphological diversity arises through the evolution of developmental mechanisms. Two of the best-studied developmental model organisms, the arthropod Drosophila and the nematode Caenorhabditis elegans, have been found to belong to a single protostome superclade, the Ecdysozoa. This finding suggests that a closely related ecdysozoan phylum could serve as a valuable model for studying how developmental mechanisms evolve in ways that can produce diverse body plans. Tardigrades, also called water bears, make up a phylum of microscopic ecdysozoan animals. Tardigrades share many characteristics with C. elegans and Drosophila that could make them useful laboratory models, but long-term culturing of tardigrades historically has been a challenge, and there have been few studies of tardigrade development. Here, we show that the tardigrade Hypsibius dujardini can be cultured continuously for decades and can be cryopreserved. We report that H. dujardini has a compact genome, a little smaller than that of C. elegans or Drosophila, and that sequence evolution has occurred at a typical rate. H. dujardini has a short generation time, 13-14 days at room temperature. We have found that the embryos of H. dujardini have a stereotyped cleavage pattern with asymmetric cell divisions, nuclear migrations, and cell migrations occurring in reproducible patterns. We present a cell lineage of the early embryo and an embryonic staging series. We expect that these data can serve as a platform for using H. dujardini as a model for studying the evolution of developmental mechanisms.

The Hedgehog gene family of the cnidarian, Nematostella vectensis, and implications for understanding metazoan Hedgehog pathway evolution

Hedgehog signaling is an important component of cell-cell communication during bilaterian development, and abnormal Hedgehog signaling contributes to disease and birth defects. Hedgehog genes are composed of a ligand ("hedge") domain and an autocatalytic intein ("hog") domain. Hedgehog (hh) ligands bind to a conserved set of receptors and activate downstream signal transduction pathways terminating with Gli/Ci transcription factors. We have identified five intein-containing genes in the anthozoan cnidarian Nematostella vectensis, two of which (NvHh1 and NvHh2) contain definitive hedgehog ligand domains, suggesting that to date, cnidarians are the earliest branching metazoan phylum to possess definitive Hh orthologs. Expression analysis of NvHh1 and NvHh2, the receptor NvPatched, and a downstream transcription factor NvGli (a Gli3/Ci ortholog) indicate that these genes may have conserved roles in planar and trans-epithelial signaling during gut and germline development, while the three remaining intein-containing genes (NvHint1,2,3) are expressed in a cell-type-specific manner in putative neural precursors. Metazoan intein-containing genes that lack a hh ligand domain have previously only been identified within nematodes. However, we have identified intein-containing genes from both Nematostella and in two newly annotated lophotrochozoan genomes. Phylogenetic analyses suggest that while nematode inteins may be derived from an ancestral true hedgehog gene, the newly identified cnidarian and lophotrochozoan inteins may be orthologous, suggesting that both true hedgehog and hint genes may have been present in the cnidarian-bilaterian ancestor. Genomic surveys of N. vectensis suggest that most of the components of both protostome and deuterostome Hh signaling pathways are present in anthozoans and that some appear to have been lost in ecdysozoan lineages. Cnidarians possess many bilaterian cell-cell signaling pathways (Wnt, TGFbeta, FGF, and Hh) that appear to act in concert to pattern tissues along the oral-aboral axis of the polyp. Cnidarians represent a diverse group of animals with a predominantly epithelial body plan, and perhaps selective pressures to pattern epithelia resulted in the ontogeny of the hedgehog pathway in the common ancestor of the Cnidaria and Bilateria.

Support for the Coelomata clade of animals from a rigorous analysis of the pattern of intron conservation

Many intron positions are conserved in varying subsets of eukaryotic genomes and, consequently, comprise a potentially informative class of phylogenetic characters. Roy and Gilbert developed a method of phylogenetic reconstruction using the patterns of intron presence-absence in eukaryotic genes and, applying this method to the analysis of animal phylogeny, obtained support for an Ecdysozoa clade (Roy SW, Gilbert W. 2005. Resolution of a deep animal divergence by the pattern of intron conservation. Proc Natl Acad Sci USA. 102:4403-4408). The critical assumption in the method was the independence of intron loss in different branches of the phylogenetic tree. Here, this assumption is refuted by showing that the branch-specific intron loss rates are strongly correlated. We show that different tree topologies are obtained, in each case with a significant statistical support, when different subsets of intron positions are analyzed. The analysis of the conserved intron positions supports the Coelomata topology, that is, a clade comprised of arthropods and chordates, whereas the analysis of more variable intron positions favors the Ecdysozoa topology, that is, a clade of arthropods and nematodes. We show, however, that the support for Ecdysozoa is fully explained by parallel loss of introns in nematodes and arthropods, a factor that does not contribute to the analysis of the conserved introns. The developed procedure for the identification and analysis of conserved introns and other characters with minimal or no homoplasy is expected to be useful for resolving many hard phylogenetic problems.

Comprehensive analysis of animal TALE homeobox genes: new conserved motifs and cases of accelerated evolution

TALE homeodomain proteins are an ancient subgroup within the group of homeodomain transcription factors that play important roles in animal, plant, and fungal development. We have extracted the full complement of TALE superclass homeobox genes from the genome projects of seven protostomes, seven deuterostomes, and Nematostella. This was supplemented with TALE homeobox genes from additional species and phylogenetic analyses were carried out with 276 sequences. We found 20 homeobox genes and 4 pseudogenes in humans, 21 genes in mouse, 8 genes in Drosophila, and 5 genes plus one truncated gene in Caenorhabditis elegans. Apart from the previously identified TALE classes MEIS, PBC, IRO, and TGIF, a novel class is identified, termed MOHAWK (MKX). Further, we show that the MEIS class can be divided into two families, PREP and MEIS. Prep genes have previously only been described in vertebrates but are lacking in Drosophila. Here we identify orthologues in other insect taxa as well as in the cnidarian Nematostella. In C. elegans, a divergent Prep protein has lost the homeodomain. Full-length multiple sequence alignment of the protostome and deuterostome sequences allowed us to identify several novel conserved motifs within the MKX, TGIF, and MEIS classes. Phylogenetic analyses revealed fast-evolving PBC class genes; in particular, some X-linked PBC genes in nematodes are subject to rapid evolution. In addition, several instances of gene loss were identified. In conclusion, our comprehensive analysis provides a defining framework for the classification of animal TALE homeobox genes and the understanding of their evolution.

Evolution of a novel subfamily of nuclear receptors with members that each contain two DNA binding domains

Background

Nuclear receptors (NRs) are important transcriptional modulators in metazoans which regulate transcription through binding to the promoter region of their target gene by the DNA binding domain (DBD) and activation or repression of mRNA synthesis through co-regulators bound to the ligand binding domain (LBD). NRs typically have a single DBD with a LBD.

Results

Three nuclear receptors named 2DBD-NRs, were identified from the flatworm Schistosoma mansoni that each possess a novel set of two DBDs in tandem with a LBD. They represent a novel NR modular structure: A/B-DBD-DBD-hinge-LBD. The 2DBD-NRs form a new subfamily of NRs, VII. By database mining, 2DBD-NR genes from other flatworm species (Schmidtea mediterranea and Dugesia japonica), from Mollusks (Lottia gigantean) and from arthropods (Daphnia pulex) were also identified. All 2DBD-NRs possess a P-box sequence of CEACKK in the first DBD, which is unique to 2DBD-NRs, and a P-box sequence of CEGCKG in the second DBD. Phylogenetic analyses of both DBD and ligand binding domain sequences showed that 2DBD-NR genes originate from a common two DBD-containing ancestor gene. A single 2DBD-NR orthologue was found in Arthropoda, Platyhelminths and Mollusca. Subsequent 2DBD-NR gene evolution in Mollusks and Platyhelminths involved gene duplication. Chromosome localization of S. mansoni 2DBD-NR genes by Fluorescent in situ hybridization (FISH) suggests that 2DBD-NR genes duplicated on different chromosomes in the Platyhelminths. Dimerization of Sm2DBDα indicates that 2DBD-NRs may act as homodimers, suggesting either that two repeats of a half-site are necessary for each DBD of 2DBD-NRs to bind to its target gene, or that each 2DBD-NR can recognize multiple sites.

Conclusion

2DBD-NRs share a common ancestor gene which possessed an extra DBD that likely resulted from a recombination event. After the split of the Arthropods, Mollusks and Platyhelminths, 2DBD-NR underwent a recent duplication in a common ancestor of Mollusks, while two rounds of duplication occurred in a common ancestor of the Platyhelminths. This demonstrates that certain NR gene underwent recent duplication in Prostostome lineages after the split of the Prostostomia and Deuterostomia.

Mitogenomics and phylogenomics reveal priapulid worms as extant models of the ancestral Ecdysozoan

Research into arthropod evolution is hampered by the derived nature and rapid evolution of the best-studied out-group: the nematodes. We consider priapulids as an alternative out-group. Priapulids are a small phylum of bottom-dwelling marine worms; their tubular body with spiny proboscis or introvert has changed little over 520 million years and recognizable priapulids are common among exceptionally preserved Cambrian fossils. Using the complete mitochondrial genome and 42 nuclear genes from Priapulus caudatus, we show that priapulids are slowly evolving ecdysozoans; almost all these priapulid genes have evolved more slowly than nematode orthologs and the priapulid mitochondrial gene order may be unchanged since the Cambrian. Considering their primitive bodyplan and embryology and the great conservation of both nuclear and mitochondrial genomes, priapulids may deserve the popular epithet of "living fossil." Their study is likely to yield significant new insights into the early evolution of the Ecdysozoa and the origins of the arthropods and their kin as well as aiding inference of the morphology of ancestral Ecdysozoa and Bilateria and their genomes.

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.

Darwin's dilemma: the realities of the Cambrian 'explosion'

The Cambrian 'explosion' is widely regarded as one of the fulcrum points in the history of life, yet its origins and causes remain deeply controversial. New data from the fossil record, especially of Burgess Shale-type Lagerstätten, indicate, however, that the assembly of bodyplans is not only largely a Cambrian phenomenon, but can already be documented in fair detail. This speaks against a much more ancient origin of the metazoans, and current work is doing much to reconcile the apparent discrepancies between the fossil record, including the Ediacaran assemblages of latest Neoproterozoic age and molecular 'clocks'. Hypotheses to explain the Cambrian 'explosion' continue to be generated, but the recurrent confusion of cause and effect suggests that the wrong sort of question is being asked. Here I propose that despite its step-like function this evolutionary event is the inevitable consequence of Earth and biospheric change.

Expansion of the SOX gene family predated the emergence of the Bilateria

Members of the SOX gene family are involved in regulating many developmental processes including neuronal determination and differentiation, and in carcinogenesis. So far they have only been identified in species from the Bilateria (deuterostomes and protostomes). To understand the origins of the SOX family, we used a PCR-based strategy to obtain 28 new sequences of SOX gene HMG domains from four non-bilaterian Metazoa: two sponge species, one ctenophore and one cnidarian. One additional SOX sequence was retrieved from EST sequences of the cnidarian species Clytia hemisphaerica. Unexpected SOX gene diversity was found in these species, especially in the cnidarian and the ctenophore. The topology of gene relationships deduced by Maximum Likelihood analysis, although not supported by bootstrap values, suggested that the SOX family started to diversify in the metazoan stem branch prior to the divergence of demosponges, and that further diversification occurred in the eumetazoan branch, as well as later in calcisponges, ctenophores, cnidarians and vertebrates. In contrast, gene loss appears to have occurred in the nematode and probably in other protostome lineages, explaining their lower number of SOX genes.

Cannabinoid receptors in invertebrates

Two cannabinoid receptors, CB1 and CB2, are expressed in mammals, birds, reptiles, and fish. The presence of cannabinoid receptors in invertebrates has been controversial, due to conflicting evidence. We conducted a systematic review of the literature, using expanded search parameters. Evidence presented in the literature varied in validity, ranging from crude in vivo behavioural assays to robust in silico ortholog discovery. No research existed for several clades of invertebrates; we therefore tested for cannabinoid receptors in seven representative species, using tritiated ligand binding assays with [3H]CP55,940 displaced by the CB1-selective antagonist SR141716A. Specific binding of [3H]CP55,940 was found in neural membranes of Ciona intestinalis (Deuterstoma, a positive control), Lumbricusterrestris (Lophotrochozoa), and three ecdysozoans: Peripatoides novae-zealandiae (Onychophora), Jasus edwardi (Crustacea) and Panagrellus redivivus (Nematoda); the potency of displacement by SR141716A was comparable to measurements on rat cerebellum. No specific binding was observed in Actinothoe albocincta (Cnidaria) or Tethya aurantium (Porifera). The phylogenetic distribution of cannabinoid receptors may address taxonomic questions; previous studies suggested that the loss of CB1 was a synapomorphy shared by ecdysozoans. Our discovery of cannabinoid receptors in some nematodes, onychophorans, and crustaceans does not contradict the Ecdysozoa hypothesis, but gives it no support. We hypothesize that cannabinoid receptors evolved in the last common ancestor of bilaterians, with secondary loss occurring in insects and other clades. Conflicting data regarding Cnidarians precludes hypotheses regarding the last common ancestor of eumetazoans. No cannabinoid receptors are expressed in sponges, which probably diverged before the origin of the eumetazoan ancestor.

Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements

Phylogenetic relationships among the metazoan phyla are the subject of an ongoing controversy. Analysis of mitochondrial gene arrangements is a powerful tool to investigate these relationships; however, its previous application outside of individual animal phyla has been hampered by the lack of informative out-group data. To address this shortcoming, we determined complete mitochondrial DNA sequences for the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. With sponges as an outgroup, we investigated phylogenetic relationships of nine bilaterian phyla using both breakpoint analysis of global mitochondrial gene arrangements and maximum parsimony analysis of mitochondrial gene adjacencies. Our results provide strong support for a group that includes protostome (but not deuterostome) coelomate, pseudocoelomate, and acoelomate animals, thus clearly rejecting the Coelomata hypothesis. Two other groups of bilaterian animals, Lophotrochozoa and Ambulacraria, are also supported by our analyses. However, due to the remarkable stability of mitochondrial gene arrangements in Deuterostomia and the Ecdysozoa, conclusions on their evolutionary history cannot be drawn.

The anatomy, affinity, and phylogenetic significance of Markuelia

The fossil record provides a paucity of data on the development of extinct organisms, particularly for their embryology. The recovery of fossilized embryos heralds new insight into the evolution of development but advances are limited by an almost complete absence of phylogenetic constraint. Markuelia is an exception to this, known from cleavage and pre-hatchling stages as a vermiform and profusely annulated direct-developing bilaterian with terminal circumoral and posterior radial arrays of spines. Phylogenetic analyses have hitherto suggested assignment to stem-Scalidophora (phyla Kinorhyncha, Loricifera, Priapulida). We test this assumption with additional data and through the inclusion of additional taxa. The available evidence supports stem-Scalidophora affinity, leading to the conclusion that scalidophorans, cyclonerualians, and ecdysozoans are primitive direct developers, and the likelihood that scalidophorans are primitively metameric.

The evolution of metazoan axial properties

Renewed interest in the developmental basis of organismal complexity, and the emergence of new molecular tools, is improving our ability to study the evolution of metazoan body plans. The most substantial changes in body-plan organization occurred early in metazoan evolution; new model systems for studying basal metazoans are now being developed, and total-genome-sequencing initiatives are underway for at least three of the four most important taxa. The elucidation of how the gene networks that are involved in axial organization, germ-layer formation and cell differentiation are used differently during development is generating a more detailed understanding of the events that have led to the current diversity of multicellular life.

Animal phylogeny: fatal attraction

MPhylogenetic analyses of hundreds of genes from model animals have placed flies closer to vertebrates than to nematodes; recent work suggests this may be due to an artefact known as long branch attraction.

First Sequenced Mitochondrial Genome from the Phylum Acanthocephala(Leptorhynchoides thecatus) and Its Phylogenetic Position Within Metazoa

The complete sequence of the mitochondrial genome of Leptorhynchoides thecatus (Acanthocephala) was determined, and a phylogenetic analysis was carried out to determine its placement within Metazoa. The genome is circular, 13,888 bp, and contains at least 36 of the 37 genes typically found in animal mitochondrial genomes. The genes for the large and small ribosomal RNA subunits are shorter than those of most metazoans, and the structures of most of the tRNA genes are atypical. There are two significant noncoding regions (377 and 294 bp), which are the best candidates for a control region; however, these regions do not appear similar to any of the control regions of other animals studied to date. The amino acid and nucleotide sequences of the protein coding genes of L. thecatus and 25 other metazoan taxa were used in both maximum likelihood and maximum parsimony phylogenetic analyses. Results indicate that among taxa with available mitochondrial genome sequences, Platyhelminthes is the closest relative to L. thecatus, which together are the sister taxon of Nematoda; however, long branches and/or base composition bias could be responsible for this result. The monophyly of Ecdysozoa, molting organisms, was not supported by any of the analyses. This study represents the first mitochondrial genome of an acanthocephalan to be sequenced and will allow further studies of systematics, population genetics, and genome evolution.

Genome-scale evidence of the nematode-arthropod clade

The most extensive phylogenetic analysis carried out to date, including 11 complete genomes, is shown to support the Ecdysozoa hypothesis in the open-ended debate of the Coelomata-Ecdysozoa evolutionary problem.

Background

The issue of whether coelomates form a single clade, the Coelomata, or whether all animals that moult an exoskeleton (such as the coelomate arthropods and the pseudocoelomate nematodes) form a distinct clade, the Ecdysozoa, is the most puzzling issue in animal systematics and a major open-ended subject in evolutionary biology. Previous single-gene and genome-scale analyses designed to resolve the issue have produced contradictory results. Here we present the first genome-scale phylogenetic evidence that strongly supports the Ecdysozoa hypothesis.

Results

Through the most extensive phylogenetic analysis carried out to date, the complete genomes of 11 eukaryotic species have been analyzed in order to find homologous sequences derived from 18 human chromosomes. Phylogenetic analysis of datasets showing an increased adjustment to equal evolutionary rates between nematode and arthropod sequences produced a gradual change from support for Coelomata to support for Ecdysozoa. Transition between topologies occurred when fast-evolving sequences of Caenorhabditis elegans were removed. When chordate, nematode and arthropod sequences were constrained to fit equal evolutionary rates, the Ecdysozoa topology was statistically accepted whereas Coelomata was rejected.

Conclusions

The reliability of a monophyletic group clustering arthropods and nematodes was unequivocally accepted in datasets where traces of the long-branch attraction effect were removed. This is the first phylogenomic evidence to strongly support the 'moulting clade' hypothesis.

Consistency of genome-based methods in measuring Metazoan evolution

Seven distinct genome-wide divergence measures were applied pairwise to the nine sequenced animal genomes of human, mouse, rat, chicken, pufferfish, fruit fly, mosquito, and two nematode worms (Caenorhabditis briggsae and Caenorhabditis elegans). Qualitatively, all of these divergence measures are found to correlate with the estimated time since speciation; however, marked deviations are observed in a few lineages. The distinct genome divergence measures also correlate well among themselves, indicating that most of the processes shaping genomes are dominated by neutral events. The deviations from the clock-like scenario in some lineages are observed consistently by several measures, implicitly confirming their reliability.

Lack of heme synthesis in a free-living eukaryote

In most free-living eukaryotes studied thus far, heme is synthesized from a series of intermediates through a well defined evolutionarily conserved pathway. We found that free-living worms, including the model genetic organism Caenorhabditis elegans, and parasitic helminths are unable to synthesize heme de novo, even though these animals contain hemoproteins that function in key biological processes. Radioisotope, fluorescence labeling, and heme analog studies suggest that C. elegans acquires heme from exogenous sources. Iron-deprived worms were unable to grow in the presence of adequate heme unless rescued by increasing heme levels in the growth medium. These data indicate that although worms use dietary heme for incorporation into hemoproteins, ingested heme is also used as an iron source when iron is limiting. Our results provide a biochemical basis for the dependence of worm growth and development on heme, and they suggest that pharmacologic targeting of heme transport pathways in worms could be an important control measure for helminthic infections.

Evolutionary sequence analysis of complete eukaryote genomes

Background

Gene duplication and gene loss during the evolution of eukaryotes have hindered attempts to estimate phylogenies and divergence times of species. Although current methods that identify clusters of orthologous genes in complete genomes have helped to investigate gene function and gene content, they have not been optimized for evolutionary sequence analyses requiring strict orthology and complete gene matrices. Here we adopt a relatively simple and fast genome comparison approach designed to assemble orthologs for evolutionary analysis. Our approach identifies single-copy genes representing only species divergences (panorthologs) in order to minimize potential errors caused by gene duplication. We apply this approach to complete sets of proteins from published eukaryote genomes specifically for phylogeny and time estimation.

Results

Despite the conservative criterion used, 753 panorthologs (proteins) were identified for evolutionary analysis with four genomes, resulting in a single alignment of 287,000 amino acids. With this data set, we estimate that the divergence between deuterostomes and arthropods took place in the Precambrian, approximately 400 million years before the first appearance of animals in the fossil record. Additional analyses were performed with seven, 12, and 15 eukaryote genomes resulting in similar divergence time estimates and phylogenies.

Conclusion

Our results with available eukaryote genomes agree with previous results using conventional methods of sequence data assembly from genomes. They show that large sequence data sets can be generated relatively quickly and efficiently for evolutionary analyses of complete genomes.

Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia

Almost a decade ago, a new phylogeny of bilaterian animals was inferred from small-subunit ribosomal RNA (rRNA) that claimed the monophyly of two major groups of protostome animals: Ecdysozoa (e.g., arthropods, nematodes, onychophorans, and tardigrades) and Lophotrochozoa (e.g., annelids, molluscs, platyhelminths, brachiopods, and rotifers). However, it received little additional support. In fact, several multigene analyses strongly argued against this new phylogeny. These latter studies were based on a large amount of sequence data and therefore showed an apparently strong statistical support. Yet, they covered only a few taxa (those for which complete genomes were available), making systematic artifacts of tree reconstruction more probable. Here we expand this sparse taxonomic sampling and analyze a large data set (146 genes, 35,371 positions) from a diverse sample of animals (35 species). Our study demonstrates that the incongruences observed between rRNA and multigene analyses were indeed due to long-branch attraction artifacts, illustrating the enormous impact of systematic biases on phylogenomic studies. A refined analysis of our data set excluding the most biased genes provides strong support in favor of the new animal phylogeny and in addition suggests that urochordates are more closely related to vertebrates than are cephalochordates. These findings have important implications for the interpretation of morphological and genomic data.

The Opisthokonta and the Ecdysozoa May Not Be Clades: Stronger Support for the Grouping of Plant and Animal than for Animal and Fungi and Stronger Support for the Coelomata than Ecdysozoa

In considering the best possible solutions for answering phylogenetic questions from genomic sequences, we have chosen a strategy that we suggest is superior to others that have gone previously. We have ignored multigene families and instead have used single-gene families. This minimizes the inadvertent analysis of paralogs. We have employed strict data controls and have reasoned that if a protein is not capable of recovering the uncontroversial parts of a phylogenetic tree, then why should we use it for the more controversial parts? We have sliced and diced the data in as many ways as possible in order to uncover the signals in that data. Using this strategy, we have tested two controversial hypotheses concerning eukaryotic phylogenetic relationships: the placement of arthropoda and nematodes and the relationships of animals, plants, and fungi. We have constructed phylogenetic trees from 780 single-gene families from 10 completed genomes and amalgamated these into a single supertree. We have also carried out a total evidence analysis on the only universally distributed protein families that can accurately reconstruct the uncontroversial parts of the phylogenetic tree: a total of five families. In doing so, we ignore the majority of single-gene families that are universally distributed as they do not have the appropriate signals to recover the uncontroversial parts of the tree. We have also ignored every protein that has ever been used previously to address this issue, simply because none of them meet our strict criteria. Using these data controls, site stripping, and multiple analyses, 24 out of 26 analyses strongly support the grouping of vertebrates with arthropods (Coelomata hypothesis) and plants with animals. In the other two analyses, the data were ambivalent. The latter finding overturns an 11-year theory of Eukaryotic evolution; the first confirms what has already been said by others. In the light of this new tree, we re-analyze the evolution of intron gain and loss in the rpL14 gene and find that it is much more compatible with the hypothesis presented here than with the Opisthokonta hypothesis.

Functional evolution of ADAMTS genes: evidence from analyses of phylogeny and gene organization

Background

The ADAMTS (A Disintegrin-like and Metalloprotease with Thrombospondin motifs) proteins are a family of metalloproteases with sequence similarity to the ADAM proteases, that contain the thrombospondin type 1 sequence repeat motifs (TSRs) common to extracellular matrix proteins. ADAMTS proteins have recently gained attention with the discovery of their role in a variety of diseases, including tissue and blood disorders, cancer, osteoarthritis, Alzheimer's and the genetic syndromes Weill-Marchesani syndrome (ADAMTS10), thrombotic thrombocytopenic purpura (ADAMTS13), and Ehlers-Danlos syndrome type VIIC (ADAMTS2) in humans and belted white-spotting mutation in mice (ADAMTS20).

Results

Phylogenetic analysis and comparison of the exon/intron organization of vertebrate (Homo, Mus, Fugu), chordate (Ciona) and invertebrate (Drosophila and Caenorhabditis) ADAMTS homologs has elucidated the evolutionary relationships of this important gene family, which comprises 19 members in humans.

Conclusions

The evolutionary history of ADAMTS genes in vertebrate genomes has been marked by rampant gene duplication, including a retrotransposition that gave rise to a distinct ADAMTS subfamily (ADAMTS1, -4, -5, -8, -15) that may have distinct aggrecanase and angiogenesis functions.

Convergent evolution of domain architectures (is rare)

MOTIVATION

In this paper, we shall examine the evolution of domain architectures across 62 genomes of known phylogeny including all kingdoms of life. We look in particular at the possibility of convergent evolution, with a view to determining the extent to which the architectures observed in the genomes are due to functional necessity or evolutionary descent. We used domains of known structure, because from this and other information we know their evolutionary relationships. We use a range of methods including phylogenetic grouping, sequence similarity/alignment, mutation rates and comparative genomics to approach this difficult problem from several angles.

RESULTS

Although we do not claim an exhaustive analysis, we conclude that between 0.4 and 4% of sequences are involved in convergent evolution of domain architectures, and expect the actual number to be close to the lower bound. We also made two incidental observations, albeit on a small sample: the events leading to convergent evolution appear to be random with no functional or structural preferences, and changes in the number of tandem repeat domains occur more readily than changes which alter the domain composition.

CONCLUSION

The principal conclusion is that the observed domain architectures of the sequences in the genomes are driven by evolutionary descent rather than functional necessity.

CONTACT

gough@supfam.org.

Phylogenetic analysis of the formin homology 2 domain

Formin proteins are key regulators of eukaryotic actin filament assembly and elongation, and many species possess multiple formin isoforms. A nomenclature system based on fundamental features would be desirable, to aid the rapid identification and characterization of novel formins. In this article, we attempt to systematize the formin family by performing phylogenetic analyses of the formin homology 2 (FH2) domain, an independently folding region common to all formins, which alone can influence actin dynamics. Through database searches, we identify 101 FH2 domains from 26 eukaryotic species, including 15 in mice. Sequence alignments reveal a highly conserved yeast-specific insert in the "knob loop" region of the FH2 domain, with unknown functional consequences. Phylogenetic analysis using minimum evolution (ME), maximum parsimony (MP), and maximum likelihood (ML) algorithms strongly supports the existence of seven metazoan groups. Yeast FH2 domains segregate from all other eukaryotes, including metazoans, other fungi, plants, and protists. Sequence comparisons of non-FH2 regions support relationships between three metazoan groups (Dia, DAAM, and FRL) and examine previously identified coiled-coil and Diaphanous auto-regulatory domain sequences. This analysis allows for a formin nomenclature system based on sequence relationships, as well as suggesting strategies for the determination of biochemical and cellular activities of these proteins.