Measuring Stratigraphic Congruence Across Trees, Higher Taxa, and Time

Relative performance of Bayesian morphological clock and parsimony methods for phylogenetic reconstructions: Insights from the case of Myomiminae and Dryomyinae glirid rodents

Extinct organisms provide crucial information about the origin and time of origination of extant groups. The importance of morphological phylogenetics for rigorously dating the tree of life is now widely recognized and has been revitalized by methodological developments such as the application of tip-dating Bayesian approaches. Traditionally, molecular clocks have been node calibrated. However, node calibrations are often unsatisfactory because they do not allow the fossil age to inform about phylogenetic hypothesis. The introduction of tip calibrations allow fossil species to be included alongside their living relatives, and the absence of molecular sequence data for these taxa to be remedied by supplementing the sequence alignments for living taxa with phenotype character matrices for both living and fossil taxa. Therefore, only phylogenetic analyses that take into account morphological characters can incorporate both fossil and extant species. Herein we present an unprecedented morphological dataset for a vast group of glirid rodents, to which different phylogenetic methodologies have been applied. We have compared the tree topologies resulting from traditional parsimony and Bayesian phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods seem to outperform parsimony with our dataset. The strict consensus tree recovered by tip dating invalidates the classic classification and allows dates to be proposed for the divergence and origin of the different clades.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

Molecular phylogenies map to biogeography better than morphological ones

Phylogenetic relationships are inferred principally from two classes of data: morphological and molecular. Currently, most phylogenies of extant taxa are inferred from molecules and when morphological and molecular trees conflict the latter are often preferred. Although supported by simulations, the superiority of molecular trees has rarely been assessed empirically. Here we test phylogenetic accuracy using two independent data sources: biogeographic distributions and fossil first occurrences. For 48 pairs of morphological and molecular trees we show that, on average, molecular trees provide a better fit to biogeographic data than their morphological counterparts and that biogeographic congruence increases over research time. We find no significant differences in stratigraphic congruence between morphological and molecular trees. These results have implications for understanding the distribution of homoplasy in morphological data sets, the utility of morphology as a test of molecular hypotheses and the implications of analysing fossil groups for which molecular data are unavailable.

Using biogeographical and phylogenetic data, it is shown that molecular trees fit species geographical data better than trees inferred from morphology, and that these differences are not simply due to better tree resolution.

Bayesian Morphological Clock versus Parsimony: An insight into the relationships and dispersal events of postvacuum Cricetidae (Rodentia, Mammalia)

Establishing an evolutionary timeline is fundamental for tackling a great variety of topics in evolutionary biology, including the reconstruction of patterns of historical biogeography, coevolution and diversification. However, the tree of life is pruned by extinction and molecular data cannot be gathered for extinct lineages. Until recently methodological challenges have prevented the application of tip-dating Bayesian approaches in morphology-based fossil-only datasets. Herein, we present a morphological dataset for a group of cricetid rodents to which we apply an array of methods fairly new in palaeontology that can be used by palaeontologists for the analysis of entirely extinct clades. We compare the tree topologies obtained by traditional parsimony, time-calibrated and non-calibrated Bayesian inference phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods outperform parsimony in the case of our dataset, which includes highly homoplastic morphological characters. Regardless, all three topologies support the monophyly of Megacricetodontinae, Democricetodontinae and Cricetodontinae. Dispersal and speciation events inferred through Bayesian Binary Markov chain Monte Carlo and biodiversity analyses provide evidence for a correlation between biogeographic events, climatic changes and diversification in cricetids.

Evolution of gigantism in right and bowhead whales (Cetacea: Mysticeti: Balaenidae)

The evolution of gigantic body size represents a key to understand the ecological role of baleen whales in oceanic ecosystems. Many efforts have been devoted to the formulation of equations relating different body parts to total body length and mass in living and fossil mysticetes, mainly focusing on balaenopterid and balaenopterid-like mysticetes. Right whales (family Balaenidae) have a unique head-to-body length ratio, suggesting that their body proportions cannot be predicted effectively using equations based primarily on non-balaenid mysticetes. A new morphometric dataset of living and fossil balaenids is provided herein, and new regression equations allow one to predict the body length and mass of extinct species based on the expected head-to-body length ratio of extant balaenids. The reconstructed values are mapped on a new phylogenetic analysis of the Balaenidae, inferring body size and mass at ancestral nodes. The variations of body size and mass in Balaenidae since the early Miocene are reconstructed, revealing that: (1) a reduction in total body length occurred in the early Pliocene; (2) the origin of the gigantic body size in the bowhead whale (Balaena mysticetus) is probably related to invasion of the Arctic Ocean in the last 3 Myr; and (3) the origin of the gigantic body size in the right whales (genus Eubalaena) occurred since the latest Miocene, probably concomitant with pulses of nutrients sustaining large zooplankton populations. We suggest that the evolution of gigantism in Balaenidae occurred independently in two lineages and, probably, in response to different palaeoenvironmental drivers.

Bayesian Tip-Dated Phylogenetics in Paleontology: Topological Effects and Stratigraphic Fit

