Identifying unstable taxa: Efficient implementation of triplet-based measures of stability, and comparison with Phyutility and RogueNaRok

A new abelisaurid dinosaur from the end Cretaceous of Patagonia and evolutionary rates among the Ceratosauria

Gondwanan dinosaur faunae during the 20 Myr preceding the Cretaceous-Palaeogene (K/Pg) extinction included several lineages that were absent or poorly represented in Laurasian landmasses. Among these, the South American fossil record contains diverse abelisaurids, arguably the most successful groups of carnivorous dinosaurs from Gondwana in the Cretaceous, reaching their highest diversity towards the end of this period. Here we describe Koleken inakayali gen. et sp. n., a new abelisaurid from the La Colonia Formation (Maastrichtian, Upper Cretaceous) of Patagonia. Koleken inakayali is known from several skull bones, an almost complete dorsal series, complete sacrum, several caudal vertebrae, pelvic girdle and almost complete hind limbs. The new abelisaurid shows a unique set of features in the skull and several anatomical differences from Carnotaurus sastrei (the only other abelisaurid known from the La Colonia Formation). Koleken inakayali is retrieved as a brachyrostran abelisaurid, clustered with other South American abelisaurids from the latest Cretaceous (Campanian-Maastrichtian), such as Aucasaurus, Niebla and Carnotaurus. Leveraging our phylogeny estimates, we explore rates of morphological evolution across ceratosaurian lineages, finding them to be particularly high for elaphrosaurine noasaurids and around the base of Abelisauridae, before the Early Cretaceous radiation of the latter clade. The Noasauridae and their sister clade show contrasting patterns of morphological evolution, with noasaurids undergoing an early phase of accelerated evolution of the axial and hind limb skeleton in the Jurassic, and the abelisaurids exhibiting sustained high rates of cranial evolution during the Early Cretaceous. These results provide much needed context for the evolutionary dynamics of ceratosaurian theropods, contributing to broader understanding of macroevolutionary patterns across dinosaurs.

Dating the evolution of the complex thalloid liverworts (Marchantiopsida): total-evidence dating analysis supports a Late Silurian-Early Devonian origin and post-Mesozoic morphological stasis

Divergence times based on molecular clock analyses often differ from those derived from total-evidence dating (TED) approaches. For bryophytes, fossils have been excluded from previous assessments of divergence times, and thus, their utility in dating analyses remains unexplored. Here, we conduct the first TED analyses of the complex thalloid liverworts (Marchantiopsida) that include fossils and evaluate macroevolutionary trends in morphological 'diversity' (disparity) and rates. Phylogenetic analyses were performed on a combined dataset of 130 discrete characters and 11 molecular markers (sampled from nuclear, plastid and mitochondrial genomes). Taxon sampling spanned 56 extant species - representing all the orders within Marchantiophyta and extant genera within Marchantiales - and eight fossil taxa. Total-evidence dating analyses support the radiation of Marchantiopsida during Late Silurian-Early Devonian (or Middle Ordovician when the outgroup is excluded) and that of Ricciaceae in the Middle Jurassic. Morphological change rate was high early in the history of the group, but it barely increased after Late Cretaceous. Disparity-through-time analyses support a fast increase in diversity until the Middle Triassic (c. 250 Ma), after which phenotypic evolution slows down considerably. Incorporating fossils in analyses challenges previous assumptions on the affinities of extinct taxa and indicates that complex thalloid liverworts radiated c. 125 Ma earlier than previously inferred.

Identifying and addressing methodological incongruence in phylogenomics: A review

