Evolution of phenotypic disparity in the plant kingdom.

The plant kingdom exhibits diverse bodyplans, from single-celled algae to complex multicellular land plants, but it is unclear how this phenotypic disparity was achieved. Here we show that the living divisions comprise discrete clusters within morphospace, separated largely by reproductive innovations, the extinction of evolutionary intermediates and lineage-specific evolution. Phenotypic complexity correlates not with disparity but with ploidy history, reflecting the role of genome duplication in plant macroevolution. Overall, the plant kingdom exhibits a pattern of episodically increasing disparity throughout its evolutionary history that mirrors the evolutionary floras and reflects ecological expansion facilitated by reproductive innovations. This pattern also parallels that seen in the animal and fungal kingdoms, suggesting a general pattern for the evolution of multicellular bodyplans.

Coevolution of enamel, ganoin, enameloid, and their matrix SCPP genes in osteichthyans.

We resolve debate over the evolution of vertebrate hypermineralized tissues through analyses of matrix protein-encoding secretory calcium-binding phosphoprotein (SCPP) genes and phylogenetic inference of hypermineralized tissues. Among these genes, AMBN and ENAM are found in both sarcopterygians and actinopterygians, whereas AMEL and SCPP5 are found only in sarcopterygians and actinopterygians, respectively. Actinopterygian AMBN, ENAM, and SCPP5 are expressed during the formation of hypermineralized tissues on scales and teeth: ganoin, acrodin, and collar enamel in gar, and acrodin and collar enameloid in zebrafish. Our phylogenetic analyses indicate the emergence of an ancestral enamel in stem-osteichthyans, whereas ganoin emerged in stem-actinopterygians and true enamel in stem-sarcopterygians. Thus, AMBN and ENAM originated in concert with ancestral enamel, SCPP5 evolved in association with ganoin, and AMEL evolved with true enamel. Shifts in gene expression domain and timing explain the evolution of different hypermineralized tissues. We propose that hypermineralized tissues in osteichthyans coevolved with matrix SCPP genes.

Arachnid monophyly: Morphological, palaeontological and molecular support for a single terrestrialization within Chelicerata.

The majority of extant arachnids are terrestrial, but other chelicerates are generally aquatic, including horseshoe crabs, sea spiders, and the extinct eurypterids. It is necessary to determine whether arachnids are exclusively descended from a single common ancestor (monophyly), because only that relationship is compatible with one land colonisation in chelicerate evolutionary history. Some studies have cast doubt on arachnid monophyly and recast the origins of their terrestrialization. These include some phylogenomic analyses placing horseshoe crabs within Arachnida, and from aquatic Palaeozoic stem-group scorpions. Here, we evaluate the possibility of arachnid monophyly by considering morphology, fossils and molecules holistically. We argue arachnid monophyly obviates the need to posit reacquisition/retention of aquatic characters such as gnathobasic feeding and book gills without trabeculae from terrestrial ancestors in horseshoe crabs, and that the scorpion total-group contains few aquatic taxa. We built a matrix composed of 200 slowly-evolving genes and re-analysed two published molecular datasets. We retrieved arachnid monophyly where other studies did not - highlighting the difficulty of resolving chelicerate relationships from current molecular data. As such, we consider arachnid monophyly the best-supported hypothesis. Finally, we inferred that arachnids terrestrialized during the Cambrian-Ordovician using the slow-evolving molecular matrix, in agreement with recent analyses.

Shifting spaces: Which disparity or dissimilarity measurement best summarize occupancy in multidimensional spaces?

Multidimensional analysis of traits are now common in ecology and evolution and are based on trait spaces in which each dimension summarizes the observed trait combination (a morphospace or an ecospace). Observations of interest will typically occupy a subset of this space, and researchers will calculate one or more measures to quantify how organisms inhabit that space. In macroevolution and ecology, these measures called disparity or dissimilarity metrics are generalized as space occupancy measures. Researchers use these measures to investigate how space occupancy changes through time, in relation to other groups of organisms, or in response to global environmental changes. However, the mathematical and biological meaning of most space occupancy measures is vague with the majority of widely used measures lacking formal description. Here, we propose a broad classification of space occupancy measures into three categories that capture changes in size, density, or position. We study the behavior of 25 measures to changes in trait space size, density, and position on simulated and empirical datasets. We find that no measure describes all of trait space aspects but that some are better at capturing certain aspects. Our results confirm the three broad categories (size, density, and position) and allow us to relate changes in any of these categories to biological phenomena. Because the choice of space occupancy measures is specific to the data and question, we introduced https://tguillerme.shinyapps.io/moms/moms, a tool to both visualize and capture changes in space occupancy for any measurement. https://tguillerme.shinyapps.io/moms/moms is designed to help workers choose the right space occupancy measures, given the properties of their trait space and their biological question. By providing guidelines and common vocabulary for space occupancy analysis, we hope to help bridging the gap in multidimensional research between ecology and evolution.

