Early-branching or fast-evolving eukaryotes? An answer based on slowly evolving positions

Sources of signal in 62 protein-coding nuclear genes for higher-level phylogenetics of arthropods

BACKGROUND

This study aims to investigate the strength of various sources of phylogenetic information that led to recent seemingly robust conclusions about higher-level arthropod phylogeny and to assess the role of excluding or downweighting synonymous change for arriving at those conclusions.

METHODOLOGY/PRINCIPAL FINDINGS

The current study analyzes DNA sequences from 68 gene segments of 62 distinct protein-coding nuclear genes for 80 species. Gene segments analyzed individually support numerous nodes recovered in combined-gene analyses, but few of the higher-level nodes of greatest current interest. However, neither is there support for conflicting alternatives to these higher-level nodes. Gene segments with higher rates of nonsynonymous change tend to be more informative overall, but those with lower rates tend to provide stronger support for deeper nodes. Higher-level nodes with bootstrap values in the 80% - 99% range for the complete data matrix are markedly more sensitive to substantial drops in their bootstrap percentages after character subsampling than those with 100% bootstrap, suggesting that these nodes are likely not to have been strongly supported with many fewer data than in the full matrix. Data set partitioning of total data by (mostly) synonymous and (mostly) nonsynonymous change improves overall node support, but the result remains much inferior to analysis of (unpartitioned) nonsynonymous change alone. Clusters of genes with similar nonsynonymous rate properties (e.g., faster vs. slower) show some distinct patterns of node support but few conflicts. Synonymous change is shown to contribute little, if any, phylogenetic signal to the support of higher-level nodes, but it does contribute nonphylogenetic signal, probably through its underlying heterogeneous nucleotide composition. Analysis of seemingly conservative indels does not prove useful.

CONCLUSIONS

Generating a robust molecular higher-level phylogeny of Arthropoda is currently possible with large amounts of data and an exclusive reliance on nonsynonymous change.

Evolution and diversity of the Golgi

The Golgi is an ancient and fundamental eukaryotic organelle. Evolutionary cell biological studies have begun establishing the repertoire, processes, and level of complexity of membrane-trafficking machinery present in early eukaryotic cells. This article serves as a review of the literature on the topic of Golgi evolution and diversity and reports a novel comparative genomic survey addressing Golgi machinery in the widest taxonomic diversity of eukaryotes sampled to date. Finally, the article is meant to serve as a primer on the rationale and design of evolutionary cell biological studies, hopefully encouraging readers to consider this approach as an addition to their cell biological toolbox. It is clear that the major machinery involved in vesicle trafficking to and from the Golgi was already in place by the time of the divergence of the major eukaryotic lineages, nearly 2 billion years ago. Much of this complexity was likely generated by an evolutionary process involving gene duplication and coevolution of specificity encoding membrane-trafficking proteins. There have also been clear cases of loss of Golgi machinery in some lineages as well as innovation of novel machinery. The Golgi is a wonderfully complex and diverse organelle and its continued exploration promises insight into the evolutionary history of the eukaryotic cell.

Phylogenomics-based reconstruction of protozoan species tree

We have developed a semi-automatic methodology to reconstruct the phylogenetic species tree in Protozoa, integrating different phylogenetic algorithms and programs, and demonstrating the utility of a supermatrix approach to construct phylogenomics-based trees using 31 universal orthologs (UO). The species tree obtained was formed by three major clades that were related to three groups of data: i) Species containing at least 80% of UO (25/31) in the concatenated multiple alignment or supermatrix, this clade was called C1, ii) Species containing between 50%–79% (15–24/31) of UO called C2, and iii) Species containing less than 50% (1–14/31) of UO called C3. C1 was composed by only protozoan species, C2 was composed by species related to Protozoa, and C3 was composed by some species of C1 (Protozoa) and C2 (related to Protozoa). Our phylogenomics-based methodology using a supermatrix approach proved to be reliable with protozoan genome data and using at least 25 UO, suggesting that (a) the more UO used the better, (b) using the entire UO sequence or just a conserved block of it for the supermatrix produced similar phylogenomic trees.

Splintrons in Giardia intestinalis: Spliceosomal introns in a split form

The divergent eukaryotic unicellular organism Giardia intestinalis is an intestinal parasite in humans and various animals. An analysis of a draft genome sequence suggested that G. intestinalis has a much simpler genome organization and gene repertoire than those of other model eukaryotic organisms (e.g., Arabidopsis and human). This general picture of the G. intestinalis genome seemingly agrees with the fact that only four spliceosomal (cis-spliced) introns have been identified in this organism to date. We have recently shown that G. intestinalis possesses a unique gene expression system incorporating spliceosome-mediated trans-splicing. Some protein-coding genes in G. intestinalis are split into multiple pieces in the genome and each gene fragment is independently transcribed. Two particular pre-mRNAs directly interact with each other by forming an intermolecular-stem structure and are then trans-spliced into a mature mRNA by spliceosomes. We believe that this trans-splicing secondarily arose from the system that excises canonical (cis-splicing) introns. Based on these findings, we suspect that similar phenomena-split genes and post-transcriptional assemblage of their transcripts via trans-splicing-may be prevalent in more distinct eukaryotic lineages than previously known, particularly in organisms possessing "intron-poor" genomes.

Informational gene phylogenies do not support a fourth domain of life for nucleocytoplasmic large DNA viruses

Mimivirus is a nucleocytoplasmic large DNA virus (NCLDV) with a genome size (1.2 Mb) and coding capacity ( 1000 genes) comparable to that of some cellular organisms. Unlike other viruses, Mimivirus and its NCLDV relatives encode homologs of broadly conserved informational genes found in Bacteria, Archaea, and Eukaryotes, raising the possibility that they could be placed on the tree of life. A recent phylogenetic analysis of these genes showed the NCLDVs emerging as a monophyletic group branching between Eukaryotes and Archaea. These trees were interpreted as evidence for an independent "fourth domain" of life that may have contributed DNA processing genes to the ancestral eukaryote. However, the analysis of ancient evolutionary events is challenging, and tree reconstruction is susceptible to bias resulting from non-phylogenetic signals in the data. These include compositional heterogeneity and homoplasy, which can lead to the spurious grouping of compositionally-similar or fast-evolving sequences. Here, we show that these informational gene alignments contain both significant compositional heterogeneity and homoplasy, which were not adequately modelled in the original analysis. When we use more realistic evolutionary models that better fit the data, the resulting trees are unable to reject a simple null hypothesis in which these informational genes, like many other NCLDV genes, were acquired by horizontal transfer from eukaryotic hosts. Our results suggest that a fourth domain is not required to explain the available sequence data.

A late origin of the extant eukaryotic diversity: divergence time estimates using rare genomic changes

Background

