Bayesian placement of fossils on phylogenies using quantitative morphometric data

Accuracy of phylogenetic reconstructions from continuous characters analysed under parsimony and its parametric correlates

Quantitative traits are a source of evolutionary information often difficult to handle in cladistics. Tools exist to analyse this kind of data without subjective discretization, avoiding biases in the delimitation of categorical states. Nonetheless, our ability to accurately infer relationships from continuous characters is incompletely understood, particularly under parsimony analysis. This study evaluates the accuracy of phylogenetic reconstructions from simulated matrices of continuous characters evolving under alternative evolutionary processes and analysed by parsimony. We sampled 100 empirical trees to simulate 9000 matrices, each containing between 25 and 50 taxa, and 50 and 150 continuous characters evolving under three evolutionary processes: Brownian motion, Ornstein-Uhlenbeck and early burst with variable parametrizations. Our cladogram comparisons revealed that continuous character matrices, when properly coded and analysed by parsimony in TNT, carry phylogenetic signals from which species relationships can be inferred, regardless of the evolutionary models and parameterization schemes. Interestingly, implementing equal weighting or implied weighting with varying penalization strengths against homoplasies did not affect cladogram reconstructions based on continuous characters. Finally, the accuracy of continuous characters in resolving species relationships is skewed towards apical nodes of the recovered trees. Our findings provide general insights of the utility of quantitative traits in cladistics and demonstrate that their effectiveness in estimating shallower nodes is independent of the underlying evolutionary model, parameters and weighting schemes.

Transitions Into Freezing Environments Linked With Shifts in Phylogenetic Integration Between Vitaceae Leaf Traits

Understanding how the intrinsic ability of populations and species to meet shifting selective demands shapes evolutionary patterns over both short and long timescales is a major question in biology. One major axis of evolutionary flexibility can be measured by phenotypic integration and modularity. The strength, scale, and structure of integration may constrain or catalyze evolution in the face of new selective pressures. We analyze a dataset of seven leaf measurements across Vitaceae to examine how correlations in trait divergence are linked to transitions between freezing and nonfreezing habitats. We assess this by applying a custom algorithm to compare the timing of habitat shifts to changes in the structure of evolutionary trait correlation at discrete points along a phylogeny. We also explore these patterns in relation to lineage diversification rates to understand how and whether patterns in the evolvability of complex multivariate phenotypes are linked to higher-level macroevolutionary dynamics. We found that shifts in the structure, but not the overall strength, of phylogenetic integration of leaves precipitate colonization of freezing climates. Lineages that underwent associated shifts in leaf trait integration and subsequent movement into freezing habitats also displayed lower turnover and higher net diversification, suggesting a link among shifting vectors of selection, internal constraint, and lineage persistence in the face of changing environments.

Do morphometric data improve phylogenetic reconstruction? A systematic review and assessment

BACKGROUND

Isolating phylogenetic signal from morphological data is crucial for accurately merging fossils into the tree of life and for calibrating molecular dating. However, subjective character definition is a major limitation which can introduce biases that mislead phylogenetic inferences and divergence time estimation. The use of quantitative data, e.g., geometric morphometric (GMM; shape) data can allow for more objective integration of morphological data into phylogenetic inference. This systematic review describes the current state of the field in using continuous morphometric data (e.g., GMM data) for phylogenetic reconstruction and assesses the efficacy of these data compared to discrete characters using the PRISMA-EcoEvo v1.0. reporting guideline, and offers some pathways for approaching this task with GMM data. A comprehensive search string yielded 11,123 phylogenetic studies published in English up to Oct 2023 in the Web of Science database. Title and abstract screening removed 10,975 articles, and full-text screening was performed for 132 articles. Of these, a total of twelve articles met final inclusion criteria and were used for downstream analyses.

RESULTS

Phylogenetic performance was compared between approaches that employed continuous morphometric and discrete morphological data. Overall, the reconstructed phylogenies did not show increased resolution or accuracy (i.e., benchmarked against molecular phylogenies) as continuous data alone or combined with discrete morphological datasets.