The availability of phylogenetic data has greatly expanded in recent years. As a result, a new era in phylogenetic analysis is dawning-one in which the methods we use to analyse and assess our data are the bottleneck to producing valuable phylogenetic hypotheses, rather than the need to acquire more data. This makes the ability to accurately appraise and evaluate new methods of phylogenetic analysis and phylogenetic artefact identification more important than ever. Incongruence in phylogenetic reconstructions based on different datasets may be due to two major sources: biological and methodological. Biological sources comprise processes like horizontal gene transfer, hybridization and incomplete lineage sorting, while methodological ones contain falsely assigned data or violations of the assumptions of the underlying model. While the former provides interesting insights into the evolutionary history of the investigated groups, the latter should be avoided or minimized as best as possible. However, errors introduced by methodology must first be excluded or minimized to be able to conclude that biological sources are the cause. Fortunately, a variety of useful tools exist to help detect such misassignments and model violations and to apply ameliorating measurements. Still, the number of methods and their theoretical underpinning can be overwhelming and opaque. Here, we present a practical and comprehensive review of recent developments in techniques to detect artefacts arising from model violations and poorly assigned data. The advantages and disadvantages of the different methods to detect such misleading signals in phylogenetic reconstructions are also discussed. As there is no one-size-fits-all solution, this review can serve as a guide in choosing the most appropriate detection methods depending on both the actual dataset and the computational power available to the researcher. Ultimately, this informed selection will have a positive impact on the broader field, allowing us to better understand the evolutionary history of the group of interest.

Isolated tooth reveals hidden spinosaurid dinosaur diversity in the British Wealden Supergroup (Lower Cretaceous)

Isolated spinosaurid teeth are relatively well represented in the Lower Cretaceous Wealden Supergroup of southern England, UK. Until recently it was assumed that these teeth were referable to Baryonyx, the type species (B. walkeri) and specimen of which is from the Barremian Upper Weald Clay Formation of Surrey. British spinosaurid teeth are known from formations that span much of the c. 25 Ma depositional history of the Wealden Supergroup, and recent works suggest that British spinosaurids were more taxonomically diverse than previously thought. On the basis of both arguments, it is appropriate to doubt the hypothesis that isolated teeth from outside the Upper Weald Clay Formation are referable to Baryonyx. Here, we use phylogenetic, discriminant and cluster analyses to test whether an isolated spinosaurid tooth (HASMG G369a, consisting of a crown and part of the root) from a non-Weald Clay Formation unit can be referred to Baryonyx. HASMG G369a was recovered from an uncertain Lower Cretaceous locality in East Sussex but is probably from a Valanginian exposure of the Hastings Group and among the oldest spinosaurid material known from the UK. Spinosaurid affinities are both quantitatively and qualitatively supported, and HASMG G369a does not associate with Baryonyx in any analysis. This supports recent reinterpretations of the diversity of spinosaurid in the Early Cretaceous of Britain, which appears to have been populated by multiple spinosaurid lineages in a manner comparable to coeval Iberian deposits. This work also reviews the British and global records of early spinosaurids (known mainly from dental specimens), and revisits evidence for post-Cenomanian spinosaurid persistence.

Experimental and computational analysis of the ancestry of an evolutionary young enzyme from histidine biosynthesis

The conservation of fold and chemistry of the enzymes associated with histidine biosynthesis suggests that this pathway evolved prior to the diversification of Bacteria, Archaea, and Eukaryotes. The only exception is the histidinol phosphate phosphatase (HolPase). So far, non-homologous HolPases that possess distinct folds and belong to three different protein superfamilies have been identified in various phylogenetic clades. However, their evolution has remained unknown to date. Here, we analyzed the evolutionary history of the HolPase from γ-Proteobacteria (HisB-N). It has been argued that HisB-N and its closest homologue d-glycero-d-manno-heptose-1,7-bisphosphate 7-phosphatase (GmhB) have emerged from the same promiscuous ancestral phosphatase. GmhB variants catalyze the hydrolysis of the anomeric d-glycero-d-manno-heptose-1,7-bisphosphate (αHBP or βHBP) with a strong preference for one anomer (αGmhB or βGmhB). We found that HisB-N from Escherichia coli shows promiscuous activity for βHBP but not αHBP, while βGmhB from Crassaminicella sp. shows promiscuous activity for HolP. Accordingly, a combined phylogenetic tree of αGmhBs, βGmhBs, and HisB-N sequences revealed that HisB-Ns form a compact subcluster derived from βGmhBs. Ancestral sequence reconstruction and in vitro analysis revealed a promiscuous HolPase activity in the resurrected enzymes prior to functional divergence of the successors. The following increase in catalytic efficiency of the HolP turnover is reflected in the shape and electrostatics of the active site predicted by AlphaFold. An analysis of the phylogenetic tree led to a revised evolutionary model that proposes the horizontal gene transfer of a promiscuous βGmhB from δ- to γ-Proteobacteria where it evolved to the modern HisB-N.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Africa's oldest dinosaurs reveal early suppression of dinosaur distribution

