Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences.

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Fossilization can mislead analyses of phenotypic disparity.

Analyses of morphological disparity can incorporate living and fossil taxa to facilitate the exploration of how phenotypic variation changes through time. However, taphonomic processes introduce non-random patterns of data loss in fossil data and their impact on perceptions of disparity is unclear. To address this, we characterize how measures of disparity change when simulated and empirical data are degraded through random and structured data loss. We demonstrate that both types of data loss can distort the disparity of clades, and that the magnitude and direction of these changes varies between the most commonly employed distance metrics and disparity indices. The inclusion of extant taxa and exceptionally preserved fossils mitigates these distortions and clarifies the full extent of the data lost, most of which would otherwise go uncharacterized. This facilitates the use of ancestral state estimation and evolutionary simulations to further control for the effects of data loss. Where the addition of such reference taxa is not possible, we urge caution in the extrapolation of general patterns in disparity from datasets that characterize subsets of phenotype, which may represent no more than the traits that they sample.

Phylogenetic congruence, conflict and consilience between molecular and morphological data.

Morphology and molecules are important data sources for estimating evolutionary relationships. Modern studies often utilise morphological and molecular partitions alongside each other in combined analyses. However, the effect of combining phenomic and genomic partitions is unclear. This is exacerbated by their size imbalance, and conflict over the efficacy of different inference methods when using morphological characters. To systematically address the effect of topological incongruence, size imbalance, and tree inference methods, we conduct a meta-analysis of 32 combined (molecular + morphology) datasets across metazoa. Our results reveal that morphological-molecular topological incongruence is pervasive: these data partitions yield very different trees, irrespective of which method is used for morphology inference. Analysis of the combined data often yields unique trees that are not sampled by either partition individually, even with the inclusion of relatively small quantities of morphological characters. Differences between morphology inference methods in terms of resolution and congruence largely relate to consensus methods. Furthermore, stepping stone Bayes factor analyses reveal that morphological and molecular partitions are not consistently combinable, i.e. data partitions are not always best explained under a single evolutionary process. In light of these results, we advise that the congruence between morphological and molecular data partitions needs to be considered in combined analyses. Nonetheless, our results reveal that, for most datasets, morphology and molecules can, and should, be combined in order to best estimate evolutionary history and reveal hidden support for novel relationships. Studies that analyse only phenomic or genomic data in isolation are unlikely to provide the full evolutionary picture.

Multivariate mapping of ontogeny, taphonomy and phylogeny to reconstruct problematic fossil taxa.

Exceptionally preserved fossils of soft-bodied organisms provide unique evidence of evolutionary history, but they are often contentious; different approaches frequently produce radically different reconstructions of taxa and their affinities. Conflict arises due to difficulties in disentangling the three non-independent factors that underlie all morphological variation within and between fossils: ontogeny, taphonomy and phylogeny. Comparative data from extant organisms can be extremely powerful in this context, but is often difficult to apply given the multi-dimensionality of anatomical variation. Here, we present a multivariate ordination method using discrete morphological character data from modern taxa at different ontogenetic and taphonomic stages (semaphoront and 'semataphonts'). Analysing multiple axes of morphological variation simultaneously allows us to visualize character combinations that are likely to exist in fossil specimens at intersecting stages of growth and decay, and thus constrain interpretation of fossils. Application to early vertebrates finds variation in fossil specimens to be accounted for by all three axes: primarily decay in Mayomyzon, ontogeny in Priscomyzon and phylogeny in 'euphaneropoids' and Palaeospondylus. Our demonstration of empirical multi-factorial variation underscores the power of multivariate approaches to fossil interpretation, especially non-biomineralized taxa. As such, this conceptual approach provides a new method for resolving enigmatic taxa throughout the fossil record.

A well-preserved 'placoderm' (stem-group Gnathostomata) upper jaw from the Early Devonian of Mongolia clarifies jaw evolution.