CONCLUSIONS

An exhaustive search of the literature for existing empirical continuous data resulted in a total of twelve articles for final inclusion following title/abstract, and full-text screening. Our study was performed under a rigorous framework for systematic reviews, which showed that the lack of available comparisons between discrete and continuous data hinders our understanding of the performance of continuous data. Our study demonstrates the problem surrounding the efficacy of continuous data as remaining relatively intractable despite an exhaustive search, due in part to the difficulty in obtaining relevant comparisons from the literature. Thus, we implore researchers to address this issue with studies that collect discrete and continuous data sets with directly comparable properties (i.e., describing shape, or size).

A macroevolutionary analysis of European Late Upper Palaeolithic stone tool shape using a Bayesian phylodynamic framework

Phylogenetic models are commonly used in palaeobiology to study the patterns and processes of organismal evolution. In the human sciences, phylogenetic methods have been deployed for reconstructing ancestor-descendant relationships using linguistic and material culture data. Within evolutionary archaeology specifically, phylogenetic analyses based on maximum parsimony and discrete traits dominate, which sets limitations for the downstream role cultural phylogenies, once derived, can play in more elaborate analytical pipelines. Recent methodological advances in Bayesian phylogenetics, however, now allow us to infer evolutionary dynamics using continuous characters. Capitalizing on these developments, we here present an exploratory analysis of cultural macroevolution of projectile point shape evolution in the European Final Palaeolithic and earliest Mesolithic (approx. 15 000-11 000 BP) using a Bayesian phylodynamic approach and the fossilized birth-death process model. This model-based approach leaps far beyond the application of parsimony, in that it not only produces a tree, but also divergence times, and diversification rates while incorporating uncertainties. This allows us to compare rates to the pronounced climatic changes that occurred during our time frame. While common in cultural evolutionary analyses of language, the extension of Bayesian phylodynamic models to archaeology arguably represents a major methodological breakthrough.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Cell type evolution reconstruction across species through cell phylogenies of single-cell RNA sequencing data

The origin and evolution of cell types has emerged as a key topic in evolutionary biology. Driven by rapidly accumulating single-cell datasets, recent attempts to infer cell type evolution have largely been limited to pairwise comparisons because we lack approaches to build cell phylogenies using model-based approaches. Here we approach the challenges of applying explicit phylogenetic methods to single-cell data by using principal components as phylogenetic characters. We infer a cell phylogeny from a large, comparative single-cell dataset of eye cells from five distantly related mammals. Robust cell type clades enable us to provide a phylogenetic, rather than phenetic, definition of cell type, allowing us to forgo marker genes and phylogenetically classify cells by topology. We further observe evolutionary relationships between diverse vessel endothelia and identify the myelinating and non-myelinating Schwann cells as sister cell types. Finally, we examine principal component loadings and describe the gene expression dynamics underlying the function and identity of cell type clades that have been conserved across the five species. A cell phylogeny provides a rigorous framework towards investigating the evolutionary history of cells and will be critical to interpret comparative single-cell datasets that aim to ask fundamental evolutionary questions.

Multiple modes of inference reveal less phylogenetic signal in marsupial basicranial shape compared with the rest of the cranium

Incorporating morphological data into modern phylogenies allows integration of fossil evidence, facilitating divergence dating and macroevolutionary inferences. Improvements in the phylogenetic utility of morphological data have been sought via Procrustes-based geometric morphometrics (GMM), but with mixed success and little clarity over what anatomical areas are most suitable. Here, we assess GMM-based phylogenetic reconstructions in a heavily sampled source of discrete characters for mammalian phylogenetics-the basicranium-in 57 species of marsupial mammals, compared with the remainder of the cranium. We show less phylogenetic signal in the basicranium compared with a 'Rest of Cranium' partition, using diverse metrics of phylogenetic signal (Kmult, phylogenetically aligned principal components analysis, comparisons of UPGMA/neighbour-joining/parsimony trees and cophenetic distances to a reference phylogeny) for scaled, Procrustes-aligned landmarks and allometry-corrected residuals. Surprisingly, a similar pattern emerged from parsimony-based analyses of discrete cranial characters. The consistent results across methods suggest that easily computed metrics such as Kmult can provide good guidance on phylogenetic information in a landmarking configuration. In addition, GMM data may be less informative for intricate but conservative anatomical regions such as the basicranium, while better-but not necessarily novel-phylogenetic information can be expected for broadly characterized shapes such as entire bones. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Fossil berries reveal global radiation of the nightshade family by the early Cenozoic