The vertebrate lineages that would shape Mesozoic and Cenozoic terrestrial ecosystems originated across Triassic Pangaea^1-11. By the Late Triassic (Carnian stage, ~235 million years ago), cosmopolitan 'disaster faunas' (refs. ^12-14) had given way to highly endemic assemblages^12,13 on the supercontinent. Testing the tempo and mode of the establishment of this endemism is challenging-there were few geographic barriers to dispersal across Pangaea during the Late Triassic. Instead, palaeolatitudinal climate belts, and not continental boundaries, are proposed to have controlled distribution^15-18. During this time of high endemism, dinosaurs began to disperse and thus offer an opportunity to test the timing and drivers of this biogeographic pattern. Increased sampling can test this prediction: if dinosaurs initially dispersed under palaeolatitudinal-driven endemism, then an assemblage similar to those of South America^4,19-21 and India^19,22-including the earliest dinosaurs-should be present in Carnian deposits in south-central Africa. Here we report a new Carnian assemblage from Zimbabwe that includes Africa's oldest definitive dinosaurs, including a nearly complete skeleton of the sauropodomorph Mbiresaurus raathi gen. et sp. nov. This assemblage resembles other dinosaur-bearing Carnian assemblages, suggesting that a similar vertebrate fauna ranged high-latitude austral Pangaea. The distribution of the first dinosaurs is correlated with palaeolatitude-linked climatic barriers, and dinosaurian dispersal to the rest of the supercontinent was delayed until these barriers relaxed, suggesting that climatic controls influenced the initial composition of the terrestrial faunas that persist to this day.

Using information theory to detect rogue taxa and improve consensus trees

“Rogue” taxa of uncertain affinity can confound attempts to summarize the results of phylogenetic analyses. Rogues reduce resolution and support values in consensus trees, potentially obscuring strong evidence for relationships between other taxa. Information theory provides a principled means of assessing the congruence between a set of trees and their consensus, allowing rogue taxa to be identified more effectively than when using ad hoc measures of tree quality. A basic implementation of this approach in R recovers reduced consensus trees that are better resolved, more accurate, and more informative than those generated by existing methods. [Consensus trees; information theory; phylogenetic software; Rogue taxa.]

New spinosaurids from the Wessex Formation (Early Cretaceous, UK) and the European origins of Spinosauridae

Spinosaurids are among the most distinctive and yet poorly-known of large-bodied theropod dinosaurs, a situation exacerbated by their mostly fragmentary fossil record and competing views regarding their palaeobiology. Here, we report two new Early Cretaceous spinosaurid specimens from the Wessex Formation (Barremian) of the Isle of Wight. Large-scale phylogenetic analyses using parsimony and Bayesian techniques recover the pair in a new clade within Baryonychinae that also includes the hypodigm of the African spinosaurid Suchomimus. Both specimens represent distinct and novel taxa, herein named Ceratosuchops inferodios gen. et sp. nov. and Riparovenator milnerae gen. et sp. nov. A palaeogeographic reconstruction suggests a European origin for Spinosauridae, with at least two dispersal events into Africa. These new finds provide welcome information on poorly sampled areas of spinosaurid anatomy, suggest that sympatry was present and potentially common in baryonychines and spinosaurids as a whole, and contribute to updated palaeobiogeographic reconstructions for the clade.

Morphological phylogeny and review of the generic classification of Colletinae (Hymenoptera: Colletidae)

