Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.     

Is cladogenesis heritable?

The heritability of speciation rates and extinction risks is a crucial parameter in models of macroevolution, but little direct evidence has been available to assess the occurrence, strength, or generality of this heritability. We tested for heritability using correlations between ancestral and descendent branch lengths in phylogenetic trees, an approach first applied to a bird phylogeny by Harvey et al. (1991, pages 123-137 in Genes in ecology [R. J. Berry et al., eds.], Blackwell Scientific, Oxford). We applied Harvey et al.'s test to some of the largest DNA sequence-based phylogenetic analyses published to date for plants, insects, fungi, and bacteria. If one of two parent lineages splits first and if this is the case for any heritable reason, then on average we expect its daughter lineages to also split first. We also used a randomization procedure to assess significance of branch length heritability. Using maximum parsimony and maximum likelihood branch lengths and trees made ultrametric after nonparametric rate smoothing or by enforcing a molecular clock, we found a pattern for most clades consistent with heritable net cladogenesis. Heritability of cladogenesis may be a general phenomenon, detectable across a large number of lineages and a broad range of taxa.     

Trends and barriers to lateral gene transfer in prokaryotes

Gene acquisition by lateral gene transfer (LGT) is an important mechanism for natural variation among prokaryotes. Laboratory experiments show that protein-coding genes can be laterally transferred extremely fast among microbial cells, inherited to most of their descendants, and adapt to a new regulatory regime within a short time. Recent advance in the phylogenetic analysis of microbial genomes using networks approach reveals a substantial impact of LGT during microbial genome evolution. Phylogenomic networks of LGT among prokaryotes reconstructed from completely sequenced genomes uncover barriers to LGT in multiple levels. Here we discuss the kinds of barriers to gene acquisition in nature including physical barriers for gene transfer between cells, genomic barriers for the integration of acquired DNA, and functional barriers for the acquisition of new genes.     

New Findings on Evolution of Metal Homeostasis Genes: Evidence from Comparative Genome Analysis of Bacteria and Archaea

In order to examine the natural history of metal homeostasis genes in prokaryotes, open reading frames with homology to characterized P_IB-type ATPases from the genomes of 188 bacteria and 22 archaea were investigated. Major findings were as follows. First, a high diversity in N-terminal metal binding motifs was observed. These motifs were distributed throughout bacterial and archaeal lineages, suggesting multiple loss and acquisition events. Second, the CopA locus separated into two distinct phylogenetic clusters, CopA1, which contained ATPases with documented Cu(I) influx activity, and CopA2, which contained both efflux and influx transporters and spanned the entire diversity of the bacterial domain, suggesting that CopA2 is the ancestral locus. Finally, phylogentic incongruences between 16S rRNA and P_IB-type ATPase gene trees identified at least 14 instances of lateral gene transfer (LGT) that had occurred among diverse microbes. Results from bootstrapped supported nodes indicated that (i) a majority of the transfers occurred among proteobacteria, most likely due to the phylogenetic relatedness of these organisms, and (ii) gram-positive bacteria with low moles percent G+C were often involved in instances of LGT. These results, together with our earlier work on the occurrence of LGT in subsurface bacteria (J. M. Coombs and T. Barkay, Appl. Environ. Microbiol. 70:1698-1707, 2004), indicate that LGT has had a minor role in the evolution of P_IB-type ATPases, unlike other genes that specify survival in metal-stressed environments. This study demonstrates how examination of a specific locus across microbial genomes can contribute to the understanding of phenotypes that are critical to the interactions of microbes with their environment.     

New findings on evolution of metal homeostasis genes: evidence from comparative genome analysis of bacteria and archaea