Accurate estimation of the divergence time of the extant eukaryotes is a fundamentally important but extremely difficult problem owing primarily to gross violations of the molecular clock at long evolutionary distances and the lack of appropriate calibration points close to the date of interest. These difficulties are intrinsic to the dating of ancient divergence events and are reflected in the large discrepancies between estimates obtained with different approaches. Estimates of the age of Last Eukaryotic Common Ancestor (LECA) vary approximately twofold, from ~1,100 million years ago (Mya) to ~2,300 Mya.

Results

We applied the genome-wide analysis of rare genomic changes associated with conserved amino acids (RGC_CAs) and used several independent techniques to obtain date estimates for the divergence of the major lineages of eukaryotes with calibration intervals for insects, land plants and vertebrates. The results suggest an early divergence of monocot and dicot plants, approximately 340 Mya, raising the possibility of plant-insect coevolution. The divergence of bilaterian animal phyla is estimated at ~400-700 Mya, a range of dates that is consistent with cladogenesis immediately preceding the Cambrian explosion. The origin of opisthokonts (the supergroup of eukaryotes that includes metazoa and fungi) is estimated at ~700-1,000 Mya, and the age of LECA at ~1,000-1,300 Mya. We separately analyzed the red algal calibration interval which is based on single fossil. This analysis produced time estimates that were systematically older compared to the other estimates. Nevertheless, the majority of the estimates for the age of the LECA using the red algal data fell within the 1,200-1,400 Mya interval.

Conclusion

The inference of a "young LECA" is compatible with the latest of previously estimated dates and has substantial biological implications. If these estimates are valid, the approximately 1 to 1.4 billion years of evolution of eukaryotes that is open to comparative-genomic study probably was preceded by hundreds of millions years of evolution that might have included extinct diversity inaccessible to comparative approaches.

Reviewers

This article was reviewed by William Martin, Herve Philippe (nominated by I. King Jordan), and Romain Derelle.

Eukaryotic systematics: a user's guide for cell biologists and parasitologists

Single-celled parasites like Entamoeba, Trypanosoma, Phytophthora and Plasmodium wreak untold havoc on human habitat and health. Understanding the position of the various protistan pathogens in the larger context of eukaryotic diversity informs our study of how these parasites operate on a cellular level, as well as how they have evolved. Here, we review the literature that has brought our understanding of eukaryotic relationships from an idea of parasites as primitive cells to a crystallized view of diversity that encompasses 6 major divisions, or supergroups, of eukaryotes. We provide an updated taxonomic scheme (for 2011), based on extensive genomic, ultrastructural and phylogenetic evidence, with three differing levels of taxonomic detail for ease of referencing and accessibility (see supplementary material at Cambridge Journals On-line). Two of the most pressing issues in cellular evolution, the root of the eukaryotic tree and the evolution of photosynthesis in complex algae, are also discussed along with ideas about what the new generation of genome sequencing technologies may contribute to the field of eukaryotic systematics. We hope that, armed with this user's guide, cell biologists and parasitologists will be encouraged about taking an increasingly evolutionary point of view in the battle against parasites representing real dangers to our livelihoods and lives.

Detecting lateral gene transfers by statistical reconciliation of phylogenetic forests

Background

To understand the evolutionary role of Lateral Gene Transfer (LGT), accurate methods are needed to identify transferred genes and infer their timing of acquisition. Phylogenetic methods are particularly promising for this purpose, but the reconciliation of a gene tree with a reference (species) tree is computationally hard. In addition, the application of these methods to real data raises the problem of sorting out real and artifactual phylogenetic conflict.

Results

We present Prunier, a new method for phylogenetic detection of LGT based on the search for a maximum statistical agreement forest (MSAF) between a gene tree and a reference tree. The program is flexible as it can use any definition of "agreement" among trees. We evaluate the performance of Prunier and two other programs (EEEP and RIATA-HGT) for their ability to detect transferred genes in realistic simulations where gene trees are reconstructed from sequences. Prunier proposes a single scenario that compares to the other methods in terms of sensitivity, but shows higher specificity. We show that LGT scenarios carry a strong signal about the position of the root of the species tree and could be used to identify the direction of evolutionary time on the species tree. We use Prunier on a biological dataset of 23 universal proteins and discuss their suitability for inferring the tree of life.

Conclusions

The ability of Prunier to take into account branch support in the process of reconciliation allows a gain in complexity, in comparison to EEEP, and in accuracy in comparison to RIATA-HGT. Prunier's greedy algorithm proposes a single scenario of LGT for a gene family, but its quality always compares to the best solutions provided by the other algorithms. When the root position is uncertain in the species tree, Prunier is able to infer a scenario per root at a limited additional computational cost and can easily run on large datasets.

Prunier is implemented in C++, using the Bio++ library and the phylogeny program Treefinder. It is available at: http://pbil.univ-lyon1.fr/software/prunier

Identifying cliques of convergent characters: concerted evolution in the cormorants and shags

A phylogenetic tree comprising clades with high bootstrap values or other strong measures of statistical support is usually interpreted as providing a good estimate of the true phylogeny. Convergent evolution acting on groups of characters in concert, however, can lead to highly supported but erroneous phylogenies. Identifying such groups of phylogenetically misleading characters is obviously desirable. Here we present a procedure that uses an independent data source to identify sets of characters that have undergone concerted convergent evolution. We examine the problematic case of the cormorants and shags, for which trees constructed using osteological and molecular characters both have strong statistical support and yet are fundamentally incongruent. We find that the osteological characters can be separated into those that fit the phylogenetic history implied by the molecular data set and those that do not. Moreover, these latter nonfitting osteological characters are internally consistent and form groups of mutually compatible characters or "cliques," which are significantly larger than cliques of shuffled characters. We suggest, therefore, that these cliques of characters are the result of similar selective pressures and are a signature of concerted convergence.

Phylogenomic analyses predict sistergroup relationship of nucleariids and fungi and paraphyly of zygomycetes with significant support

BACKGROUND

Resolving the evolutionary relationships among Fungi remains challenging because of their highly variable evolutionary rates, and lack of a close phylogenetic outgroup. Nucleariida, an enigmatic group of amoeboids, have been proposed to emerge close to the fungal-metazoan divergence and might fulfill this role. Yet, published phylogenies with up to five genes are without compelling statistical support, and genome-level data should be used to resolve this question with confidence.

RESULTS

Our analyses with nuclear (118 proteins) and mitochondrial (13 proteins) data now robustly associate Nucleariida and Fungi as neighbors, an assemblage that we term 'Holomycota'. With Nucleariida as an outgroup, we revisit unresolved deep fungal relationships.

CONCLUSION

Our phylogenomic analysis provides significant support for the paraphyly of the traditional taxon Zygomycota, and contradicts a recent proposal to include Mortierella in a phylum Mucoromycotina. We further question the introduction of separate phyla for Glomeromycota and Blastocladiomycota, whose phylogenetic positions relative to other phyla remain unresolved even with genome-level datasets. Our results motivate broad sampling of additional genome sequences from these phyla.