Fossil discoveries can transform our understanding of plant diversification over time and space. Recently described fossils in many plant families have pushed their known records farther back in time, pointing to alternative scenarios for their origin and spread. Here, we describe two new Eocene fossil berries of the nightshade family (Solanaceae) from the Esmeraldas Formation in Colombia and the Green River Formation in Colorado (USA). The placement of the fossils was assessed using clustering and parsimony analyses based on 10 discrete and five continuous characters, which were also scored in 291 extant taxa. The Colombian fossil grouped with members of the tomatillo subtribe, and the Coloradan fossil aligned with the chili pepper tribe. Along with two previously reported early Eocene fossils from the tomatillo genus, these findings indicate that Solanaceae were distributed at least from southern South America to northwestern North America by the early Eocene. Together with two other recently discovered Eocene berries, these fossils demonstrate that the diverse berry clade and, in turn, the entire nightshade family, is much older and was much more widespread in the past than previously thought.

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.

Morphometric analyses of the vertebrae of Ambystoma (Tschudi, 1838) and the implications for identification of fossil salamanders

Ambystoma (Tschudi, 1838) represents a speciose clade of salamanders that are found across much of North America. Fossils referred to Ambystoma are reported from early Cenozoic deposits and are common in Quaternary fossil deposits. Most fossils identified as Ambystoma are isolated vertebrae. Both quantitative and qualitative characters were reported as being useful for identifying fossilized vertebrae of Ambystoma below the genus level. However, there is limited information on intraspecific variation in those characters and previous studies noted intracolumnar variation which affects the utility of those characters for fossil identification. A lack of understanding of variation in modern species of Ambystoma casts uncertainty on our ability to identify fossil vertebrae confidently. We aimed to document intraspecific and intracolumnar variation in vertebral morphology among species of Ambystoma and examine the implications for fossil identification. We assembled one of the largest skeletal datasets for Ambystoma and took linear measurements on 15 species. We used 2D geometric morphometric analyses to characterize atlantal shape variation in Ambystoma. We apply those morphometric data in a case study where we identify fossil vertebrae from Hall's Cave, a Quaternary fossil locality in central Texas. We found patterns of intraspecific and intracolumnar variation that have substantial implications for fossil identification. Classification accuracies for species and clades within Ambystoma varied considerably. Overall classification accuracies based on size adjusted measurements and 2D geometric morphometric landmarks were lower compared to classifications from non-size adjusted linear measurements and were similar to accuracies based on size adjusted linear measurements. We identified fossil vertebrae from our case study as likely belonging to the tiger salamander clade within Ambystoma, but found that some fossils with lower classification probabilities are of uncertain identity. We discuss biogeographic implications for our fossil identifications and comment on challenges and next steps for advancing our understanding of morphological variation in Ambystoma.

SPASOS 1.1: a program for the inference of ancestral shape ontogenies

We recently published a method to infer ancestral landmark-based shape ontogenies that takes into account the possible existence of changes in developmental timing. Here we describe SPASOS, a software to perform that analysis. SPASOS is an open-source Windows program written in C. Input data include landmark coordinates for each specimen -with the corresponding information about developmental timing- and a phylogenetic tree showing the relationships among the species sampled. As output, the program produces image files for an easy visualization of the results and data files useful for post-processing. The program incorporates an interpolating function, based on weighting moving averages, which allows analysis of data with scarce information along the ontogenetic trajectory. An empirical evaluation of this function showed its suitability to fill in incomplete ontogenetic trajectories. Finally, we present the results of a reanalysis in SPASOS of a published dataset, where changes in developmental timing were originally inferred by considering PCA scores as shape variables. Both approaches retrieved the same four largest changes in developmental timing, but differed in the ancestral shapes inferred.

