Molecular ecology studies of species radiations: current research gaps,opportunities and challenges

Understanding the drivers and limits of species radiations is a crucial goal of evolutionary genetics and molecular ecology, yet research on this topic has been hampered by the notorious difficulty of connecting micro‐ and macroevolutionary approaches to studying the drivers of diversification. To chart the current research gaps, opportunities and challenges of molecular ecology approaches to studying radiations, we examine the literature in the journal Molecular Ecology and revisit recent high‐profile examples of evolutionary genomic research on radiations. We find that available studies of radiations are highly unevenly distributed among taxa, with many ecologically important and species‐rich organismal groups remaining severely understudied, including arthropods, plants and fungi. Most studies employed molecular methods suitable over either short or long evolutionary time scales, such as microsatellites or restriction site‐associated DNA sequencing (RAD‐seq) in the former case and conventional amplicon sequencing of organellar DNA in the latter. The potential of molecular ecology studies to address and resolve patterns and processes around the species level in radiating groups of taxa is currently limited primarily by sample size and a dearth of information on radiating nuclear genomes as opposed to organellar ones. Based on our literature survey and personal experience, we suggest possible ways forward in the coming years. We touch on the potential and current limitations of whole‐genome sequencing (WGS) in studies of radiations. We suggest that WGS and targeted (‘capture’) resequencing emerge as the methods of choice for scaling up the sampling of populations, species and genomes, including currently understudied organismal groups and the genes or regulatory elements expected to matter most to species radiations.     

Spatiotemporal sampling patterns in the 230 million year fossil record of terrestrial crocodylomorphs and their impact on diversity

The 24 extant crocodylian species are the remnants of a once much more diverse and widespread clade. Crocodylomorpha has an approximately 230 million year evolutionary history, punctuated by a series of radiations and extinctions. However, the group's fossil record is biased. Previous studies have reconstructed temporal patterns in subsampled crocodylomorph palaeobiodiversity, but have not explicitly examined variation in spatial sampling, nor the quality of this record. We compiled a dataset of all taxonomically diagnosable non‐marine crocodylomorph species (393). Based on the number of phylogenetic characters that can be scored for all published fossils of each species, we calculated a completeness value for each taxon. Mean average species completeness (56%) is largely consistent within subgroups and for different body size classes, suggesting no significant biases across the crocodylomorph tree. In general, average completeness values are highest in the Mesozoic, with an overall trend of decreasing completeness through time. Many extant taxa are identified in the fossil record from very incomplete remains, but this might be because their provenance closely matches the species’ present‐day distribution, rather than through autapomorphies. Our understanding of nearly all crocodylomorph macroevolutionary ‘events’ is essentially driven by regional patterns, with no global sampling signal. Palaeotropical sampling is especially poor for most of the group's history. Spatiotemporal sampling bias impedes our understanding of several Mesozoic radiations, whereas molecular divergence times for Crocodylia are generally in close agreement with the fossil record. However, the latter might merely be fortuitous, i.e. divergences happened to occur during our ephemeral spatiotemporal sampling windows.     

Why coelacanths are not ‘living fossils’

A series of recent studies on extant coelacanths has emphasised the slow rate of molecular and morphological evolution in these species. These studies were based on the assumption that a coelacanth is a ‘living fossil’ that has shown little morphological change since the Devonian, and they proposed a causal link between low molecular evolutionary rate and morphological stasis. Here, we have examined the available molecular and morphological data and show that: (i) low intra‐specific molecular diversity does not imply low mutation rate, (ii) studies not showing low substitution rates in coelacanth are often neglected, (iii) the morphological stability of coelacanths is not supported by paleontological evidence. We recall that intra‐species levels of molecular diversity, inter‐species genome divergence rates and morphological divergence rates are under different constraints and they are not necessarily correlated. Finally, we emphasise that concepts such as ‘living fossil’, ‘basal lineage’, or ‘primitive extant species’ do not make sense from a tree‐thinking perspective. Editor's suggested further reading in BioEssays Tree thinking for all biology: the problem with reading phylogenies as ladders of progress Abstract     

The continuity of microevolution and macroevolution