The complex effects of mass extinctions on morphological disparity.

Studies of biodiversity through deep time have been a staple for biologists and paleontologists for over 60 years. Investigations of species richness (diversity) revealed that at least five mass extinctions punctuated the last half billion years, each seeing the rapid demise of a large proportion of contemporary taxa. In contrast to diversity, the response of morphological diversity (disparity) to mass extinctions is unclear. Generally, diversity and disparity are decoupled, such that diversity may decline as morphological disparity increases, and vice versa. Here, we develop simulations to model disparity changes across mass extinctions using continuous traits and birth-death trees. We find no simple null for disparity change following a mass extinction but do observe general patterns. The range of trait values decreases following either random or trait-selective mass extinctions, whereas variance and the density of morphospace occupation only decline following trait-selective events. General trends may differentiate random and trait-selective mass extinctions, but methods struggle to identify trait selectivity. Long-term effects of mass extinction trait selectivity change support for phylogenetic comparative methods away from the simulated Brownian motion toward Ornstein-Uhlenbeck and Early Burst models. We find that morphological change over mass extinction is best studied by quantifying multiple aspects of morphospace occupation.

Disparities in the analysis of morphological disparity.

Analyses of morphological disparity have been used to characterize and investigate the evolution of variation in the anatomy, function and ecology of organisms since the 1980s. While a diversity of methods have been employed, it is unclear whether they provide equivalent insights. Here, we review the most commonly used approaches for characterizing and analysing morphological disparity, all of which have associated limitations that, if ignored, can lead to misinterpretation. We propose best practice guidelines for disparity analyses, while noting that there can be no 'one-size-fits-all' approach. The available tools should always be used in the context of a specific biological question that will determine data and method selection at every stage of the analysis.

A Cambrian-Ordovician Terrestrialization of Arachnids.

Understanding the temporal context of terrestrialization in chelicerates depends on whether terrestrial groups, the traditional Arachnida, have a single origin and whether or not horseshoe crabs are primitively or secondarily marine. Molecular dating on a phylogenomic tree that recovers arachnid monophyly, constrained by 27 rigorously vetted fossil calibrations, estimates that Arachnida originated during the Cambrian or Ordovician. After the common ancestor colonized the land, the main lineages appear to have rapidly radiated in the Cambrian-Ordovician boundary interval, coinciding with high rates of molecular evolution. The highest rates of arachnid diversification are detected between the Permian and Early Cretaceous. A pattern of ancient divergence estimates for terrestrial arthropod groups in the Cambrian while the oldest fossils are Silurian (seen in both myriapods and arachnids) is mirrored in the molecular and fossil records of land plants. We suggest the discrepancy between molecular and fossil evidence for terrestrialization is likely driven by the extreme sparseness of terrestrial sediments in the rock record before the late Silurian.

Origin of horsetails and the role of whole-genome duplication in plant macroevolution.

