The Ediacaran macrofossil Charnia masoni Ford is perhaps the most iconic member of the Rangeomorpha: a group of seemingly sessile, frondose organisms that dominates late Ediacaran benthic, deep‐marine fossil assemblages. Despite C. masoni exhibiting broad palaeogeographical and stratigraphical ranges, there have been few morphological studies that consider the variation observed among populations of specimens derived from multiple global localities. We present an analysis of C. masoni that evaluates specimens from the UK, Canada and Russia, representing the largest morphological study of this taxon to date. We describe substantial morphological variation within C. masoni and present a new morphological model for this species that has significant implications both for interpretation of rangeomorph architecture, and potentially for existing taxonomic schemes. Previous reconstructions of Charnia include assumptions regarding the presence of structures seen in other rangeomorphs (e.g. an internal stalk) and of homogeneity in higher order branch morphology; observations that are not borne out by our investigations. We describe variation in the morphology of third and fourth order branches, as well as variation in gross structure near the base of the frond. The diagnosis of Charnia masoni is emended to take account of these new features. These findings highlight the need for large‐scale analyses of rangeomorph morphology in order to better understand the biology of this long‐enigmatic group.

In order to understand patterns and processes of the diversification of life we require an accurate understanding of taxa interrelationships. Recent studies have suggested that analyses of morphological character data using the Bayesian and Maximum likelihood Mk model provide phylogenies of higher accuracy compared to parsimony methods. These studies have proved controversial, particularly simulating morphology-data under Markov models that assume shared branch lengths for characters, as it is claimed this leads to bias favouring the Bayesian or Maximum likelihood Mk model over parsimony models which do not explicitly make this assumption. We avoid these potential issues by employing a simulation protocol in which character states are randomly assigned to tips, but datasets are constrained to an empirically-realistic distribution of homoplasy as measured by the Consistency Index. Datasets were analysed with equal-weights and implied weights parsimony, and the Maximum Likelihood and Bayesian Mk model. We find that consistent (low homoplasy) datasets render method choice largely irrelevant, as all methods perform well with high consistency (low homoplasy) datasets, but the largest discrepancies in accuracy occur with low consistency datasets (high homoplasy). In such cases, the Bayesian Mk model is significantly more accurate than alternative models, and Implied weights parsimony never significantly out-performs the Bayesian Mk model. When poorly-supported branches are collapsed, the Bayesian Mk model recovers trees with higher resolution compared to other methods. Since it is not possible to assess homoplasy independently of a tree estimate, the Bayesian Mk model emerges as the most reliable method for categorical morphological analyses.

Timescales are of fundamental importance to evolutionary biology as they facilitate hypothesis tests of historical evolutionary processes. Through the incorporation of fossil occurrence data, the fossilised birth-death (FBD) process provides a framework for estimating divergence times using more palaeontological data than traditional node calibration approaches have allowed. The inclusion of more data can refine evolutionary timescale estimates, but for many taxonomic groups it is computationally infeasible to include all fossil occurrence data. Here, we utilise both empirical data and a simulation framework to identify approaches to subsampling fossil occurrence data that result in the most accurate estimates of divergence times. To achieve this we assess the performance of the FBD-Skyline model when implementing multiple approaches to incorporating subsampled fossil occurrences. Our results demonstrate that it is necessary to account for all available fossil occurrence data to achieve the most accurate estimates of clade age. We show that this can be achieved if an empirical Bayes approach to account for fossil sampling through time is applied to the FBD process. Random subsampling of occurrence data can lead to estimates of clade age that are incompatible with fossil evidence if no control over the affinities of fossil occurrences is enforced. Our results call into question the accuracy of previous divergence time studies incorporating the FBD process that have used only a subsample of all available fossil occurrence data.