A persistent debate in evolutionary biology is one over the continuity of microevolution and macroevolution – whether macroevolutionary trends are governed by the principles of microevolution. The opposition of evolutionary trends over different time scales is taken as evidence that selection is uncoupled over these scales. I argue that the paradox inferred by trend opposition is eliminated by a hierarchical application of the ‘geometric‐mean fitness’ principle, a principle that has been invoked only within the limited context of microevolution in response to environmental variance. This principle implies the elimination of well adapted genotypes – even those with the highest arithmetic mean fitness over a shorter time scale. Contingent on premises concerning the temporal structure of environmental variance, selectivity of extinction, and clade‐level heritability, the evolutionary outcome of major environmental change may be viewed as identical in principle to the outcome of minor environmental fluctuations over the short‐term. Trend reversals are thus recognized as a fundamental property of selection operating at any phylogenetic level that occur in response to event severities of any magnitude over all time scales. This ‘bet‐hedging’ perspective differs from others in that a specified, single hierarchical selective process is proposed to explain observed hierarchical patterns of extinction.     

Extant‐only comparative methods fail to recover the disparity preserved in the bird fossil record

Most extant species are in clades with poor fossil records, and recent studies of comparative methods show they have low power to infer even highly simplified models of trait evolution without fossil data. Birds are a well‐studied radiation, yet their early evolutionary patterns are still contentious. The fossil record suggests that birds underwent a rapid ecological radiation after the end‐Cretaceous mass extinction, and several smaller, subsequent radiations. This hypothesized series of repeated radiations from fossil data is difficult to test using extant data alone. By uniting morphological and phylogenetic data on 604 extant genera of birds with morphological data on 58 species of extinct birds from 50 million years ago, the “halfway point” of avian evolution, I have been able to test how well extant‐only methods predict the diversity of fossil forms. All extant‐only methods underestimate the disparity, although the ratio of within‐ to between‐clade disparity does suggest high early rates. The failure of standard models to predict high early disparity suggests that recent radiations are obscuring deep time patterns in the evolution of birds. Metrics from different models can be used in conjunction to provide more valuable insights than simply finding the model with the highest relative fit.     

Understanding the effect of competition during evolutionary radiations: an integrated model of phenotypic and species diversification

Competition can drive macroevolutionary change, for example during adaptive radiations. However, we still lack a clear understanding of how it shapes diversification processes and patterns. To better understand the macroevolutionary consequences of competition, as well as the signal left on phylogenetic data, we developed a model linking trait evolution and species diversification in an ecological context. We find four main results: first, competition spurs trait diversity but not necessarily species richness; second, competition produces slowdowns in species diversification even in the absence of explicit ecological limits, but not in phenotypic diversification even in the presence of such limits; third, early burst patterns do not provide a reliable way of testing for adaptive radiations; and fourth, looking for phylogenetic signal in trait data and support for phenotypic models incorporating competition is a better alternative. Our results clarify the macroevolutionary consequences of competition and could help design more powerful tests of adaptive radiations in nature.     

Body Size Evolution in Extant Oryzomyini Rodents: Cope's Rule or Miniaturization?