Whole-genome duplication (WGD) has occurred commonly in land plant evolution and it is often invoked as a causal agent in diversification, phenotypic and developmental innovation, as well as conferring extinction resistance. The ancient and iconic lineage of Equisetum is no exception, where WGD has been inferred to have occurred prior to the Cretaceous-Palaeogene (K-Pg) boundary, coincident with WGD events in angiosperms. In the absence of high species diversity, WGD in Equisetum is interpreted to have facilitated the long-term survival of the lineage. However, this characterization remains uncertain as these analyses of the Equisetum WGD event have not accounted for fossil diversity. Here, we analyse additional available transcriptomes and summarize the fossil record. Our results confirm support for at least one WGD event shared among the majority of extant Equisetum species. Furthermore, we use improved dating methods to constrain the age of gene duplication in geological time and identify two successive Equisetum WGD events. The two WGD events occurred during the Carboniferous and Triassic, respectively, rather than in association with the K-Pg boundary. WGD events are believed to drive high rates of trait evolution and innovations, but analysed trends of morphological evolution across the historical diversity of Equisetum provide little evidence for further macroevolutionary consequences following WGD. WGD events cannot have conferred extinction resistance to the Equisetum lineage through the K-Pg boundary since the ploidy events occurred hundreds of millions of years before this mass extinction and we find evidence of extinction among fossil polyploid Equisetum lineages. Our findings precipitate the need for a review of the proposed roles of WGDs in biological innovation and extinction survival in angiosperm and non-angiosperm lineages alike.

MCMCtreeR: functions to prepare MCMCtree analyses and visualize posterior ages on trees.