The phylogenetic relationships of geoemydid turtles from the Eocene Messel Pit Quarry: a first assessment using methods for continuous and discrete characters

The geoemydid turtles of the Eocoene Messel Pit Quarry of Hesse, Germany, are part of a rich Western European fossil record of testudinoids. Originally referred to as "Ocadia" kehreri and "Ocadia" messeliana, their systematic relationships remain unclear. A previous study proposed that a majority of the Western European geoemydids, including the Messel geoemydids, are closely related to the Recent European representatives of the clade Mauremys. Another study hypothesised that the Western European geoemydid fauna is more phylogenetically diverse, and that the Messel geoemydids are closely related to the East Asian turtles Orlitia and Malayemys. Here we present the first quantitative analyses to date that investigate this question. We use continuous characters in the form of ratios to estimate the placement of the Messel geoemydids in a reference tree that was estimated from molecular data. We explore the placement error obtained from that data with maximum likelihood and Bayesian methods, as well as linear parsimony in combination with discrete characters. We find good overall performance with Bayesian and parsimony analyses. Parsimony performs even better when we also incorporated discrete characters. Yet, we cannot pin down the position of the Messel geoemydids with high confidence. Depending on how intraspecific variation of the ratio characters is treated, parsimony favours a placement of the Messel fossils sister to Orlitia borneensis or sister to Geoemyda spengleri, with weak bootstrap support. The latter placement is suspect because G. spengleri is a phylogenetically problematic species with molecular and morphological data. There is even less support for placements within the Mauremys clade.

Evolution: Morphological saturation and release in mammals

Mammal evolution in the Mesozoic was thought to be heavily constrained by competition and predation by dinosaurs. A new study suggests that placental mammals remained constrained for several million years after non-avian dinosaurs perished, perhaps due to competition from archaic mammals.

Relative tooth size, Bayesian inference, and Homo naledi

OBJECTIVES

Size-corrected tooth crown measurements were used to estimate phenetic affinities among Homo naledi (~335-236 ka) and 11 other Plio-Pleistocene and recent species. To assess further their efficacy, and identify dental evolutionary trends, the data were then quantitatively coded for phylogenetic analyses. Results from both methods contribute additional characterization of H. naledi relative to other hominins.

MATERIALS AND METHODS

After division by their geometric mean, scaled mesiodistal and buccolingual dimensions were used in tooth size apportionment analysis to compare H. naledi with Australopithecus africanus, A. afarensis, Paranthropus robustus, P. boisei, H. habilis, H. ergaster, H. erectus, H. heidelbergensis, H. neanderthalensis, H. sapiens, and Pan troglodytes. These data produce equivalently scaled samples unaffected by interspecific size differences. The data were then gap-weighted for Bayesian inference.

RESULTS

Congruence in interspecific relationships is evident between methods, and with many inferred from earlier systematic studies. However, the present results place H. naledi as a sister taxon to H. habilis, based on a symplesiomorphic pattern of relative tooth size. In the preferred Bayesian phylogram, H. naledi is nested within a clade comprising all Homo species, but it shares some characteristics with australopiths and, particularly, early Homo.

DISCUSSION

Phylogenetic analyses of relative tooth size yield information about evolutionary dental trends not previously reported in H. naledi and the other hominins. Moreover, with an appropriate model these data recovered plausible evolutionary relationships. Together, the findings support recent study suggesting H. naledi originated long before the geological date of the Dinaledi Chamber, from which the specimens under study were recovered.

Investigating Morphological Complexes Using Informational Dissonance and Bayes Factors: A Case Study in Corbiculate Bees