The incorporation of stratigraphic data into phylogenetic analysis has a long history of debate but is not currently standard practice for paleontologists. Bayesian tip-dated (or morphological clock) phylogenetic methods have returned these arguments to the spotlight, but how tip dating affects the recovery of evolutionary relationships has yet to be fully explored. Here I show, through analysis of several data sets with multiple phylogenetic methods, that topologies produced by tip dating are outliers as compared to topologies produced by parsimony and undated Bayesian methods, which retrieve broadly similar trees. Unsurprisingly, trees recovered by tip dating have better fit to stratigraphy than trees recovered by other methods under both the Gap Excess Ratio (GER) and the Stratigraphic Completeness Index (SCI). This is because trees with better stratigraphic fit are assigned a higher likelihood by the fossilized birth-death tree model. However, the degree to which the tree model favors tree topologies with high stratigraphic fit metrics is modulated by the diversification dynamics of the group under investigation. In particular, when net diversification rate is low, the tree model favors trees with a higher GER compared to when net diversification rate is high. Differences in stratigraphic fit and tree topology between tip dating and other methods are concentrated in parts of the tree with weaker character signal, as shown by successive deletion of the most incomplete taxa from two data sets. These results show that tip dating incorporates stratigraphic data in an intuitive way, with good stratigraphic fit an expectation that can be overturned by strong evidence from character data. [fossilized birth-death; fossils; missing data; morphological clock; morphology; parsimony; phylogenetics.].

Craniodental and Postcranial Characters of Non-Avian Dinosauria Often Imply Different Trees

Despite the increasing importance of molecular sequence data, morphology still makes an important contribution to resolving the phylogeny of many groups, and is the only source of data for most fossils. Most systematists sample morphological characters as broadly as possible on the principle of total evidence. However, it is not uncommon for sampling to be focused on particular aspects of anatomy, either because characters therein are believed to be more informative, or because preservation biases restrict what is available. Empirically, the optimal trees from partitions of morphological data sets often represent significantly different hypotheses of relationships. Previous work on hard-part versus soft-part characters across animal phyla revealed significant differences in about a half of sampled studies. Similarly, studies of the craniodental versus postcranial characters of vertebrates revealed significantly different trees in about one-third of cases, with the highest rates observed in non-avian dinosaurs. We test whether this is a generality here with a much larger sample of 81 published data matrices across all major dinosaur groups. Using the incongruence length difference test and two variants of the incongruence relationship difference test, we found significant incongruence in about 50% of cases. Incongruence is not uniformly distributed across major dinosaur clades, being highest (63%) in Theropoda and lowest (25%) in Thyreophora. As in previous studies, our partition tests show some sensitivity to matrix dimensions and the amount and distribution of missing entries. Levels of homoplasy and retained synapomorphy are similar between partitions, such that incongruence must partly reflect differences in patterns of homoplasy between partitions, which may itself be a function of modularity and mosaic evolution. Finally, we implement new tests to determine which partition yields trees most similar to those from the entire matrix. Despite no bias across dinosaurs overall, there are striking differences between major groups. The craniodental characters of Ornithischia and the postcranial characters of Saurischia yield trees most similar to the "total evidence" trees derived from the entire matrix. Trees from these same character partitions also tend to be most stratigraphically congruent: a mutual consilience suggesting that those partitions yield more accurate trees. [Dinosauria; homoplasy; partition homogeneity.].

Phylogenetic incongruence and homoplasy in the appendages and bodies of arthropods: why broad character sampling is best

Notwithstanding the rapidly increasing sampling density of molecular sequence data, morphological characters still make an important contribution to our understanding of the evolutionary relationships of arthropod groups. In many clades, characters relating to the number and morphological specialization of appendages are ascribed particular phylogenetic significance and may be preferentially sampled. However, previous studies have shown that partitions of morphological character matrices often imply significantly different phylogenies. Here, we ask whether a similar incongruence is observed in the appendage and non-appendage characters of arthropods. We apply tree length (incongruence length difference, ILD) and tree distance (incongruence relationship difference, IRD) tests to these partitions in an empirical sample of 53 published neontological datasets for arthropods. We find significant incongruence about one time in five: more often than expected, but markedly less often than in previous partition studies. We also find similar levels of homoplasy in limb and non-limb characters, both in terms of internal consistency and consistency relative to molecular trees. Taken together, these findings imply that sampled limb and non-limb characters are of similar phylogenetic utility and quality, and that a total evidence approach to their analysis is preferable.

Parsimony, not Bayesian analysis, recovers more stratigraphically congruent phylogenetic trees

Reconstructing evolutionary histories requires accurate phylogenetic trees. Recent simulation studies suggest that probabilistic phylogenetic analyses of morphological data are more accurate than traditional parsimony techniques. Here, we use empirical data to compare Bayesian and parsimony phylogenies in terms of their congruence with the distribution of age ranges of the component taxa. Analysis of 167 independent morphological data matrices of fossil tetrapods finds that Bayesian trees exhibit significantly lower stratigraphic congruence than the equivalent parsimony trees. As such, taking stratigraphic data as an independent benchmark indicates that parsimony analyses are more accurate for phylogenetic reconstruction of morphological data. The discrepancy between simulated and empirical studies may result from historic data peaking practices or some complexities of empirical data as yet unaccounted for.