AIR: A batch-oriented web program package for construction of supermatrices ready for phylogenomic analyses

BACKGROUND

Large multigene sequence alignments have over recent years been increasingly employed for phylogenomic reconstruction of the eukaryote tree of life. Such supermatrices of sequence data are preferred over single gene alignments as they contain vastly more information about ancient sequence characteristics, and are thus more suitable for resolving deeply diverging relationships. However, as alignments are expanded, increasingly numbers of sites with misleading phylogenetic information are also added. Therefore, a major goal in phylogenomic analyses is to maximize the ratio of information to noise; this can be achieved by the reduction of fast evolving sites.

RESULTS

Here we present a batch-oriented web-based program package, named AIR that allows 1) transformation of several single genes to one multigene alignment, 2) identification of evolutionary rates in multigene alignments and 3) removal of fast evolving sites. These three processes can be done with the programs AIR-Appender, AIR-Identifier, and AIR-Remover (AIR), which can be used independently or in a semi-automated pipeline. AIR produces user-friendly output files with filtered and non-filtered alignments where residues are colored according to their evolutionary rates. Other bioinformatics applications linked to the AIR package are available at the Bioportal http://www.bioportal.uio.no, University of Oslo; together these greatly improve the flexibility, efficiency and quality of phylogenomic analyses.

CONCLUSION

The AIR program package allows for efficient creation of multigene alignments and better assessment of evolutionary rates in sequence alignments. Removing fast evolving sites with the AIR programs has been employed in several recent phylogenomic analyses resulting in improved phylogenetic resolution and increased statistical support for branching patterns among the early diverging eukaryotes.

Analysis of rare genomic changes does not support the unikont-bikont phylogeny and suggests cyanobacterial symbiosis as the point of primary radiation of eukaryotes

The deep phylogeny of eukaryotes is an important but extremely difficult problem of evolutionary biology. Five eukaryotic supergroups are relatively well established but the relationship between these supergroups remains elusive, and their divergence seems to best fit a “Big Bang” model. Attempts were made to root the tree of eukaryotes by using potential derived shared characters such as unique fusions of conserved genes. One popular model of eukaryotic evolution that emerged from this type of analysis is the unikont–bikont phylogeny: The unikont branch consists of Metazoa, Choanozoa, Fungi, and Amoebozoa, whereas bikonts include the rest of eukaryotes, namely, Plantae (green plants, Chlorophyta, and Rhodophyta), Chromalveolata, excavates, and Rhizaria. We reexamine the relationships between the eukaryotic supergroups using a genome-wide analysis of rare genomic changes (RGCs) associated with multiple, conserved amino acids (RGC_CAMs and RGC_CAs), to resolve trifurcations of major eukaryotic lineages. The results do not support the basal position of Chromalveolata with respect to Plantae and unikonts or the monophyly of the bikont group and appear to be best compatible with the monophyly of unikonts and Chromalveolata. Chromalveolata show a distinct, additional signal of affinity with Plantae, conceivably, owing to genes transferred from the secondary, red algal symbiont. Excavates are derived forms, with extremely long branches that complicate phylogenetic inference; nevertheless, the RGC analysis suggests that they are significantly more likely to cluster with the unikont–Chromalveolata assemblage than with the Plantae. Thus, the first split in eukaryotic evolution might lie between photosynthetic and nonphotosynthetic forms and so could have been triggered by the endosymbiosis between an ancestral unicellular eukaryote and a cyanobacterium that gave rise to the chloroplast.

Molecular phylogeny: reconstructing the forest

Phylogeny, be it morphological or molecular, has long tried to explain the extant biodiversity by the Tree of Species, which is a logical consequence of strict Darwinian evolutionary principles. Through constant improvement of both methods and data sets, some parts of this diversity have actually been demonstrated to be the result of a tree-like process. For some other parts, and especially for prokaryotes, different molecular markers have, however, produced different evolutionary trees, preventing the reconstruction of such a Tree. While technical artifacts could be blamed for these discrepancies, Lateral Gene Transfers are now largely held for responsible, and their existence requires an extension of the Darwinian framework, since genetic material is not always vertically inherited from parents to offspring. Through a variety of biological processes, sometimes large parts of DNA are exchanged between phylogenetically distant contemporary organisms, especially between those sharing the same environment. While mainly concerning prokaryotes, Lateral Gene Transfers have been also demonstrated to affect eukaryotes, and even multicellular ones, like plants or animals. Most of the time, these transfers allow important adaptations and the colonisation of new niches. The quantitative and qualitative importance of genetic transfers has thus severely challenged the very existence of a universal Tree of Species, since genetic connections, at least for microbes, seem more reticulated than tree-like. Even traditional biological concepts, like the concept of species, need to be re-evaluated in the light of recent discoveries. In short, instead of focusing on a elusive universal tree, biologists are now considering the whole forest corresponding to the multiple processes of inheritance, both vertical and horizontal. This constitutes the major challenge of evolutionary biology for the years to come.

Toward resolving deep neoaves phylogeny: data, signal enhancement, and priors

We report three developments toward resolving the challenge of the apparent basal polytomy of neoavian birds. First, we describe improved conditional down-weighting techniques to reduce noise relative to signal for deeper divergences and find increased agreement between data sets. Second, we present formulae for calculating the probabilities of finding predefined groupings in the optimal tree. Finally, we report a significant increase in data: nine new mitochondrial (mt) genomes (the dollarbird, New Zealand kingfisher, great potoo, Australian owlet-nightjar, white-tailed trogon, barn owl, a roadrunner [a ground cuckoo], New Zealand long-tailed cuckoo, and the peach-faced lovebird) and together they provide data for each of the six main groups of Neoaves proposed by Cracraft J (2001). We use his six main groups of modern birds as priors for evaluation of results. These include passerines, cuckoos, parrots, and three other groups termed "WoodKing" (woodpeckers/rollers/kingfishers), "SCA" (owls/potoos/owlet-nightjars/hummingbirds/swifts), and "Conglomerati." In general, the support is highly significant with just two exceptions, the owls move from the "SCA" group to the raptors, particularly accipitrids (buzzards/eagles) and the osprey, and the shorebirds may be an independent group from the rest of the "Conglomerati". Molecular dating mt genomes support a major diversification of at least 12 neoavian lineages in the Late Cretaceous. Our results form a basis for further testing with both nuclear-coding sequences and rare genomic changes.

Characterizing positive and negative selection and their phylogenetic effects

Total evidence and the use of large datasets to overcome uncertainty are the state of the art in systematic analysis. This assumes that the only true phylogenetic signal is ancestry and that functional, structural, and other factors will not add an alternative signal. Using gene families, where individual codon positions were sorted into bins based upon average-pairwise dN/dS ratio, we show that standard, common phylogenetic methods that were designed for stochastic, neutral, site-independent processes, generate less robust phylogenetic signal for bins with strong negative or positive selection. This was true for phylogenetic reconstruction with parsimony, distance, and likelihood methods. Further, we present a case for the potential existence of systematic functional or structural signal that competes with ancestral signal. For the example of positive selection, we simulate the evolution of sequences through three dimensional lattice constructs with folding constraint and changing binding functionality and show that total evidence for these lattice genes presents trees with functional signal, but that the neutral synonymous sites in these genes show the true ancestral signal. In this case, sequence convergence is promoted by functional convergence.