It is widely recognized that different regions of a genome often have different evolutionary histories and that ignoring this variation when estimating phylogenies can be misleading. However, the extent to which this is also true for morphological data is still largely unknown. Discordance among morphological traits might plausibly arise due to either variable convergent selection pressures or else phenomena such as hemiplasy. Here, we investigate patterns of discordance among 282 morphological characters, which we scored for 50 bee species particularly targeting corbiculate bees, a group that includes the well-known eusocial honeybees and bumblebees. As a starting point for selecting the most meaningful partitions in the data, we grouped characters as morphological modules, highly integrated trait complexes that as a result of developmental constraints or coordinated selection we expect to share an evolutionary history and trajectory. In order to assess conflict and coherence across and within these morphological modules, we used recently developed approaches for computing Bayesian phylogenetic information allied with model comparisons using Bayes factors. We found that despite considerable conflict among morphological complexes, accounting for among-character and among-partition rate variation with individual gamma distributions, rate multipliers, and linked branch lengths can lead to coherent phylogenetic inference using morphological data. We suggest that evaluating information content and dissonance among partitions is a useful step in estimating phylogenies from morphological data, just as it is with molecular data. Furthermore, we argue that adopting emerging approaches for investigating dissonance in genomic datasets may provide new insights into the integration and evolution of anatomical complexes. [Apidae; entropy; morphological modules; phenotypic integration; phylogenetic information.].

Estimating Phylogenies from Shape and Similar Multidimensional Data: Why It Is Not Reliable

In recent years, there has been controversy whether multidimensional data such as geometric morphometric data or information on gene expression can be used for estimating phylogenies. This study uses simulations of evolution in multidimensional phenotype spaces to address this question and to identify specific factors that are important for answering it. Most of the simulations use phylogenies with four taxa, so that there are just three possible unrooted trees and the effect of different combinations of branch lengths can be studied systematically. In a comparison of methods, squared-change parsimony performed similarly well as maximum likelihood, and both methods outperformed Wagner and Euclidean parsimony, neighbor-joining and UPGMA. Under an evolutionary model of isotropic Brownian motion, phylogeny can be estimated reliably if dimensionality is high, even with relatively unfavorable combinations of branch lengths. By contrast, if there is phenotypic integration such that most variation is concentrated in one or a few dimensions, the reliability of phylogenetic estimates is severely reduced. Evolutionary models with stabilizing selection also produce highly unreliable estimates, which are little better than picking a phylogenetic tree at random. To examine how these results apply to phylogenies with more than four taxa, we conducted further simulations with up to eight taxa, which indicated that the effects of dimensionality and phenotypic integration extend to more than four taxa, and that convergence among internal nodes may produce additional complications specifically for greater numbers of taxa. Overall, the simulations suggest that multidimensional data, under evolutionary models that are plausible for biological data, do not produce reliable estimates of phylogeny. [Brownian motion; gene expression data; geometric morphometrics; morphological integration; squared-change parsimony; phylogeny; shape; stabilizing selection.].

Morphological and phylogeographic evidence for budding speciation: an example in hominins

Parametric phylogenetic approaches that attempt to delineate between distinct 'modes' of speciation (splitting cladogenesis, budding cladogenesis and anagenesis) between fossil taxa have become increasingly popular among comparative biologists. But it is not yet well understood how clearly morphological data from fossil taxa speak to detailed questions of speciation mode when compared with the lineage diversification models that serve as their basis. In addition, the congruence of inferences made using these approaches with geographical patterns has not been explored. Here, I extend a previously introduced maximum-likelihood approach for the examination of ancestor-descendant relationships to accommodate budding speciation and apply it to a dataset of fossil hominins. I place these results in a phylogeographic context to better understand spatial dynamics underlying the hypothesized speciation patterns. The spatial patterns implied by the phylogeny hint at the complex demographic processes underlying the spread and diversification of hominins throughout the Pleistocene. I also find that inferences of budding are driven primarily by stratigraphic, versus morphological, data and discuss the ramifications for interpretations of speciation process in hominins specifically and from phylogenetic data in general.

Linking morphological and molecular sources to disentangle the case of Xylodon australis