At the macroevolutionary level, one of the first and most important hypotheses that proposes an evolutionary tendency in the evolution of body sizes is "Cope's rule". This rule has considerable empirical support in the fossil record and predicts that the size of species within a lineage increases over evolutionary time. Nevertheless, there is also a large amount of evidence indicating the opposite pattern of miniaturization over evolutionary time. A recent analysis using a single phylogenetic tree approach and a bayesian based model of evolution found no evidence for Cope's rule in extant mammal species. Here we utilize a likelihood-based phylogenetic method, to test the evolutionary trend in body size, which considers phylogenetic uncertainty, to discern between Cope's rule and miniaturization, using extant Oryzomyini rodents as a study model. We evaluated body size trends using two principal predictions: (a) phylogenetically related species are more similar in their body size, than expected by chance; (b) body size increased (Cope's rule)/decreased (miniaturization) over time. Consequently the distribution of forces and/or constraints that affect the tendency are homogenous and generate this directional process from a small/large sized ancestor. Results showed that body size in the Oryzomyini tribe evolved according to phylogenetic relationships, with a positive trend, from a small sized ancestor. Our results support that the high diversity and specialization currently observed in the Oryzomyini tribe is a consequence of the evolutionary trend of increased body size, following and supporting Cope's rule.     

Tempo and mode of mandibular shape and size evolution reveal mixed support for incumbency effects in two clades of island‐endemic rodents (Muridae: Murinae)*

Existing radiations in a spatially limited system such as an oceanic island may limit the ecological opportunity experienced by later colonists, resulting in lower macroevolutionary rates for secondary radiations. Additionally, potential colonists may be competitively excluded by these incumbent (resident) species, unless they are biologically distinct (biotic filtering). The extant phenotypic diversity of secondary colonists may thus be impacted by lower rates of phenotypic evolution, exclusion from certain phenotypes, and transitions to new morphotypes to escape competition from incumbent lineages. We used geometric morphometric methods to test whether the rates and patterns of mandibular evolution of the Luzon “old endemic” rodent clades, Phloeomyini and Chrotomyini, are consistent with these predictions. Each clade occupied nearly completely separate shape space and partially separate size space. We detected limited support for decelerating and clade‐specific evolutionary rates for both shape and size, with strong evidence for a shift in evolutionary mode within Chrotomyini. Our results suggest that decelerating phenotypic evolutionary rates are not a necessary result of incumbency interactions; rather, incumbency effects may be more likely to determine which clades can become established in the system. Nonincumbent clades that pass a biotic filter can potentially exhibit relatively unfettered evolution.     

Cultural and climatic changes shape the evolutionary history of the Uralic languages

Quantitative phylogenetic methods have been used to study the evolutionary relationships and divergence times of biological species, and recently, these have also been applied to linguistic data to elucidate the evolutionary history of language families. In biology, the factors driving macroevolutionary processes are assumed to be either mainly biotic (the Red Queen model) or mainly abiotic (the Court Jester model) or a combination of both. The applicability of these models is assumed to depend on the temporal and spatial scale observed as biotic factors act on species divergence faster and in smaller spatial scale than the abiotic factors. Here, we used the Uralic language family to investigate whether both ‘biotic’ interactions (i.e. cultural interactions) and abiotic changes (i.e. climatic fluctuations) are also connected to language diversification. We estimated the times of divergence using Bayesian phylogenetics with a relaxed‐clock method and related our results to climatic, historical and archaeological information. Our timing results paralleled the previous linguistic studies but suggested a later divergence of Finno‐Ugric, Finnic and Saami languages. Some of the divergences co‐occurred with climatic fluctuation and some with cultural interaction and migrations of populations. Thus, we suggest that both ‘biotic’ and abiotic factors contribute either directly or indirectly to the diversification of languages and that both models can be applied when studying language evolution.     

Discrete and morphometric traits reveal contrasting patterns and processes in the macroevolutionary history of a clade of scorpions

Many palaeontological studies have investigated the evolution of entire body plans, generally relying on discrete character‐taxon matrices. In contrast, macroevolutionary studies performed by neontologists have mostly focused on morphometric traits. Although these data types are very different, some studies have suggested that they capture common patterns. Nonetheless, the tests employed to support this claim have not explicitly incorporated a phylogenetic framework and may therefore be susceptible to confounding effects due to the presence of common phylogenetic structure. We address this question using the scorpion genus Brachistosternus Pocock 1893 as case study. We make use of a time‐calibrated multilocus molecular phylogeny, and compile discrete and traditional morphometric data sets, both capturing the overall morphology of the organisms. We find that morphospaces derived from these matrices are significantly different, and that the degree of discordance cannot be replicated by simulations of random character evolution. Moreover, we find strong support for contrasting modes of evolution, with discrete characters being congruent with an ‘early burst’ scenario whereas morphometric traits suggest species‐specific adaptations to have driven morphological evolution. The inferred macroevolutionary dynamics are therefore contingent on the choice of character type. Finally, we confirm that metrics of correlation fail to detect these profound differences given common phylogenetic structure in both data sets, and that methods incorporating a phylogenetic framework and accounting for expected covariance should be favoured.     

POTENTIAL PITFALLS OF RECONSTRUCTING DEEP TIME EVOLUTIONARY HISTORY WITH ONLY EXTANT DATA,A CASE STUDY USING THE CANIDAE (MAMMALIA,CARNIVORA)

Reconstructing evolutionary patterns and their underlying processes is a central goal in biology. Yet many analyses of deep evolutionary histories assume that data from the fossil record is too incomplete to include, and rely solely on databases of extant taxa. Excluding fossil taxa assumes that character state distributions across living taxa are faithful representations of a clade's entire evolutionary history. Many factors can make this assumption problematic. Fossil taxa do not simply lead‐up to extant taxa; they represent now‐extinct lineages that can substantially impact interpretations of character evolution for extant groups. Here, we analyze body mass data for extant and fossil canids (dogs, foxes, and relatives) for changes in mean and variance through time. AIC‐based model selection recovered distinct models for each of eight canid subgroups. We compared model fit of parameter estimates for (1) extant data alone and (2) extant and fossil data, demonstrating that the latter performs significantly better. Moreover, extant‐only analyses result in unrealistically low estimates of ancestral mass. Although fossil data are not always available, reconstructions of deep‐time organismal evolution in the absence of deep‐time data can be highly inaccurate, and we argue that every effort should be made to include fossil data in macroevolutionary studies.     

Evolutionary analyses of non‐family genes in plants

There are a large number of ‘non‐family’ (NF) genes that do not cluster into families with three or more members per genome. While gene families have been extensively studied, a systematic analysis of NF genes has not been reported. We performed comparative studies on NF genes in 14 plant species. Based on the clustering of protein sequences, we identified ~94 000 NF genes across these species that were divided into five evolutionary groups: Viridiplantae wide, angiosperm specific, monocot specific, dicot specific, and those that were species specific. Our analysis revealed that the NF genes resulted largely from less frequent gene duplications and/or a higher rate of gene loss after segmental duplication relative to genes in both low‐copy‐number families (LF; 3–10 copies per genome) and high‐copy‐number families (HF; >10 copies). Furthermore, we identified functions enriched in the NF gene set as compared with the HF genes. We found that NF genes were involved in essential biological processes shared by all plant lineages (e.g. photosynthesis and translation), as well as gene regulation and stress responses associated with phylogenetic diversification. In particular, our analysis of an Arabidopsis protein–protein interaction network revealed that hub proteins with the top 10% most connections were over‐represented in the NF set relative to the HF set. This research highlights the roles that NF genes may play in evolutionary and functional genomics research.     

Catching the phylogenic history through the ontogenic hourglass: a phylogenomic analysis of Drosophila body segmentation genes

SUMMARY The phylogenetic information content of different developmental stages is a long‐standing issue in the study of development and evolution. We performed phylogenetic analyses of 51 body segmentation genes in 12 species of Drosophila in order to investigate the impact of the mode of evolution of development on phylogeny inference. Previous studies of these genes in Drosophila using pairwise phenetic comparisons at the species group level revealed the presence of an “hourglass model” (HG), wherein mid‐embryonic stages are the most evolutionarily constrained. We utilized two character‐based approaches: taxonomic congruence using the relative consensus fork index (RCFI), in which phylogenies are inferred from each gene separately and compared with a total evidence tree (TET), and partitioned simultaneous analysis using several indices such as branch support (BS) and localized incongruence length difference (LILD) test. We also proposed a new index, the recapitulatory index (R), which divides the number of synapomorphies on the total number of informative characters in a data set. Polynomial adjustment of both BS and R indices showed strong support for the hourglass model regardless of the taxonomic level (species subgroup vs. subgenera), showing less phylogenetic information content for mid‐developmental stages (mainly the zygotic segment polarity stage). Significant LILD scores were randomly distributed among developmental stages revealing the absence of differential selective constraints, but were significantly related to chromosomal location showing physical (linkage) impact on phylogenetic incongruence. RCFI was the most sensitive measure to taxonomic level, having a convex parabola at the species subgroup level in support of the hourglass model and a concave parabola at the subgeneric level in support of the adaptive penetrance model. This time‐dependent discrepancy of best fit developmental model parallels previous conflicting results from the vertebrates. Because of the quasi‐phenetic nature of this index, we argue that the discrepancy is due to the evolutionary rate heterogeneity of developmental genes rather than to fundamental differences among organisms. We suggest that simultaneous character‐based analyses give better macroevolutionary support to the hourglass model of the developmental constraints on genome evolution than pairwise phenetic comparisons.     

Convergence and divergence in the evolution of cat skulls: temporal and spatial patterns of morphological diversity

Background

Studies of biological shape evolution are greatly enhanced when framed in a phylogenetic perspective. Inclusion of fossils amplifies the scope of macroevolutionary research, offers a deep-time perspective on tempo and mode of radiations, and elucidates life-trait changes. We explore the evolution of skull shape in felids (cats) through morphometric analyses of linear variables, phylogenetic comparative methods, and a new cladistic study of saber-toothed cats.

Methodology/Principal Findings

A new phylogenetic analysis supports the monophyly of saber-toothed cats (Machairodontinae) exclusive of Felinae and some basal felids, but does not support the monophyly of various saber-toothed tribes and genera. We quantified skull shape variation in 34 extant and 18 extinct species using size-adjusted linear variables. These distinguish taxonomic group membership with high accuracy. Patterns of morphospace occupation are consistent with previous analyses, for example, in showing a size gradient along the primary axis of shape variation and a separation between large and small-medium cats. By combining the new phylogeny with a molecular tree of extant Felinae, we built a chronophylomorphospace (a phylogeny superimposed onto a two-dimensional morphospace through time). The evolutionary history of cats was characterized by two major episodes of morphological divergence, one marking the separation between saber-toothed and modern cats, the other marking the split between large and small-medium cats.

Conclusions/Significance

Ancestors of large cats in the ‘Panthera’ lineage tend to occupy, at a much later stage, morphospace regions previously occupied by saber-toothed cats. The latter radiated out into new morphospace regions peripheral to those of extant large cats. The separation between large and small-medium cats was marked by considerable morphologically divergent trajectories early in feline evolution. A chronophylomorphospace has wider applications in reconstructing temporal transitions across two-dimensional trait spaces, can be used in ecophenotypical and functional diversity studies, and may reveal novel patterns of morphospace occupation.     

Reconstructing the evolution of umbonal sculptures in the Unionida

‘Umbonal sculptures’ of freshwater mussels (Unionida), which ornament the early ontogenetic shell, have long been used for species identification. Specificity of these sculptures to higher taxonomic levels and their value for phylogenetic reconstruction are still under considerable scientific debate. In particular, the distribution of beak sculpture morphotypes across the unionoid phylogeny and, consequently, evolution of this character remain poorly understood. Based on an examination of 187 taxa, covering five of the six extant unionoid families, this study presents a new model of character evolution of umbonal sculptures in the order. Ten morphotypes were recognized and conceptualized into the cladistic characters sculpture presence and category. Optimization of sculpture presence on two recent hypotheses of palaeoheterodont phylogenetic relationships using the program Mesquite indicates a sculptured common ancestor of the extant Unionida, with multiple losses of the umbonal ornament occurring subsequently within the clade. Reconstruction of changes in sculpture category is ambiguous and demonstrates the need for further research into the evolutionary relationships of freshwater mussels in general and of their early ontogenetic sculptures in particular. Ambiguity is reduced in analyses applying a model with unequal costs of transformation between character states, which was derived from observations on intermediate forms and polymorphisms. These analyses suggest ‘V‐shaped’ or ‘nodulous’ sculpture as the plesiomorphic sculptural category for Unionida. The relatively low levels of homoplasy inferred for V‐shaped, pseudo‐radial and double‐looped sculptures suggest that these types may comprise useful guides to relationships within Unionida. The high degree of homoplasy of W‐shaped, pseudo‐concentric, wrinkled and single‐looped sculptures, on the other hand, renders these sculpture types less fit for such purposes.     

INTEGRATING FOSSILS WITH MOLECULAR PHYLOGENIES IMPROVES INFERENCE OF TRAIT EVOLUTION

Comparative biologists often attempt to draw inferences about tempo and mode in evolution by comparing the fit of evolutionary models to phylogenetic comparative data consisting of a molecular phylogeny with branch lengths and trait measurements from extant taxa. These kinds of approaches ignore historical evidence for evolutionary pattern and process contained in the fossil record. In this article, we show through simulation that incorporation of fossil information dramatically improves our ability to distinguish among models of quantitative trait evolution using comparative data. We further suggest a novel Bayesian approach that allows fossil information to be integrated even when explicit phylogenetic hypotheses are lacking for extinct representatives of extant clades. By applying this approach to a comparative dataset comprising body sizes for caniform carnivorans, we show that incorporation of fossil information not only improves ancestral state estimates relative to those derived from extant taxa alone, but also results in preference of a model of evolution with trend toward large body size over alternative models such as Brownian motion or Ornstein–Uhlenbeck processes. Our approach highlights the importance of considering fossil information when making macroevolutionary inference, and provides a way to integrate the kind of sparse fossil information that is available to most evolutionary biologists.     

Mitochondrial rate variation among lineages of passerine birds

The order Passeriformes comprises the majority of extant avian species. Analyses of molecular data have provided important insights into the evolution of this diverse order. However, molecular estimates of the evolutionary and demographic timescales of passerine species have been hindered by a lack of reliable calibrations. This has led to a reliance on the application of standard substitution rates to mitochondrial DNA data, particularly rates estimated from analyses of the gene encoding cytochrome b (CYTB). To investigate patterns of rate variation across passerine lineages, we used a Bayesian phylogenetic approach to analyse the protein‐coding genes of 183 mitochondrial genomes. We found that the most commonly used mitochondrial marker, CYTB, has low variation in rates across passerine lineages. This lends support to its widespread use as a molecular clock in birds. However, we also found that the patterns of among‐lineage rate variation in CYTB are only weakly related to the evolutionary rate of the mitochondrial genome as a whole. Our analyses confirmed the presence of mutational saturation at third codon positions across the protein‐coding genes of the mitochondrial genome, reinforcing the view that these sites should be excluded in studies of deep passerine relationships. The results of our analyses have provided information that will be useful for molecular‐clock studies of passerine evolution.     

‘Fish’ (Actinopterygii and Elasmobranchii) diversification patterns through deep time

Actinopterygii (ray‐finned fishes) and Elasmobranchii (sharks, skates and rays) represent more than half of today's vertebrate taxic diversity (approximately 33000 species) and form the largest component of vertebrate diversity in extant aquatic ecosystems. Yet, patterns of ‘fish’ evolutionary history remain insufficiently understood and previous studies generally treated each group independently mainly because of their contrasting fossil record composition and corresponding sampling strategies. Because direct reading of palaeodiversity curves is affected by several biases affecting the fossil record, analytical approaches are needed to correct for these biases. In this review, we propose a comprehensive analysis based on comparison of large data sets related to competing phylogenies (including all Recent and fossil taxa) and the fossil record for both groups during the Mesozoic–Cainozoic interval. This approach provides information on the ‘fish’ fossil record quality and on the corrected ‘fish’ deep‐time phylogenetic palaeodiversity signals, with special emphasis on diversification events. Because taxonomic information is preserved after analytical treatment, identified palaeodiversity events are considered both quantitatively and qualitatively and put within corresponding palaeoenvironmental and biological settings. Results indicate a better fossil record quality for elasmobranchs due to their microfossil‐like fossil distribution and their very low diversity in freshwater systems, whereas freshwater actinopterygians are diverse in this realm with lower preservation potential. Several important diversification events are identified at familial and generic levels for elasmobranchs, and marine and freshwater actinopterygians, namely in the Early–Middle Jurassic (elasmobranchs), Late Jurassic (actinopterygians), Early Cretaceous (elasmobranchs, freshwater actinopterygians), Cenomanian (all groups) and the Paleocene–Eocene interval (all groups), the latter two representing the two most exceptional radiations among vertebrates. For each of these events along with the Cretaceous‐Paleogene extinction, we provide an in‐depth review of the taxa involved and factors that may have influenced the diversity patterns observed. Among these, palaeotemperatures, sea‐levels, ocean circulation and productivity as well as continent fragmentation and environment heterogeneity (reef environments) are parameters that largely impacted on ‘fish’ evolutionary history, along with other biotic constraints.     

A phylogeny of Cenozoic macroperforate planktonic foraminifera from fossil data

We present a complete phylogeny of macroperforate planktonic foraminifer species of the Cenozoic Era (∼65 million years ago to present). The phylogeny is developed from a large body of palaeontological work that details the evolutionary relationships and stratigraphic (time) distributions of species‐level taxa identified from morphology (‘morphospecies’). Morphospecies are assigned to morphogroups and ecogroups depending on test morphology and inferred habitat, respectively. Because gradual evolution is well documented in this clade, we have identified many instances of morphospecies intergrading over time, allowing us to eliminate ‘pseudospeciation’ and ‘pseudoextinction’ from the record and thereby permit the construction of a more natural phylogeny based on inferred biological lineages. Each cladogenetic event is determined as either budding or bifurcating depending on the pattern of morphological change at the time of branching. This lineage phylogeny provides palaeontologically calibrated ages for each divergence that are entirely independent of molecular data. The tree provides a model system for macroevolutionary studies in the fossil record addressing questions of speciation, extinction, and rates and patterns of evolution.