Phylogeny of endocytic components yields insight into the process of nonendosymbiotic organelle evolution

The process by which some eukaryotic organelles, for example the endomembrane system, evolved without endosymbiotic input remains poorly understood. This problem largely arises because many major cellular systems predate the last common eukaryotic ancestor (LCEA) and thus do not provide examples of organellogenesis in progress. A model is emerging whereby gene duplication and divergence of multiple "specificity-" or "identity-" encoding proteins for the various endomembranous organelles produced the diversity of nonendosymbiotically derived cellular compartments present in modern eukaryotes. To address this possibility, we analyzed three molecular components of the endocytic membrane-trafficking machinery. Phylogenetic analyses of the endocytic syntaxins, Rab 5, and the beta-adaptins each reveal a pattern of ancestral, undifferentiated endocytic homologues in the LCEA. Subsequently, these undifferentiated progenitors independently duplicated in widely divergent lineages, convergently producing components with similar endocytic roles, e.g., beta1 and beta2-adaptin. In contrast, beta3, beta4, and all other adaptin complex subunits, as well as paralogues of the syntaxins and Rabs specific for the other membrane-trafficking organelles, all evolved before the LCEA. Thus, the process giving rise to the differentiated organelles of the endocytic system appears to have been interrupted by the major speciation event that produced the extant eukaryotic lineages. These results suggest that although many endocytic components evolved before the LCEA, other major features evolved independently and convergently after diversification into the primary eukaryotic supergroups. This finding provides an example of a basic cellular system that was simpler in the LCEA than in many extant eukaryotes and yields insight into nonendosymbiotic organelle evolution.

Spectral partitioning of phylogenetic data sets based on compatibility

We describe two new methods to partition phylogenetic data sets of discrete characters based on pairwise compatibility. The partitioning methods make no assumptions regarding the phylogeny, model of evolution, or characteristics of the data. The methods first build a compatibility graph, in which each node represents a character in the data set. Edges in the compatibility graph may represent strict compatibility of characters or they may be weighted based on a fractional compatibility scoring procedure that measures how close the characters are to being compatible. Given the desired number of partitions, the partitioning methods then seek to cluster the characters with the highest average pairwise compatibility, so that characters in each cluster are more compatible with each other than they are with characters in the other cluster(s). Partitioning according to these criteria is computationally intractable (NP-hard); however, spectral methods can quickly provide high-quality solutions. We demonstrate that the spectral partitioning effectively identifies characters with different evolutionary histories in simulated data sets, and it is better at highlighting phylogenetic conflict within empirical data sets than previously used partitioning methods.

Implications of the new eukaryotic systematics for parasitologists

An accurate understanding of evolutionary relationships is central in biology. For parasitologists, understanding the relationships among eukaryotic organisms allows the prediction of virulence mechanisms, reconstruction of metabolic pathways, identification of potential drug targets, elucidation of parasite-specific cellular processes and understanding of interactions with the host or vector. Here we consider the impact of major recent revisions of eukaryotic systematics and taxonomy on parasitology. The previous, ladder-like model placed some protists as early diverging, with the remaining eukaryotes "progressing" towards a "crown radiation" of animals, plants, Fungi and some additional protistan lineages. This model has been robustly disproven. The new model is based on vastly increased amounts of molecular sequence data, integration with morphological information and the rigorous application of phylogenetic methods to those data. It now divides eukaryotes into six major supergroups; the relationships between those groups and the order of branching remain unknown. This new eukaryotic phylogeny emphasizes that organisms including Giardia, Trypanosoma and Trichomonas are not primitive, but instead highly evolved and specialised for their specific environments. The wealth of newly available comparative genomic data has also allowed the reconstruction of ancient suites of characteristics and mapping of character evolution in diverse parasites. For example, the last common eukaryotic ancestor was apparently complex, suggesting that lineage-specific adaptations and secondary losses have been important in the evolution of protistan parasites. Referring to the best evidence-based models for eukaryotic evolution will allow parasitologists to make more accurate and reliable inferences about pathogens that cause significant morbidity and mortality.

An elaborate classification of SNARE proteins sheds light on the conservation of the eukaryotic endomembrane system

Proteins of the SNARE (soluble N-ethylmalemide-sensitive factor attachment protein receptor) family are essential for the fusion of transport vesicles with an acceptor membrane. Despite considerable sequence divergence, their mechanism of action is conserved: heterologous sets assemble into membrane-bridging SNARE complexes, in effect driving membrane fusion. Within the cell, distinct functional SNARE units are involved in different trafficking steps. These functional units are conserved across species and probably reflect the conservation of the particular transport step. Here, we have systematically analyzed SNARE sequences from 145 different species and have established a highly accurate classification for all SNARE proteins. Principally, all SNAREs split into four basic types, reflecting their position in the four-helix bundle complex. Among these four basic types, we established 20 SNARE subclasses that probably represent the original repertoire of a eukaryotic cenancestor. This repertoire has been modulated independently in different lines of organisms. Our data are in line with the notion that the ur-eukaryotic cell was already equipped with the various compartments found in contemporary cells. Possibly, the development of these compartments is closely intertwined with episodes of duplication and divergence of a prototypic SNARE unit.

Evolution of the eukaryotic membrane-trafficking system: origin, tempo and mode

The emergence of an endomembrane system was a crucial stage in the prokaryote-to-eukaryote evolutionary transition. Recent genomic and molecular evolutionary analyses have provided insight into how this critical system arrived at its modern configuration. The apparent relative absence of prokaryotic antecedents for the endomembrane machinery contrasts with the situation for mitochondria, plastids and the nucleus. Overall, the evidence suggests an autogenous origin for the eukaryotic membrane-trafficking machinery. The emerging picture is that early eukaryotic ancestors had a complex endomembrane system, which implies that this cellular system evolved relatively rapidly after the proto-eukaryote diverged away from the other prokaryotic lines. Many of the components of the trafficking system are the result of gene duplications that have produced proteins that have similar functions but differ in their subcellular location. A proto-eukaryote possessing a very simple trafficking system could thus have evolved to near modern complexity in the last common eukaryotic ancestor (LCEA) via paralogous gene family expansion of the proteins encoding organelle identity. The descendents of this common ancestor have undergone further modification of the trafficking machinery; unicellular simplicity and multicellular complexity are the prevailing trend, but there are some remarkable counter-examples.

Toward Resolving the Eukaryotic Tree: The Phylogenetic Positions of Jakobids and Cercozoans