The use of different sources of evidence has been recommended in order to conduct species delimitation analyses to solve taxonomic issues. In this study, we use a maximum likelihood framework to combine morphological and molecular traits to study the case of Xylodon australis (Hymenochaetales, Basidiomycota) using the locate.yeti function from the phytools R package. Xylodon australis has been considered a single species distributed across Australia, New Zealand and Patagonia. Multi-locus phylogenetic analyses were conducted to unmask the actual diversity under X. australis as well as the kinship relations respect their relatives. To assess the taxonomic position of each clade, locate.yeti function was used to locate in a molecular phylogeny the X. australis type material for which no molecular data was available using morphological continuous traits. Two different species were distinguished under the X. australis name, one from Australia–New Zealand and other from Patagonia. In addition, a close relationship with Xylodon lenis, a species from the South East of Asia, was confirmed for the Patagonian clade. We discuss the implications of our results for the biogeographical history of this genus and we evaluate the potential of this method to be used with historical collections for which molecular data is not available.

Morphological Phylogenetics Evaluated Using Novel Evolutionary Simulations

Evolutionary inferences require reliable phylogenies. Morphological data have traditionally been analyzed using maximum parsimony, but recent simulation studies have suggested that Bayesian analyses yield more accurate trees. This debate is ongoing, in part, because of ambiguity over modes of morphological evolution and a lack of appropriate models. Here, we investigate phylogenetic methods using two novel simulation models-one in which morphological characters evolve stochastically along lineages and another in which individuals undergo selection. Both models generate character data and lineage splitting simultaneously: the resulting trees are an emergent property, rather than a fixed parameter. Standard consensus methods for Bayesian searches (Mki) yield fewer incorrect nodes and quartets than the standard consensus trees recovered using equal weighting and implied weighting parsimony searches. Distances between the pool of derived trees (most parsimonious or posterior distribution) and the true trees-measured using Robinson-Foulds (RF), subtree prune and regraft (SPR), and tree bisection reconnection (TBR) metrics-demonstrate that this is related to the search strategy and consensus method of each technique. The amount and structure of homoplasy in character data differ between models. Morphological coherence, which has previously not been considered in this context, proves to be a more important factor for phylogenetic accuracy than homoplasy. Selection-based models exhibit relatively lower homoplasy, lower morphological coherence, and higher inaccuracy in inferred trees. Selection is a dominant driver of morphological evolution, but we demonstrate that it has a confounding effect on numerous character properties which are fundamental to phylogenetic inference. We suggest that the current debate should move beyond considerations of parsimony versus Bayesian, toward identifying modes of morphological evolution and using these to build models for probabilistic search methods. [Bayesian; evolution; morphology; parsimony; phylogenetics; selection; simulation.].

Mosaic evolution, preadaptation, and the evolution of evolvability in apes

A major goal in postsynthesis evolutionary biology has been to better understand how complex interactions between traits drive movement along and facilitate the formation of distinct evolutionary pathways. I present analyses of a character matrix sampled across the haplorrhine skeleton that revealed several modules of characters displaying distinct patterns in macroevolutionary disparity. Comparison of these patterns to those in neurological development showed that early ape evolution was characterized by an intense regime of evolutionary and developmental flexibility. Shifting and reduced constraint in apes was met with episodic bursts in phenotypic innovation that built a wide array of functional diversity over a foundation of shared developmental and anatomical structure. Shifts in modularity drove dramatic evolutionary changes across the ape body plan in two distinct ways: (1) an episode of relaxed integration early in hominoid evolution coincided with bursts in evolutionary rate across multiple character suites; (2) the formation of two new trait modules along the branch leading to chimps and humans preceded rapid and dramatic evolutionary shifts in the carpus and pelvis. Changes to the structure of evolutionary mosaicism may correspond to enhanced evolvability that has a "preadaptive" effect by catalyzing later episodes of dramatic morphological remodeling.

Detecting Mosaic Patterns in Macroevolutionary Disparity

