The timing of divergences among metazoan lineages is integral to understanding the processes of animal evolution, placing the biological events of species divergences into the correct geological timeframe. Recent fossil discoveries and molecular clock dating studies have suggested a divergence of bilaterian phyla >100 million years before the Cambrian, when the first definite crown-bilaterian fossils occur. Most previous molecular clock dating studies, however, have suffered from limited data and biases in methodologies, and virtually all have failed to acknowledge the large uncertainties associated with the fossil record of early animals, leading to inconsistent estimates among studies. Here we use an unprecedented amount of molecular data, combined with four fossil calibration strategies (reflecting disparate and controversial interpretations of the metazoan fossil record) to obtain Bayesian estimates of metazoan divergence times. Our results indicate that the uncertain nature of ancient fossils and violations of the molecular clock impose a limit on the precision that can be achieved in estimates of ancient molecular timescales. For example, although we can assert that crown Metazoa originated during the Cryogenian (with most crown-bilaterian phyla diversifying during the Ediacaran), it is not possible with current data to pinpoint the divergence events with sufficient accuracy to test for correlations between geological and biological events in the history of animals. Although a Cryogenian origin of crown Metazoa agrees with current geological interpretations, the divergence dates of the bilaterians remain controversial. Thus, attempts to build evolutionary narratives of early animal evolution based on molecular clock timescales appear to be premature.

Here we show that the precision of molecular clock estimates of times has been grossly over-estimated. Using a Bayesian method to estimate the timescale of metazoan diversification, we performed sensitivity analyses to explore the impact of the different sources of uncertainties. We used a large amino acid alignment (38,577 sites) of 203 nuclear encoded proteins for 71 species (based on []). We employed four fossil calibration strategies that accommodate different interpretations of the fossil record and show that these have a dramatic impact on the estimated times. We also explored the use of different relaxed clock models and show that at this level of divergence the molecular clock is significantly violated. We tested for the effects of different data partitioning strategies and show that this, too, has a significant impact on divergence time estimates. Finally, we show that competing phylogenetic hypotheses yield very different divergence time estimates. An evolutionary timescale for metazoan diversification that accommodates these uncertainties has precision that is insufficient to discriminate among causal hypotheses. Though some of this uncertainty can be reduced through increased precision of calibrations afforded by statistical modeling of fossil occurrence, more sequence data, reduced topological uncertainty, etc., the limitations of the fossil record and the confounding effect of times and rates will remain, making it difficult to achieve the precision required to test competing hypotheses on the causes and consequences of metazoan diversification.

It has been hoped that these questions may be answered, and a timescale for animal evolutionary history established, using molecular clock methodology. Indeed, there is a long history of attempts to estimate the timing of animal diversification [], yielding ages for crown Metazoa that range between 1,298 Ma [] and 615 Ma []. The disparity between molecular clock estimates and fossil evidence of clade age minima has diminished in association with the development of molecular clock methodology, particularly in accommodating rate variation. Molecular clock timescales are considered good enough by some to begin to synthesize evolutionary narratives integrating evidence of extrinsic environmental and ecological evolution from the geological record with intrinsic biological evolution []. Most recent divergence time analyses have been undertaken within a framework of Bayesian inference because it is capable of integrating much of the uncertainty associated with divergence time estimation, viz. the relationships between fossil evidence and clade age, rate variation among lineages (the relaxed clock), branch length estimation, tree topology, and parameters such as data partitioning. Yet, few studies have considered the cumulative impact of these uncertainties on the precision of evolutionary timescales (e.g., []).

Unequivocal fossil evidence of animals is limited to the Phanerozoic. Older records of animals are controversial: organic biomarkers indicative of demosponges [] are apparently derived ultimately from now symbiotic bacteria []; putative animal embryo fossils [] are alternately interpreted as protists []; and contested reports of sponges [], molluscs [], and innumerable cnidarians [], as well as putative traces of eumetazoan or bilaterian grade animals [], all from the Ediacaran. Certainly, there are no unequivocal records of crown-group bilaterians prior to the Cambrian [], and robust evidence for bilaterian phyla does not occur until some 20 million years into the Cambrian []. There is, nevertheless, increasingly general acceptance of a Precambrian history to animal evolution, and it is only its extent that remains open to debate. Was there an explosive radiation of bilaterian phyla close to the base of Cambrian []? Or is there an extensive Precambrian bilaterian history that extends deep into the Cryogenian [], the absence of a fossil record merely reflecting preservation, collection, and/or interpretation biases?