The bee subfamily Colletinae includes 542 species, the vast majority of which (518 spp.) belong to Colletes. The generic placement of the remaining 24 species has been controversial, resulting in several classifications being proposed. Despite several recently published molecular phylogenies of Colletinae, it remains unknown (1) what morphological synapomorphies support the recognized genera, (2) in which direction some relevant functional traits (e.g. basitibial plate) have evolved and (3) whether morphology supports the available molecular data. Herein, we provide a morphological phylogeny of Colletinae, which was constructed through parsimony analyses of 186 characters. In total, 50 ingroup species were included representing all major lineages of Colletes (29 spp.), plus all but three of the non-Colletes species of Colletinae (21 spp.). Trees were estimated through equal weights and extended implied weighting. Both provide strong support for the monophyly of Colletinae and indicate that the subfamily is defined by four unique synapomorphies. Our results also confirm recent phylogenetic hypotheses showing that Colletinae can be subdivided into two major clades: one comprising the reciprocally monophyletic Mourecotelles and Xanthocotelles; the other includes Colletes plus Hemicotelles, which are also reciprocally monophyletic. We also provide a fully illustrated key to facilitate generic identification of the Colletinae.

A phylogenetic C interpreter for TNT

MOTIVATION

TNT (a widely used program for phylogenetic analysis) includes an interpreter for a scripting language, but that implementation is nonstandard and uses several conventions of its own. This article describes the implementation and basic usage of a C interpreter (with all the ISO essentials) now included in TNT. A phylogenetic library includes functions that can be used for manipulating trees and data, as well as other phylogeny-specific tasks. This greatly extends the capabilities of TNT.

AVAILABILITY AND IMPLEMENTATION

Versions of TNT including the C interpreter for scripts can be downloaded from http://www.lillo.org.ar/phylogeny/tnt/.

The phylogenetic relationships of neosuchian crocodiles and their implications for the convergent evolution of the longirostrine condition

Since their origin in the Late Triassic, crocodylomorphs have had a long history of evolutionary change. Numerous studies examined their phylogeny, but none have attempted to unify their morphological characters into a single, combined dataset. Following a comprehensive review of published character sets, we present a new dataset for the crocodylomorph clade Neosuchia consisting of 569 morphological characters for 112 taxa. For the first time in crocodylian phylogenetic studies, quantitative variation was treated as continuous data (82 characters). To provide the best estimate of neosuchian relationships, and to investigate the origins of longirostry, these data were analysed using a variety of approaches. Our results show that equally weighted parsimony and Bayesian methods cluster unrelated longirostrine forms together, producing a topology that conflicts strongly with their stratigraphic distributions. By contrast, applying extended implied weighting improves stratigraphic congruence and removes longirostrine clustering. The resulting topologies resolve the major neosuchian clades, confirming several recent hypotheses regarding the phylogenetic placements of particular species (e.g. Baryphracta deponiae as a member of Diplocynodontinae) and groups (e.g. Tethysuchia as non-eusuchian neosuchians). The longirostrine condition arose at least three times independently by modification of the maxilla and premaxilla, accompanied by skull roof changes unique to each longirostrine clade.

Taxonomic and phylogenetic review of Necrosuchus ionensis (Alligatoroidea: Caimaninae) and the early evolution and radiation of caimanines

Alligatoroidea is the most species-rich crocodylomorph clade of the Cenozoic of South America, with nearly all species belonging to the Caimaninae clade. However, the earliest records of Caimaninae in South America, which are from the Palaeocene, are based mostly on incomplete specimens, which increases the importance of detailed taxonomic and phylogenetic studies on these taxa. This paper offers a taxonomic and phylogenetic review of Necrosuchus ionensis, a caimanine species from the Salamanca Formation of the Palaeocene of Argentina. Necrosuchus ionensis is considered a valid species, albeit with a different diagnosis from that proposed by previous authors. The phylogenetic analysis shows, for the first time, that N. ionensis belongs to the derived Caimaninae clade Jacarea. However, a better understanding of the Jacarea clade is needed, and alternative placements for N. ionensis might be considered. Nevertheless, the placement of N. ionensis as a derived caimanine raises interesting perspectives on the early evolution and radiation of caimanines, which are thoroughly discussed in this paper together with other results obtained in this study, such as the recovery of the North American caimanines Bottosaurus and Tsoabichi as a clade.

Bayesian and parsimony approaches reconstruct informative trees from simulated morphological datasets