Resolving the global phylogeny of eukaryotes has proven to be challenging. Among the eukaryotic groups of uncertain phylogenetic position are jakobids, a group of bacterivorous flagellates that possess the most bacteria-like mitochondrial genomes known. Jakobids share several ultrastructural features with malawimonads and an assemblage of anaerobic protists (e.g., diplomonads and oxymonads). These lineages together with Euglenozoa and Heterolobosea have collectively been designated "excavates". However, published molecular phylogenies based on the sequences of nuclear rRNAs and up to six nucleus-encoded proteins do not provide convincing support for the monophyly of excavates, nor do they uncover their relationship to other major eukaryotic groups. Here, we report the first large-scale eukaryotic phylogeny, inferred from 143 nucleus-encoded proteins comprising 31,604 amino acid positions, that includes jakobids, malawimonads and cercozoans. We obtain compelling support for the monophyly of jakobids, Euglenozoa plus Heterolobosea (JEH group), and for the association of cercozoans with stramenopiles plus alveolates. Furthermore, we observe a sister-group relationship between the JEH group and malawimonads after removing fast-evolving species from the dataset. We discuss the implications of these results for the concept of "excavates" and for the elucidation of eukaryotic phylogeny in general.

Reconstructing the evolution of the mitochondrial ribosomal proteome

For production of proteins that are encoded by the mitochondrial genome, mitochondria rely on their own mitochondrial translation system, with the mitoribosome as its central component. Using extensive homology searches, we have reconstructed the evolutionary history of the mitoribosomal proteome that is encoded by a diverse subset of eukaryotic genomes, revealing an ancestral ribosome of alpha-proteobacterial descent that more than doubled its protein content in most eukaryotic lineages. We observe large variations in the protein content of mitoribosomes between different eukaryotes, with mammalian mitoribosomes sharing only 74 and 43% of its proteins with yeast and Leishmania mitoribosomes, respectively. We detected many previously unidentified mitochondrial ribosomal proteins (MRPs) and found that several have increased in size compared to their bacterial ancestral counterparts by addition of functional domains. Several new MRPs have originated via duplication of existing MRPs as well as by recruitment from outside of the mitoribosomal proteome. Using sensitive profile-profile homology searches, we found hitherto undetected homology between bacterial and eukaryotic ribosomal proteins, as well as between fungal and mammalian ribosomal proteins, detecting two novel human MRPs. These newly detected MRPs constitute, along with evolutionary conserved MRPs, excellent new screening targets for human patients with unresolved mitochondrial oxidative phosphorylation disorders.

Multiple losses of sex within a single genus of Microsporidia

Background

Most asexual eukaryotic lineages have arisen recently from sexual ancestors and contain few ecologically distinct species, providing evidence for long-term advantages of sex. Ancient asexual lineages provide rare exceptions to this rule and so can yield valuable information relating to the evolutionary forces underlying the maintenance of sex. Microsporidia are parasitic, unicellular fungi. They include many asexual species which have traditionally been grouped together into large, presumably ancient taxonomic groups. However, these putative ancient asexual lineages have been identified on the basis of morphology, life cycles and small subunit ribosomal RNA (16S rRNA) gene sequences, all of which hold questionable value in accurately inferring phylogenetic relationships among microsporidia.

Results

The hypothesis of a single, ancient loss of sex within the Nosema/Vairimorpha group of microsporidia was tested using phylogenetic analyses based on alignments of rRNA and RPB1 gene sequences from sexual and asexual species. Neither set of gene trees supported ancient asexuality, instead indicating at least two, recent losses of sex.

Conclusion

Sex has been lost on multiple, independent occasions within the Nosema/Vairimorpha group of microsporidia and there is no evidence for ancient asexual lineages. It appears therefore that sex confers important long-term advantages even upon highly simplified eukaryotes such as microsporidia. The rapid evolution of microsporidian life cycles indicated by this study also suggests that even closely related microsporidia cannot be assumed to have similar life cycles and the life cycle of each newly discovered species must therefore be completely described. These findings are relevant to the use of microsporidia as biological control agents, since several species under consideration as potential agents have life cycles that have been incompletely described.

Suppression of long-branch attraction artefacts in the animal phylogeny using a site-heterogeneous model

BACKGROUND

Thanks to the large amount of signal contained in genome-wide sequence alignments, phylogenomic analyses are converging towards highly supported trees. However, high statistical support does not imply that the tree is accurate. Systematic errors, such as the Long Branch Attraction (LBA) artefact, can be misleading, in particular when the taxon sampling is poor, or the outgroup is distant. In an otherwise consistent probabilistic framework, systematic errors in genome-wide analyses can be traced back to model mis-specification problems, which suggests that better models of sequence evolution should be devised, that would be more robust to tree reconstruction artefacts, even under the most challenging conditions.

METHODS

We focus on a well characterized LBA artefact analyzed in a previous phylogenomic study of the metazoan tree, in which two fast-evolving animal phyla, nematodes and platyhelminths, emerge either at the base of all other Bilateria, or within protostomes, depending on the outgroup. We use this artefactual result as a case study for comparing the robustness of two alternative models: a standard, site-homogeneous model, based on an empirical matrix of amino-acid replacement (WAG), and a site-heterogeneous mixture model (CAT). In parallel, we propose a posterior predictive test, allowing one to measure how well a model acknowledges sequence saturation.

RESULTS

Adopting a Bayesian framework, we show that the LBA artefact observed under WAG disappears when the site-heterogeneous model CAT is used. Using cross-validation, we further demonstrate that CAT has a better statistical fit than WAG on this data set. Finally, using our statistical goodness-of-fit test, we show that CAT, but not WAG, correctly accounts for the overall level of saturation, and that this is due to a better estimation of site-specific amino-acid preferences.

CONCLUSION

The CAT model appears to be more robust than WAG against LBA artefacts, essentially because it correctly anticipates the high probability of convergences and reversions implied by the small effective size of the amino-acid alphabet at each site of the alignment. More generally, our results provide strong evidence that site-specificities in the substitution process need be accounted for in order to obtain more reliable phylogenetic trees.

Evaluating support for the current classification of eukaryotic diversity

Perspectives on the classification of eukaryotic diversity have changed rapidly in recent years, as the four eukaryotic groups within the five-kingdom classification—plants, animals, fungi, and protists—have been transformed through numerous permutations into the current system of six “supergroups.” The intent of the supergroup classification system is to unite microbial and macroscopic eukaryotes based on phylogenetic inference. This supergroup approach is increasing in popularity in the literature and is appearing in introductory biology textbooks. We evaluate the stability and support for the current six-supergroup classification of eukaryotes based on molecular genealogies. We assess three aspects of each supergroup: (1) the stability of its taxonomy, (2) the support for monophyly (single evolutionary origin) in molecular analyses targeting a supergroup, and (3) the support for monophyly when a supergroup is included as an out-group in phylogenetic studies targeting other taxa. Our analysis demonstrates that supergroup taxonomies are unstable and that support for groups varies tremendously, indicating that the current classification scheme of eukaryotes is likely premature. We highlight several trends contributing to the instability and discuss the requirements for establishing robust clades within the eukaryotic tree of life.

