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

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.

Metalloproteinases in Restorative Dentistry: An In Silico Study toward an Ideal Animal Model

In dentistry, various animal models are used to evaluate adhesive systems, dental caries and periodontal diseases. Metalloproteinases (MMPs) are enzymes that degrade collagen in the dentin matrix and are categorized in over 20 different classes. Collagenases and gelatinases are intrinsic constituents of the human dentin organic matrix fibrillar network and are the most abundant MMPs in this tissue. Understanding such enzymes' action on dentin is important in the development of approaches that could reduce dentin degradation and provide restorative procedures with extended longevity. This in silico study is based on dentistry's most used animal models and intends to search for the most suitable, evolutionarily close to Homo sapiens. We were able to retrieve 176,077 mammalian MMP sequences from the UniProt database. These sequences were manually curated through a three-step process. After such, the remaining 3178 sequences were aligned in a multifasta file and phylogenetically reconstructed using the maximum likelihood method. Our study inferred that the animal models most evolutionarily related to Homo sapiens were Orcytolagus cuniculus (MMP-1 and MMP-8), Canis lupus (MMP-13), Rattus norvegicus (MMP-2) and Orcytolagus cuniculus (MMP-9). Further research will be needed for the biological validation of our findings.

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'.

Species delimitation, molecular phylogeny and historical biogeography of the sweetlips fish (Perciformes, Haemulidae)

The subfamily Plectorhinchinae (sweetlips) is composed of poorly-known species with high commercially and ecologically values that exhibit phenotypic plasticity and various morphologies. Few studies have assessed the validity of sweetlips, intergeneric relationships and evolutionary survey in this subfamily, which have not yet been resolved. This study investigated the DNA sequences of (1) the mitochondrial COI gene to delimit species, and (2) two mitochondrial (COI and Cyt b), and one nuclear (RAG1) markers to infer phylogenetic relationships and evolutionary and biogeographic history. The molecular results could differentiate Diagramma punctatum from the other species, but failed to distinguish D. labiosum as a distinct species with considerably lower genetic distances for the COI (0.53%) and Cyt b (0.51%) markers. However, additional taxonomic investigations are required to shed light on this issue. All previously described nominal species of sweetlips in the northwest Indian Ocean were found to be well supported. The monophyly of Plectorhinchus is not supported and Diagramma pictum and D. punctatum should be assigned to the genus Plectorhinchus. The biogeographic history of Plectorhinchinae likely originated in the Indo-Pacific ca. 34 Ma (30–39 Ma; late Eocene/ middle Oligocene) and subsequently colonised the Western Indian Ocean and the Central Indo-Pacific. Maximum diversification within the subfamily occurred from the middle Miocene to Pliocene, coinciding with dispersal and vicariance events. Diversification was probably driven by both biological and geographical factors.

Methods for studying allometry in geometric morphometrics: a comparison of performance

Allometry has been the focus of growing interest in studies using geometric morphometric methods to address a wide range of research questions at the interface of ecology and evolution. This study uses computer simulations to compare four methods for estimating allometric vectors from landmark data: the multivariate regression of shape on a measure of size, the first principal component (PC1) of shape, the PC1 in conformation space, and a recently proposed method, the PC1 of Boas coordinates. Simulations with no residual variation around the allometric relationship showed that all four methods are logically consistent with one another, up to minor nonlinearities in the mapping between conformation space and shape tangent space. In simulations that included residual variation, either isotropic or with a pattern independent of allometry, regression of shape on size performed consistently better than the PC1 of shape. The PC1s of conformation and of Boas coordinates were very similar and very close to the simulated allometric vectors under all conditions. An extra series of simulations to elucidate the relation between conformation and Boas coordinates indicated that they are almost identical, with a marginal advantage for conformation. Empirical examples of ontogenetic allometry in rat skulls and rockfish body shape illustrate simple biological applications of the methods. The paper concludes with recommendations how these methods for estimating allometry can be used in studies of evolution and ecology.

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.

Covariation between the cranium and the cervical vertebrae in hominids

The analysis of patterns of integration is crucial for the reconstruction and understanding of how morphological changes occur in a taxonomic group throughout evolution. These patterns are relatively constant; however, both patterns and the magnitudes of integration may vary across species. These differences may indicate morphological diversification, in some cases related to functional adaptations to the biomechanics of organisms. In this study, we analyze patterns of integration between two functional and developmental structures, the cranium and the cervical spine in hominids, and we quantify the amount of divergence of each anatomical element through phylogeny. We applied these methods to three-dimensional data from 168 adult hominid individuals, summing a total of more than 1000 cervical vertebrae. We found the atlas (C1) and axis (C2) display the lowest covariation with the cranium in hominids (Homo sapiens, Pan troglodytes, Pan paniscus, Gorilla gorilla, Gorilla beringei, Pongo pygmaeus). H. sapiens show a relatively different pattern of craniocervical correlation compared with chimpanzees and gorillas, especially in variables implicated in maintaining the balance of the head. Finally, the atlas and axis show lower magnitude of shape change during evolution than the rest of the cervical vertebrae, especially those located in the middle of the subaxial cervical spine. Overall, results suggest that differences in the pattern of craniocervical correlation between humans and gorillas and chimpanzees could reflect the postural differences between these groups. Also, the stronger craniocervical integration and larger magnitude of shape change during evolution shown by the middle cervical vertebrae suggests that they have been selected to play an active role in maintaining head balance.

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.

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.

How Exactly Did the Nose Get That Long? A Critical Rethinking of the Pinocchio Effect and How Shape Changes Relate to Landmarks

The Pinocchio effect has long been discussed in the literature on geometric morphometrics. It denotes the observation that Procrustes superimposition tends to distribute shape changes over many landmarks, even though a different superimposition may exist for the same landmark configurations that concentrates changes in just one or a few landmarks. This is widely seen as a flaw of Procrustes methods. Visualizations illustrating the Pinocchio effect use a comparison of the same pair of shapes superimposed in two different ways: in a resistant-fit superimposition that concentrates the shape difference in just one or a few landmarks, and in Procrustes superimposition, which distributes differences over most or all landmarks. Because these superimpositions differ only in the non-shape aspects of size, position and orientation, they are equivalent from the perspective of shape analysis. Simulation studies of the Pinocchio effect usually generate data, either single pairs or larger samples of landmark configurations, in a particular superimposition so that differences occur mostly or exclusively at just one or a few landmarks, but no steps are taken to remove variation from size, position and orientation. When these configurations are then compared with Procrustes-superimposed data, differences appear and are attributed to the Pinocchio effect. Overall, it is ironic that all manifestations of the Pinocchio effect in one way or another rely on differences in the non-shape properties of position and orientation. Rigorous thinking about shape variation and careful choice of visualization methods can prevent confusion over this issue.