Our understanding of the phylogenetic relationships among eukaryotic lineages has improved dramatically over the few past decades thanks to the development of sophisticated phylogenetic methods and models of evolution, in combination with the increasing availability of sequence data for a variety of eukaryotic lineages. Concurrently, efforts have been made to infer the age of major evolutionary events along the tree of eukaryotes using fossil-calibrated molecular clock-based methods. Here, we review the progress and pitfalls in estimating the age of the last eukaryotic common ancestor (LECA) and major lineages. After reviewing previous attempts to date deep eukaryote divergences, we present the results of a Bayesian relaxed-molecular clock analysis of a large dataset (159 proteins, 85 taxa) using 19 fossil calibrations. We show that for major eukaryote groups estimated dates of divergence, as well as their credible intervals, are heavily influenced by the relaxed molecular clock models and methods used, and by the nature and treatment of fossil calibrations. Whereas the estimated age of LECA varied widely, ranging from 1007 (943–1102) Ma to 1898 (1655–2094) Ma, all analyses suggested that the eukaryotic supergroups subsequently diverged rapidly (i.e., within 300 Ma of LECA). The extreme variability of these and previously published analyses preclude definitive conclusions regarding the age of major eukaryote clades at this time. As more reliable fossil data on eukaryotes from the Proterozoic become available and improvements are made in relaxed molecular clock modeling, we may be able to date the age of extant eukaryotes more precisely.
Fossil-calibrated molecular clock-based approaches are being used to estimate the age of major evolutionary events in eukaryotes. The estimates are heavily influenced by the methods and models used, but are steadily improving.
Our conception of the tree of eukaryotes has changed dramatically over the last few decades. In the 1980s and early 1990s, prevailing views were based on small subunit ribosomal RNA (SSU rRNA) gene phylogenies (e.g., Sogin 1991). However, as multiple protein-coding gene datasets were developed and more sophisticated phylogenetic methods were used, it became clear that the deep structure of the rRNA tree was the result of a methodological artifact known as long branch attraction (LBA) (Budin and Philippe 1998; Roger et al. 1999; Philippe et al. 2000a,b). Analyses based on multiple protein genes instead hinted at the existence of higher-level eukaryotic “supergroups” that encompassed both protistan and multicellular eukaryotic lineages (Baldauf et al. 2000). More recently, a better understanding of protistan ultrastructural diversity and the development of phylogenomic approaches have refined this picture and further delineated these groups (see also Fig. 1) (Bapteste et al. 2002; Burki et al. 2007; Hampl et al. 2009; Brown et al. 2012; Zhao et al. 2012).
Figure 1.Maximum likelihood phylogenetic tree of eukaryotes based on a phylogenomic dataset. Additional taxa were added to the original 159-gene Brown et al. (2013) dataset to maximize available fossil calibrations (total of 85 taxa, 43,099 sites). Black dots represent nodes on which fossil calibration constraints were imposed; yellow stars indicate the various positions of the root of the eukaryote tree considered; pink dots indicate the origin of major eukaryotic groups discussed here. A maximum likelihood (ML) phylogenetic tree was obtained from 60 heuristic searches using RAxML version 7.2.6 (Stamatakis 2006) under the Le and Gascuel (LG) + Γ +F amino acid substitution model (Le and Gascuel 2008). Numbers at nodes indicate bootstrap support (BS) for splits estimated from 500 bootstrap replicates. Most splits received maximum support and only BS < 100% is reported. Tree is shown rooted at the base of Amorphea, although roots at the base of either Obazoa or Excavates were also explored. Bayesian inference was also conducted using PhyloBayes 3.2 (Lartillot et al. 2009) by running four chains under either the CAT-GTR, CAT-Poisson, or the catfix C60-Poisson models of evolution, all combined with a gamma rates model. Bayesian calculations were not completed because of lack of convergence between chains, although the postburn-in consensus phylogeny from all runs was identical to the ML tree except for an unresolved multifurcation at the base of Excavata. Relaxed molecular clock (RMC) analyses were conducted with Phylobayes using the ML tree as a fixed topology. For all analyses, a birth–death tree prior was applied. Two chains were run until diagnostic statistics indicated convergence or estimated dates on nodes of interest for the two chains were <5% different. Fossil calibrations were taken from Parfrey et al. (2011) with the following modifications: four calibrations (Gonyaulacales, Spirotrichs, Foraminifera, Euglenids) were removed because of insufficient gene coverage within the clade of interest; the “Ciliate” calibration based on the tetrahymenol biomarker was removed (see text); as insufficient gene data was available from the haptophyte Isochrysis galbana, the upper bound on the coccolithophorid calibration was adjusted to an uninformative maximum (3000 Ma); the oldest cestode fossil (tapeworm) (Dentzien-Dias et al. 2013) was added as a calibration for Platyhelminths. The minimum age (250 Ma) was taken from the youngest possible age of the fossil and the upper boundary was set equal to the next-oldest calibrated node (Bilateria).
As our understanding of eukaryote phylogeny improved, fossil-calibrated molecular clock-based methods were beginning to be applied to date the major diversification events in this domain (Hedges et al. 2001; Douzery et al. 2004; Hedges and Kumar 2004; Berney and Pawlowski 2006; Parfrey et al. 2011). Molecular clock analyses were first introduced by Zuckerkandl and Pauling (1965). They showed that the differences between homologous proteins of different species are approximately proportional to their divergence time. Since then, sophisticated RMC methods have been developed that combine fossil data with molecular phylogenies for the inference of divergence times. However, attempts to estimate the age of deep divisions within eukaryotes using these methods have yielded vastly different estimates (e.g., see Douzery et al. 2004 vs. Hedges et al. 2004). These discrepancies can be explained by a myriad of sources of variability and error including (1) the assumed phylogeny of eukaryotes, (2) the sparse fossil record of protists and other organisms lacking hard structures for fossilization, (3) how fossil constraints are applied to phylogenetic trees, (4) methods and models used in RMC analysis, and (5) the selection of taxa and genes included.
Here, we review the progress and pitfalls in estimating the age of the last eukaryotic common ancestor (LECA) and supergroups using molecular clock-based analyses. We first discuss recent progress in our understanding of eukaryotic phylogeny and the ancient eukaryotic fossil record, and then we review the development of molecular clock-based methods and how fossil constraints are treated. Next, we describe attempts to date ancient eukaryotic divergences using RMC methods. Finally, we present an RMC analysis of a very large dataset comprised of 159 proteins and 85 taxa, using 19 fossil calibrations.
THE EUKARYOTIC TREE OF LIFEEstimates of divergence dates from molecular clock analyses are only meaningful if the phylogeny on which they are based is correct. However, recovering deep phylogenetic relationships among eukaryotes has proven to be an extremely challenging task.
The most recent conceptions of the eukaryotic tree of life feature five or six “supergroups” (Keeling et al. 2005; Roger and Simpson 2009; Burki 2014) including Opisthokonta, Amoebozoa, Excavata, the SAR group, Archaeplastida, and Hacrobia (Haptophyta and Cryptophyta). Whereas phylogenomic analyses robustly recover the monophyly of Opisthokonta, Amoebozoa and SAR, the phylogenetic coherence of Excavata, Archaeplastida, and Hacrobia is less certain (see also Fig. 1) (Burki et al. 2008, 2012; Hampl et al. 2009; Parfrey et al. 2010; Brown et al. 2012; Zhao et al. 2012; Burki 2014).
Another challenge inherent to dating ancient events in eukaryotic evolution is the current uncertainty regarding the location of the root (Stechmann and Cavalier-Smith 2002, 2003a,b; Cavalier-Smith 2010; Derelle and Lang 2012; Katz et al. 2012). Tackling this question is made especially difficult by the absence of closely related outgroups to eukaryotes. The large phylogenetic distance between sequences from eukaryotes and their archaeal or bacterial orthologs makes their use as outgroups highly prone to phylogenetic artefacts like LBA (Felsenstein 1978; Roger and Hug 2006). Consequently, various alternative methods have been used in the last decade, yielding different results regarding the placement of the root; between Amorphea (or “unikonts”) and all other eukaryotes (i.e., “bikonts”) (Richards and Cavalier-Smith 2005; Roger and Simpson 2009; Cavalier-Smith 2010), which, despite being the current leading working hypothesis, is challenged by several lines of evidence (Arisue et al. 2005; Roger and Simpson 2009): at the base of an Excavata lineage, Euglenozoa (Cavalier-Smith 2010); between Archaeplastida and other eukaryotes (Rogozin et al. 2009; Koonin 2010); and between Opisthokonta and other eukaryotes (Katz et al. 2012). Therefore, although the unikont/bikont root position remains the most popular, the evidence supporting it is not definitive and the debate remains open.
CALIBRATING ESTIMATES OF EVOLUTIONARY RATES: BIOMARKERS AND FOSSILSThe evolutionary distance between sequences is the product of their evolutionary rate and the time that passed since they diverged. Therefore, to estimate the molecular rate, and thus to infer the precise timing of an evolutionary event, it is necessary to calibrate the tree of life with “known” dates associated with the available paleobiological data. For ancient evolutionary events, calibrations are commonly based on the fossil record and, to a lesser extent, on biomarkers (organic molecules in the rock record that are characteristic of particular organismal groups) (for a more comprehensive discussion, see Knoll 2013). When considering fossil evidence for ancient eukaryotes, it is important to distinguish between “stem” and “crown” lineages. Briefly, all lineages that descend from LECA are known as crown eukaryotes, whereas stem groups are extinct lineages that emerged before LECA, but diverged from the eukaryote nucleocytoplasmic lineage after the eukaryote/archaeal split. Any characteristic feature of eukaryotes evident in ancient microfossils or biomarkers can, in principle, be a property of either stem or crown organisms. Therefore, unless there are specific features that definitively associate the fossils/biomarkers with particular crown eukaryote groups, they cannot be used as divergence time constraints in molecular clock analyses.
The oldest proposed biomarker evidence for the existence of eukaryotes is ∼2.7-Ga-old steranes (i.e., breakdown products of sterols) (Brocks 1999). It was later shown that these sterol biomarkers were younger than the rocks in which they were found (Fischer 2008; Rasmussen et al. 2008), although there have been new strictly controlled studies reporting steranes in similar aged sediments (Waldbauer et al. 2009). A bacterial origin of these biomarkers cannot be excluded as some bacteria make sterols. However, these differ from many that are characteristic of specific eukaryotes (Summons et al. 2006; Desmond and Gribaldo 2009). In general, biomarkers need to be considered carefully as contamination is difficult to rule out, and their specificity to one lineage cannot be guaranteed, especially for microbes in which lateral gene transfer is prevalent. For example, gammacerane (a breakdown product of the triterpenoid tetrahymenol) has been suggested to indicate the presence of ciliates (e.g., Tetrahymena) that produce tetrahymenol (Summons and Walter 1990). However, the gene responsible for tetrahymenol synthesis was shown to have been transferred between a wide range of microaerophilic eukaryotes (Takishita et al. 2012).
Fossilized features that are consistent with eukaryotic affinity include a combination of a large size, the presence of complex morphological features, wall ultrastructure or ornamentation, or typical excystment structures unknown in prokaryotic organisms (Javaux et al. 2001, 2003, 2004; Knoll et al. 2006; Strother et al. 2011). Many microfossils of suggested eukaryotic origin are large (>50 μm) organic-walled structures known as acritarchs (Buick and Young 2010; Javaux et al. 2010). The oldest of these are microfossils described from the ∼3.2-billion-year-old (Ga) Moodies group (Buick and Young 2010; Javaux et al. 2010), far older than most acritarchs described to date (Javaux 2007). Javaux and colleagues were very cautious in their interpretation and argue that a prokaryote origin could not be definitively excluded. The next oldest remains of possible eukaryotes are large spiral ribbon-shaped fossils identified as the most ancient putative representatives of Grypania spiralis and found in 1800–2100-million-year-old iron formations (Han and Runnegar 1992), although their eukaryotic origin has been challenged (Samuelsson and Butterfield 2001; for review, see Javaux et al. 2003; Porter et al. 2003; Knoll et al. 2006). Additionally, large acritarchs with possible eukaryotic attributes were discovered in 1.8-Ga formations from China (Zhongying 1986), although their eukaryotic affiliation is still debated (Javaux et al. 2003). In fact, most of the ancient Proterozoic assemblages (i.e., late Paleoproterozoic and Mesoproterozoic rocks, 1800–1000 million years [Ma]) include fossils that are difficult, if not impossible, to associate with crown group eukaryotes. One of the most convincingly “eukaryotic” of these is Shuiyousphaeridium macroreticulatum from the ∼1.7-Ga Ruyang group, an ornamented acritarch with ridged walls made up of regularly packed hexagonal plates (Pang et al. 2013). Although an affinity to dinoflagellates (Alveolata) has been suggested (Leiming et al. 2005), there is little evidence for a crown group affiliation and they could represent a stem eukaryote lineage (Javaux 2006). Another example are microfossils of Tappania plana from the early Mesoproterozoic (Yin 1997) that were proposed to be of fungal affinity (Butterfield 2005), but which are still viewed as uncertain because of the limited number of systematically informative characters (Knoll et al. 2006). A notable exception is Bangiomorpha pubescens from the 1.2-Ga hunting formation (Butterfield 2000) that represent the oldest fossils confidently assigned to a crown eukaryotic lineage, the bangiophyte red algae. Many unambiguously eukaryotic fossils from the Neoproterozoic (1000–543 Ma) are also difficult to assign to specific extant eukaryotic clades (Javaux et al. 2003). For example, although some vase-shaped microfossils found in >742-Ma rocks have been relatively confidently identified as arcellinids (Amoebozoa), others from this assemblage could correspond to either amoebozoan or euglyphid (Rhizaria) amoebae (Porter and Knoll 2000; Porter et al. 2003). Based on an apparent increase in the diversity of complex organic-walled microfossils, a number of paleontologists suggest that the diversification of most eukaryotic kingdoms occurred at ∼800 Ma (Porter et al. 2003; Knoll et al. 2006; Porter 2006; Javaux 2007; Knoll 2013).
APPLYING FOSSIL CONSTRAINTSThere has been much debate in the last decade over how fossil dates should be treated in the context of molecular clock analyses (Graur and Martin 2004; Hedges and Kumar 2004; Reisz and Müller 2004; Blair and Hedges 2005; Glazko et al. 2005; reviewed in detail by Parham et al. 2012; Ronquist et al. 2012; Warnock et al. 2012).
First, there is inherent uncertainty associated with the dating of the rocks in which the fossils are found. Second, a systematic bias is introduced by the fact that the true divergence date must be older than the age of the fossil itself, and the time gap between the two is often unclear (Hedges and Kumar 2004). Because genetic divergence precedes detectable morphological variation, genetic divergence times are commonly underestimated by paleontological evidence, leading to overestimates of molecular rates (van Tuinen and Hadly 2004; Near et al. 2005; Roger and Hug 2006; Ho et al. 2011). Consequently, how fossil calibrations are applied on a phylogenetic tree has a significant impact on age estimates, affecting the age of the deepest nodes of the eukaryotic tree by hundreds of millions of years (Hug and Roger 2007). In some cases, fossil calibration methods and the presence of particular fossil calibrations has a greater impact on age estimates than taxonomic sampling (Hug and Roger 2007).
There are several approaches to applying fossil constraints. “Hard constraints” (or “hard bounds”) treat calibration points as fixed and accurate time intervals (specified by the paleontological evidence) assigned to particular nodes on the tree (Kishino et al. 2001). Hard bounds disallow the estimated age of a constrained node to be outside the specified interval—an assumption that is hardly justified given the uncertainties discussed above. If the fossil of interest is accurately dated, it can provide a realistic lower bound on the age of a divergence, but rarely a good upper bound. To overcome this problem, extremely large hard upper bounds can be used, although this is likely to bias time estimates to be too old (Yang 2006). In contrast, the “soft bound” approach allows for a smoothly decreasing probability of the node age falling outside the interval (Drummond et al. 2006; Yang and Rannala 2006; Rannala and Yang 2007; Inoue et al. 2010). However, the nature of these probability distributions and how they are applied (i.e., with or without maximum constraints) can significantly alter the resulting estimates (Warnock et al. 2012). In the best-case scenario, several fossils are available that display apomorphies corresponding to different clades within a group, allowing for sequential constraints on minimum and maximum ages of the nodes in that part of the phylogeny. Unfortunately, the patchy and difficult-to-interpret ancient fossil record corresponding to the deeper eukaryote divergences is extremely limiting for the accuracy of this type of upper bound implementation. Some progress has been made in the recent years through the greater use of the taxonomically diverse older microfossil record that documents the appearance of a variety of protistan groups (Berney and Pawlowski 2006; Parfrey et al. 2011).
RMC METHODSMolecular dating approaches theoretically allow divergence times to be estimated from genetic distances. To do this, the phylogenetic tree is calibrated with one or more known dates, usually based on fossil or biomarker records, and divergence times are extrapolated throughout the tree. Originally, molecular dating relied on the assumption of a strict molecular clock postulating a constant rate of evolution over the whole tree (Zuckerkandl and Pauling 1965). However, variation in substitution rates has been widely documented (Smith and Peterson 2002; Bromham and Penny 2003; Davies et al. 2004). Consequently, RMC methods were developed that allow the rate of sequence evolution to vary across different branches (for reviews, see Welch and Bromham 2005; Lepage et al. 2007; Ho and Phillips 2009). A number of different methods have been developed and debate continues as to which best captures biological reality (Drummond et al. 2006; Lepage et al. 2007; Ho and Phillips 2009; Linder et al. 2011).
Local clock implementations estimate a separate molecular rate for each user-defined part of the tree (Yoder and Yang 2000). This requires arbitrary choices to be made regarding the number of different rate classes and how they should be assigned to branches; it is unclear how this can be achieved in a rigorous way in the absence of substantial prior information.
Rate-smoothing algorithms (Sanderson 1997, 2002) assume that the rate of evolution itself is evolving and is correlated across adjacent branches on the phylogeny such that related lineages have similar rates. For nonparametric rate smoothing (NPRS), rates are optimized to minimize a smoothing function that summarizes differences between rates on adjacent branches, taking into account fossil constraints (Sanderson 1997). The penalized-likelihood approach (Sanderson 2002), in contrast, combines a probabilistic model of sequence evolution with a penalty function (similar to the “smoothing function” of NPRS methods). These rate-smoothing methods require arbitrary decisions to be made regarding the parameters of the smoothing/penalty functions, and doubts persist regarding their statistical properties (Yang 2006).
The most sophisticated RMC methods are Bayesian approaches that probabilistically model every feature of molecular evolution over the tree, including the substitution process, tree generation processes, and substitution rate changes across branches (i.e., the RMC process). The various RMC process models can be divided in two main classes differing in whether temporal autocorrelation among adjacent lineages is assumed (Kishino et al. 2001; Aris-Brosou and Yang 2003; Rannala and Yang 2007) or not (Drummond et al. 2006; Akerborg et al. 2008). When rates are not autocorrelated, each rate is independent from the neighboring ones and randomly sampled from a probability distribution (e.g., the uncorrelated models introduced by Drummond et al. 2006). In autocorrelated models, rates follow a diffusion process along lineages, with the rate in each branch being drawn a priori from a parametric distribution whose mean is a function of the rate on the parent branch (e.g., the log-normal [LogN] distribution) (Thorne et al. 1998; Kishino et al. 2001), in which the evolutionary rate varies according to a lognormal distribution; or the Cox–Ingersoll–Rand distribution (CIR) (Lepage et al. 2007), which possesses a stationary distribution, contrasting with the linearly increasing variance of the LogN model.
RECENT ATTEMPTS TO DATE LECA AND MAJOR EUKARYOTIC GROUPSIn 2004, Hedges and colleagues analyzed multiprotein datasets with both constant rate and RMC methods to estimate the age of a number of deep eukaryote divergences (Hedges et al. 2004). After screening out proteins with detectable departures from the strict molecular clock, they obtained similar results with all methods, estimating LECA to be 2309 (2115–2503) Ma old, although several methodological aspects of this work have been debated (Graur and Martin 2004; Hedges and Kumar 2004; Roger and Hug 2006). Another potential problem with these analyses is that inferences regarding LECA were based on rooting the eukaryote tree on the lineage leading to the diplomonad Giardia. As discussed earlier, this rooting position is likely to be an artefactual result of LBA.
Douzery et al. (2004) estimated the age of LECA and other major groups using a Bayesian RMC approach calibrated with, for the first time, six paleontological constraints from the Phanerozoic. To reduce the impact of stochastic error resulting from the variation in evolutionary rates among individual genes, they analyzed a large supermatrix of 129 proteins and a tree rooted on the branch leading to Dictyostelium (close to the unikont/bikont root). They estimated a substantially younger age for LECA, dated at 1085 Ma (950–1259 Ma), and inferred that the subsequent diversification of all major groups occurred within ∼200 Myr.
To avoid some sources of error associated with the fossil record, Berney and Pawlowski (2006) used many more fossil constraints including Phanerozoic protist microfossils as a source of 26 calibration points (four maximum and 22 minimum time constraints). Their phylogeny, inferred from a single gene (SSU rRNA), was analyzed with fossil constraints in a Bayesian RMC framework. Based on a unikont/bikonts rooting, they estimated the age of LECA to be 1126 (948–1357) Ma. They concluded that most Proterozoic fossils suggested to be eukaryotes should not be assumed to belong to extant groups or used as calibration points for molecular dating analyses.
Most recently, Parfrey et al. (2011) performed Bayesian RMC analyses of a taxon-rich multigene dataset (91–109 taxa and 15 proteins) combined with 23 calibration points derived from Proterozoic and Phanerozoic eukaryotic fossils and biomarkers. They concluded that LECA was 1679–1866 Ma old. In contrast to Berney and Pawlowski (2006), Parfrey et al. (2011) argue that estimated ages of major groups are consistent with the tentative taxonomic assignments of many of the putative eukaryotic fossils from the Proterozoic.
ON THE DIFFICULTY OF DATING DEEP EVENTS IN EUKARYOTIC EVOLUTIONAs discussed in the foregoing sections, the discrepancies between date estimates from these previous RMC analyses can be attributed to a myriad of different causes. Roger and Hug (Roger and Hug 2006; Hug and Roger 2007) examined a number of these sources of variation for dates estimated from two published datasets (Douzery et al. 2004; Peterson and Butterfield 2005). They found that both the estimated ages of nodes and size of the confidence (or credible) intervals associated with them were extremely sensitive to the RMC method used, the number and nature of fossils chosen for constraints, and how these constraints were applied.
Here, we further explore sources of variation in molecular clock analyses focusing on the chosen root position for eukaryotes, the amino acid substitution model, the RMC model, and the manner in which calibration constraints are applied. In these analyses, we rely on a sophisticated Bayesian implementation of RMC models (Lartillot et al. 2009) and analyze a large phylogenomic dataset (159 proteins, 85 taxa). Nineteen calibration points were used, the majority of them taken from Parfrey et al. (2011) (for details, see Fig. 1).
Impact of Calibration ConstraintsFirst, we investigated the impact of the way of applying fossil calibrations by using both “hard-” and “soft-” bounds constraints. Strikingly, hard bounds gave notably older estimates for all nodes over a range of RMC models and root positions (9%–41% older) (compare filled and open shapes in Figs. 2 and 3). This suggests that at least one of the fossil calibrations is not consistent with the others (i.e., the estimated date of the fossil itself or its taxonomic affiliation/node assignment is incorrect). Closer inspection revealed that, under soft bounds, the fossil Bangiomorpha pubescens (assigned to the basal red algal split) was estimated to be 690 (639–771)–1026 (835–1206) Ma old, departing hugely from the age constraints applied (1174–3000 Myr). This was the only Proterozoic fossil in our analyses that was estimated to be younger than its lower age bound when soft-bounds constraints were applied. That the age or taxonomic assignment of Bangiomorpha is not consistent with other fossil dates in our analyses is similar to the observations of Berney and Pawlowski (2006). Similarly, Parfrey et al. (2011) found that excluding Proterozoic fossil constraints shifted the age of major groups to 200–300 Ma younger and also noted a marked change in the estimated age of red algae. Therefore, all recent molecular clock analyses seem to concur that the Bangiomorpha fossil date or taxonomic assignment is problematic. An alternative explanation is that the rates of evolution have dramatically shifted along the red algal lineage and this phenomenon is not well captured by relaxed clock models, in general.
Figure 2.Impact of the root position and calibration constraints on estimated divergence times. Estimated age (in Myr) of LECA and major eukaryotic groups when the root was placed at the base of Amorphea (squares), Obazoa (diamonds), and Excavata (circles), using hard- (filled shapes) and soft- (open shapes) bounds constraints. Error bars represent 95% credible intervals. All estimations were performed using the C60 substitution model and the uncorrelated γ distribution of rates (UGam) clock model.
Figure 3.Impact of the molecular clock models and calibration constraints on estimated divergence times. Estimated age (in Myr) of LECA and major eukaryotic groups using the UGam (squares), LogN (diamonds), and CIR (circles) clock models, applying hard- (filled shapes) and soft- (open shapes) bounds constraints. Error bars represent 95% credible intervals. All analyses are based on the Amorphea root, using the C60 substitution model.
Impact of the Root PositionDuring our analyses, we considered three possible positions for the eukaryote root. We first placed the root between Amorphea and the remainder of eukaryotes (i.e., unikont-bikont rooting). We tested a second rooting position between the Obazoa (Opisthokonta + apusomonads + breviates) and all other taxa—a variation of the hypothesis by Katz et al. (2012) in which they located the root at the base of Opisthokonta. This modification is based on recent work presenting evidence for the monophyly of Obazoa (Brown et al. 2013a). Finally, we tested a root at the base of Excavata as many excavate protists (e.g., diplomonads [Hedges et al. 2004], Euglenozoa [Cavalier-Smith 2010], or jakobids [Brinkmann et al. 2007]) have variously been suggested to be “basal” eukaryote lineages.
Interestingly, these three alternatives for the position of the root did not have much impact on the estimated age of the root itself and the supergroups under the uncorrelated gamma (UGam) RMC model (this model setting is discussed further below) (Fig. 2). For fixed model settings, the three tested positions gave similar date estimates for a given specific timepoint and varied by less than 8%, with the 95% credible intervals largely overlapping (Fig. 2) (compare the variation among filled shapes and open shapes). Berney and Pawlowski also found that four different placements of the root (at the base of unikonts, Opisthokonta, Amoebozoa, and Excavata, respectively) yielded similar time estimates (Berney and Pawlowski 2006). Parfrey et al. (2011) also reached similar conclusions for the unikont, Opisthokonta, Excavata, and Excavata + unikont roots.
For all sets of parameters, our results show that Opisthokonta is estimated as being the youngest of all major groups. Nevertheless, the estimated age difference between LECA and the last common ancestor of Opisthokonta is relatively small, and ranges between 23 and 334 Ma with LECA estimated to be between 1007 and 1898 Ma old, and Opisthokonta between 904 and 1579 Ma old (Fig. 3, filled circles and open diamonds). This suggests that the six groups considered here diversified over a relatively short period of time after LECA and explains why the various placements of the root in this region of the tree do not have a major impact on its estimated age. This observation is in agreement with the “big-bang hypothesis” for eukaryotic evolution, according to which major eukaryotic groups emerged rapidly, virtually leading to a massive multifurcation (Philippe et al. 2000a,b). However, caution is warranted as problems with saturation of sequence changes, model misspecification, and other sources of conflict within the data could also lead to the estimation of extremely short internal branches deep in the eukaryote tree, artefactually generating a multifurcating topology (see Roger and Hug 2006). It should be noted that the roots considered here and in the other aforementioned studies (Berney and Pawlowski 2006; Parfrey et al. 2011) are all topologically “close”; we expect that more distantly placed putative roots (i.e., well within one of the supergroups) would give significantly different estimates for the age of LECA, although we did not test this in our analyses.
Impact of the Substitution Model and Relaxed Clock ModelTo test the impact of the substitution and the relaxed clock models, we estimated the age of each major group by applying six model combinations (Fig. 4).
Figure 4.Impact of the substitution and molecular clock models on estimated divergence times. Estimated age (in Myr) of LECA and major eukaryotic groups using the UGam (squares), LogN (diamonds), and CIR (circles) clock models. Substitution model was either C60 (filled shapes) or LG (open shapes). Error bars represent 95% credible intervals. All estimations were calculated with the Amorphea root, using soft bounds calibration constraints.
The two substitution models used were either the site-heterogeneous empirical profile mixture model C60-Poisson (C60) (Le et al. 2008) or the more classical site-homogeneous LG substitution matrix (Le and Gascuel 2008). The three RMC models included the UGam and the autocorrelated LogN and CIR models.
Isolating the impact of the substitution model and the RMC model turns out to be a difficult task. C60 + UGam and LG + UGam both gave similar date estimates (difference <7%) (Fig. 4, filled and empty squares) as did C60 + LogN and LG + LogN (<2%) (Fig. 4, filled and empty diamonds). In contrast, C60 + CIR yielded very different results from LG + CIR (27%–32% difference) (Fig. 4, filled and empty squares).
However, comparing the estimated dates under the LG substitution model showed that they were similar when the LogN or CIR relaxed clock models were used (1%–6% difference between LG + LogN and LG + CIR) (Fig. 4, open diamonds and circles), whereas LG + UGam gave much older estimates (difference up to 30% with LG + LogN and LG + CIR) (Fig. 4, compare open squares with open diamonds and circles, respectively). Finally, C60 + UGam and C60 + CIR models led to similar date estimates (1%–8% difference) (Fig. 4, filled squares and circles), whereas C60 + LogN gave much younger ages (22%–36% younger) (Fig. 4, compare filled diamonds with filled squares and circles).
The complexity of the relationship between these parameters is further amplified by the choice of hard- versus “soft-bounds calibration constraints. All previous comparisons used soft bounds and application of hard bounds led to very different conclusions. For example, with hard bounds, the C60 + UGam and C60 + LogN models gave similar estimates (<5% difference) (Fig. 3, filled squares and diamonds), that were drastically younger than estimates of the C60 + CIR model (19%–23%) (Fig. 3, compare filled squares and diamonds with filled circles). In contrast, with soft bounds, the C60 + CIR model gave estimates similar to the C60 + UGam model (<8% variation) (Fig. 3, open circles and squares), whereas C60 + LogN yielded much younger dates (24%–36%) (Fig. 3, compare open circles and squares to open diamonds).
Overall, the largest differences among estimates is observed between the CIR clock model with hard bounds calibrations combination (which gave the oldest dates) (Fig. 3, filled circles), and the “LogN + soft bounds” combination (estimating the youngest ones) (Fig. 3, open diamonds). Using these two settings, the eukaryotic root was estimated to be 1898 versus 1007 Ma old, respectively. Because the soft-bounds approach is a more realistic way to treat fossil constraints, the younger date estimates obtained under these constraints are likely more reasonable. Restricting attention to soft-bounds analyses, combinations of RMC models and calibration constraints gave estimates for LECA between 1007–1456 Ma (Fig. 3, open shapes). It is important to note, however, that credible intervals on the age of LECA and most other nodes can often be quite large, sometimes spanning up to ∼650 Ma (e.g., CIR model) (Fig. 3).
MODEL COMPARISON AND ASSESSMENTAs suggested by previous studies and the foregoing analysis, the various features of RMC analyses interact in a complex manner in influencing the age estimates obtained. This inevitably leads to the question of which models or settings should be trusted over others. Generally, this involves model selection via the assessment of the relative fit of various models used in the analysis (e.g., the substitution or molecular clock models); better fitting models should, ceteris paribus, provide better estimates. Therefore, we suggest that model selection techniques be used to choose among alternative RMC and substitution models. Below, we briefly review the various model selection techniques in the Bayesian context, the statistical inference paradigm in which most RMC models have been implemented to date.
The Bayes factor is the most widely used Bayesian model selection criterion and is defined as the ratio of the marginal likelihoods under the two models of interest (Linder et al. 2011). Bayes factors measure how strongly the data support a given model and have the advantages of implicitly penalizing more complex models as well as allowing for a general comparison among models that are not necessarily nested (Lepage et al. 2007). Unfortunately, in the phylogenetic context, exact calculation of Bayes factors is computationally infeasible and approximation techniques must be used. Although the harmonic mean estimator has been used for this purpose, it turns out to be extremely unstable (Lartillot and Philippe 2006). More reliable estimators are complex approximation methods such as thermodynamic integration (path-sampling) (Ogata 1989; Lartillot and Philippe 2004, 2006; Blanquart and Lartillot 2006; Rodrigue et al. 2006) and the stepping-stone method (Xie et al. 2011) that have been implemented in programs like Phylobayes (Lartillot et al. 2009) and BEAST (Baele et al. 2012), respectively. Note that these methods are generally used for RMC model selection but can be difficult to implement for selecting among substitution models.
Several other methods have been developed for substitution model selection in a Bayesian context. Lartillot and colleagues have developed cross-validation (CV) (Lartillot et al. 2007) methods that calculate the ability of one part of the data (the learning set) to predict the remaining data (the test set) after the dataset is randomly split into two uneven parts. In Phylobayes, Markov chain Monte Carlo (MCMC) is run on the learning set and, after burn-in is discarded, the average likelihood of the topology/parameters/model visited in the chain is evaluated in the test set. The procedure is repeated 10 times for each model, and the model with the greatest average likelihood over the test sets is the one selected as better fitting. CV analyses virtually always select site-heterogeneous substitution models (e.g., CAT or C60) over the simpler single-matrix models (e.g., LG) (Lartillot et al. 2009).
Posterior predictive simulation is another way of evaluating relative fit of substitution models (Bollback 2002; Nielsen and Huelsenbeck 2002; Lartillot and Philippe 2004; Blanquart and Lartillot 2008). From each postburn-in MCMC sample, a dataset is simulated from the model parameters. A statistic measuring some property (e.g., mean number of distinct residues observed at each column of the alignment) can be calculated for the real dataset and compared with the distribution of this statistic from the simulated datasets (e.g., Lartillot et al. 2009). A significant deviation between the statistics evaluated on real and simulated data implies that some of the model assumptions are unrealistic, indicating misspecification. The model in which the statistic evaluated on the simulated distribution falls closest to that of the real dataset is then considered the better fitting model.
FOSSIL ACCURACY ASSESSMENTCV procedures were also developed to assess consistency between fossil and molecular age estimates (Near and Sanderson 2004; Near et al. 2005; Rutschmann et al. 2007). They aim to detect point calibrations that poorly predict other calibration dates in the dataset, and are therefore suspected to be erroneous. The main pitfall of these approaches stems from their tendency to discard “outliers” to obtain an internally consistent set without taking into account the relative credibility of the corresponding fossils. Another weakness comes from the fact that they consider point calibrations and not minimum age constraints (Parham and Irmis 2008).
As point calibrations were progressively being replaced by distributions that better represent palaeontological uncertainty (Drummond et al. 2006; Yang and Rannala 2006), Sanders and Lee (2007) proposed to test for calibration accuracy in a Bayesian context through the comparison of the posterior age distribution of a calibration node with the prior specified for it.
CONCLUSIONMolecular clock estimates of ancient divergence times in the tree of life are affected by numerous sources of errors and uncertainties. Although resolution in the tree of eukaryotes appears to be steadily improving, the location of the root (i.e., LECA) remains uncertain. Moreover, controversy in assigning Proterozoic fossils to extant eukaryote groups means that molecular clock analyses must rely heavily on extrapolation from the younger, but richer, Phanerozoic fossil record. There are also inherent biases and uncertainties associated with assigning fossil calibrations to nodes in molecular phylogenies. These factors, combined with variability in estimates and credible intervals yielded by different molecular clock model assumptions, have led to the wide ranges of estimated ages of LECA and major eukaryote supergroups published in the last decade. The analyses of a taxon-rich supermatrix with 19 fossil calibrations presented here provide estimates for the age of LECA in the range of ∼1000–1900 Ma depending on the methods used, with the credibility intervals on estimates in some cases spanning up to ∼650 Myr. Despite this uncertainty about precise ages, both our and other recent molecular clock analyses recover a relatively short time interval (<300 Myr) between the age of LECA and emergence of all of the various supergroups, consistent with a rapid big-bang diversification of eukaryotes (Philippe et al. 2000a).
It is unclear how much more precision we will be able to achieve in dating the extant eukaryote clade. For molecular phylogenetics, difficult to resolve parts of trees can be clarified by the addition of more taxa and more genes (sites) because, in theory, likelihood-based methods are statistically consistent (see discussion in Yang 2006). However, this is not the case for molecular clock analyses because inherent uncertainties in fossil ages and biases in their node assignments imply that molecular clock age estimates will always be associated with error. Nevertheless, the addition of more fossils and genes may help narrow confidence or credible intervals, and steady progress is being made in the development of better RMC methods, model selection, and fossil assignment validation. These methods, combined with novel sources of time constraints such as those provided by horizontal or endosymbiotic gene transfer (see Shih and Matzke 2013), may ultimately help in the quest to determine the age of extant eukaryotes.
ACKNOWLEDGMENTSThis work and M.W.B.'s postdoctoral fellowship is supported by Discovery grant 227085-11 and an Accelerator grant from the Natural Sciences and Engineering Research Council of Canada (NSERC) awarded to A.J.R. L.E. is supported by a Centre for Comparative Genomics and Evolutionary Bioinformatics postdoctoral fellowship from the Tula Foundation; S.C.S. is supported by graduate scholarships from NSERC and Killam trusts. We thank Nicolas Lartillot, Laura Wegener Parfrey, and Andrew H. Knoll for very stimulating discussions and insights. Computations were partially performed on the supercomputers at the SciNet HPC Consortium. SciNet is funded by the Canada Foundation for Innovation under the auspices of Compute Canada, the Government of Ontario, Ontario Research Fund—Research Excellence, and the University of Toronto (Loken et al. 2010).
FootnotesEditors: Patrick J. Keeling and Eugene V. Koonin
Additional Perspectives on The Origin and Evolution of Eukaryotes available at www.cshperspectives.org
REFERENCES*Reference is also in this collection.
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4