Evolutionary perspectives, including the classification of living organisms, provide the unifying scaffold on which biological knowledge is assembled. Researchers in many areas of biology use evolutionary classifications (taxonomy) in many ways, including as a means for interpreting the origin of evolutionary innovations, as a framework for comparative genetics/genomics, and as the basis for drawing broad conclusions about the diversity of living organisms. Thus, it is essential that taxonomy be robust. Here the authors evaluate the stability of and support for the current classification system of eukaryotic cells (cells with nuclei) in which eukaryotes are divided into six kingdom level categories, or supergroups. These six supergroups unite diverse microbial and macrobial eukaryotic lineages, including the well-known groups of plants, animals, and fungi. The authors assess the stability of supergroup classifications through time and reveal a rapidly changing taxonomic landscape that is difficult to navigate for the specialist and generalist alike. Additionally, the authors find variable support for each of the supergroups in published analyses based on DNA sequence variation. The support for supergroups differs according to the taxonomic area under study and the origin of the genes (e.g., nuclear, plastid) used in the analysis. Encouragingly, combining a conservative approach to taxonomy with increased sampling of microbial eukaryotes and the use of multiple types of data is likely to produce a robust scaffold for the eukaryotic tree of life.

Evolutionary Relationships of Apusomonads Inferred from Taxon-Rich Analyses of 6 Nuclear Encoded Genes

The phylogenetic relationships of the biflagellate protist group Apusomonadidae have been unclear despite the availability of some molecular data. We analyzed sequences from 6 nuclear encoded genes-small-subunit rRNA, large-subunit rRNA, alpha-tubulin, beta-tubulin, actin, and heat shock protein 90-to infer the phylogenetic position of Apusomonas proboscidea Aléxéieff 1924. To increase the taxon richness of the study, we also obtained new sequences from representatives of several other major eukaryotic groups: Chrysochromulina sp. National Institute for Environmental Studies 1333 (Haptophyta), Cyanophora paradoxa (Glaucophyta), Goniomonas truncata (Cryptophyceae), Leucocryptos marina (Kathablepharidae), Mesostigma viride (Streptophyta, Viridiplantae), Peridinium limbatum (Alveolata), Pterosperma cristatum (Prasinophytae, Viridiplantae), Synura sphagnicola (Stramenopiles), and Thaumatomonas sp. (Rhizaria). In most individual gene phylogenies, Apusomonas branched close to either of the 2 related taxa-Opisthokonta (including animals, fungi, and choanoflagellates) or Amoebozoa. Combined analyses of all 4 protein-coding genes or all 6 studied genes strongly supported the hypothesis that Apusomonadidae is closely related to Opisthokonta (or to all other eukaryotic groups except Opisthokonta, depending on the position of the eukaryotic root). Alternative hypotheses were rejected in approximately unbiased tests at the 5% level. However, the strong phylogenetic signal supporting a specific affiliation between Apusomonadidae and Opisthokonta largely originated from the alpha-tubulin data. If alpha-tubulin is not considered, topologies in which Apusomonadidae is sister to Opisthokonta or is sister to Amoebozoa were more or less equally supported. One current model for deep eukaryotic evolution holds that eukaryotes are divided into primary "unikont" and "bikont" clades and are descended from a "uniflagellate" common ancestor. Together with other information, our data suggest instead that unikonts (=Opisthokonta and Amoebozoa) are not strictly monophyletic and are descended from biflagellate ancestors.

Multiple secondary origins of the anaerobic lifestyle in eukaryotes

Classical ideas for early eukaryotic evolution often posited a period of anaerobic evolution producing a nucleated phagocytic cell to engulf the mitochondrial endosymbiont, whose presence allowed the host to colonize emerging aerobic environments. This idea was given credence by the existence of contemporary anaerobic eukaryotes that were thought to primitively lack mitochondria, thus providing examples of the type of host cell needed. However, the groups key to this hypothesis have now been shown to contain previously overlooked mitochondrial homologues called hydrogenosomes or mitosomes; organelles that share common ancestry with mitochondria but which do not carry out aerobic respiration. Mapping these data on the unfolding eukaryotic tree reveals that secondary adaptation to anaerobic habitats is a reoccurring theme among eukaryotes. The apparent ubiquity of mitochondrial homologues bears testament to the importance of the mitochondrial endosymbiosis, perhaps as a founding event, in eukaryotic evolution. Comparative study of different mitochondrial homologues is needed to determine their fundamental importance for contemporary eukaryotic cells.

The origin and diversification of eukaryotes: problems with molecular phylogenetics and molecular clock estimation

Determining the relationships among and divergence times for the major eukaryotic lineages remains one of the most important and controversial outstanding problems in evolutionary biology. The sequencing and phylogenetic analyses of ribosomal RNA (rRNA) genes led to the first nearly comprehensive phylogenies of eukaryotes in the late 1980s, and supported a view where cellular complexity was acquired during the divergence of extant unicellular eukaryote lineages. More recently, however, refinements in analytical methods coupled with the availability of many additional genes for phylogenetic analysis showed that much of the deep structure of early rRNA trees was artefactual. Recent phylogenetic analyses of a multiple genes and the discovery of important molecular and ultrastructural phylogenetic characters have resolved eukaryotic diversity into six major hypothetical groups. Yet relationships among these groups remain poorly understood because of saturation of sequence changes on the billion-year time-scale, possible rapid radiations of major lineages, phylogenetic artefacts and endosymbiotic or lateral gene transfer among eukaryotes. Estimating the divergence dates between the major eukaryote lineages using molecular analyses is even more difficult than phylogenetic estimation. Error in such analyses comes from a myriad of sources including: (i) calibration fossil dates, (ii) the assumed phylogenetic tree, (iii) the nucleotide or amino acid substitution model, (iv) substitution number (branch length) estimates, (v) the model of how rates of evolution change over the tree, (vi) error inherent in the time estimates for a given model and (vii) how multiple gene data are treated. By reanalysing datasets from recently published molecular clock studies, we show that when errors from these various sources are properly accounted for, the confidence intervals on inferred dates can be very large. Furthermore, estimated dates of divergence vary hugely depending on the methods used and their assumptions. Accurate dating of divergence times among the major eukaryote lineages will require a robust tree of eukaryotes, a much richer Proterozoic fossil record of microbial eukaryotes assignable to extant groups for calibration, more sophisticated relaxed molecular clock methods and many more genes sampled from the full diversity of microbial eukaryotes.

Monophyly of Rhizaria and multigene phylogeny of unicellular bikonts