Evolutionary biologists have long sought to understand the full complexity in pattern and process that shapes organismal diversity. Although phylogenetic comparative methods are often used to reconstruct complex evolutionary dynamics, they are typically limited to a single phenotypic trait. Extensions that accommodate multiple traits lack the ability to partition multidimensional data sets into a set of mosaic suites of evolutionarily linked characters. I introduce a comparative framework that identifies heterogeneity in evolutionary patterns across large data sets of continuous traits. Using a model of continuous trait evolution based on the differential accumulation of disparity across lineages in a phylogeny, the approach algorithmically partitions traits into a set of character suites that best explains the data, where each suite displays a distinct pattern in phylogenetic morphological disparity. When applied to empirical data, the approach revealed a mosaic pattern predicted by developmental biology. The evolutionary distinctiveness of individual suites can be investigated in more detail either by fitting conventional comparative models or by directly studying the phylogenetic patterns in disparity recovered during the analysis. This framework can supplement existing comparative approaches by inferring the complex, integrated patterns that shape evolution across the body plan from disparate developmental, morphometric, and environmental sources of phenotypic data.

Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation

Big, time-scaled phylogenies are fundamental to connecting evolutionary processes to modern biodiversity patterns. Yet inferring reliable phylogenetic trees for thousands of species involves numerous trade-offs that have limited their utility to comparative biologists. To establish a robust evolutionary timescale for all approximately 6,000 living species of mammals, we developed credible sets of trees that capture root-to-tip uncertainty in topology and divergence times. Our "backbone-and-patch" approach to tree building applies a newly assembled 31-gene supermatrix to two levels of Bayesian inference: (1) backbone relationships and ages among major lineages, using fossil node or tip dating, and (2) species-level "patch" phylogenies with nonoverlapping in-groups that each correspond to one representative lineage in the backbone. Species unsampled for DNA are either excluded ("DNA-only" trees) or imputed within taxonomic constraints using branch lengths drawn from local birth-death models ("completed" trees). Joining time-scaled patches to backbones results in species-level trees of extant Mammalia with all branches estimated under the same modeling framework, thereby facilitating rate comparisons among lineages as disparate as marsupials and placentals. We compare our phylogenetic trees to previous estimates of mammal-wide phylogeny and divergence times, finding that (1) node ages are broadly concordant among studies, and (2) recent (tip-level) rates of speciation are estimated more accurately in our study than in previous "supertree" approaches, in which unresolved nodes led to branch-length artifacts. Credible sets of mammalian phylogenetic history are now available for download at http://vertlife.org/phylosubsets, enabling investigations of long-standing questions in comparative biology.

Estimating the phylogeny of geoemydid turtles (Cryptodira) from landmark data: an assessment of different methods

BACKGROUND

In the last 20 years, a general picture of the evolutionary relationships between geoemydid turtles (ca. 70 species distributed over the Northern hemisphere) has emerged from the analysis of molecular data. However, there is a paucity of good traditional morphological characters that correlate with the phylogeny, which are essential for the robust integration of fossil and molecular data. Part of this problem might be due to intrinsic limitations of traditional discrete characters. Here, we explore the use of continuous data in the form of 3D coordinates of homologous landmarks on the turtle shell for phylogenetic inference and the phylogenetic placement of single species on a scaffold molecular tree. We focus on the performance yielded by sampling the carapace and/or plastral lobes and using various phylogenetic methods.

METHODS

We digitised the landmark coordinates of the carapace and plastron of 42 and 46 extant geoemydid species, respectively. The configurations were superimposed and we estimated the phylogenetic tree of geoemydids with landmark analysis under parsimony, traditional Farris parsimony, unweighted squared-change parsimony, maximum likelihood with a Brownian motion model, and neighbour-joining on a matrix of pairwise Procrustes distances. We assessed the performance of those analyses by comparing the trees against a reference phylogeny obtained from seven molecular markers. For comparisons between trees we used difference measures based on quartets and splits. We used the same reference tree to evaluate phylogenetic placement performance by a leave-one-out validation procedure.

RESULTS

Whatever method we used, similarity to the reference phylogeny was low. The carapace alone gave slightly better results than the plastron or the complete shell. Assessment of the potential for placement of single species on the reference tree with landmark data gave much better results, with similar accuracy and higher precision compared to the performance of discrete characters with parsimony.