In order to examine the natural history of metal homeostasis genes in prokaryotes, open reading frames with homology to characterized P(IB)-type ATPases from the genomes of 188 bacteria and 22 archaea were investigated. Major findings were as follows. First, a high diversity in N-terminal metal binding motifs was observed. These motifs were distributed throughout bacterial and archaeal lineages, suggesting multiple loss and acquisition events. Second, the CopA locus separated into two distinct phylogenetic clusters, CopA1, which contained ATPases with documented Cu(I) influx activity, and CopA2, which contained both efflux and influx transporters and spanned the entire diversity of the bacterial domain, suggesting that CopA2 is the ancestral locus. Finally, phylogentic incongruences between 16S rRNA and P(IB)-type ATPase gene trees identified at least 14 instances of lateral gene transfer (LGT) that had occurred among diverse microbes. Results from bootstrapped supported nodes indicated that (i) a majority of the transfers occurred among proteobacteria, most likely due to the phylogenetic relatedness of these organisms, and (ii) gram-positive bacteria with low moles percent G+C were often involved in instances of LGT. These results, together with our earlier work on the occurrence of LGT in subsurface bacteria (J. M. Coombs and T. Barkay, Appl. Environ. Microbiol. 70:1698-1707, 2004), indicate that LGT has had a minor role in the evolution of P(IB)-type ATPases, unlike other genes that specify survival in metal-stressed environments. This study demonstrates how examination of a specific locus across microbial genomes can contribute to the understanding of phenotypes that are critical to the interactions of microbes with their environment.     

Recent events dominate interdomain lateral gene transfers between prokaryotes and eukaryotes and,with the exception of endosymbiotic gene transfers,few ancient transfer events persist

While there is compelling evidence for the impact of endosymbiotic gene transfer (EGT; transfer from either mitochondrion or chloroplast to the nucleus) on genome evolution in eukaryotes, the role of interdomain transfer from bacteria and/or archaea (i.e. prokaryotes) is less clear. Lateral gene transfers (LGTs) have been argued to be potential sources of phylogenetic information, particularly for reconstructing deep nodes that are difficult to recover with traditional phylogenetic methods. We sought to identify interdomain LGTs by using a phylogenomic pipeline that generated 13 465 single gene trees and included up to 487 eukaryotes, 303 bacteria and 118 archaea. Our goals include searching for LGTs that unite major eukaryotic clades, and describing the relative contributions of LGT and EGT across the eukaryotic tree of life. Given the difficulties in interpreting single gene trees that aim to capture the approximately 1.8 billion years of eukaryotic evolution, we focus on presence–absence data to identify interdomain transfer events. Specifically, we identify 1138 genes found only in prokaryotes and representatives of three or fewer major clades of eukaryotes (e.g. Amoebozoa, Archaeplastida, Excavata, Opisthokonta, SAR and orphan lineages). The majority of these genes have phylogenetic patterns that are consistent with recent interdomain LGTs and, with the notable exception of EGTs involving photosynthetic eukaryotes, we detect few ancient interdomain LGTs. These analyses suggest that LGTs have probably occurred throughout the history of eukaryotes, but that ancient events are not maintained unless they are associated with endosymbiotic gene transfer among photosynthetic lineages.     

Phylogenetic analyses of diplomonad genes reveal frequent lateral gene transfers affecting eukaryotes

BACKGROUND: Lateral gene transfer (LGT) is an important evolutionary mechanism among prokaryotes. The situation in eukaryotes is less clear; the human genome sequence failed to give strong support for any recent transfers from prokaryotes to vertebrates, yet a number of LGTs from prokaryotes to protists (unicellular eukaryotes) have been documented. Here, we perform a systematic analysis to investigate the impact of LGT on the evolution of diplomonads, a group of anaerobic protists.RESULTS: Phylogenetic analyses of 15 genes present in the genome of the Atlantic Salmon parasite Spironucleus barkhanus and/or the intestinal parasite Giardia lamblia show that most of these genes originated via LGT. Half of the genes are putatively involved in processes related to an anaerobic lifestyle, and this finding suggests that a common ancestor, which most probably was aerobic, of Spironucleus and Giardia adapted to an anaerobic environment in part by acquiring genes via LGT from prokaryotes. The sources of the transferred diplomonad genes are found among all three domains of life, including other eukaryotes. Many of the phylogenetic reconstructions show eukaryotes emerging in several distinct regions of the tree, strongly suggesting that LGT not only involved diplomonads, but also involved other eukaryotic groups.CONCLUSIONS: Our study shows that LGT is a significant evolutionary mechanism among diplomonads in particular and protists in general. These findings provide insights into the evolution of biochemical pathways in early eukaryote evolution and have important implications for studies of eukaryotic genome evolution and organismal relationships. Furthermore, "fusion" hypotheses for the origin of eukaryotes need to be rigorously reexamined in the light of these results.     