When life got smart: the evolution of behavioral complexity through the Ediacaran and Early Cambrian of NW Canada.

The timing and tempo of the evolutionary emergence of animal biodiversity has been among the most enduring problems in evolutionary biology. Innumerable hypotheses have been proposed to explain how the transition to multicellularity was effected, why it occurred when it did, and why it did not occur much earlier in Earth history []. Much of the molecular genetic toolkit required for animal development originated deep in eukaryote evolutionary history [], and it has been widely held that the emergence of complex multicellular organisms was precluded until the oxygenation of the biosphere []. Other potential, but not necessarily mutually exclusive, triggers for animal diversification include the release of their forebears from the environmental strictures of the Cryogenian or Ediacaran Snowball Earth [] and the effects of cosmic radiation [], polar wander [], continental fragmentation [], HS toxicity [], salinity [], a scarcity of trace metal micronutrients [], a pulse of continental weathering yielding nutrients to the oceans [], global warming [], or an escalatory predator-prey arms race []. Although these hypotheses propose more or less proximal causal mechanistic relationships with metazoan diversification, they rely ultimately on presumed temporal coincidence. This is challenging since the timing and the threshold of extrinsic environmental effects is invariably as unclear as the timing of the phenomena that they have been invoked to explain (e.g., []), varying from the origin of animals, eumetazoans, or bilaterians, to the origin of the animal phyla themselves, their crown radiations, or the sum total of this diversity. Indeed, it has been argued that the oxygenation of global oceans was a consequence, not a cause, of metazoan evolution []. Reconciling these competing hypotheses requires calibration to a common absolute timescale.

Temperature and salinity history of the Precambrian ocean: implications for the course of microbial evolution.

The Nebula Winter: the united view of the snowball Earth, mass extinctions, and explosive evolution in the late Neoproterozoic and Cambrian periods.

All of the preceding analyses employed a fixed tree topology ( Figure 1 ), yet the phylogenetic position of some metazoan taxa remains the subject of debate []. To account for this uncertainty, we analyzed 161 alternative binary trees, accounting for uncertainties in the positioning of Bilateria, chaetognaths, molluscs, nematodes, and xenacoelomorphs. The results of these analyses show that nodes are affected differently depending on the tree topology. For example, some nodes are characterized by time estimates that remain similar across all topologies ( Figure 5 ). These nodes are usually well calibrated and/or the local phylogeny well accepted, such as in crown deuterostomes and arthropods ( Figure 5 ). In contrast, nodes with uncertain phylogenetic relationships exhibit considerable variation in estimated ages. These include the nodes close to the root of the tree, such as Metazoa, Bilateria, and Cnidaria; this variation increases with proximity to the root. For example, moving the position of Placozoa around the eumetazoan node has a profound impact on the estimated age of the root ( Figure 5 ).

Four nodes (A–D) can be rearranged in three different ways (1–3), and a fifth node (E) can be rearranged in two ways, resulting in a total of 162 tree topologies reflecting the uncertain relationships around these five nodes. Divergence times were estimated using strategy 1, the IR model, and a single partition using each tree (bottom panel). Some phylogenetic hypothesis had a strong effect on posterior mean times; for example, placing the Placozoa as the most basal with respect to Cnidaria and Bilateria (A) leads to substantially older divergence times for the Metazoa (bottom panel), whereas placing Cnidaria as the most basal leads to substantially older times for the divergence of Eumetazoa.

The 95% HPD width is plotted against the mean of the divergence times estimated without molecular data (prior) and with the 203 gene alignments divided into one, two, four, five, and ten partitions. This plot indicates how much of the uncertainty in the posterior time estimates is due to the uncertain fossil calibrations and how much is due to the limited amount of sequence data. Thus, the low correlations indicate that the limited amount of sequence data contributes substantially to posterior uncertainty and the regression coefficients also indicate that the fossil calibrations involve much uncertainty.

Partitioning of the molecular sequence alignment may impact on divergence time estimates []. To explore this, the protein alignment was divided into two, four, five, and ten partitions, according to the relative amino acid substitution rates among genes (see the Supplemental Experimental Procedures ). The posterior mean times for the most ancient nodes tended to increase as the number of partitions increases ( Figure 2 C). For example, divergence time estimates for crown Metazoa vary from 833–681 Ma (single partition) to 834–787 Ma (ten partitions; Table 1 ). The discrepancy between age estimates increases with proximity to the root, regardless of whether or not the nodes are calibrated ( Figure 2 C). Age estimates on intermediate nodes (e.g., all vertebrates and most arthropod nodes) do not vary significantly with partition strategy; for a small number of nodes, younger date estimates were obtained when more partitions were used ( Figure 2 C and Table 1 ). Overall, nodes with highly variable time estimates among different partitions are those without calibration or are close to the root, where the calibrations are least informative ( Table S4 and Figure S2 ).

