Placing cryptic, recently extinct, or hypothesized taxa into an ultrametric phylogeny using continuous character data: A case study with the lizardAnolis roosevelti

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.

phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things)

Phylogenetic comparative methods comprise the general endeavor of using an estimated phylogenetic tree (or set of trees) to make secondary inferences: about trait evolution, diversification dynamics, biogeography, community ecology, and a wide range of other phenomena or processes. Over the past ten years or so, the phytools R package has grown to become an important research tool for phylogenetic comparative analysis. phytools is a diverse contributed R library now consisting of hundreds of different functions covering a variety of methods and purposes in phylogenetic biology. As of the time of writing, phytools included functionality for fitting models of trait evolution, for reconstructing ancestral states, for studying diversification on trees, and for visualizing phylogenies, comparative data, and fitted models, as well numerous other tasks related to phylogenetic biology. Here, I describe some significant features of and recent updates to phytools, while also illustrating several popular workflows of the phytools computational software.

Mosaic adaptive peak shifts underlie body shape diversification in pelagiarian fishes (Acanthomorpha: Percomorpha)

Extreme body elongation in fishes is a major evolutionary transformation that extends the boundaries of morphological diversity and alters aspects of function, behaviour and ecology. Prior studies have identified features of the cranial and axial skeleton that characterize elongate fishes, but a lack of detailed reconstructions of anatomical evolution has limited inferences about factors that underlie major shifts in body shape. In this study, we fitted multi-peak adaptive (Ornstein–Uhlenbeck) evolutionary models to species body shape and anatomical dimensions in Pelagiaria, a radiation of open-ocean fishes whose species span a continuum from deep bodied to highly elongate. We inferred an ancestral fusiform adaptive peak that is retained by several major pelagiarian lineages (e.g. Scombridae) and found robust support for multiple transitions to deep-bodied optima (in the families Stromateidae, Bramidae and Caristiidae) and elongate-bodied optima (within Trichiuroidei), including two instances of sequential shifts towards increasingly elongate optima that followed distinct paths of anatomical evolution. Within Trichiuridae, initial increases in head length and the number of vertebrae were followed by changes in head and vertebral shape. Within an elongate-bodied subclade of taxa traditionally identified as ‘gempylids’, changes in head and vertebral shape and in the number of precaudal vertebrae preceded an increase in the number of caudal vertebrae. Altogether, this mosaic of anatomical peak shifts suggests that body shape transformations were associated with differing selective demands and developmental changes.

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.

Body shape transformations by alternate anatomical adaptive peak shifts in blenniiform fishes

Extreme body elongation has occurred repeatedly in the evolutionary history of ray-finned fishes. Lengthening of the anterior-posterior body axis relative to depth and width can involve changes in the cranial skeleton and vertebral column, but to what extent is anatomical evolution determined by selective factors and intrinsic constraints that are shared broadly among closely related lineages? In this study, we fit adaptive (Ornstein-Uhlenbeck) evolutionary models to body shape and its anatomical determinants and identified two instances of extreme elongation by divergent anatomical peak shifts in the Blenniiformes, a radiation of small-bodied substrate-associated marine teleost fishes. Species in the genus Xiphasia (hairtail blennies) evolved toward a peak defined by a highly elongated caudal vertebral region but ancestral cranial and precaudal vertebral morphology. In contrast, a clade that includes the genera Chaenopsis and Lucayablennius (pike and arrow blennies) evolved toward a peak with a long slender skull but ancestral axial skeletal anatomy. Neither set of anatomical peak shifts aligns closely with the major axis of anatomical diversification in other blenniiform fishes. These results provide little evidence that ancestral constraints have affected body shape transformation, and instead suggest that extreme elongation arose with distinct shifts in selective factors and development.

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.

Estimating Diversification Rates on Incompletely Sampled Phylogenies: Theoretical Concerns and Practical Solutions

Molecular phylogenies are a key source of information about the tempo and mode of species diversification. However, most empirical phylogenies do not contain representatives of all species, such that diversification rates are typically estimated from incompletely sampled data. Most researchers recognize that incomplete sampling can lead to biased rate estimates, but the statistical properties of methods for accommodating incomplete sampling remain poorly known. In this point of view, we demonstrate theoretical concerns with the widespread use of analytical sampling corrections for sparsely sampled phylogenies of higher taxonomic groups. In particular, corrections based on "sampling fractions" can lead to low statistical power to infer rate variation when it is present, depending on the likelihood function used for inference. In the extreme, the sampling fraction correction can lead to spurious patterns of diversification that are driven solely by unbalanced sampling across the tree in concert with low overall power to infer shifts. Stochastic polytomy resolution provides an alternative to sampling fraction approaches that avoids some of these biases. We show that stochastic polytomy resolvers can greatly improve the power of common analyses to estimate shifts in diversification rates. We introduce a new stochastic polytomy resolution method (Taxonomic Addition for Complete Trees [TACT]) that uses birth-death-sampling estimators across an ultrametric phylogeny to estimate branching times for unsampled taxa, with taxonomic information to compatibly place new taxa onto a backbone phylogeny. We close with practical recommendations for diversification inference under several common scenarios of incomplete sampling. [Birth-death process; diversification; incomplete sampling; phylogenetic uncertainty; rate heterogeneity; rate shifts; stochastic polytomy resolution.].

Bayesian placement of fossils on phylogenies using quantitative morphometric data

Jointly developing a comprehensive tree of life from living and fossil taxa has long been a fundamental goal in evolutionary biology. One major challenge has stemmed from difficulties in merging evidence from extant and extinct organisms. While these efforts have resulted in varying stages of synthesis, they have been hindered by their dependence on qualitative descriptions of morphology. Though rarely applied to phylogenetic inference, traditional and geometric morphometric data can improve these issues by generating more rigorous ways to quantify variation in morphological structures. They may also facilitate the rapid and objective aggregation of large morphological datasets. I describe a new Bayesian method that leverages quantitative trait data to reconstruct the positions of fossil taxa on fixed reference trees composed of extant taxa. Unlike most formulations of phylogenetic Brownian motion models, this method expresses branch lengths in units of morphological disparity, suggesting a new framework through which to construct Bayesian node calibration priors for molecular dating and explore comparative patterns in morphological disparity. I am hopeful that the approach described here will help to facilitate a deeper integration of neo- and paleontological data to move morphological phylogenetics further into the genomic era.

Brain size evolution in pipefishes and seahorses: the role of feeding ecology, life history and sexual selection

Brain size varies greatly at all taxonomic levels. Feeding ecology, life history and sexual selection have been proposed as key components in generating contemporary diversity in brain size across vertebrates. Analyses of brain size evolution have, however, been limited to lineages where males predominantly compete for mating and females choose mates. Here, we present the first original data set of brain sizes in pipefishes and seahorses (Syngnathidae) a group in which intense female mating competition occurs in many species. After controlling for the effect of shared ancestry and overall body size, brain size was positively correlated with relative snout length. Moreover, we found that females, on average, had 4.3% heavier brains than males and that polyandrous species demonstrated more pronounced (11.7%) female-biased brain size dimorphism. Our results suggest that adaptations for feeding on mobile prey items and sexual selection in females are important factors in brain size evolution of pipefishes and seahorses. Most importantly, our study supports the idea that sexual selection plays a major role in brain size evolution, regardless of on which sex sexual selection acts stronger.