A Total-Evidence Dated Phylogeny of Echinoidea Combining Phylogenomic and Paleontological Data

Next-generation bioinformatics: An ultrafast and user-friendly tool for phylogenomic data exploration

With increasingly large genomic datasets, even routine bioinformatic tasks can be arduous, computationally demanding, and time-consuming. Additionally, most bioinformatic programs do not have a graphical user interface (GUI) and thus, require users to be minimally competent in command-line. These impediments present significant economic and technological barriers for practitioners who do not have access to advanced computational resources and support. In this issue of Molecular Ecology Resources, Handika and Esselstyn (2024) developed an ultrafast and memory-efficient bioinformatic tool, SEGUL, that performs common manipulations and calculations of summary statistics on phylogenomic datasets. SEGUL has two main features that distinguish it from other bioinformatic programs: (1) it is based on the recently emergent, high-performance programming language Rust, and (2) it has a GUI written using Flutter, a cross-platform programming framework that also supports mobile operating systems (mobile iOS, iPadOS and Android). By leveraging and combining the power of Rust and Flutter, SEGUL achieves significantly faster computation times and lower memory usage across different platforms and CPU architectures compared to similar programs. The unique and innovative features of SEGUL pave the way for a new era of bioinformatics that can be more energy-efficient, inclusive, and available to a broader swathe of users.

Phylogenomic analyses of ochrophytes (stramenopiles) with an emphasis on neglected lineages

Ochrophyta is a photosynthetic lineage that crowns the phylogenetic tree of stramenopiles, one of the major eukaryotic supergroups. Due to their ecological impact as a major primary producer, ochrophytes are relatively well-studied compared to the rest of the stramenopiles, yet their evolutionary relationships remain poorly understood. This is in part due to a number of missing lineages in large-scale multigene analyses, and an apparently rapid radiation leading to many short internodes between ochrophyte subgroups in the tree. These short internodes are also found across deep-branching lineages of stramenopiles with limited phylogenetic signal, leaving many relationships controversial overall. We have addressed this issue with other deep-branching stramenopiles recently, and now examine whether contentious relationships within the ochrophytes may be resolved with the help of filling in missing lineages in an updated phylogenomic dataset of ochrophytes, along with exploring various gene filtering criteria to identify the most phylogenetically informative genes. We generated ten new transcriptomes from various culture collections and a single-cell isolation from an environmental sample, added these to an existing phylogenomic dataset, and examined the effects of selecting genes with high phylogenetic signal or low phylogenetic noise. For some previously contentious relationships, we find a variety of analyses and gene filtering criteria consistently unite previously unstable groupings with strong statistical support. For example, we recovered a robust grouping of Eustigmatophyceae with Raphidophyceae-Phaeophyceae-Xanthophyceae while Olisthodiscophyceae formed a sister-lineage to Pinguiophyceae. Selecting genes with high phylogenetic signal or data quality recovered more stable topologies. Overall, we find that adding under-represented groups across different lineages is still crucial in resolving phylogenetic relationships, and discrete gene properties affect lineages of stramenopiles differently. This is something which may be explored to further our understanding of the molecular evolution of stramenopiles.

Scoutknife: A naïve, whole genome informed phylogenetic robusticity metric

Background: The phylogenetic bootstrap, first proposed by Felsenstein in 1985, is a critically important statistical method in assessing the robusticity of phylogenetic datasets. Core to its concept was the use of pseudo sampling - assessing the data by generating new replicates derived from the initial dataset that was used to generate the phylogeny. In this way, phylogenetic support metrics could overcome the lack of perfect, infinite data. With infinite data, however, it is possible to sample smaller replicates directly from the data to obtain both the phylogeny and its statistical robusticity in the same analysis. Due to the growth of whole genome sequencing, the depth and breadth of our datasets have greatly expanded and are set to only expand further. With genome-scale datasets comprising thousands of genes, we can now obtain a proxy for infinite data. Accordingly, we can potentially abandon the notion of pseudo sampling and instead randomly sample small subsets of genes from the thousands of genes in our analyses. Methods: We introduce Scoutknife, a jackknife-style subsampling implementation that generates 100 datasets by randomly sampling a small number of genes from an initial large-gene dataset to jointly establish both a phylogenetic hypothesis and assess its robusticity. We assess its effectiveness by using 18 previously published datasets and 100 simulation studies. Results: We show that Scoutknife is conservative and informative as to conflicts and incongruence across the whole genome, without the need for subsampling based on traditional model selection criteria. Conclusions: Scoutknife reliably achieves comparable results to selecting the best genes on both real and simulation datasets, while being resistant to the potential biases caused by selecting for model fit. As the amount of genome data grows, it becomes an even more exciting option to assess the robusticity of phylogenetic hypotheses.