Reconstructing a global phylogeny of eukaryotes is an ongoing challenge of molecular phylogenetics. The availability of genomic data from a broad range of eukaryotic phyla helped in resolving the eukaryotic tree into a topology with a rather small number of large assemblages, but the relationships between these "supergroups" are yet to be confirmed. Rhizaria is the most recently recognized "supergroup," but, in spite of this important position within the tree of life, their representatives are still missing in global phylogenies of eukaryotes. Here, we report the first large-scale analysis of eukaryote phylogeny including data for 2 rhizarian species, the foraminiferan Reticulomyxa filosa and the chlorarachniophyte Bigelowiella natans. Our results confirm the monophyly of Rhizaria (Foraminifera + Cercozoa), with very high bootstrap supports in all analyses. The overall topology of our trees is in agreement with the current view of eukaryote phylogeny with basal division into "unikonts" (Opisthokonts and Ameobozoa) and "bikonts" (Plantae, alveolates, stramenopiles, and excavates). As expected, Rhizaria branch among bikonts; however, their phylogenetic position is uncertain. Depending on the data set and the type of analysis, Rhizaria branch as sister group to either stramenopiles or excavates. Overall, the relationships between the major groups of unicellular bikonts are poorly resolved, despite the use of 85 proteins and the largest taxonomic sampling for this part of the tree available to date. This may be due to an acceleration of evolutionary rates in some bikont phyla or be related to their rapid diversification in the early evolution of eukaryotes.

On the origin and functions of RNA-mediated silencing: from protists to man

Double-stranded RNA has been shown to induce gene silencing in diverse eukaryotes and by a variety of pathways. We have examined the taxonomic distribution and the phylogenetic relationship of key components of the RNA interference (RNAi) machinery in members of five eukaryotic supergroups. On the basis of the parsimony principle, our analyses suggest that a relatively complex RNAi machinery was already present in the last common ancestor of eukaryotes and consisted, at a minimum, of one Argonaute-like polypeptide, one Piwi-like protein, one Dicer, and one RNA-dependent RNA polymerase. As proposed before, the ancestral (but non-essential) role of these components may have been in defense responses against genomic parasites such as transposable elements and viruses. From a mechanistic perspective, the RNAi machinery in the eukaryotic ancestor may have been capable of both small-RNA-guided transcript degradation as well as transcriptional repression, most likely through histone modifications. Both roles appear to be widespread among living eukaryotes and this diversification of function could account for the evolutionary conservation of duplicated Argonaute-Piwi proteins. In contrast, additional RNAi-mediated pathways such as RNA-directed DNA methylation, programmed genome rearrangements, meiotic silencing by unpaired DNA, and miRNA-mediated gene regulation may have evolved independently in specific lineages.

How much data are needed to resolve a difficult phylogeny?: case study in Lamiales

Reconstructing phylogeny is a crucial target of contemporary biology, now commonly approached through computerized analysis of genetic sequence data. In angiosperms, despite recent progress at the ordinal level, many relationships between families remain unclear. Here we take a case study from Lamiales, an angiosperm order in which interfamilial relationships have so far proved particularly problematic. We examine the effect of changing one factor-the quantity of sequence data analyzed-on phylogeny reconstruction in this group. We use simulation to estimate a priori the sequence data that would be needed to resolve an accurate, supported phylogeny of Lamiales. We investigate the effect of increasing the length of sequence data analyzed, the rate of substitution in the sequences used, and of combining gene partitions. This method could be a valuable technique for planning systematic investigations in other problematic groups. Our results suggest that increasing sequence length is a better way to improve support, resolution, and accuracy than employing sequences with a faster substitution rate. Indeed, the latter may in some cases have detrimental effects on phylogeny reconstruction. Further molecular sequencing-of at least 10,000 bp-should result in a fully resolved and supported phylogeny of Lamiales, but at present the problematic aspects of this tree model remain.

Eukaryotic evolution, changes and challenges

The idea that some eukaryotes primitively lacked mitochondria and were true intermediates in the prokaryote-to-eukaryote transition was an exciting prospect. It spawned major advances in understanding anaerobic and parasitic eukaryotes and those with previously overlooked mitochondria. But the evolutionary gap between prokaryotes and eukaryotes is now deeper, and the nature of the host that acquired the mitochondrion more obscure, than ever before.

A “Green” Phosphoribulokinase in Complex Algae with Red Plastids: Evidence for a Single Secondary Endosymbiosis Leading to Haptophytes, Cryptophytes, Heterokonts, and Dinoflagellates

Phosphoribulokinase (PRK) is an essential enzyme of photosynthetic eukaryotes which is active in the plastid-located Calvin cycle and regenerates the substrate for ribulose-bisphosphate carboxylase/oxygenase (Rubisco). Rhodophytes and chlorophytes (red and green algae) recruited their nuclear-encoded PRK from the cyanobacterial ancestor of plastids. The plastids of these organisms can be traced back to a single primary endosymbiosis, whereas, for example, haptophytes, dinoflagellates, and euglenophytes obtained their "complex" plastids through secondary endosymbioses, comprising the engulfment of a unicellular red or green alga by a eukaryotic host cell. We have cloned eight new PRK sequences from complex algae as well as a rhodophyte in order to investigate their evolutionary origin. All available PRK sequences were used for phylogenetic analyses and the significance of alternative topologies was estimated by the approximately unbiased test. Our analyses led to several astonishing findings. First, the close relationship of PRK genes of haptophytes, heterokontophytes, cryptophytes, and dinophytes (complex red lineage) supports a monophyletic origin of their sequences and hence their plastids. Second, based on PRK genes the complex red lineage forms a highly supported assemblage together with chlorophytes and land plants, to the exclusion of the rhodophytes. This green affinity is in striking contrast to the expected red algal origin and our analyses suggest that the PRK gene was acquired once via lateral transfer from a green alga. Third, surprisingly the complex green lineages leading to Bigelowiella and Euglena probably also obtained their PRK genes via lateral gene transfers from a red alga and a complex alga with red plastids, respectively.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

Local origin and diversification among Lampornis hummingbirds: A Mesoamerican taxon