The morphology of the vertebrate lower jaw has been used to infer feeding ecology.; transformations in mandibular shape and structure likely to have facilitated the emergence of different feeding behaviours in vertebrate evolution. Here we present elliptical Fourier shape and principal component analyses, characterizing and comparing the disparity of jaw shape in early gnathostomes and their modern primitively aquatic counterparts. 83% of shape variation is summarized on the first three principal component axes and all component clades of early gnathostomes exhibit overlapping morphological variation. Non-tetrapodomorph Palaeozoic sarcopterygians are more disparate than their extant counterparts whereas extant chondrichthyans are more disparate than their Palaeozoic counterparts. More generally, extant jawed fishes are more disparate than their Palaeozoic relatives largely because of the extensive shape variation exhibited by mandibles of extant actinopterygians. Only some areas of shape space vacated by Palaeozoic gnathostomes have been convergently refilled by living taxa. Characterization of theoretical jaw morphologies demonstrates that less than half of all possible shapes are realised by the jawed fishes that comprise our empirical dataset; many of these morphologies are realised by unrepresented terrestrial tetrapods, implying environmental constraint. Our results are incompatible with the early burst model of clade evolution and contradict the hypothesis that maximum disparity is reached early in the evolutionary history of jawed fishes.

Fossil taxa are critical to inferences of historical diversity and the origins of modern biodiversity, but realizing their evolutionary significance is contingent on restoring fossil species to their correct position within the tree of life. For most fossil species, morphology is the only source of data for phylogenetic inference; this has traditionally been analysed using parsimony, the predominance of which is currently challenged by the development of probabilistic models that achieve greater phylogenetic accuracy. Here, based on simulated and empirical datasets, we explore the relative efficacy of competing phylogenetic methods in terms of clade support. We characterize clade support using bootstrapping for parsimony and Maximum Likelihood, and intrinsic Bayesian posterior probabilities, collapsing branches that exhibit less than 50% support. Ignoring node support, Bayesian inference is the most accurate method in estimating the tree used to simulate the data. After assessing clade support, Bayesian and Maximum Likelihood exhibit comparable levels of accuracy, and parsimony remains the least accurate method. However, Maximum Likelihood is less precise than Bayesian phylogeny estimation, and Bayesian inference recaptures more correct nodes with higher support compared to all other methods, including Maximum Likelihood. We assess the effects of these findings on empirical phylogenies. Our results indicate probabilistic methods should be favoured over parsimony.

Consensus trees are required to summarise trees obtained through MCMC sampling of a posterior distribution, providing an overview of the distribution of estimated parameters such as topology, branch lengths and divergence times. Numerous consensus tree construction methods are available, each presenting a different interpretation of the tree sample. The rise of morphological clock and sampled-ancestor methods of divergence time estimation, in which times and topology are co-estimated, has increased the popularity of the maximum clade credibility (MCC) consensus tree method. The MCC method assumes that the sampled, fully resolved topology with the highest clade credibility contains an adequate summary of the most probable clades, with parameter estimates from compatible sampled trees used to obtain the marginal distributions of parameters such as clade ages and branch lengths. Using both simulated and empirical data, we demonstrate that MCC trees, and trees constructed using the similar maximum a posteriori (MAP) method, often include poorly supported and incorrect clades when summarising diffuse posterior samples of trees. We demonstrate that the paucity of information in morphological datasets contributes to the inability of MCC and MAP trees to present an accurate summary of the posterior distribution. Conversely, majority-rule consensus (MRC) trees report a lower proportion of incorrect nodes when summarising the same posterior samples of trees. Thus, we advocate the use of MRC trees, in place of MCC or MAP trees, in attempts to summarise the results of Bayesian phylogenetic analyses of morphological data.

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.