Phylogenetic analysis aims to establish the true relationships between taxa. Different analytical methods, however, can reach different conclusions. In order to establish which approach best reconstructs true relationships, previous studies have simulated datasets from known tree topologies, and identified the method that reconstructs the generative tree most accurately. On this basis, researchers have argued that morphological datasets should be analysed by Bayesian approaches, which employ an explicit probabilistic model of evolution, rather than parsimony methods-with implied weights parsimony sometimes identified as particularly inaccurate. Accuracy alone, however, is an inadequate measure of a tree's utility: a fully unresolved tree is perfectly accurate, yet contains no phylogenetic information. The highly resolved trees recovered by implied weights parsimony in fact contain as much useful information as the more accurate, but less resolved, trees recovered by Bayesian methods. By collapsing poorly supported groups, this superior resolution can be traded for accuracy, resulting in trees as accurate as those obtained by a Bayesian approach. By contrast, equally weighted parsimony analysis produces trees that are less resolved and less accurate, leading to less reliable evolutionary conclusions.

TreeShrink: fast and accurate detection of outlier long branches in collections of phylogenetic trees

BACKGROUND

Sequence data used in reconstructing phylogenetic trees may include various sources of error. Typically errors are detected at the sequence level, but when missed, the erroneous sequences often appear as unexpectedly long branches in the inferred phylogeny.

RESULTS

We propose an automatic method to detect such errors. We build a phylogeny including all the data then detect sequences that artificially inflate the tree diameter. We formulate an optimization problem, called the k-shrink problem, that seeks to find k leaves that could be removed to maximally reduce the tree diameter. We present an algorithm to find the exact solution for this problem in polynomial time. We then use several statistical tests to find outlier species that have an unexpectedly high impact on the tree diameter. These tests can use a single tree or a set of related gene trees and can also adjust to species-specific patterns of branch length. The resulting method is called TreeShrink. We test our method on six phylogenomic biological datasets and an HIV dataset and show that the method successfully detects and removes long branches. TreeShrink removes sequences more conservatively than rogue taxon removal and often reduces gene tree discordance more than rogue taxon removal once the amount of filtering is controlled.

CONCLUSIONS

TreeShrink is an effective method for detecting sequences that lead to unrealistically long branch lengths in phylogenetic trees. The tool is publicly available at https://github.com/uym2/TreeShrink .

Measuring inferential importance of taxa using taxon influence indices

Assessing the importance of different taxa for inferring evolutionary history is a critical, but underutilized, aspect of systematics. Quantifying the importance of all taxa within a dataset provides an empirical measurement that can establish a ranking of extant taxa for ecological study and/or quantify the relative importance of newly announced or redescribed specimens to enable the disentangling of novelty and inferential influence. Here, we illustrate the use of taxon influence indices through analysis of both molecular and morphological datasets, introducing a modified Bayesian approach to the taxon influence index that accounts for model and topological uncertainty. Quantification of taxon influence using the Bayesian approach produced clear rankings for both dataset types. Bayesian taxon rankings differed from maximum likelihood (ML)‐derived rankings from a mitogenomic dataset, and the highest ranking taxa exhibited the largest interquartile range in influence estimate, suggesting variance in the estimate must be taken into account when the ranking of taxa is the feature of interest. Application of the Bayesian taxon influence index to a recent morphological analysis of the Tully Monster (Tullimonstrum) reveals that it exhibits consistently low inferential importance across two recent treatments of the taxon with alternative character codings. These results lend support to the idea that taxon influence indices may be robust to character coding and therefore effective for morphological analyses. These results underscore a need for the development of approaches to, and application of, taxon influence analyses both for the purpose of establishing robust rankings for future inquiry and for explicitly quantifying the importance of individual taxa. Quantifying the importance of individual taxa refocuses debates in morphological studies from questions of character choice/significance and taxon sampling to explicitly analytical techniques, and guides discussion of the context of new discoveries.

Inferring 'weak spots' in phylogenetic trees: application to mosasauroid nomenclature