The origin of jaws and teeth remains contentious in vertebrate evolution. 'Placoderms' (Silurian-Devonian armoured jawed fishes) are central to debates on the origins of these anatomical structures. 'Acanthothoracids' are generally considered the most primitive 'placoderms'. However, they are so far known mainly from disarticulated skeletal elements that are typically incomplete. The structure of the jaws-particularly the jaw hinge-is poorly known, leaving open questions about their jaw function and comparison with other placoderms and modern gnathostomes. Here we describe a near-complete 'acanthothoracid' upper jaw, allowing us to reconstruct the likely orientation and angle of the bite and compare its morphology with that of other known 'placoderm' groups. We clarify that the bite position is located on the upper jaw cartilage rather than on the dermal cheek and thus show that there is a highly conserved bite morphology among most groups of 'placoderms', regardless of their overall cranial geometry. Incorporation of the dermal skeleton appears to provide a sound biomechanical basis for jaw origins. It appears that 'acanthothoracid' dentitions were fundamentally similar in location to that of arthrodire 'placoderms', rather than resembling bony fishes. Irrespective of current phylogenetic uncertainty, the new data here resolve the likely general condition for 'placoderms' as a whole, and as such, ancestral morphology of known jawed vertebrates.

Endochondral bone in an Early Devonian 'placoderm' from Mongolia.

Endochondral bone is the main internal skeletal tissue of nearly all osteichthyans-the group comprising more than 60,000 living species of bony fishes and tetrapods. Chondrichthyans (sharks and their kin) are the living sister group of osteichthyans and have primarily cartilaginous endoskeletons, long considered the ancestral condition for all jawed vertebrates (gnathostomes). The absence of bone in modern jawless fishes and the absence of endochondral ossification in early fossil gnathostomes appear to lend support to this conclusion. Here we report the discovery of extensive endochondral bone in Minjinia turgenensis, a new genus and species of 'placoderm'-like fish from the Early Devonian (Pragian) of western Mongolia described using X-ray computed microtomography. The fossil consists of a partial skull roof and braincase with anatomical details providing strong evidence of placement in the gnathostome stem group. However, its endochondral space is filled with an extensive network of fine trabeculae resembling the endochondral bone of osteichthyans. Phylogenetic analyses place this new taxon as a proximate sister group of the gnathostome crown. These results provide direct support for theories of generalized bone loss in chondrichthyans. Furthermore, they revive theories of a phylogenetically deeper origin of endochondral bone and its absence in chondrichthyans as a secondary condition.

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

Bite marks and predation of fossil jawless fish during the rise of jawed vertebrates.

Although modern vertebrate diversity is dominated by jawed vertebrates, early vertebrate assemblages were predominantly composed of jawless fishes. Hypotheses for this faunal shift and the Devonian decline of jawless vertebrates include predation and competitive replacement. The nature and prevalence of ecological interactions between jawed and jawless vertebrates are highly relevant to both hypotheses, but direct evidence is limited. Here, we use the occurrence and distribution of bite mark type traces in fossil jawless armoured heterostracans to infer predation interactions. A total of 41 predated specimens are recorded; their prevalence increases through time, reaching a maximum towards the end of the Devonian. The bite mark type traces significantly co-occur with jawed vertebrates, and their distribution through time is correlated with jawed vertebrate diversity patterns, particularly placoderms and sarcopterygians. Environmental and ecological turnover in the Devonian, especially relating to the nekton revolution, have been inferred as causes of the faunal shift from jawless to jawed vertebrates. Here, we provide direct evidence of escalating predation from jawed vertebrates as a potential contributing factor to the demise and extinction of ostracoderms.

Variation in appendages in early Cambrian bradoriids reveals a wide range of body plans in stem-euarthropods.

Traditionally, the origin and evolution of modern arthropod body plans has been revealed through increasing levels of appendage specialisation exhibited by Cambrian euarthropods. Here we show significant variation in limb morphologies and patterns of limb-tagmosis among three early Cambrian arthropod species conventionally assigned to the Bradoriida. These arthropods are recovered as a monophyletic stem-euarthropod group (and sister taxon to crown-group euarthropods, i.e. Chelicerata, Mandibulata and their extinct relatives), thus implying a radiation of stem-euarthropods where trends towards increasing appendage specialisation were explored convergently with other euarthropod groups. The alternative solution, where bradoriids are polyphyletic, representing several independent origins of a small, bivalved body plan in lineages from diverse regions of the euarthropod and mandibulate stems, is only marginally less parsimonious. The new data reveal a previously unknown disparity of body plans in stem-euarthropods and both solutions support remarkable evolutionary convergence, either of fundamental body plans or appendage specialization patterns.