A Systematist's Guide to Estimating Bayesian Phylogenies From Morphological Data

Phylogenetic trees are crucial to many aspects of taxonomic and comparative biology. Many researchers have adopted Bayesian methods to estimate their phylogenetic trees. In this family of methods, a model of morphological evolution is assumed to have generated the data observed by the researcher. These models make a variety of assumptions about the evolution of morphological characters, and these assumptions are translated into mathematics as parameters. The incorporation of prior distributions further allows researchers to quantify their prior beliefs about the value any one parameter can take. How to translate biological knowledge into mathematical language is difficult, and can be confusing to many biologists. This review aims to help systematics researchers understand the biological meaning of common models and assumptions. Using examples from the insect fossil record, I will demonstrate empirically what assumptions mean in concrete terms, and discuss how researchers can use and understand Bayesian methods for phylogenetic estimation.

Bayesian Estimation of Species Divergence Times Using Correlated Quantitative Characters

Discrete morphological data have been widely used to study species evolution, but the use of quantitative (or continuous) morphological characters is less common. Here, we implement a Bayesian method to estimate species divergence times using quantitative characters. Quantitative character evolution is modeled using Brownian diffusion with character correlation and character variation within populations. Through simulations, we demonstrate that ignoring the population variation (or population “noise”) and the correlation among characters leads to biased estimates of divergence times and rate, especially if the correlation and population noise are high. We apply our new method to the analysis of quantitative characters (cranium landmarks) and molecular data from carnivoran mammals. Our results show that time estimates are affected by whether the correlations and population noise are accounted for or ignored in the analysis. The estimates are also affected by the type of data analyzed, with analyses of morphological characters only, molecular data only, or a combination of both; showing noticeable differences among the time estimates. Rate variation of morphological characters among the carnivoran species appears to be very high, with Bayesian model selection indicating that the independent-rates model fits the morphological data better than the autocorrelated-rates model. We suggest that using morphological continuous characters, together with molecular data, can bring a new perspective to the study of species evolution. Our new model is implemented in the MCMCtree computer program for Bayesian inference of divergence times.

A new phylogenetic analysis of Phytosauria (Archosauria: Pseudosuchia) with the application of continuous and geometric morphometric character coding

Phytosauria is a clade of large, carnivorous, semi-aquatic archosauromorphs which reached its peak diversity and an almost global distribution in the Late Triassic (c. 230-201 Mya). Previous phylogenetic analyses of Phytosauria have either focused primarily on the relationships of specific subclades, or were limited in taxonomic scope, and no taxonomically comprehensive dataset is currently available. We here present the most taxonomically comprehensive cladistic dataset of phytosaurs to date, based on extensive first-hand study, identification of novel characters and synthesis of previous matrices. This results in an almost twofold increase in phylogenetic information scored per taxon over previous analyses. Alongside a traditional discrete character matrix, three variant matrices were analysed in which selected characters were coded using continuous and landmarking methods, to more rigorously explore phytosaur relationships. Based on these four data matrices, four tree topologies were recovered. Relationships among non-leptosuchomorph phytosaurs are largely consistent between these four topologies, whereas those of more derived taxa are more variable. Rutiodon carolinensis consistently forms a sister relationship with Angistorhinus. In three topologies Nicrosaurus nests deeply within a group of traditionally non-Mystriosuchini taxa, leading us to redefine Mystriosuchini by excluding Nicrosaurus as an internal specifier. Two distinct patterns of relationships within Mystriosuchini are present in the four topologies, distinguished largely by the variable position of Mystriosuchus. In two topologies Mystriosuchus forms the most basal clade in Mystriosuchini, whilst in the others it occupies a highly derived position within the Machaeroprosopus clade. 'Redondasaurus' is consistently recovered as monophyletic; however, it also nests within the Machaeroprosopus clade. The greatest impact on tree topology was associated with the incorporation of continuous data into our matrices, with landmark characters exerting a relatively modest influence. All topologies correlated significantly with stratigraphic range estimates. Topological variability in our results highlights clades in which further investigation may better elucidate phytosaur relationships.