Mosasauroid squamates represented the apex predators within the Late Cretaceous marine and occasionally also freshwater ecosystems. Proper understanding of the origin of their ecological adaptations or paleobiogeographic dispersals requires adequate knowledge of their phylogeny. The studies assessing the position of mosasauroids on the squamate evolutionary tree and their origins have long given conflicting results. The phylogenetic relationships within Mosasauroidea, however, have experienced only little changes throughout the last decades. Considering the substantial improvements in the development of phylogenetic methodology that have undergone in recent years, resulting, among others, in numerous alterations in the phylogenetic hypotheses of other fossil amniotes, we test the robustness in our understanding of mosasauroid beginnings and their evolutionary history. We re-examined a data set that results from modifications assembled in the course of the last 20 years and performed multiple parsimony analyses and Bayesian tip-dating analysis. Following the inferred topologies and the 'weak spots' in the phylogeny of mosasauroids, we revise the nomenclature of the 'traditionally' recognized mosasauroid clades, to acknowledge the overall weakness among branches and the alternative topologies suggested previously, and discuss several factors that might have an impact on the differing phylogenetic hypotheses and their statistical support.

Phylogenetic evaluation of Amynthas earthworms from South China reveals the initial ancestral state of spermathecae

Our knowledge of the phylogeny of the earthworm genus Amynthas under the family Megascolecidae, which is comprised of a huge number of species, is very limited compared to the better-known and much smaller family Lumbricidae. In order to investigate the phylogenetic relationships among the species within the genus Amynthas, which is the largest genus of the Megascolecidae family, nuclear and mitochondrial DNA sequences of 77 species, including 76 in-group Amynthas species collected from South China and 1 out-group species, were analyzed. A 5402bp segments composed of whole nuclear 18S rDNA and the mitochondrial genes COI, COII, ND1, 12S, and 16S was assembled from 77 species. Maximum Likelihood and Bayesian analyses of the concatenated sequences were performed. The results revealed evolution of two geographically independent lineages, both showing the ancestral state of two pairs of spermatheca (Sp.p 7/8/9). We found the species groups described by Sims and Easton (1972) to be non-monophyletic, and the origin of the parthenogenetic species group to likely be a quadthecal ancestor. These results provide modest evidence in support of an Indochinese peninsula origin of the Chinese Amynthas species and divergence of the genus once it had spread to mainland China. The findings of this study are consistent with a divergence scenario that resulted in at least one branch spreading to the Southeast of China and another branch spreading to the Southwest of China, but further research is required to confirm this interpretation of the Amynthas phylogeny.

Plasticity and Convergence in the Evolution of Short-Necked Plesiosaurs

Plesiosaurs were the longest-surviving group of secondarily marine tetrapods, comparable in diversity to today's cetaceans. During their long evolutionary history, which spanned the Jurassic and the Cretaceous (201 to 66 Ma), plesiosaurs repeatedly evolved long- and short-necked body plans [1, 2]. Despite this postcranial plasticity, short-necked plesiosaur clades have traditionally been regarded as being highly constrained to persistent and clearly distinct ecological niches: advanced members of Pliosauridae (ranging from the Middle Jurassic to the early Late Cretaceous) have been characterized as apex predators [2-5], whereas members of the distantly related clade Polycotylidae (middle to Late Cretaceous) were thought to have been fast-swimming piscivores [1, 5-7]. We report a new, highly unusual pliosaurid from the Early Cretaceous of Russia that shows close convergence with the cranial structure of polycotylids: Luskhan itilensis gen. et sp. nov. Using novel cladistic and ecomorphological data, we show that pliosaurids iteratively evolved polycotylid-like cranial morphologies from the Early Jurassic until the Early Cretaceous. This underscores the ecological diversity of derived pliosaurids and reveals a more complex evolutionary history than their iconic representation as gigantic apex predators of Mesozoic marine ecosystems suggests. Collectively, these data demonstrate an even higher degree of morphological plasticity and convergence in the evolution of plesiosaurs than previously thought and suggest the existence of an optimal ecomorphology for short-necked piscivorous plesiosaurs through time and across phylogeny.

Exploring Phylogenetic Relationships within Myriapoda and the Effects of Matrix Composition and Occupancy on Phylogenomic Reconstruction