Detecting the historical signature of key innovations using stochastic models of character evolution and cladogenesis

Phylogenetic evidence for biological traits that increase the net diversification rate of lineages (key innovations) is most commonly drawn from comparisons of clade size. This can work well for ancient, unreversed traits and for correlating multiple trait origins with higher diversification rates, but it is less suitable for unique events, recently evolved innovations, and traits that exhibit homoplasy. Here I present a new method for detecting the phylogenetic signature of key innovations that tests whether the evolutionary history of the candidate trait is associated with shorter waiting times between cladogenesis events. The method employs stochastic models of character evolution and cladogenesis and integrates well into a Bayesian framework in which uncertainty in historical inferences (such as phylogenetic relationships) is allowed. Applied to a well-known example in plants, nectar spurs in columbines, the method gives much stronger support to the key innovation hypothesis than previous tests.     

Evolution of the cutinase gene family: evidence for lateral gene transfer of a candidate Phytophthora virulence factor

Lateral gene transfer (LGT) can facilitate the acquisition of new functions in recipient lineages, which may enable them to colonize new environments. Several recent publications have shown that gene transfer between prokaryotes and eukaryotes occurs with appreciable frequency. Here we present a study of interdomain gene transfer of cutinases -- well documented virulence factors in fungi -- between eukaryotic plant pathogens Phytophthora species and prokaryotic bacterial lineages. Two putative cutinase genes were cloned from Phytophthora brassicae and Northern blotting experiments showed that these genes are expressed early during the infection of the host Arabidopsis thaliana and induced during cyst germination of the pathogen. Analysis of the gene organisation of this gene family in Phytophthora ramorum and P. sojae showed three and ten copies in tight succession within a region of 5 and 25 kb, respectively, probably indicating a recent expansion in Phytophthora lineages by gene duplications. Bioinformatic analyses identified orthologues only in three genera of Actinobacteria, and in two distantly related eukaryotic groups: oomycetes and fungi. Together with phylogenetic analyses this limited distribution of the gene in the tree of life strongly support a scenario where cutinase genes originated after the origin of land plants in a microbial lineage living in proximity of plants and subsequently were transferred between distantly related plant-degrading microbes. More precisely, a cutinase gene was likely acquired by an ancestor of P. brassicae, P. sojae, P. infestans and P. ramorum, possibly from an actinobacterial source, suggesting that gene transfer might be an important mechanism in the evolution of their virulence. These findings could indeed provide an interesting model system to study acquisition of virulence factors in these important plant pathogens.     

Whole-tree methods for detecting differential diversification rates

Prolific cladogenesis, adaptive radiation, species selection, key innovations, and mass extinctions are a few examples of biological phenomena that lead to differential diversification among lineages. Central to the study of differential diversification rates is the ability to distinguish chance variation from that which requires deterministic explanation. To detect diversification rate variation among lineages, we propose a number of methods that incorporate information on the topological distribution of species diversity from all internal nodes of a phylogenetic tree. These whole-tree methods (M(Pi), M(Sigma), and M(R)) are explicitly connected to a null model of random diversification--the equal-rates Markov (ERM) random branching model--and an alternative model of differential diversification: M(Pi) is based on the product of individual nodal ERM probabilities; M(Sigma) is based on the sum of individual nodal ERM probabilities, and M(R) is based on a transformation of ERM probabilities that corresponds to a formalized system that orders trees by their relative symmetry. These methods have been implemented in a freely available computer program, SYMMETREE, to detect clades with variable diversification rates, thereby allowing the study of biological processes correlated with and possibly causal to shifts in diversification rate. Application of these methods to several published phylogenies demonstrates their ability to contend with relatively large, incompletely resolved trees. These topology-based methods do not require estimates of relative branch lengths, which should facilitate the analysis of phylogenies, such as supertrees, for which such data are unreliable or unavailable.     