SUMMARY: The fossil record is incomplete, so molecular divergence time analysis is a crucial tool in estimating evolutionary timescales. MCMCtree contained in the PAML software provides Bayesian methods to estimate divergence times of genomic-sized sequences. Here, I present MCMCtreeR, a flexible R package to prepare time priors for MCMCtree analysis and plot time-scaled phylogenies. The package provides functions to refine parameters and visualize time-calibrated node prior distributions so that these priors accurately reflect confidence in known, usually fossil, time information. After the parameters have been chosen, the package produces output files ready for MCMCtree analysis. Following analysis, the package has tools to compare prior and posterior calibrated node age distributions and produce plots of the time-scaled phylogenies. The plotting functions allow for the inclusion of age uncertainty on time-scaled phylogenies, including the display of full posterior distributions on nodes. Options also allow for the inclusion of the geological timescale, and these plotting functions are applicable with posterior age estimates from any Bayesian divergence time estimation software. AVAILABILITY AND IMPLEMENTATION: MCMCtreeR is an R package available on CRAN (https://CRAN.R-project.org/package=MCMCtreeR). MCMCtreeR depends on the R packages ape, sn and stats4.

Characterization of melanosomes involved in the production of non-iridescent structural feather colours and their detection in the fossil record.

Non-iridescent structural colour in avian feathers is produced by coherent light scattering through quasi-ordered nanocavities in the keratin cortex of the barbs. To absorb unscattered light, melanosomes form a basal layer underneath the nanocavities. It has been shown that throughout Aves, melanosome morphology reflects broad categories of melanin-based coloration, as well as iridescence, allowing identification of palaeocolours in exceptionally preserved fossils. However, no studies have yet investigated the morphology of melanosomes in non-iridescent structural colour. Here, we analyse a wide sample of melanosomes from feathers that express non-iridescent structural colour from a phylogenetically broad range of extant avians to describe their morphology and compare them with other avian melanosome categories. We find that investigated melanosomes are typically wide (approx. 300 nm) and long (approx. 1400 nm), distinct from melanosomes found in black, brown and iridescent feathers, but overlapping significantly with melanosomes from grey feathers. This may suggest a developmental, and perhaps evolutionary, relationship between grey coloration and non-iridescent structural colours. We show that through analyses of fossil melanosomes, melanosomes indicative of non-iridescent structural colour can be predicted in an Eocene stem group roller ( Eocoracias: Coraciiformes) and with phylogenetic comparative methods the likely hue can be surmised. The overlap between melanosomes from grey and non-iridescent structurally coloured feathers complicates their distinction in fossil samples where keratin does not preserve. However, the abundance of grey coloration relative to non-iridescent structural coloration makes the former a more likely occurrence except in phylogenetically bracketed specimens like the specimen of Eocoracias studied here.

Evolution of metazoan morphological disparity.

The animal kingdom exhibits a great diversity of organismal form (i.e., disparity). Whether the extremes of disparity were achieved early in animal evolutionary history or clades continually explore the limits of possible morphospace is subject to continuing debate. Here we show, through analysis of the disparity of the animal kingdom, that, even though many clades exhibit maximal initial disparity, arthropods, chordates, annelids, echinoderms, and mollusks have continued to explore and expand the limits of morphospace throughout the Phanerozoic, expanding dramatically the envelope of disparity occupied in the Cambrian. The "clumpiness" of morphospace occupation by living clades is a consequence of the extinction of phylogenetic intermediates, indicating that the original distribution of morphologies was more homogeneous. The morphological distances between phyla mirror differences in complexity, body size, and species-level diversity across the animal kingdom. Causal hypotheses of morphologic expansion include time since origination, increases in genome size, protein repertoire, gene family expansion, and gene regulation. We find a strong correlation between increasing morphological disparity, genome size, and microRNA repertoire, but no correlation to protein domain diversity. Our results are compatible with the view that the evolution of gene regulation has been influential in shaping metazoan disparity whereas the invasion of terrestrial ecospace appears to represent an additional gestalt, underpinning the post-Cambrian expansion of metazoan disparity.

Integrated genomic and fossil evidence illuminates life's early evolution and eukaryote origin.

Establishing a unified timescale for the early evolution of Earth and life is challenging and mired in controversy because of the paucity of fossil evidence, the difficulty of interpreting it and dispute over the deepest branching relationships in the tree of life. Surprisingly, it remains perhaps the only episode in the history of life where literal interpretations of the fossil record hold sway, revised with every new discovery and reinterpretation. We derive a timescale of life, combining a reappraisal of the fossil material with new molecular clock analyses. We find the last universal common ancestor of cellular life to have predated the end of late heavy bombardment (>3.9 billion years ago (Ga)). The crown clades of the two primary divisions of life, Eubacteria and Archaebacteria, emerged much later (<3.4 Ga), relegating the oldest fossil evidence for life to their stem lineages. The Great Oxidation Event significantly predates the origin of modern Cyanobacteria, indicating that oxygenic photosynthesis evolved within the cyanobacterial stem lineage. Modern eukaryotes do not constitute a primary lineage of life and emerged late in Earth's history (<1.84 Ga), falsifying the hypothesis that the Great Oxidation Event facilitated their radiation. The symbiotic origin of mitochondria at 2.053-1.21 Ga reflects a late origin of the total-group Alphaproteobacteria to which the free living ancestor of mitochondria belonged.

Well-Annotated microRNAomes Do Not Evidence Pervasive miRNA Loss.

microRNAs are conserved noncoding regulatory factors implicated in diverse physiological and developmental processes in multicellular organisms, as causal macroevolutionary agents and for phylogeny inference. However, the conservation and phylogenetic utility of microRNAs has been questioned on evidence of pervasive loss. Here, we show that apparent widespread losses are, largely, an artefact of poorly sampled and annotated microRNAomes. Using a curated data set of animal microRNAomes, we reject the view that miRNA families are never lost, but they are rarely lost (92% are never lost). A small number of families account for a majority of losses (1.7% of families account for >45% losses), and losses are associated with lineages exhibiting phenotypic simplification. Phylogenetic analyses based on the presence/absence of microRNA families among animal lineages, and based on microRNA sequences among Osteichthyes, demonstrate the power of these small data sets in phylogenetic inference. Perceptions of widespread evolutionary loss of microRNA families are due to the uncritical use of public archives corrupted by spurious microRNA annotations, and failure to discriminate false absences that occur because of incomplete microRNAome annotation.

Probabilistic methods surpass parsimony when assessing clade support in phylogenetic analyses of discrete morphological data.

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.

The timescale of early land plant evolution.

Establishing the timescale of early land plant evolution is essential for testing hypotheses on the coevolution of land plants and Earth's System. The sparseness of early land plant megafossils and stratigraphic controls on their distribution make the fossil record an unreliable guide, leaving only the molecular clock. However, the application of molecular clock methodology is challenged by the current impasse in attempts to resolve the evolutionary relationships among the living bryophytes and tracheophytes. Here, we establish a timescale for early land plant evolution that integrates over topological uncertainty by exploring the impact of competing hypotheses on bryophyte-tracheophyte relationships, among other variables, on divergence time estimation. We codify 37 fossil calibrations for Viridiplantae following best practice. We apply these calibrations in a Bayesian relaxed molecular clock analysis of a phylogenomic dataset encompassing the diversity of Embryophyta and their relatives within Viridiplantae. Topology and dataset sizes have little impact on age estimates, with greater differences among alternative clock models and calibration strategies. For all analyses, a Cambrian origin of Embryophyta is recovered with highest probability. The estimated ages for crown tracheophytes range from Late Ordovician to late Silurian. This timescale implies an early establishment of terrestrial ecosystems by land plants that is in close accord with recent estimates for the origin of terrestrial animal lineages. Biogeochemical models that are constrained by the fossil record of early land plants, or attempt to explain their impact, must consider the implications of a much earlier, middle Cambrian-Early Ordovician, origin.

The Interrelationships of Land Plants and the Nature of the Ancestral Embryophyte.

The evolutionary emergence of land plant body plans transformed the planet. However, our understanding of this formative episode is mired in the uncertainty associated with the phylogenetic relationships among bryophytes (hornworts, liverworts, and mosses) and tracheophytes (vascular plants). Here we attempt to clarify this problem by analyzing a large transcriptomic dataset with models that allow for compositional heterogeneity between sites. Zygnematophyceae is resolved as sister to land plants, but we obtain several distinct relationships between bryophytes and tracheophytes. Concatenated sequence analyses that can explicitly accommodate site-specific compositional heterogeneity give more support for a mosses-liverworts clade, "Setaphyta," as the sister to all other land plants, and weak support for hornworts as the sister to all other land plants. Bryophyte monophyly is supported by gene concatenation analyses using models explicitly accommodating lineage-specific compositional heterogeneity and analyses of gene trees. Both maximum-likelihood analyses that compare the fit of each gene tree to proposed species trees and Bayesian supertree estimation based on gene trees support bryophyte monophyly. Of the 15 distinct rooted relationships for embryophytes, we reject all but three hypotheses, which differ only in the position of hornworts. Our results imply that the ancestral embryophyte was more complex than has been envisaged based on topologies recognizing liverworts as the sister lineage to all other embryophytes. This requires many phenotypic character losses and transformations in the liverwort lineage, diminishes inconsistency between phylogeny and the fossil record, and prompts re-evaluation of the phylogenetic affinity of early land plant fossils, the majority of which are considered stem tracheophytes.

Uncertain-tree: discriminating among competing approaches to the phylogenetic analysis of phenotype data.

Morphological data provide the only means of classifying the majority of life's history, but the choice between competing phylogenetic methods for the analysis of morphology is unclear. Traditionally, parsimony methods have been favoured but recent studies have shown that these approaches are less accurate than the Bayesian implementation of the Mk model. Here we expand on these findings in several ways: we assess the impact of tree shape and maximum-likelihood estimation using the Mk model, as well as analysing data composed of both binary and multistate characters. We find that all methods struggle to correctly resolve deep clades within asymmetric trees, and when analysing small character matrices. The Bayesian Mk model is the most accurate method for estimating topology, but with lower resolution than other methods. Equal weights parsimony is more accurate than implied weights parsimony, and maximum-likelihood estimation using the Mk model is the least accurate method. We conclude that the Bayesian implementation of the Mk model should be the default method for phylogenetic estimation from phenotype datasets, and we explore the implications of our simulations in reanalysing several empirical morphological character matrices. A consequence of our finding is that high levels of resolution or the ability to classify species or groups with much confidence should not be expected when using small datasets. It is now necessary to depart from the traditional parsimony paradigms of constructing character matrices, towards datasets constructed explicitly for Bayesian methods.

Partially incorrect fossil data augment analyses of discrete trait evolution in living species.

Ancestral state reconstruction of discrete character traits is often vital when attempting to understand the origins and homology of traits in living species. The addition of fossils has been shown to alter our understanding of trait evolution in extant taxa, but researchers may avoid using fossils alongside extant species if only few are known, or if the designation of the trait of interest is uncertain. Here, I investigate the impacts of fossils and incorrectly coded fossils in the ancestral state reconstruction of discrete morphological characters under a likelihood model. Under simulated phylogenies and data, likelihood-based models are generally accurate when estimating ancestral node values. Analyses with combined fossil and extant data always outperform analyses with extant species alone, even when around one quarter of the fossil information is incorrect. These results are especially pronounced when model assumptions are violated, such as when there is a trend away from the root value. Fossil data are of particular importance when attempting to estimate the root node character state. Attempts should be made to include fossils in analysis of discrete traits under likelihood, even if there is uncertainty in the fossil trait data.