Which clock model should be used? Bayes factors have been used to decide between competing clock models (such as the IR and AR models, e.g., []) in phylogenetic analysis. MCMCTree does not yet implement Bayes factors, and so we did not calculate them here. Further work will be required to assess the suitability of the various clock models to describe rate evolution in the metazoan phylogeny. Thus, in attempting to encompass the uncertainty in the rate drift model, we consider here the spread of node age estimates that arise from both rate models.

The AR model penalizes extreme rate variation over short time intervals and effectively imposes local clocks for closely related species while allowing large rate variation among distant clades. This contrasts with the IR model, which assumes that the variance of the rate is independent of the divergence time, so that the variance is the same whether the species are closely or distantly related. Figure 3 shows the change in the shape of the log-normal distribution of rates under the AR model across 500 million years of evolution and highlights the extreme level of rate variation in Metazoan phylogeny. At short timescales, the distribution is more symmetrical and has a smaller variance than at longer timescales. In the case of the IR model with μ = 0.089/100 million years and σ= 0.468/100 million years, the log-normal distribution has the same shape as that for 100 million years for the AR ( Figure 3 , third plot).

In the AR model, the rates at the tips of a star phylogeny are log-normally distributed with mean r A (the ancestral rate at the root) and log-variance of the rate σ 2 = tν. For the metazoan phylogeny, the posterior mean of r A is 0.089 s/s/100 million years and of v is 0.468/100 million years. In (A)–(F), the evolution of the rate of molecular evolution is shown through 500 million years of metazoan history assuming the AR model to be correct. The numbers in brackets are the 95% equal-tail range of the distribution of the rate for the given time. As the star phylogeny evolves, the variance of the rates increases exponentially. After 500 million years of evolution, the 95% equal-tail range encompasses two orders of magnitude. Note that in case of the IR model with μ = 0.089/100 million years and σ 2 = 0.468/100 million years, the shape of the log-normal distribution is the same as that for 100 million years for the AR at any time point. Note that here the timescale is given in million years from the root (i.e., 0 million years is the root, and 500 million years is present time).

To examine the impact of this uncertainty, we re-estimated the divergence times of metazoans assuming an autocorrelated rates (AR) model [] under calibration strategy 1. This relaxed-clock model imposes a correlation of rates between ancestral and descendant branches by modeling rate change on the tree as a geometric Brownian diffusion process []. We found that the choice between AR versus IR relaxed-clock models has a strong impact on the estimated divergences ( Table 1 and Figure 2 B). Our results show that many posterior time estimates for young nodes using the AR model are older than those derived using the IR model, whereas a few nodes, especially the deep nodes, are younger ( Table 1 and Figure 2 B). In particular, the divergences of crown Metazoa (764–650 Ma), crown Bilateria (619–596 Ma), crown Deuterostomia (611–587 Ma) and crown Protostomia (599–578 Ma) are substantially younger.

The posterior mean times estimated under calibration strategy 1, independent rates (IR) model, and a single partition are plotted against (A) estimates using strategy 2, (B) estimates under the autocorrelated rates (AR) model, and (C) estimates obtained when the 203 gene alignments are divided into ten partitions according to substitution rate. The bars indicate the 95% HPDs.

When rate variation across a phylogeny is extreme (that is, when the molecular clock is seriously violated), the rates calculated on one part of the phylogeny will serve as a poor proxy for estimating divergence times in other parts of the tree. In such instances, divergence time estimation is challenging and the analysis becomes sensitive to the rate model used.

Calibration strategies 1–4 are based on a protostome interpretation of the Ediacaran Kimberella (552.85 Ma), to constrain the minimum time of divergence of Protostomia, Bilateria, Eumetazoa, and Metazoa ( Table 1 ). However, to some, there is no unequivocal fossil evidence of metazoans prior to the Cambrian. In this view, interpreting Kimberella as a protostome leads to unduly ancient estimates for the origin of all the more universal clades encompassing Protostomia. To assess the impact of using Kimberella as a minimum constraint on the age of the protostome clade, we employed a variation of calibration strategy 1 in which the next-oldest record of Protostomia and oldest unequivocal total-group mollusc, the Cambrian Aldanella yanjiahensis (532 Ma), was used in place of Kimberella. The resulting divergence time estimates are effectively the same as those derived using strategy 1 ( Figure S3 ). Thus, even under the assumption that the fossil record of metazoan is limited to the Cambrian, our estimates require an Ediacaran origin for most crown-bilaterian phyla, a late Cryogenian-early Ediaracan origin of crown Bilateria, and an early Cryogenian origin of crown Metazoa.