The chastity of amoebae: re-evaluating evidence for sex in amoeboid organisms

Amoebae are generally assumed to be asexual. We argue that this view is a relict of early classification schemes that lumped all amoebae together inside the 'lower' protozoa, separated from the 'higher' plants, animals and fungi. This artificial classification allowed microbial eukaryotes, including amoebae, to be dismissed as primitive, and implied that the biological rules and theories developed for macro-organisms need not apply to microbes. Eukaryotic diversity is made up of 70+ lineages, most of which are microbial. Plants, animals and fungi are nested among these microbial lineages. Thus, theories on the prevalence and maintenance of sex developed for macro-organisms should in fact apply to microbial eukaryotes, though the theories may need to be refined and generalized (e.g. to account for the variation in sexual strategies and prevalence of facultative sex in natural populations of many microbial eukaryotes). We use a revised phylogenetic framework to assess evidence for sex in several amoeboid lineages that are traditionally considered asexual, and we interpret this evidence in light of theories on the evolution of sex developed for macro-organisms. We emphasize that the limited data available for many lineages coupled with natural variation in microbial life cycles overestimate the extent of asexuality. Mapping sexuality onto the eukaryotic tree of life demonstrates that the majority of amoeboid lineages are, contrary to popular belief, anciently sexual, and that most asexual groups have probably arisen recently and independently. Additionally, several unusual genomic traits are prevalent in amoeboid lineages, including cyclic polyploidy, which may serve as alternative mechanisms to minimize the deleterious effects of asexuality.     

Assessing the effects of a sequestered germline on interdomain lateral gene transfer in Metazoa

A sequestered germline in Metazoa has been argued to be an obstacle to lateral gene transfer (LGT), though few studies have specifically assessed this claim. Here, we test the hypothesis that the origin of a sequestered germline reduced LGT events in Bilateria (i.e., triploblast lineages) as compared to early‐diverging Metazoa (i.e., Ctenophora, Cnidaria, Porifera, and Placozoa). We analyze single‐gene phylogenies generated with over 900 species sampled from among Bacteria, Archaea, and Eukaryota to identify well‐supported interdomain LGTs. We focus on ancient interdomain LGT (i.e., those between prokaryotes and multiple lineages of Metazoa) as systematic errors in single‐gene tree reconstruction create uncertainties for interpreting eukaryote‐to‐eukaryote transfer. The breadth of the sampled Metazoa enables us to estimate the timing of LGTs, and to examine the pattern before versus after the evolution of a sequestered germline. We identified 58 LGTs found only in Metazoa and prokaryotes (i.e., bacteria and/or archaea), and seven genes transferred from prokaryotes into Metazoa plus one other eukaryotic clade. Our analyses indicate that more interdomain transfers occurred before the development of a sequestered germline, consistent with the hypothesis that this feature is an obstacle to LGT.     

Evolution of glutamate dehydrogenase genes: evidence for lateral gene transfer within and between prokaryotes and eukaryotes

Background

Lateral gene transfer can introduce genes with novel functions into genomes or replace genes with functionally similar orthologs or paralogs. Here we present a study of the occurrence of the latter gene replacement phenomenon in the four gene families encoding different classes of glutamate dehydrogenase (GDH), to evaluate and compare the patterns and rates of lateral gene transfer (LGT) in prokaryotes and eukaryotes.

Results