Selective concentration of iron, titanium, and zirconium substrate minerals within Gregory's diverticulum, an organ unique to derived sand dollars (Echinoidea: Scutelliformes)

Gregory's diverticulum, a digestive tract structure unique to a derived group of sand dollars (Echinoidea: Scutelliformes), is filled with sand grains obtained from the substrate the animals inhabit. The simple methods of shining a bright light through a specimen or testing response to a magnet can reveal the presence of a mineral-filled diverticulum. Heavy minerals with a specific gravity of >2.9 g/cm3 are selectively concentrated inside the organ, usually at concentrations one order of magnitude, or more, greater than found in the substrate. Analyses of diverticulum content for thirteen species from nine genera, using optical mineralogy, powder X-ray diffraction, scanning electron microscopy and energy dispersive X-ray spectroscopy, as well as micro-computed tomography shows the preference for selection of five major heavy minerals: magnetite (Fe3O4), hematite (Fe2O3), ilmenite (FeTiO3), rutile (TiO2), and zircon (ZrSiO4). Minor amounts of heavy or marginally heavy amphibole, pyroxene and garnet mineral grains may also be incorporated. In general, the animals exhibit a preference for mineral grains with a specific gravity of >4.0 g/cm3, although the choice is opportunistic and the actual mix of mineral species depends on the mineral composition of the substrate. The animals also select for grain size, with mineral grains generally in the range of 50 to 150 μm, and do not appear to alter this preference during ontogeny. A comparison of analytical methods demonstrates that X-ray attenuation measured using micro-computed tomography is a reliable non-destructive method for heavy mineral quantification when supported by associated analyses of mineral grains extracted destructively from specimens or from substrate collected together with the specimens. Commonalities in the electro-chemical surface properties of the ingested minerals suggest that such characteristics play an important role in the selection process.

A 50-million-year-old, three-dimensionally preserved bat skull supports an early origin for modern echolocation

Bats are among the most recognizable, numerous, and widespread of all mammals. But much of their fossil record is missing, and bat origins remain poorly understood, as do the relationships of early to modern bats. Here, we describe a new early Eocene bat that helps bridge the gap between archaic stem bats and the hyperdiverse modern bat radiation of more than 1,460 living species. Recovered from ∼50 million-year-old cave sediments in the Quercy Phosphorites of southwestern France, Vielasia sigei's remains include a near-complete, three-dimensionally preserved skull-the oldest uncrushed bat cranium yet found. Phylogenetic analyses of a 2,665 craniodental character matrix, with and without 36.8 kb of DNA sequence data, place Vielasia outside modern bats, with total evidence tip-dating placing it sister to the crown clade. Vielasia retains the archaic dentition and skeletal features typical of early Eocene bats, but its inner ear shows specializations found in modern echolocating bats. These features, which include a petrosal only loosely attached to the basicranium, an expanded cochlea representing ∼25% basicranial width, and a long basilar membrane, collectively suggest that the kind of laryngeal echolocation used by most modern bats predates the crown radiation. At least 23 individuals of V. sigei are preserved together in a limestone cave deposit, indicating that cave roosting behavior had evolved in bats by the end of the early Eocene; this period saw the beginning of significant global climate cooling that may have been an evolutionary driver for bats to first congregate in caves.

Confusion will be my epitaph: genome-scale discordance stifles phylogenetic resolution of Holothuroidea