Age estimates for the younger nodes are similar under all four calibration strategies (e.g., nodes 68, 86, and 92; Table S3 ). However, the posterior age estimates of nodes close to the root exhibit dramatic differences among the different calibration strategies (e.g., Figure 2 A). This appears to reflect a paucity of palaeontological evidence, requiring very different scenarios for the timing and tempo of metazoan diversification. Strategies 3 and 4 yield timescales that strongly favor an early Cryogenian (834–780 Ma) diversification, evidently constrained by the root age, while the age estimates arising from calibration strategies 1 and 2 are compatible with metazoans diversifying at any time within the Cryogenian, though these analyses are not otherwise very informative ( Figure 1 D).

Calibration strategy has a large impact on estimated divergence times ( Figure 1 A, Table S3 , and Figure S1 ). Estimates under strategy 1 indicate that Metazoa originated 833–681 Ma, Bilateria 638–615 Ma, Deuterostomia 628–594 Ma, and Protostomia 626–598 Ma ( Table S3 ). When the skew-normal distribution is employed (which places the majority of its probability mass near the minimum age bound; strategy 2), the resulting posterior time estimates agree largely with those obtained using the uniform prior time distribution of strategy 1 ( Table 1 ). In contrast, calibration densities modeled with the Cauchy distribution (strategies 3 and 4) exhibit strong truncation effects in the time priors ( Figures 1 B and 1C), resulting in substantially older time estimates ( Figure 1 D). This can be seen, for example, in association with crown Bilateria, Deuterostomia, and Protostomia, where truncation caused the effective priors to place considerable probability mass beyond the maximum bound of 636.1 Ma ( Figure 1 C). This differs significantly from the specified calibration densities (cf. Figure 1 B), resulting in posterior time estimates that are substantially older than those derived using strategies 1 and 2 ( Figure 1 D). For example, estimates under strategy 3 indicate Metazoa originated 834–795 Ma, Bilateria 759–685 Ma, Deuterostomia 722–644 Ma, and Protostomia 712–644 Ma ( Table S3 , cf. strategy 1) Thus, truncation can have dramatic and perhaps surprising effects. These effects may be hard to predict, highlighting the challenges in constructing fossil calibrations, as calibrations based on the same fossil information can unintentionally lead to dramatically different priors and posterior estimates of divergence times.

The program MCMCTree [] was used to obtain posterior time estimates under these four strategies and on the fixed tree topology of Figure 1 . The evolutionary rates on branches of the tree were assumed to vary independently among lineages (the independent rates [IR] model []). All gene alignments were concatenated and analyzed as a single partition (1P) under the LG + Γ amino acid substitution model. In all instances, we first ran the analyses without sequence data to establish the effective time prior. This allowed us to evaluate the impact of truncation, which can yield marginal time priors that differ considerably from the original fossil evidence [].

(B–D) Calibration, prior, and posterior densities for four ancient nodes in the metazoan phylogeny; coloring relates to the calibration strategy employed as in (A). (The phylogeny with species names is provided in Figure 6 .)

(A) Time trees showing posterior divergence time estimates for major metazoan groups. Nodes are drawn at the posterior means obtained and horizontal bars represent 95% HPD intervals. Estimates were obtained with MCMCTree using the LG + Γ 4 + F model, IR, and with the 203 proteins concatenated into a super alignment.

To assess the robustness of estimated Metazoan divergences to calibration choice, we established temporal constraints on the ages of 34 clades based on fossil evidence ( Table 1 ). These were used as the basis for four competing sets of calibration densities, reflecting different interpretations of the fossil evidence ( Table S1 ). Under strategies 1 and 2, the age of crown Metazoa has the minimum constraint based on a protostome interpretation of the Ediacaran Kimberella, whereas in strategies 3 and 4 it is based on the disputed biogeochemical evidence of Cryogenian demosponges []. In strategy 1, all 34 calibrations were modeled as uniform distributions with soft bounds []. In strategies 2–4, we used different calibration densities for 14 phylum and superphylum crown nodes. In strategy 2, these 14 nodes are modeled using a skew-normal distribution with the mode of the distribution near the minimum bound and the tail extending into the past. These calibrations represent an optimistic interpretation of the fossil minima as a close approximation of the true clade age. In strategies 3 and 4, the 14 nodes use truncated Cauchy distributions [] with either a long tail (strategy 3) or a short tail (strategy 4), extending back in time. This represents a pessimistic interpretation of palaeontological evidence in which the first fossil records of clades are a poor approximation of their antiquity. Note that the Cauchy is a heavy-tailed distribution, that is, it places considerable probability mass on its tail (contrary to the skew normal, which is light tailed). The calibration-based time prior is shown in Table S2