The nearshore cradle of early vertebrate diversification.

Ancestral vertebrate habitats are subject to controversy and obscured by limited, often contradictory paleontological data. We assembled fossil vertebrate occurrence and habitat datasets spanning the middle Paleozoic (480 million to 360 million years ago) and found that early vertebrate clades, both jawed and jawless, originated in restricted, shallow intertidal-subtidal environments. Nearshore divergences gave rise to body plans with different dispersal abilities: Robust fishes shifted shoreward, whereas gracile groups moved seaward. Fresh waters were invaded repeatedly, but movement to deeper waters was contingent upon form and short-lived until the later Devonian. Our results contrast with the onshore-offshore trends, reef-centered diversification, and mid-shelf clustering observed for benthic invertebrates. Nearshore origins for vertebrates may be linked to the demands of their mobility and may have influenced the structure of their early fossil record and diversification.

Morphospace saturation in the stem-gnathostomes pteraspidiformes heterostracans: an early radiation of a 'bottom' heavy clade.

Ostracoderms (fossil armoured jawless fishes) shed light on early vertebrate evolution by revealing the step-wise acquisition of jawed vertebrate characters, and were important constituents of Middle Palaeozoic vertebrate faunas. A wide variety of head shield shapes are observed within and between the ostracoderm groups, but the timing of these diversifications and the consistency between different measures of their morphospace are unclear. Here, we present the first disparity (explored morphospace) versus diversity (number of taxa) analysis of Pteraspidiformes heterostracans using continuous and discrete characters. Patterns of taxic diversity and morphological disparity are in accordance: they both show a rise to a peak in the Lochkovian followed by a gradual decline in the Middle-Late Devonian. Patterns are largely consistent for disparity measures using sum of ranges or total variance, and when using continuous or discrete characters. Pteraspidiformes heterostracans can be classified as a "bottom-heavy clade", i.e., a group where a high initial disparity decreasing over time is detected. In fact, the group explored morphospace early in its evolutionary history, with much of the subsequent variation in dermal armour occurring as variation in the proportions of already evolved anatomical features. This Early Devonian radiation is also in agreement with the paleobiogeographic distribution of the group, with a maximum of dispersal and explored morphospace during the Lochkovian and Pragian time bins.

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

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

Differences between hard and soft phylogenetic data.

When building the tree of life, variability of phylogenetic signal is often accounted for by partitioning gene sequences and testing for differences. The same considerations, however, are rarely applied to morphological data, potentially undermining its use in evolutionary contexts. Here, we apply partition heterogeneity tests to 59 animal datasets to demonstrate that significant differences exist between the phylogenetic signal conveyed by 'hard' and 'soft' characters (bones, teeth and shells versus myology, integument etc). Furthermore, the morphological partitions differ significantly in their consistency relative to independent molecular trees. The observed morphological differences correspond with missing data biases, and as such their existence presents a problem not only for phylogeny reconstruction, but also for interpretations of fossil data. Evolutionary inferences drawn from clades in which hard, readily fossilizable characters are relatively less consistent and different from other morphology (mammals, bivalves) may be less secure. More secure inferences might be drawn from the fossil record of clades that exhibit fewer differences, or exhibit more consistent hard characters (fishes, birds). In all cases, it will be necessary to consider the impact of missing data on empirical data, and the differences that exist between morphological modules.

Dental Data Perform Relatively Poorly in Reconstructing Mammal Phylogenies: Morphological Partitions Evaluated with Molecular Benchmarks.