Molecular sequence data provide information about relative times only, and fossil-based age constraints are the ultimate source of information about absolute times in molecular clock dating analyses. Thus, fossil calibrations are critical to molecular clock dating, but competing methods are difficult to evaluate empirically because the true evolutionary time scale is never known. Here, we combine mechanistic models of fossil preservation and sequence evolution in simulations to evaluate different approaches to constructing fossil calibrations and their impact on Bayesian molecular clock dating, and the relative impact of fossil versus molecular sampling. We show that divergence time estimation is impacted by the model of fossil preservation, sampling intensity and tree shape. The addition of sequence data may improve molecular clock estimates, but accuracy and precision is dominated by the quality of the fossil calibrations. Posterior means and medians are poor representatives of true divergence times; posterior intervals provide a much more accurate estimate of divergence times, though they may be wide and often do not have high coverage probability. Our results highlight the importance of increased fossil sampling and improved statistical approaches to generating calibrations, which should incorporate the non-uniform nature of ecological and temporal fossil species distributions.

Fossilized embryos afford direct insight into the pattern of development in extinct organisms, providing unique tests of hypotheses of developmental evolution based in comparative embryology. However, these fossils can only be effective in this role if their embryology and phylogenetic affinities are well constrained. We elucidate and interpret the development of Olivooides from embryonic and adult stages and use these data to discriminate among competing interpretations of their anatomy and affinity. The embryology of Olivooides is principally characterized by the development of an ornamented periderm that initially forms externally and is subsequently formed internally, released at the aperture, facilitating the direct development of the embryo into an adult theca. Internal anatomy is known only from embryonic stages, revealing two internal tissue layers, the innermost of which is developed into three transversally arranged walls that partly divide the lumen into an abapertural region, interpreted as the gut of a polyp, and an adapertural region that includes structures that resemble the peridermal teeth of coronate scyphozoans. The anatomy and pattern of development exhibited by Olivooides appears common to the other known genus of olivooid, Quadrapyrgites, which differs in its tetraradial, as opposed to pentaradial symmetry. We reject previous interpretations of the olivooids as cycloneuralians, principally on the grounds that they lack a through gut and introvert, in embryo and adult. Instead we consider the affinities of the olivooids among medusozoan cnidarians; our phylogenetic analysis supports their classification as total-group Coronata, within crown-Scyphozoa. Olivooides and Quadrapyrgites evidence a broader range of life history strategies and bodyplan symmetry than is otherwise commonly represented in extant Scyphozoa specifically, and Cnidaria more generally.

NOTE: PLEASE ALSO SEE THE CORRIGENDUM TO THE ORIGINAL ARTICLE, PUBLISHED AT http://dx.doi.org/10.1111/pala.12193. Cyanobacteria are among the most ancient of evolutionary lineages, oxygenic photosynthesizers that may have originated before 3.0 Ga, as evidenced by free oxygen levels. Throughout the Precambrian, cyanobacteria were one of the most important drivers of biological innovations, strongly impacting early Earth's environments. At the end of the Archean Eon, they were responsible for the rapid oxygenation of Earth's atmosphere during an episode referred to as the Great Oxidation Event (GOE). However, little is known about the origin and diversity of early cyanobacterial taxa, due to: (1) the scarceness of Precambrian fossil deposits; (2) limited characteristics for the identification of taxa; and (3) the poor preservation of ancient microfossils. Previous studies based on 16S rRNA have suggested that the origin of multicellularity within cyanobacteria might have been associated with the GOE. However, single-gene analyses have limitations, particularly for deep branches. We reconstructed the evolutionary history of cyanobacteria using genome scale data and re-evaluated the Precambrian fossil record to get more precise calibrations for a relaxed clock analysis. For the phylogenomic reconstructions, we identified 756 conserved gene sequences in 65 cyanobacterial taxa, of which eight genomes have been sequenced in this study. Character state reconstructions based on maximum likelihood and Bayesian phylogenetic inference confirm previous findings, of an ancient multicellular cyanobacterial lineage ancestral to the majority of modern cyanobacteria. Relaxed clock analyses provide firm support for an origin of cyanobacteria in the Archean and a transition to multicellularity before the GOE. It is likely that multicellularity had a greater impact on cyanobacterial fitness and thus abundance, than previously assumed. Multicellularity, as a major evolutionary innovation, forming a novel unit for selection to act upon, may have served to overcome evolutionary constraints and enabled diversification of the variety of morphotypes seen in cyanobacteria today.