Nodes are numbered as in Figure 6 . Note: posterior times are the 95% highest probability density (HPD) interval, estimated with MCMCTree v4.8 under the LG + Γ+ F model. S1, calibration strategy 1; S2, strategy 2; IR, independent rates model; AR, autocorrelated rates model; 1P, the 203 proteins analyzed as a single partition; 10P, the proteins are grouped into ten partitions according to their evolutionary rates. Nodes in bold have calibrations that differ in strategy 1 and strategy 2. Composite: 95% confidence interval (CI) is a composite of the 95% CI across all analysis, except those under strategy 3 and strategy 4 and under alternative topologies.

Minimum and Maximum Fossil Constraints and 95% HPD Limits of Posterior Divergence Times for Various Metazoan Clades, in Millions of Years before Present

Table 1 Minimum and Maximum Fossil Constraints and 95% HPD Limits of Posterior Divergence Times for Various Metazoan Clades, in Millions of Years before Present

Fossil calibrations are affected by numerous sources of uncertainty, including phylogenetic interpretation, dating of the rocks, and estimates of the time gap between the fossil minimum and the true clade ages []. This uncertainty is accommodated by statistical distributions describing the prior paleontological estimates of the true node ages within a phylogeny. Since a node cannot be older than its ancestors, the specified calibration densities are truncated to accommodate this intrinsic constraint from tree topology, generating the effective prior used by the dating program. The resulting marginal prior densities on clade ages can be quite different from the fossil calibration densities specified originally [].

Discussion

59 Darwin C. On the Origin of Species. 60 Davidson E.H.

Erwin D.H. Gene regulatory networks and the evolution of animal body plans. 61 Erwin D.H. Evolutionary uniformitarianism. 38 Lee M.S.

Soubrier J.

Edgecombe G.D. Rates of phenotypic and genomic evolution during the Cambrian explosion. The timing of the emergence of animals has troubled evolutionary biologists at least since Darwin, who was sufficiently incredulous that he considered the abrupt appearance of animal fossils in the Cambrian as a challenge to his theory of evolution by natural selection []. There has been, as a result, a long history of attempts to rationalize a rapid radiation of animals through theories of non-uniform evolutionary processes, such as homeotic mutations, removal of environmental restrictions on larger body sizes, through to the assembly of gene regulation kernels—proposed both as an explanation for rapid rates of innovation followed by subsequent constraint against fundamental innovation of new body plans after the Cambrian []. Indeed, there have been explicit attempts to accommodate rapid rates of phenotypic evolution in the early Cambrian, compatible with these hypotheses and a semi-literal (albeit phylogenetically constrained) reading of the fossil record [].

20 Love G.D.

Grosjean E.

Stalvies C.

Fike D.A.

Grotzinger J.P.

Bradley A.S.

Kelly A.E.

Bhatia M.

Meredith W.

Snape C.E.

et al. Fossil steroids record the appearance of Demospongiae during the Cryogenian period. 62 Maloof A.C.

Rose C.V.

Beach R.

Samuels B.M.

Calmet C.C.

Erwin D.H.

Poirer G.R.

Yao N.

Simons F.J. Possible animal-body fossils in pre-Marinoan limestones from South Australia. 63 Budd G.E.