Phylogenetic trees underpin reconstructions of evolutionary history and tests of evolutionary hypotheses. They are inferred from both molecular and morphological data, yet the relative value of morphology has been questioned in this context due to perceived homoplasy, developmental linkage, and nonindependence of characters. Nevertheless, fossil data are limited to incomplete subsets of preserved morphology, and different regions are treated as equivalent. Through meta-analysis of 40 data sets, we show here that the dental and osteological characters of mammals convey significantly different phylogenetic signals, and that osteological characters are significantly more compatible with molecular trees. Furthermore, the application of simplified paleontological filters (retaining only dental data) results in significantly greater loss of phylogenetic signal than random character ablation. Although the mammal fossil record is largely comprised of teeth, dental data alone are generally found to be less reliable for phylogenetic reconstruction given their incongruence with osteological and molecular data. These findings highlight the need for rigorous meta-analyses of distributions of homoplasy in morphological data. These tests, and consequent refinements to phylogenetic analyses that they permit, promise to improve the quality of all macroevolutionary studies that hinge on accurate trees. [Homoplasy; Mammalia; morphology; osteology; phylogeny; teeth.

Preservation and phylogeny of Cambrian ecdysozoans tested by experimental decay of Priapulus.

The exceptionally preserved Cambrian fossil record provides unique insight into the early evolutionary history of animals. Understanding of the mechanisms of exceptional soft tissue preservation frames all interpretations of the fauna and its evolutionary significance. This is especially true for recent interpretations of preserved nervous tissues in fossil ecdysozoans. However, models of soft tissue preservation lack empirical support from actualistic studies. Here experimental decay of the priapulid Priapulus reveal consistent bias towards rapid loss of internal non-cuticular anatomy compared with recalcitrant cuticular anatomy. This is consistent with models of Burgess Shale-type preservation and indicates that internal tissues are unlikely to be preserved with fidelity if organically preserved. This pattern, along with extreme body margin distortion, is consistent with onychophoran decay, and is therefore resolved as general for early ecdysozoans. Application of these patterns to phylogenetic data finds scalidophoran taxa to be very sensitive to taphonomically informed character coding, but not panarthropodan taxa. Priapulid decay also have unexpected relevance for interpretation of myomeres in fossil chordates. The decay data presented serve not only as a test of models of preservation but also a framework with which to interpret ecdysozoan fossil anatomies, and the subsequent evolutionary inferences drawn from them.

Pigmented anatomy in Carboniferous cyclostomes and the evolution of the vertebrate eye.

The success of vertebrates is linked to the evolution of a camera-style eye and sophisticated visual system. In the absence of useful data from fossils, scenarios for evolutionary assembly of the vertebrate eye have been based necessarily on evidence from development, molecular genetics and comparative anatomy in living vertebrates. Unfortunately, steps in the transition from a light-sensitive 'eye spot' in invertebrate chordates to an image-forming camera-style eye in jawed vertebrates are constrained only by hagfish and lampreys (cyclostomes), which are interpreted to reflect either an intermediate or degenerate condition. Here, we report-based on evidence of size, shape, preservation mode and localized occurrence-the presence of melanosomes (pigment-bearing organelles) in fossil cyclostome eyes. Time of flight secondary ion mass spectrometry analyses reveal secondary ions with a relative intensity characteristic of melanin as revealed through principal components analyses. Our data support the hypotheses that extant hagfish eyes are degenerate, not rudimentary, that cyclostomes are monophyletic, and that the ancestral vertebrate had a functional visual system. We also demonstrate integument pigmentation in fossil lampreys, opening up the exciting possibility of investigating colour patterning in Palaeozoic vertebrates. The examples we report add to the record of melanosome preservation in Carboniferous fossils and attest to surprising durability of melanosomes and biomolecular melanin.

Sampling diverse characters improves phylogenies: Craniodental and postcranial characters of vertebrates often imply different trees.

Morphological cladograms of vertebrates are often inferred from greater numbers of characters describing the skull and teeth than from postcranial characters. This is either because the skull is believed to yield characters with a stronger phylogenetic signal (i.e., contain less homoplasy), because morphological variation therein is more readily atomized, or because craniodental material is more widely available (particularly in the palaeontological case). An analysis of 85 vertebrate datasets published between 2000 and 2013 confirms that craniodental characters are significantly more numerous than postcranial characters, but finds no evidence that levels of homoplasy differ in the two partitions. However, a new partition test, based on tree-to-tree distances (as measured by the Robinson Foulds metric) rather than tree length, reveals that relationships inferred from the partitions are significantly different about one time in three, much more often than expected. Such differences may reflect divergent selective pressures in different body regions, resulting in different localized patterns of homoplasy. Most systematists attempt to sample characters broadly across body regions, but this is not always possible. We conclude that trees inferred largely from either craniodental or postcranial characters in isolation may differ significantly from those that would result from a more holistic approach. We urge the latter.

Bias and sensitivity in the placement of fossil taxa resulting from interpretations of missing data.

The utility of fossils in evolutionary contexts is dependent on their accurate placement in phylogenetic frameworks, yet intrinsic and widespread missing data make this problematic. The complex taphonomic processes occurring during fossilization can make it difficult to distinguish absence from non-preservation, especially in the case of exceptionally preserved soft-tissue fossils: is a particular morphological character (e.g., appendage, tentacle, or nerve) missing from a fossil because it was never there (phylogenetic absence), or just happened to not be preserved (taphonomic loss)? Missing data have not been tested in the context of interpretation of non-present anatomy nor in the context of directional shifts and biases in affinity. Here, complete taxa, both simulated and empirical, are subjected to data loss through the replacement of present entries (1s) with either missing (?s) or absent (0s) entries. Both cause taxa to drift down trees, from their original position, toward the root. Absolute thresholds at which downshift is significant are extremely low for introduced absences (two entries replaced, 6% of present characters). The opposite threshold in empirical fossil taxa is also found to be low; two absent entries replaced with presences causes fossil taxa to drift up trees. As such, only a few instances of non-preserved characters interpreted as absences will cause fossil organisms to be erroneously interpreted as more primitive than they were in life. This observed sensitivity to coding non-present morphology presents a problem for all evolutionary studies that attempt to use fossils to reconstruct rates of evolution or unlock sequences of morphological change. Stem-ward slippage, whereby fossilization processes cause organisms to appear artificially primitive, appears to be a ubiquitous and problematic phenomenon inherent to missing data, even when no decay biases exist. Absent characters therefore require explicit justification and taphonomic frameworks to support their interpretation.

Discriminating signal from noise in the fossil record of early vertebrates reveals cryptic evolutionary history.

The fossil record of early vertebrates has been influential in elucidating the evolutionary assembly of the gnathostome bodyplan. Understanding of the timing and tempo of vertebrate innovations remains, however, mired in a literal reading of the fossil record. Early jawless vertebrates (ostracoderms) exhibit restriction to shallow-water environments. The distribution of their stratigraphic occurrences therefore reflects not only flux in diversity, but also secular variation in facies representation of the rock record. Using stratigraphic, phylogenetic and palaeoenvironmental data, we assessed the veracity of the fossil records of the jawless relatives of jawed vertebrates (Osteostraci, Galeaspida, Thelodonti, Heterostraci). Non-random models of fossil recovery potential using Palaeozoic sea-level changes were used to calculate confidence intervals of clade origins. These intervals extend the timescale for possible origins into the Upper Ordovician; these estimates ameliorate the long ghost lineages inferred for Osteostraci, Galeaspida and Heterostraci, given their known stratigraphic occurrences and stem-gnathostome phylogeny. Diversity changes through the Silurian and Devonian were found to lie within the expected limits predicted from estimates of fossil record quality indicating that it is geological, rather than biological factors, that are responsible for shifts in diversity. Environmental restriction also appears to belie ostracoderm extinction and demise rather than competition with jawed vertebrates.

Fossilization causes organisms to appear erroneously primitive by distorting evolutionary trees.

Fossils are vital for calibrating rates of molecular and morphological change through geological time, and are the only direct source of data documenting macroevolutionary transitions. Many evolutionary studies therefore require the robust phylogenetic placement of extinct organisms. Here, we demonstrate that the inevitable bias of the fossil record to preserve just hard, skeletal morphology systemically distorts phylogeny. Removal of soft part characters from 78 modern vertebrate and invertebrate morphological datasets resulted in significant changes to phylogenetic signal; it caused individual taxa to drift from their original position, predominately downward toward the root of their respective trees. This last bias could systematically inflate evolutionary rates inferred from molecular data because first fossil occurrences will not be recognised as such. Stem-ward slippage, whereby fundamental taphonomic biases cause fossils to be interpreted as erroneously primitive, is therefore a ubiquitous problem for all biologists attempting to infer macroevolutionary rates or sequences.