Myriapods, including the diverse and familiar centipedes and millipedes, are one of the dominant terrestrial arthropod groups. Although molecular evidence has shown that Myriapoda is monophyletic, its internal phylogeny remains contentious and understudied, especially when compared to those of Chelicerata and Hexapoda. Until now, efforts have focused on taxon sampling (e.g., by including a handful of genes from many species) or on maximizing matrix size (e.g., by including hundreds or thousands of genes in just a few species), but a phylogeny maximizing sampling at both levels remains elusive. In this study, we analyzed 40 Illumina transcriptomes representing 3 of the 4 myriapod classes (Diplopoda, Chilopoda, and Symphyla); 25 transcriptomes were newly sequenced to maximize representation at the ordinal level in Diplopoda and at the family level in Chilopoda. Ten supermatrices were constructed to explore the effect of several potential phylogenetic biases (e.g., rate of evolution, heterotachy) at 3 levels of gene occupancy per taxon (50%, 75%, and 90%). Analyses based on maximum likelihood and Bayesian mixture models retrieved monophyly of each myriapod class, and resulted in 2 alternative phylogenetic positions for Symphyla, as sister group to Diplopoda + Chilopoda, or closer to Diplopoda, the latter hypothesis having been traditionally supported by morphology. Within centipedes, all orders were well supported, but 2 deep nodes remained in conflict in the different analyses despite dense taxon sampling at the family level. Relationships among centipede orders in all analyses conducted with the most complete matrix (90% occupancy) are at odds not only with the sparser but more gene-rich supermatrices (75% and 50% supermatrices) and with the matrices optimizing phylogenetic informativeness or most conserved genes, but also with previous hypotheses based on morphology, development, or other molecular data sets. Our results indicate that a high percentage of ribosomal proteins in the most complete matrices, in conjunction with distance from the root, can act in concert to compromise the estimated relationships within the ingroup. We discuss the implications of these findings in the context of the ever more prevalent quest for completeness in phylogenomic studies.

Split-specific bootstrap measures for quantifying phylogenetic stability and the influence of taxon selection

Assessing the robustness of an inferred phylogeny is an important element of phylogenetics. This is typically done with measures of stabilities at the internal branches and the variation of the positions of the leaf nodes. The bootstrap support for branches in maximum parsimony, distance and maximum likelihood estimation, or posterior probabilities in Bayesian inference, measure the uncertainty about a branch due to the sampling of the sites from genes or sampling genes from genomes. However, these measures do not reveal how taxon sampling affects branch support and the effects of taxon sampling on the estimated phylogeny. An internal branch in a phylogenetic tree can be viewed as a split that separates the taxa into two nonempty complementary subsets. We develop several split-specific measures of stability determined from bootstrap support for quartets. These include BPtaxon_split (average bootstrap percentage [BP] for all quartets involving a taxon within a split), BPsplit (BPtaxon_split averaged over taxa), BPtaxon (BPtaxon_split averaged over splits) and RBIC-taxon (average BP over all splits after removing a taxon). We also develop a pruned-tree distance metric. Application of our measures to empirical and simulated data illustrate that existing measures of overall stability can fail to detect taxa that are the primary source of a split-specific instability. Moreover, we show that the use of many reduced sets of quartets is important in being able to detect the influence of joint sets of taxa rather than individual taxa. These new measures are valuable diagnostic tools to guide taxon sampling in phylogenetic experimental design.

Human population history revealed by a supertree approach

Over the past two decades numerous new trees of modern human populations have been published extensively but little attention has been paid to formal phylogenetic synthesis. We utilized the "matrix representation with parsimony" (MRP) method to infer a composite phylogeny (supertree) of modern human populations, based on 257 genetic/genomic, as well as linguistic, phylogenetic trees and 44 admixture plots from 200 published studies (1990-2014). The resulting supertree topology includes the most basal position of S African Khoisan followed by C African Pygmies, and the paraphyletic section of all other sub-Saharan peoples. The sub-Saharan African section is basal to the monophyletic clade consisting of the N African-W Eurasian assemblage and the consistently monophyletic Eastern superclade (Sahul-Oceanian, E Asian, and Beringian-American peoples). This topology, dominated by genetic data, is well-resolved and robust to parameter set changes, with a few unstable areas (e.g., West Eurasia, Sahul-Melanesia) reflecting the existing phylogenetic controversies. A few populations were identified as highly unstable "wildcard taxa" (e.g. Andamanese, Malagasy). The linguistic classification fits rather poorly on the supertree topology, supporting a view that direct coevolution between genes and languages is far from universal.

TNT version 1.5, including a full implementation of phylogenetic morphometrics