Jensen S. A critical reappraisal of the fossil record of the bilaterian phyla. 64 Runnegar B. The Cambrian Explosion - animals or fossils. Figure 6 The Timetree of the Metazoa Encompassing Major Sources of Uncertainty in Time Estimates Show full caption Node ages are plotted at the posterior mean for the calibration strategy 1, one partition, IR, and LG + Γ analysis. The node bars are composites extending from the minimum 2.5% HPD limit to the maximum 97.5% limit across all analyses (excluding results from calibration strategies 3 and 4 and from alternative topologies). Cen, Cenozoic; K, Cretaceous; Jr, Jurassic; Tr, Triassic; Pr, Permian; Carb, Carboniferous; Dev, Devonian; S, Silurian; O, Ordovician; Cam, Cambrian; Ediacar, Ediacaran. And yet our results, as have others before them, suggest that there is no justification for invoking non-uniform mechanisms to explain the emergence of animals and their phylum-level body plans. Our analysis attempts to integrate different interpretations of the animal fossil record in informing the minimum age of animal clades. Some of these identify fossil evidence of animals extending into the Cryogenian [], whereas, at the other extreme, others argue that coherent evidence of animals is limited to the Cambrian or the terminal few millions of years of the Neoproterozoic []. Although a case may be made for the restriction of animal fossils to the Phanerozoic, there is only negative evidence (an absence of uncontroversial animal fossils) supporting a Cambrian explosion of animals. This is the long-standing conundrum of the Cambrian—whether the first animal fossils faithfully reflect an explosion in animal biodiversity or merely an explosion of fossils []. The results of our study—which integrates fossil and molecular evidence to establish an evolutionary timescale—suggest that the Cambrian explosion is a phenomenon of fossilization, while biological diversity was established in the Neoproterozoic. Integrating all of the sources of uncertainty that we explore ( Figure 6 Table 1 ) allows us to conclude that crown Metazoa originated 833–650 Ma, fully within the Cryogenian, while the component clades of crown Eumetazoa (746–626 Ma), crown Bilateria (688–596) Ma, crown Deuterostomia (662–587 Ma), and crown Protostomia (653–578 Ma) all diverged within a Cryogenian to early- or mid-Ediacaran interval.

65 Höhna S.

Stadler T.

Ronquist F.

Britton T. Inferring speciation and extinction rates under different sampling schemes. 66 Lartillot N.

Philippe H. A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process. 67 Rota-Stabelli O.

Daley A.C.

Pisani D. Molecular timetrees reveal a Cambrian colonization of land and a new scenario for ecdysozoan evolution. The results of our analyses leads us to reject the hypothesis that metazoans, eumetazoans, bilaterians, protostomes, deuterostomes, ecdysozoans, lophotrochozoans, or, for that matter, any of the component phylum-level total groups, originated in the Cambrian. The uncertainties from competing interpretations of the fossil record, through the choice of rate models and sequence partition strategies, to competing phylogenetic hypotheses all contribute to an evolutionary timescale that lacks sufficient precision to rule out many hypotheses. The situation is compounded by at least two additional sources of uncertainty that we did not study here: uncertainty introduced by the birth-death time prior, and its failure to accommodate diversified sampling of species in phylogenies [], and uncertainty due to the substitution model, which may have an important effect when estimating branch lengths in ancient phylogenies []. Some of the uncertainty in metazoan divergence times can be reduced, for example, by the addition of more sequence data, constraining local rate variation rate through the addition of more taxa. However, the improvements in precision possible even with genome-scale sequence data will be limited by the confounding effects of time and rate, which is the crux of the problem.

68 Brasier M.D.

Antcliffe J.B.

Callow R. Evolutionary trends in remarkable fossil preservation across the Ediacaran–Cambrian transition and the impact of Metazoan mixing. 14 Peters S.E.

Gaines R.R. Formation of the ‘Great Unconformity’ as a trigger for the Cambrian explosion. 9 Dalziel I.W.D. Cambrian transgression and radiation linked to an Iapetus-Pacific oceanic connection?. No matter how imprecise, our timescale for metazoan diversification still indicates a mismatch between the fossil evidence used to calibrate the molecular clock analyses and the resulting divergence time estimates. This is not altogether surprising since, by definition, minimum constraints of clade ages anticipate their antiquity. Nevertheless, it is the extent of this prehistory that is surprising, particularly since the conditions required for exceptional fossil preservation, so key to evidencing the existence of animal phyla in the early Cambrian, obtained also in the Ediacaran []. However, the early Cambrian is characterized by a global sea level rise associated with increased tectonic activity leading to the destruction of older rock sequences by erosion and subduction. Although this may have promoted the innovation and radiation of skeletonizing animals [], it will also have diminished the fossil record of their forebears []. That said, there remains a record of metazoan- and bilaterian-like fossil remains and traces in the Ediacaran that we considered insufficiently robust to substantiate a minimum constraint on metazoan clades but that invariably informed maxima. Further insights into the biology of these organisms and others like them may well explain away apparent inconsistencies between molecular clock estimates of deep metazoan clade ages and their fossil record.