Sea cucumbers (Holothuroidea) are a diverse clade of echinoderms found from intertidal waters to the bottom of the deepest oceanic trenches. Their reduced skeletons and limited number of phylogenetically informative traits have long obfuscated morphological classifications. Sanger-sequenced molecular datasets have also failed to constrain the position of major lineages. Noteworthy, topological uncertainty has hindered a resolution for Neoholothuriida, a highly diverse clade of Permo-Triassic age. We perform the first phylogenomic analysis of Holothuroidea, combining existing datasets with 13 novel transcriptomes. Using a highly curated dataset of 1100 orthologues, our efforts recapitulate previous results, struggling to resolve interrelationships among neoholothuriid clades. Three approaches to phylogenetic reconstruction (concatenation under both site-homogeneous and site-heterogeneous models, and coalescent-aware inference) result in alternative resolutions, all of which are recovered with strong support and across a range of datasets filtered for phylogenetic usefulness. We explore this intriguing result using gene-wise log-likelihood scores and attempt to correlate these with a large set of gene properties. While presenting novel ways of exploring and visualizing support for alternative trees, we are unable to discover significant predictors of topological preference, and our efforts fail to favour one topology. Neoholothuriid genomes seem to retain an amalgam of signals derived from multiple phylogenetic histories.

Phylogeny, ancestral ranges and reclassification of sand dollars

Classification of the Class Echinoidea is under significant revision in light of emerging molecular phylogenetic evidence. In particular, the sister-group relationships within the superorder Luminacea (Echinoidea: Irregularia) have been considerably updated. However, the placement of many families remains largely unresolved due to a series of incongruent evidence obtained from morphological, paleontological, and genetic data for the majority of extant representatives. In this study, we investigated the phylogenetic relationships of 25 taxa, belonging to eleven luminacean families. We proposed three new superfamilies: Astriclypeoidea, Mellitoidea, and Taiwanasteroidea (including Dendrasteridae, Taiwanasteridae, Scutellidae, and Echinarachniidae), instead of the currently recognized superfamily Scutelloidea Gray, 1825. In light of the new data obtained from ten additional species, the historical biogeography reconstructed shows that the tropical western Pacific and eastern Indian Oceans are the cradle for early sand dollar diversification. Hothouse conditions during the late Cretaceous and early Paleogene were coupled with diversification events of major clades of sand dollars. We also demonstrate that Taiwan fauna can play a key role in terms of understanding the major Cenozoic migration and dispersal events in the evolutionary history of Luminacea.

Handling Logical Character Dependency in Phylogenetic Inference: Extensive Performance Testing of Assumptions and Solutions Using Simulated and Empirical Data

Logical character dependency is a major conceptual and methodological problem in phylogenetic inference of morphological data sets, as it violates the assumption of character independence that is common to all phylogenetic methods. It is more frequently observed in higher-level phylogenies or in data sets characterizing major evolutionary transitions, as these represent parts of the tree of life where (primary) anatomical characters either originate or disappear entirely. As a result, secondary traits related to these primary characters become "inapplicable" across all sampled taxa in which that character is absent. Various solutions have been explored over the last three decades to handle character dependency, such as alternative character coding schemes and, more recently, new algorithmic implementations. However, the accuracy of the proposed solutions, or the impact of character dependency across distinct optimality criteria, has never been directly tested using standard performance measures. Here, we utilize simple and complex simulated morphological data sets analyzed under different maximum parsimony optimization procedures and Bayesian inference to test the accuracy of various coding and algorithmic solutions to character dependency. This is complemented by empirical analyses using a recoded data set on palaeognathid birds. We find that in small, simulated data sets, absent coding performs better than other popular coding strategies available (contingent and multistate), whereas in more complex simulations (larger data sets controlled for different tree structure and character distribution models) contingent coding is favored more frequently. Under contingent coding, a recently proposed weighting algorithm produces the most accurate results for maximum parsimony. However, Bayesian inference outperforms all parsimony-based solutions to handle character dependency due to fundamental differences in their optimization procedures-a simple alternative that has been long overlooked. Yet, we show that the more primary characters bearing secondary (dependent) traits there are in a data set, the harder it is to estimate the true phylogenetic tree, regardless of the optimality criterion, owing to a considerable expansion of the tree parameter space. [Bayesian inference, character dependency, character coding, distance metrics, morphological phylogenetics, maximum parsimony, performance, phylogenetic accuracy.].