We extend the taxon sampling of gdh genes with nine new eukaryotic sequences and examine the phylogenetic distribution pattern of the various GDH classes in combination with maximum likelihood phylogenetic analyses. The distribution pattern analyses indicate that LGT has played a significant role in the evolution of the four gdh gene families. Indeed, a number of gene transfer events are identified by phylogenetic analyses, including numerous prokaryotic intra-domain transfers, some prokaryotic inter-domain transfers and several inter-domain transfers between prokaryotes and microbial eukaryotes (protists).

Conclusion

LGT has apparently affected eukaryotes and prokaryotes to a similar extent within the gdh gene families. In the absence of indications that the evolution of the gdh gene families is radically different from other families, these results suggest that gene transfer might be an important evolutionary mechanism in microbial eukaryote genome evolution.

     

Turning the crown upside down: gene tree parsimony roots the eukaryotic tree of life

The first analyses of gene sequence data indicated that the eukaryotic tree of life consisted of a long stem of microbial groups "topped" by a crown-containing plants, animals, and fungi and their microbial relatives. Although more recent multigene concatenated analyses have refined the relationships among the many branches of eukaryotes, the root of the eukaryotic tree of life has remained elusive. Inferring the root of extant eukaryotes is challenging because of the age of the group (～1.7-2.1 billion years old), tremendous heterogeneity in rates of evolution among lineages, and lack of obvious outgroups for many genes. Here, we reconstruct a rooted phylogeny of extant eukaryotes based on minimizing the number of duplications and losses among a collection of gene trees. This approach does not require outgroup sequences or assumptions of orthology among sequences. We also explore the impact of taxon and gene sampling and assess support for alternative hypotheses for the root. Using 20 gene trees from 84 diverse eukaryotic lineages, this approach recovers robust eukaryotic clades and reveals evidence for a eukaryotic root that lies between the Opisthokonta (animals, fungi and their microbial relatives) and all remaining eukaryotes.     

How big is the iceberg of which organellar genes in nuclear genomes are but the tip?

As more and more complete bacterial and archaeal genome sequences become available, the role of lateral gene transfer (LGT) in shaping them becomes more and more clear. Over the long term, it may be the dominant force, affecting most genes in most prokaryotes. We review the history of LGT, suggesting reasons why its prevalence and impact were so long dismissed. We discuss various methods purporting to measure the extent of LGT, and evidence for and against the notion that there is a core of never-exchanged genes shared by all genomes, from which we can deduce the "true" organismal tree. We also consider evidence for, and implications of, LGT between prokaryotes and phagocytic eukaryotes.     

Functional replacement of a primary metabolic pathway via multiple independent eukaryote‐to‐eukaryote gene transfers and selective retention

Although lateral gene transfer (LGT) is now recognized as a major force in the evolution of prokaryotes, the contribution of LGT to the evolution and diversification of eukaryotes is less understood. Notably, transfers of complete pathways are believed to be less likely between eukaryotes, because the successful transfer of a pathway requires the physical clustering of functionally related genes. Here, we report that in one of the closest unicellular relatives of animals, the choanoflagellate, Monosiga, three genes whose products work together in the glutamate synthase cycle are of algal origin. The concerted retention of these three independently acquired genes is best explained as the consequence of a series of adaptive replacement events. More generally, this study argues that (i) eukaryote‐to‐eukaryote transfers of entire metabolic pathways are possible, (ii) adaptive functional replacements of primary pathways can occur, and (iii) functional replacements involving eukaryotic genes are likely to have also contributed to the evolution of eukaryotes. Lastly, these data underscore the potential contribution of algal genes to the evolution of nonphotosynthetic lineages.     

Concerted gene recruitment in early plant evolution

Background  

Horizontal gene transfer occurs frequently in prokaryotes and unicellular eukaryotes. Anciently acquired genes, if retained among descendants, might significantly affect the long-term evolution of the recipient lineage. However, no systematic studies on the scope of anciently acquired genes and their impact on macroevolution are currently available in eukaryotes.     