Version 1.5 of the computer program TNT completely integrates landmark data into phylogenetic analysis. Landmark data consist of coordinates (in two or three dimensions) for the terminal taxa; TNT reconstructs shapes for the internal nodes such that the difference between ancestor and descendant shapes for all tree branches sums up to a minimum; this sum is used as tree score. Landmark data can be analysed alone or in combination with standard characters; all the applicable commands and options in TNT can be used transparently after reading a landmark data set. The program continues implementing all the types of analyses in former versions, including discrete and continuous characters (which can now be read at any scale, and automatically rescaled by TNT). Using algorithms described in this paper, searches for landmark data can be made tens to hundreds of times faster than it was possible before (from T to 3T times faster, where T is the number of taxa), thus making phylogenetic analysis of landmarks feasible even on standard personal computers.

Postcranial anatomy of Pissarrachampsa sera (Crocodyliformes, Baurusuchidae) from the Late Cretaceous of Brazil: insights on lifestyle and phylogenetic significance

The postcranial anatomy of Crocodyliformes has historically been neglected, as most descriptions are based solely on skulls. Yet, the significance of the postcranium in crocodyliforms evolution is reflected in the great lifestyle diversity exhibited by the group, with members ranging from terrestrial animals to semi-aquatic and fully marine forms. Recently, studies have emphasized the importance of the postcranium. Following this trend, here we present a detailed description of the postcranial elements of Pissarrachampsa sera (Mesoeucrocodylia, Baurusuchidae), from the Adamantina Formation (Bauru Group, Late Cretaceous of Brazil). The preserved elements include dorsal vertebrae, partial forelimb, pelvic girdle, and hindlimbs. Comparisons with the postcranial anatomy of baurusuchids and other crocodyliforms, together with body-size and mass estimates, lead to a better understanding of the paleobiology of Pissarrachampsa sera, including its terrestrial lifestyle and its role as a top predator. Furthermore, the complete absence of osteoderms in P. sera, a condition previously known only in marine crocodyliforms, suggests osteoderms very likely played a minor role in locomotion of baurusuchids, unlike other groups of terrestrial crocodyliforms. Finally, a phylogenetic analysis including the newly recognized postcranial features was carried out, and exploratory analyses were performed to investigate the influence of both cranial and postcranial characters in the phylogeny of Crocodyliformes. Our results suggest that crocodyliform relationships are mainly determined by cranial characters. However, this seems to be a consequence of the great number of missing entries in the data set with only postcranial characters and not of the lack of potential (or synapomorphies) for this kind of data to reflect the evolutionary history of Crocodyliformes.

Avoiding Missing Data Biases in Phylogenomic Inference: An Empirical Study in the Landfowl (Aves: Galliformes)

Production of massive DNA sequence data sets is transforming phylogenetic inference, but best practices for analyzing such data sets are not well established. One uncertainty is robustness to missing data, particularly in coalescent frameworks. To understand the effects of increasing matrix size and loci at the cost of increasing missing data, we produced a 90 taxon, 2.2 megabase, 4,800 locus sequence matrix of landfowl using target capture of ultraconserved elements. We then compared phylogenies estimated with concatenated maximum likelihood, quartet-based methods executed on concatenated matrices and gene tree reconciliation methods, across five thresholds of missing data. Results of maximum likelihood and quartet analyses were similar, well resolved, and demonstrated increasing support with increasing matrix size and sparseness. Conversely, gene tree reconciliation produced unexpected relationships when we included all informative loci, with certain taxa placed toward the root compared with other approaches. Inspection of these taxa identified a prevalence of short average contigs, which potentially biased gene tree inference and caused erroneous results in gene tree reconciliation. This suggests that the more problematic missing data in gene tree-based analyses are partial sequences rather than entire missing sequences from locus alignments. Limiting gene tree reconciliation to the most informative loci solved this problem, producing well-supported topologies congruent with concatenation and quartet methods. Collectively, our analyses provide a well-resolved phylogeny of landfowl, including strong support for previously problematic relationships such as those among junglefowl (Gallus), and clarify the position of two enigmatic galliform genera (Lerwa, Melanoperdix) not sampled in previous molecular phylogenetic studies.