Syngnathoid Evolutionary History and the Conundrum of Fossil Misplacement

Seahorses, pipefishes, trumpetfishes, shrimpfishes, and allies are a speciose, globally distributed clade of fishes that have evolved a large number of unusual body plans. The clade that includes all these forms, Syngnathoidei, has become a model for the study of life history evolution, population biology, and biogeography. Yet, the timeline of syngnathoid evolution has remained highly contentious. This debate is largely attributable to the nature of the syngnathoid fossil record, which is both poorly described and patchy for several major lineages. Although fossil syngnathoids have been used to calibrate molecular phylogenies, the interrelationships of extinct species and their affinities to major living syngnathoid clades have scarcely been quantitatively tested. Here, I use an expanded morphological dataset to reconstruct the evolutionary relationships and clade ages of fossil and extant syngnathoids. Phylogenies generated using different analytical methodologies are largely congruent with molecular phylogenetic trees of Syngnathoidei but consistently find novel placements for several key taxa used as fossil calibrators in phylogenomic studies. Tip-dating of the syngnathoid phylogeny finds a timeline for their evolution that differs slightly from the one inferred using molecular trees but is generally congruent with a post-Cretaceous diversification event. These results emphasize the importance of quantitatively testing the relationships of fossil species, particularly when they are critical to assessing divergence times.

Worms and gills, plates and spines: the evolutionary origins and incredible disparity of deuterostomes revealed by fossils, genes, and development

Deuterostomes are the major division of animal life which includes sea stars, acorn worms, and humans, among a wide variety of ecologically and morphologically disparate taxa. However, their early evolution is poorly understood, due in part to their disparity, which makes identifying commonalities difficult, as well as their relatively poor early fossil record. Here, we review the available morphological, palaeontological, developmental, and molecular data to establish a framework for exploring the origins of this important and enigmatic group. Recent fossil discoveries strongly support a vermiform ancestor to the group Hemichordata, and a fusiform active swimmer as ancestor to Chordata. The diverse and anatomically bewildering variety of forms among the early echinoderms show evidence of both bilateral and radial symmetry. We consider four characteristics most critical for understanding the form and function of the last common ancestor to Deuterostomia: Hox gene expression patterns, larval morphology, the capacity for biomineralization, and the morphology of the pharyngeal region. We posit a deuterostome last common ancestor with a similar antero-posterior gene regulatory system to that found in modern acorn worms and cephalochordates, a simple planktonic larval form, which was later elaborated in the ambulacrarian lineage, the ability to secrete calcium minerals in a limited fashion, and a pharyngeal respiratory region composed of simple pores. This animal was likely to be motile in adult form, as opposed to the sessile origins that have been historically suggested. Recent debates regarding deuterostome monophyly as well as the wide array of deuterostome-affiliated problematica further suggest the possibility that those features were not only present in the last common ancestor of Deuterostomia, but potentially in the ur-bilaterian. The morphology and development of the early deuterostomes, therefore, underpin some of the most significant questions in the study of metazoan evolution.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Complete mitochondrial genomes of four deep-sea echinoids: conserved mitogenome organization and new insights into the phylogeny and evolution of Echinoidea