The huge biodiversity found in Mesoamerica is often explained by its geographic situation as a natural bridge between two large biogeographic regions. Often overlooked, however, are the high levels of speciation and diversification in the area. Here we assess the phylogenetic relationships within a Mesoamerican group of hummingbirds (Lampornis). We sequenced both mtDNA (1,143 bp of cyt b and 727 bp of ND5) and nuclear genes (505 bp of AK-5 intron and 567 bp of c-mos) for each of the seven recognised species and outgroups. We find two or three clades of similar age within this genus: L. clemenciae and L. amethystinus (singly or as each other's sister taxa) and a Central American clade. This Central-American clade presents a clear bipartition between northern (L. viridipallens and L. sybillae) and southern Mesoamerica, which is shared with many other Mesoamerican organisms. Our analyses suggest that L. hemileucus does not belong in the genus Lampornis. While we refrain to apply a time-scale to our data because of the lack of an appropriate calibration, our results indicate that the genus Lampornis predates the uprising of the Panama land-bridge, and that diversification among the isthmian species (L. castaneoventris and L. calolaema) is a very recent event. Our results strongly suggest a local Mesoamerican origin for this genus.

The largest subunit of RNA polymerase II from the Glaucocystophyta: functional constraint and short-branch exclusion in deep eukaryotic phylogeny

BACKGROUND

Evolutionary analyses of the largest subunit of RNA polymerase II (RPB1) have yielded important and at times provocative results. One particularly troublesome outcome is the consistent inference of independent origins of red algae and green plants, at odds with the more widely accepted view of a monophyletic Plantae comprising all eukaryotes with primary plastids. If the hypothesis of a broader kingdom Plantae is correct, then RPB1 trees likely reflect a persistent phylogenetic artifact. To gain a better understanding of RNAP II evolution, and the presumed artifact relating to green plants and red algae, we isolated and analyzed RPB1 from representatives of Glaucocystophyta, the third eukaryotic group with primary plastids.

RESULTS

Phylogenetic analyses incorporating glaucocystophytes do not recover a monophyletic Plantae; rather they result in additional conflicts with the most widely held views on eukaryotic relationships. In particular, glaucocystophytes are recovered as sister to several amoebozoans with strong support. A detailed investigation shows that this clade can be explained by what we call "short-branch exclusion," a phylogenetic artifact integrally associated with "long-branch attraction." Other systematic discrepancies observed in RPB1 trees can be explained as phylogenetic artifacts; however, these apparent artifacts also appear in regions of the tree that support widely held views of eukaryotic evolution. In fact, most of the RPB1 tree is consistent with artifacts of rate variation among sequences and co-variation due to functional constraints related to C-terminal domain based RNAP II transcription.

CONCLUSION

Our results reveal how subtle and easily overlooked biases can dominate the overall results of molecular phylogenetic analyses of ancient eukaryotic relationships. Sources of potential phylogenetic artifact should be investigated routinely, not just when obvious "long-branch attraction" is encountered.

Comprehensive multigene phylogenies of excavate protists reveal the evolutionary positions of "primitive" eukaryotes

Many of the protists thought to represent the deepest branches on the eukaryotic tree are assigned to a loose assemblage called the "excavates." This includes the mitochondrion-lacking diplomonads and parabasalids (e.g., Giardia and Trichomonas) and the jakobids (e.g., Reclinomonas). We report the first multigene phylogenetic analyses to include a comprehensive sampling of excavate groups (six nuclear-encoded protein-coding genes, nine of the 10 recognized excavate groups). Excavates coalesce into three clades with relatively strong maximum likelihood bootstrap support. Only the phylogenetic position of Malawimonas is uncertain. Diplomonads, parabasalids, and the free-living amitochondriate protist Carpediemonas are closely related to each other. Two other amitochondriate excavates, oxymonads and Trimastix, form the second monophyletic group. The third group is comprised of Euglenozoa (e.g., trypanosomes), Heterolobosea, and jakobids. Unexpectedly, jakobids appear to be specifically related to Heterolobosea. This tree topology calls into question the concept of Discicristata as a supergroup of eukaryotes united by discoidal mitochondrial cristae and makes it implausible that jakobids represent an independent early-diverging eukaryotic lineage. The close jakobids-Heterolobosea-Euglenozoa connection demands complex evolutionary scenarios to explain the transition between the presumed ancestral bacterial-type mitochondrial RNA polymerase found in jakobids and the phage-type protein in other eukaryotic lineages, including Euglenozoa and Heterolobosea.

An Empirical Assessment of Long-Branch Attraction Artefacts in Deep Eukaryotic Phylogenomics

In the context of exponential growing molecular databases, it becomes increasingly easy to assemble large multigene data sets for phylogenomic studies. The expected increase of resolution due to the reduction of the sampling (stochastic) error is becoming a reality. However, the impact of systematic biases will also become more apparent or even dominant. We have chosen to study the case of the long-branch attraction artefact (LBA) using real instead of simulated sequences. Two fast-evolving eukaryotic lineages, whose evolutionary positions are well established, microsporidia and the nucleomorph of cryptophytes, were chosen as model species. A large data set was assembled (44 species, 133 genes, and 24,294 amino acid positions) and the resulting rooted eukaryotic phylogeny (using a distant archaeal outgroup) is positively misled by an LBA artefact despite the use of a maximum likelihood-based tree reconstruction method with a complex model of sequence evolution. When the fastest evolving proteins from the fast lineages are progressively removed (up to 90%), the bootstrap support for the apparently artefactual basal placement decreases to virtually 0%, and conversely only the expected placement, among all the possible locations of the fast-evolving species, receives increasing support that eventually converges to 100%. The percentage of removal of the fastest evolving proteins constitutes a reliable estimate of the sensitivity of phylogenetic inference to LBA. This protocol confirms that both a rich species sampling (especially the presence of a species that is closely related to the fast-evolving lineage) and a probabilistic method with a complex model are important to overcome the LBA artefact. Finally, we observed that phylogenetic inference methods perform strikingly better with simulated as opposed to real data, and suggest that testing the reliability of phylogenetic inference methods with simulated data leads to overconfidence in their performance. Although phylogenomic studies can be affected by systematic biases, the possibility of discarding a large amount of data containing most of the nonphylogenetic signal allows recovering a phylogeny that is less affected by systematic biases, while maintaining a high statistical support.

Evidence from small-subunit ribosomal RNA sequences for a fungal origin of Microsporidia

The phylum Microsporidia comprises a species-rich group of minute, single-celled, and intra-cellular parasites. Lacking normal mitochondria and with unique cytology, microsporidians have sometimes been thought to be a lineage of ancient eukaryotes. Although phylogenetic analyses using small-subunit ribosomal RNA (SSU-rRNA) genes almost invariably place the Microsporidia among the earliest branches on the eukaryotic tree, many other molecules suggest instead a relationship with fungi. Using maximum likelihood methods and a diverse SSU-rRNA data set, we have re-evaluated the phylogenetic affiliations of Microsporidia. We demonstrate that tree topologies used to estimate likelihood model parameters can materially affect phylogenetic searches. We present a procedure for reducing this bias: "tree-based site partitioning," in which a comprehensive set of alternative topologies is used to estimate sequence data partitions based on inferred evolutionary rates. This hypothesis-driven approach appears to be capable of utilizing phylogenetic information that is not available to standard likelihood implementations (e.g., approximation to a gamma distribution); we have employed it in maximum likelihood and Bayesian analysis. Applying our method to a phylogenetically diverse SSU-rRNA data set revealed that the early diverging ("deep") placement of Microsporidia typically found in SSU-rRNA trees is no better than a fungal placement, and that the likeliest placement of Microsporidia among non-long-branch eukaryotic taxa is actually within fungi. These results illustrate the importance of hypothesis testing in parameter estimation, provide a way to address certain problems in difficult data sets, and support a fungal origin for the Microsporidia.