Networks of gene sharing among 329 proteobacterial genomes reveal differences in lateral gene transfer frequency at different phylogenetic depths

Lateral gene transfer (LGT) is an important mechanism of natural variation among prokaryotes. Over the full course of evolution, most or all of the genes resident in a given prokaryotic genome have been affected by LGT, yet the frequency of LGT can vary greatly across genes and across prokaryotic groups. The proteobacteria are among the most diverse of prokaryotic taxa. The prevalence of LGT in their genome evolution calls for the application of network-based methods instead of tree-based methods to investigate the relationships among these species. Here, we report networks that capture both vertical and horizontal components of evolutionary history among 1,207,272 proteins distributed across 329 sequenced proteobacterial genomes. The network of shared proteins reveals modularity structure that does not correspond to current classification schemes. On the basis of shared protein-coding genes, the five classes of proteobacteria fall into two main modules, one including the alpha-, delta-, and epsilonproteobacteria and the other including beta- and gammaproteobacteria. The first module is stable over different protein identity thresholds. The second shows more plasticity with regard to the sequence conservation of proteins sampled, with the gammaproteobacteria showing the most chameleon-like evolutionary characteristics within the present sample. Using a minimal lateral network approach, we compared LGT rates at different phylogenetic depths. In general, gene evolution by LGT within proteobacteria is very common. At least one LGT event was inferred to have occurred in at least 75% of the protein families. The average LGT rate at the species and class depth is about one LGT event per protein family, the rate doubling at the phylum level to an average of two LGT events per protein family. Hence, our results indicate that the rate of gene acquisition per protein family is similar at the level of species (by recombination) and at the level of classes (by LGT). The frequency of LGT per genome strongly depends on the species lifestyle, with endosymbionts showing far lower LGT frequencies than free-living species. Moreover, the nature of the transferred genes suggests that gene transfer in proteobacteria is frequently mediated by conjugation.     

The Vein Patterning 1 (VEP1) gene family laterally spread through an ecological network

Lateral gene transfer (LGT) is a major evolutionary mechanism in prokaryotes. Knowledge about LGT--particularly, multicellular--eukaryotes has only recently started to accumulate. A widespread assumption sees the gene as the unit of LGT, largely because little is yet known about how LGT chances are affected by structural/functional features at the subgenic level. Here we trace the evolutionary trajectory of VEin Patterning 1, a novel gene family known to be essential for plant development and defense. At the subgenic level VEP1 encodes a dinucleotide-binding Rossmann-fold domain, in common with members of the short-chain dehydrogenase/reductase (SDR) protein family. We found: i) VEP1 likely originated in an aerobic, mesophilic and chemoorganotrophic α-proteobacterium, and was laterally propagated through nets of ecological interactions, including multiple LGTs between phylogenetically distant green plant/fungi-associated bacteria, and five independent LGTs to eukaryotes. Of these latest five transfers, three are ancient LGTs, implicating an ancestral fungus, the last common ancestor of land plants and an ancestral trebouxiophyte green alga, and two are recent LGTs to modern embryophytes. ii) VEP1's rampant LGT behavior was enabled by the robustness and broad utility of the dinucleotide-binding Rossmann-fold, which provided a platform for the evolution of two unprecedented departures from the canonical SDR catalytic triad. iii) The fate of VEP1 in eukaryotes has been different in different lineages, being ubiquitous and highly conserved in land plants, whereas fungi underwent multiple losses. And iv) VEP1-harboring bacteria include non-phytopathogenic and phytopathogenic symbionts which are non-randomly distributed with respect to the type of harbored VEP1 gene. Our findings suggest that VEP1 may have been instrumental for the evolutionary transition of green plants to land, and point to a LGT-mediated 'Trojan Horse' mechanism for the evolution of bacterial pathogenesis against plants. VEP1 may serve as tool for revealing microbial interactions in plant/fungi-associated environments.     

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.