Echinoids are an important component in benthic marine environments, which occur at all depths from the shallow-water hard substrates to abyssal depths. To date, the phylogeny of the sea urchins and the macro-evolutionary processes of deep-sea and shallow water groups have not yet been fully resolved. In the present study, we sequenced the complete mitochondrial genomes (mitogenomes) of four deep-sea sea urchins (Echinoidea), which were the first representatives of the orders Aspidodiadematoida, Pedinoida and Echinothurioida, respectively. The gene content and arrangement were highly conserved in echinoid mitogenomes. The tRNA-Ser ^AGY with DHU arm was detected in the newly sequenced echinoid mitogenomes, representing an ancestral structure of tRNA-Ser ^AGY. No difference was found between deep-sea and shallow water groups in terms of base composition and codon usage. The phylogenetic analysis showed that all the orders except Spatangoida were monophyletic. The basal position of Cidaroida was supported. The closest relationship of Scutelloida and Echinolampadoida was confirmed. Our phylogenetic analysis shed new light on the position of Arbacioida, which supported that Arbacioida was most related with the irregular sea urchins instead of Stomopneustoida. The position Aspidodiadematoida (((Aspidodiadematoida + Pedinoida) + Echinothurioida) + Diadematoida) revealed by mitogenomic data discredited the hypothesis based on morphological evidences. The macro-evolutionary pattern revealed no simple onshore-offshore or an opposite hypothesis. But the basal position of the deep-sea lineages indicated the important role of deep sea in generating the current diversity of the class Echinoidea.

Phylogenomic analyses of echinoid diversification prompt a re-evaluation of their fossil record

Echinoids are key components of modern marine ecosystems. Despite a remarkable fossil record, the emergence of their crown group is documented by few specimens of unclear affinities, rendering their early history uncertain. The origin of sand dollars, one of its most distinctive clades, is also unclear due to an unstable phylogenetic context. We employ 18 novel genomes and transcriptomes to build a phylogenomic dataset with a near-complete sampling of major lineages. With it, we revise the phylogeny and divergence times of echinoids, and place their history within the broader context of echinoderm evolution. We also introduce the concept of a chronospace - a multidimensional representation of node ages - and use it to explore methodological decisions involved in time calibrating phylogenies. We find the choice of clock model to have the strongest impact on divergence times, while the use of site-heterogeneous models and alternative node prior distributions show minimal effects. The choice of loci has an intermediate impact, affecting mostly deep Paleozoic nodes, for which clock-like genes recover dates more congruent with fossil evidence. Our results reveal that crown group echinoids originated in the Permian and diversified rapidly in the Triassic, despite the relative lack of fossil evidence for this early diversification. We also clarify the relationships between sand dollars and their close relatives and confidently date their origins to the Cretaceous, implying ghost ranges spanning approximately 50 million years, a remarkable discrepancy with their rich fossil record.

Comprehensive species sampling and sophisticated algorithmic approaches refute the monophyly of Arachnida

Deciphering the evolutionary relationships of Chelicerata (arachnids, horseshoe crabs, and allied taxa) has proven notoriously difficult, due to their ancient rapid radiation and the incidence of elevated evolutionary rates in several lineages. Although conflicting hypotheses prevail in morphological and molecular data sets alike, the monophyly of Arachnida is nearly universally accepted, despite historical lack of support in molecular data sets. Some phylotranscriptomic analyses have recovered arachnid monophyly, but these did not sample all living orders, whereas analyses including all orders have failed to recover Arachnida. To understand this conflict, we assembled a data set of 506 high-quality genomes and transcriptomes, sampling all living orders of Chelicerata with high occupancy and rigorous approaches to orthology inference. Our analyses consistently recovered the nested placement of horseshoe crabs within a paraphyletic Arachnida. This result was insensitive to variation in evolutionary rates of genes, complexity of the substitution models, and alternative algorithmic approaches to species tree inference. Investigation of sources of systematic bias showed that genes and sites that recover arachnid monophyly are enriched in noise and exhibit low information content. To test the impact of morphological data, we generated a 514-taxon morphological data matrix of extant and fossil Chelicerata, analyzed in tandem with the molecular matrix. Combined analyses recovered the clade Merostomata (the marine orders Xiphosura, Eurypterida, and Chasmataspidida), but merostomates appeared nested within Arachnida. Our results suggest that morphological convergence resulting from adaptations to life in terrestrial habitats has driven the historical perception of arachnid monophyly, paralleling the history of numerous other invertebrate terrestrial groups.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Reconstructing Dipsacales phylogeny using Angiosperms353: issues and insights

PREMISE

Phylogenetic relationships within major angiosperm clades are increasingly well resolved, but largely informed by plastid data. Areas of poor resolution persist within the Dipsacales, including placement of Heptacodium and Zabelia, and relationships within the Caprifolieae and Linnaeeae, hindering our interpretation of morphological evolution. Here, we sampled a significant number of nuclear loci using a Hyb-Seq approach and used these data to infer the Dipsacales phylogeny and estimate divergence times.

METHODS

Sampling all major clades within the Dipsacales, we applied the Angiosperms353 probe set to 96 species. Data were filtered based on locus completeness and taxon recovery per locus, and trees were inferred using RAxML and ASTRAL. Plastid loci were assembled from off-target reads, and 10 fossils were used to calibrate dated trees.

RESULTS

Varying numbers of targeted loci and off-target plastomes were recovered from most taxa. Nuclear and plastid data confidently place Heptacodium with Caprifolieae, implying homoplasy in calyx morphology, ovary development, and fruit type. Placement of Zabelia, and relationships within the Caprifolieae and Linnaeeae, remain uncertain. Dipsacales diversification began earlier than suggested by previous angiosperm-wide dating analyses, but many major splitting events date to the Eocene.

CONCLUSIONS

The Angiosperms353 probe set facilitated the assembly of a large, single-copy nuclear dataset for the Dipsacales. Nevertheless, many relationships remain unresolved, and resolution was poor for woody clades with low rates of molecular evolution. We favor expanding the Angiosperms353 probe set to include more variable loci and loci of special interest, such as developmental genes, within particular clades.

Fossils improve phylogenetic analyses of morphological characters

Fossils provide our only direct window into evolutionary events in the distant past. Incorporating them into phylogenetic hypotheses of living clades can help time-calibrate divergences, as well as elucidate macroevolutionary dynamics. However, the effect fossils have on phylogenetic reconstruction from morphology remains controversial. The consequences of explicitly incorporating the stratigraphic ages of fossils using tip-dated inference are also unclear. Here, we use simulations to evaluate the performance of inference methods across different levels of fossil sampling and missing data. Our results show that fossil taxa improve phylogenetic analysis of morphological datasets, even when highly fragmentary. Irrespective of inference method, fossils improve the accuracy of phylogenies and increase the number of resolved nodes. They also induce the collapse of ancient and highly uncertain relationships that tend to be incorrectly resolved when sampling only extant taxa. Furthermore, tip-dated analyses under the fossilized birth-death process outperform undated methods of inference, demonstrating that the stratigraphic ages of fossils contain vital phylogenetic information. Fossils help to extract true phylogenetic signals from morphology, an effect that is mediated by both their distinctive morphology and their temporal information, and their incorporation in total-evidence phylogenetics is necessary to faithfully reconstruct evolutionary history.

Post-metamorphic skeletal growth in the sea urchin Paracentrotus lividus and implications for body plan evolution

Background

Understanding the molecular and cellular processes that underpin animal development are crucial for understanding the diversity of body plans found on the planet today. Because of their abundance in the fossil record, and tractability as a model system in the lab, skeletons provide an ideal experimental model to understand the origins of animal diversity. We herein use molecular and cellular markers to understand the growth and development of the juvenile sea urchin (echinoid) skeleton.

Results

We developed a detailed staging scheme based off of the first ~ 4 weeks of post-metamorphic life of the regular echinoid Paracentrotus lividus. We paired this scheme with immunohistochemical staining for neuronal, muscular, and skeletal tissues, and fluorescent assays of skeletal growth and cell proliferation to understand the molecular and cellular mechanisms underlying skeletal growth and development of the sea urchin body plan.

Conclusions

Our experiments highlight the role of skeletogenic proteins in accretionary skeletal growth and cell proliferation in the addition of new metameric tissues. Furthermore, this work provides a framework for understanding the developmental evolution of sea urchin body plans on macroevolutionary timescales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00174-1.