A phylogenomic resolution of the sea urchin tree of life

Maternal TGF-β ligand Panda breaks the radial symmetry of the sea urchin embryo by antagonizing the Nodal type II receptor ACVRII

In the highly regulative embryo of the sea urchin Paracentrotus lividus, establishment of the dorsal-ventral (D/V) axis critically depends on the zygotic expression of the TGF-β nodal in the ventral ectoderm. nodal expression is first induced ubiquitously in the 32-cell embryo and becomes progressively restricted to the presumptive ventral ectoderm by the early blastula stage. This early spatial restriction of nodal expression is independent of Lefty, and instead relies on the activity of Panda, a maternally expressed TGF-β ligand related to Lefty and Inhibins, which is required maternally for D/V axis specification. However, the mechanism by which Panda restricts the early nodal expression has remained enigmatic and it is not known if Panda works like a BMP ligand by opposing Nodal and antagonizing Smad2/3 signaling, or if it works like Lefty by sequestering an essential component of the Nodal signaling pathway. In this study, we report that Panda functions as an antagonist of the TGF-β type II receptor ACVRII (Activin receptor type II), which is the only type II receptor for Nodal signaling in the sea urchin and is also a type II receptor for BMP ligands. Inhibiting translation of acvrII mRNA disrupted D/V patterning across all 3 germ layers and caused acvrII morphants to develop with a typical Nodal loss-of-function phenotype. In contrast, embryos overexpressing acvrII displayed strong ectopic Smad1/5/8 signaling at blastula stages and developed as dorsalized larvae, a phenotype very similar to that caused by over activation of BMP signaling. Remarkably, embryos co-injected with acvrII mRNA and panda mRNA did not show ectopic Smad1/5/8 signaling and developed with a largely normal dorsal-ventral polarity. Furthermore, using an axis induction assay, we found that Panda blocks the ability of ACVRII to orient the D/V axis when overexpressed locally. Using co-immunoprecipitation, we showed that Panda physically interacts with ACVRII, as well as with the Nodal co-receptor Cripto, and with TBR3 (Betaglycan), which is a non-signaling receptor for Inhibins in mammals. At the molecular level, we have traced back the antagonistic activity of Panda to the presence of a single proline residue, conserved with all the Lefty factors, in the ACVRII binding motif of Panda, instead of a serine as in most of TGF-β ligands. Conversion of this proline to a serine converted Panda from an antagonist that opposed Nodal signaling and promoted dorsalization to an agonist that promoted Nodal signaling and triggered ventralization when overexpressed. Finally, using phylogenomics, we analyzed the emergence of the agonist and antagonist form of Panda in the course of evolution. Our data are consistent with the idea that the presence of a serine at that position, like in most TGF-β, was the ancestral condition and that the initial function of Panda was possibly in promoting and not in antagonizing Nodal signaling. These results highlight the existence of key functional and structural elements conserved between Panda and Lefty, allow to draw an intriguing parallel between sea urchin Panda and mammalian Inhibin α and raise the unexpected possibility that the original function of Panda may have been in activation of the Nodal pathway rather than in its inhibition.

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.

The mitochondrial genome of the deep-sea pyramid urchin Echinocrepis rostrata (Echinoidea: Holasteroida: Pourtalesiidae)

We present the mitochondrial genome of the deep-sea, epibenthic, irregular echinoid Echinocrepis rostrata, representing the first sequenced mitogenome of the order Holasteroida. The length of the complete E. rostrata mitochondrial genome is 15,716 base pairs, and its GC content is 34.87%. It contains 13 protein-coding genes, two rRNA genes, and 22 tRNA genes, whose order is identical to that of all other available echinoid mitogenomes. Phylogenetic analysis of available mitochondrial genomes, based on all coding loci, places E. rostrata as the sister group to spatangoids (heart urchins).

Taxonomic reassessment of Tetrapygus niger (Arbacioida, Echinoidea): molecular and morphological evidence support its placement in Arbacia

The echinoid genus Tetrapygus was initially described by L. Agassiz (1841) based on a single species, Tetrapygus niger Molina, 1782. Since the extensive work conducted by Mortensen (1935), Tetrapygus has received limited taxonomic attention over the past century. Recent discoveries of new fossil species of Arbacia Gray, 1835 from the upper Pliocene of northern Chile revealed striking morphological similarities between the two distinct Arbaciidae genera Arbacia and Tetrapygus. These findings compelled new investigations to evaluate the taxonomic status of these genera. Based on molecular mitochondrial (COI), nuclear (28S), and morphological evidence, Tetrapygus niger is here recovered as the sister species to Arbacia dufresnii, both species forming a clade within the phylogeny of South American species of Arbacia. Consequently, the diagnosis and description of Tetrapygus niger are here revised, and the species is reattributed to Arbacia, as previously proposed by A. Agassiz in Agassiz & Desor (1846) under the species name Arbacia nigra. An emended diagnosis of Arbacia is also proposed in light of these new findings.

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.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Echinoid skeleton: an insight on the species-specific pattern of the Paracentrotus lividus plate and its microstructural variability

The skeletal plates of echinoids consist of a peculiar lightweight structure, called stereom, which is organized in a porous three-dimensional lattice-like meshwork. The stereom is characterized by an extremely complex and diverse microarchitecture, largely varying not only from species to species but also among different test plates. It consists of different basic types combined in extremely different ways according to specific functional needs, creating species-specific structural patterns. These patterns can lead to specific mechanical behaviours, which can inspire biomimetic technology and design development. In this framework, the present study aimed to characterize the species-specific pattern of the Paracentrotus lividus interambulacral plate and the main microstructural features regarding its geometrical variability and mechanical responses. The results achieved quantitatively highlighted the differences between the analysed stereom types providing new insights regarding their topological configuration and isotropic and anisotropic behaviour. Interestingly, data also revealed that the galleried stereom present at the tubercle is significantly different from the one located at the suture. These analyses and findings are encouraging and provide a starting point for future research to unravel the wide range of mechanical strategies evolved in the echinoid skeletal structure.

Trioecy is maintained as a time-stable mating system in the pink sea urchin Toxopneustes roseus from the Mexican Pacific

Trioecy is a sexual system that consists of the co-occurrence of females, males and hermaphrodites in a population and is common in plants; however, in animals it is uncommon and poorly understood. In echinoderms, trioecy had never been recorded until now. Frequencies of females, males, and hermaphrodites were evaluated and gametogenic development was histologically characterized in a population of Toxopneustes roseus inhabiting the Mexican Pacific. Trioecy in this population is functional and temporally stable, since the three sexes coexisted in each sampling month. The hermaphrodites presented similar gametogenic development as the females and males and participated during the spawning season, contributing to the population's reproductive process. Trioecy is considered an evolutionarily transitory state, and it is extremely difficult to explain its presence in a species. We hypothesize that continuous ocean warming represents a threat to the survival of this population of T. roseus, since its early developmental stages, which represent a population bottleneck, are more vulnerable to high temperatures than other sea urchins inhabiting the area, while its population density is significantly lower. These conditions generate a strongly stressed environment, which is the determining factor that maintains the stability of trioecy in the species in which it has been studied.

Developmental atlas of the indirect-developing sea urchin Paracentrotus lividus: From fertilization to juvenile stages

The sea urchin Paracentrotus lividus has been used as a model system in biology for more than a century. Over the past decades, it has been at the center of a number of studies in cell, developmental, ecological, toxicological, evolutionary, and aquaculture research. Due to this previous work, a significant amount of information is already available on the development of this species. However, this information is fragmented and rather incomplete. Here, we propose a comprehensive developmental atlas for this sea urchin species, describing its ontogeny from fertilization to juvenile stages. Our staging scheme includes three periods divided into 33 stages, plus 15 independent stages focused on the development of the coeloms and the adult rudiment. For each stage, we provide a thorough description based on observations made on live specimens using light microscopy, and when needed on fixed specimens using confocal microscopy. Our descriptions include, for each stage, the main anatomical characteristics related, for instance, to cell division, tissue morphogenesis, and/or organogenesis. Altogether, this work is the first of its kind providing, in a single study, a comprehensive description of the development of P. lividus embryos, larvae, and juveniles, including details on skeletogenesis, ciliogenesis, myogenesis, coelomogenesis, and formation of the adult rudiment as well as on the process of metamorphosis in live specimens. Given the renewed interest for the use of sea urchins in ecotoxicological, developmental, and evolutionary studies as well as in using marine invertebrates as alternative model systems for biomedical investigations, this study will greatly benefit the scientific community and will serve as a reference for specialists and non-specialists interested in studying sea urchins.

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.

Dorsal-ventral axis formation in sea urchin embryos

Most sea urchin species produce planktonic feeding larvae with distinct dorsal-ventral polarity. Such morphological indicators of polarity arise after gastrulation, when several morphogenesis and cell differentiation events occur differentially along the dorsal-ventral axis. For instance, the gut bends toward the ventral side where the mouth will form, skeletogenesis occurs initially near the ventral side with the forming skeleton extending dorsally, and pigment cells differentiate and embed in the dorsal ectoderm. The patterning mechanisms and gene regulatory networks underlying these events have been extensively studied. Two opposing TGF-β signaling pathways, Nodal and BMP, play key roles in all three germ layers to respectively pattern the sea urchin ventral and dorsal sides. In this chapter, I describe our current understanding of sea urchin dorsal-ventral patterning mechanisms. Additionally, differences in the patterning mechanisms observed in lecithotrophic sea urchins (nonfeeding larvae) and in cidaroid sea urchins are also discussed, along with evolutionary insights gained from comparative analyses.

Mechanical defensive adaptations of three Mediterranean sea urchin species

In the Mediterranean, Paracentrotus lividus and Sphaerechinus granularis are important drivers of benthic ecosystems, often coexisting in sublittoral communities. However, the introduction of the invasive diadematoid Diadema setosum, which utilizes venomous spines, may affect these communities. To describe the mechanical properties of the test and spines of these three species, specimens were collected in winter of 2019 from the sublittoral zone of the Dodecanese island complex, southeastern Aegean Sea. This region serves as a gateway for invasive species to the Mediterranean Sea. Crushing test was conducted on live individuals, while 3-point bending test was used to estimate spine stiffness. Porosity and mineralogy of the test and spine, thickness of the test, and breaking length of the spine were measured and compared, while the microstructural architecture was also determined. The test of S. granularis was the most robust (194.35 ± 59.59 N), while the spines of D. setosum (4.76 ± 2.13 GPa) exhibited highest flexibility. Increased porosity and thickness of the test were related to increased robustness, whereas increased flexibility of the spine was attributed to high porosity, indicating that porosity in the skeleton plays a key role in preventing fracture. The spines of S. granularis exhibited highest length after fracture % (71.54 ± 5.5%). D. setosum exhibited higher values of Mg concentration in the test (10%) compared with the spines (4%). For the first time, the mineralogy of an invasive species is compared with its native counterpart, while a comparison of the mechanical properties of different species of the same ecosystem also takes place. This study highlights different ways, in which sea urchins utilize their skeleton and showcases the ecological significance of these adaptations, one of which is the different ways of utilization of the skeleton for defensive purposes, while the other is the ability of D. setosum to decrease the Mg % of its skeleton degrading its mechanical properties, without compromising its defense, by depending on venomous bearing spines. This enables this species to occupy not only tropical habitats, where it is indigenous, but also temperate like the eastern Mediterranean, which it has recently invaded.

Sequence Diversity, Locus Structure, and Evolutionary History of the SpTransformer Genes in the Sea Urchin Genome

The generation of large immune gene families is often driven by evolutionary pressure exerted on host genomes by their pathogens, which has been described as the immunological arms race. The SpTransformer (SpTrf) gene family from the California purple sea urchin, Strongylocentrotus purpuratus, is upregulated upon immune challenge and encodes the SpTrf proteins that interact with pathogens during an immune response. Native SpTrf proteins bind both bacteria and yeast, and augment phagocytosis of a marine Vibrio, while a recombinant SpTrf protein (rSpTrf-E1) binds a subset of pathogens and a range of pathogen associated molecular patterns. In the sequenced sea urchin genome, there are four SpTrf gene clusters for a total of 17 genes. Here, we report an in-depth analysis of these genes to understand the sequence complexities of this family, its genomic structure, and to derive a putative evolutionary history for the formation of the gene clusters. We report a detailed characterization of gene structure including the intron type and UTRs with conserved transcriptional start sites, the start codon and multiple stop codons, and locations of polyadenylation signals. Phylogenetic and percent mismatch analyses of the genes and the intergenic regions allowed us to predict the last common ancestral SpTrf gene and a theoretical evolutionary history of the gene family. The appearance of the gene clusters from the theoretical ancestral gene may have been driven by multiple duplication and deletion events of regions containing SpTrf genes. Duplications and ectopic insertion events, indels, and point mutations in the exons likely resulted in the extant genes and family structure. This theoretical evolutionary history is consistent with the involvement of these genes in the arms race in responses to pathogens and suggests that the diversification of these genes and their encoded proteins have been selected for based on the survival benefits of pathogen binding and host protection.

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.

The first complete mitochondrial genome of Stomopneustes variolaris (Lamarck, 1816) from the Stomechinidae (Echinoidea: Stomopneustoida)

Sea urchins (Echinoidea) are key components of marine benthic communities and many are commercially important fishery resources as luxury and healthy seafood. However, despite their high ecological and economic value, the mitochondrial genomes of all sea urchins have yet to be analyzed. In this study, we report the first complete mitochondrial genome of Stomopneustidae from Stomopneustes variolaris. The mitogenome has 15,767 base pairs (59.77% A + T content) and contains 37 genes (13 protein-coding, 22 transfer RNAs and 2 ribosomal RNAs), plus a putative control region. This study provides useful molecular resources for clarifying evolutionary and phylogenetic histories of sea urchins.

Exploring adaptive landscapes across deep time: A case study using echinoid body size

Adaptive landscapes are a common way of conceptualizing the phenotypic evolution of lineages across deep time. Although multiple approaches exist to implement this concept into operational models of trait evolution, inferring adaptive landscapes from comparative datasets remains challenging. Here, I explore the macroevolutionary dynamics of echinoid body size using data from over 5000 specimens and a phylogenetic framework incorporating a dense fossil sampling and spanning approximately 270 million years. Furthermore, I implement a novel approach of exploring alternative parameterizations of adaptive landscapes that succeeds in finding simpler, yet better-fitting models. Echinoid body size has been constrained to evolve within a single adaptive optimum for much of the clade's history. However, most of the morphological disparity of echinoids was generated by multiple regime shifts that drove the repeated evolution of miniaturized and gigantic forms. Events of body size innovation occurred predominantly in the Late Cretaceous and were followed by a drastic slowdown following the Cretaceous-Paleogene mass extinction. The discovery of these patterns is contingent upon directly sampling fossil taxa. The macroevolution of echinoid body size is therefore characterized by a late increase in disparity (likely linked to an expansion of ecospace), generated by active processes driving lineages toward extreme morphologies.

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.

Characterization of the complete mitochondrial genome of Macrotocinclus affinis (Siluriformes; Loricariidae) and phylogenetic studies of Siluriformes

The vertebrate mitochondrial genome is typically circular molecules made up of 14,000 to 16,000 bp, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (12 s rRNA and 16 s rRNA) and a control region. Compared with nuclear DNA, mitochondrial DNA has a higher mutation rate, so it is one of the most effective and reliable molecular markers in fish phylogeny. Macrotocinclus affinis was the only species in Macrotocinclus (it was classified as Otocinclus in the past) and currently lacks genetic information. Most of the current researches are based on the mitochondrial Cytb gene and RAG1 and RAG2 nuclear genes to study the phylogenetic analysis of Siluriformes. So, the study provides the characteristic features of the Macrotocinclus affinis mitochondrial genome and this is the first time that the phylogenetic relationship of Siluriformes has been reconstructed based on COI. We aimed to sequence the entire mitochondrial genome of Macrotocinclus affinis using conventional PCR techniques and to clarify its phylogenetic status in Siluriformes by using the COI sequence of mitochondria. In this study, we sequenced the whole mitochondrial genome of this species yielding a 16,632 bp circular assembly composed of the typical vertebrate mitochondrial features. It contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a putative control region, and one origin of replication on the light-strand. The overall base composition includes A (30.07%), T (24.43%), C (29.43%) and G (16.01%). The genome composition is A + T biased (54.5%), and exhibits AT-skew (0.1036) and GC-skew (-0.2962). Moreover, the 13 PCGs encode 3850 amino acids in total. The result of the phylogenetic tree supports Macrotocinclus affinis has a closest relationship with Otocinclus cf. hoppei far. These results will help to understand the characteristics of the mitochondrial genome of Macrotocinclus affinis and provide molecular basis for the evolutionary relationship of Loricariidae.

Revising transcriptome assemblies with phylogenetic information

A common transcriptome assembly error is to mistake different transcripts of the same gene as transcripts from multiple closely related genes. This error is difficult to identify during assembly, but in a phylogenetic analysis such errors can be diagnosed from gene phylogenies where they appear as clades of tips from the same species with improbably short branch lengths. treeinform is a method that uses phylogenetic information across species to refine transcriptome assemblies within species. It identifies transcripts of the same gene that were incorrectly assigned to multiple genes and reassign them as transcripts of the same gene. The treeinform method is implemented in Agalma, available at https://bitbucket.org/caseywdunn/agalma, and the general approach is relevant in a variety of other contexts.

Using target enrichment sequencing to study the higher-level phylogeny of the largest lichen-forming fungi family: Parmeliaceae (Ascomycota)

Parmeliaceae is the largest family of lichen-forming fungi with a worldwide distribution. We used a target enrichment data set and a qualitative selection method for 250 out of 350 genes to infer the phylogeny of the major clades in this family including 81 taxa, with both subfamilies and all seven major clades previously recognized in the subfamily Parmelioideae. The reduced genome-scale data set was analyzed using concatenated-based Bayesian inference and two different Maximum Likelihood analyses, and a coalescent-based species tree method. The resulting topology was strongly supported with the majority of nodes being fully supported in all three concatenated-based analyses. The two subfamilies and each of the seven major clades in Parmelioideae were strongly supported as monophyletic. In addition, most backbone relationships in the topology were recovered with high nodal support. The genus Parmotrema was found to be polyphyletic and consequently, it is suggested to accept the genus Crespoa to accommodate the species previously placed in Parmotrema subgen. Crespoa. This study demonstrates the power of reduced genome-scale data sets to resolve phylogenetic relationships with high support. Due to lower costs, target enrichment methods provide a promising avenue for phylogenetic studies including larger taxonomic/specimen sampling than whole genome data would allow.

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

Phylogenomic and paleontological data constitute complementary resources for unraveling the phylogenetic relationships and divergence times of lineages, yet few studies have attempted to fully integrate them. Several unique properties of echinoids (sea urchins) make them especially useful for such synthesizing approaches, including a remarkable fossil record that can be incorporated into explicit phylogenetic hypotheses. We revisit the phylogeny of crown group Echinoidea using a total-evidence dating approach that combines the largest phylogenomic data set for the clade, a large-scale morphological matrix with a dense fossil sampling, and a novel compendium of tip and node age constraints. To this end, we develop a novel method for subsampling phylogenomic data sets that selects loci with high phylogenetic signal, low systematic biases, and enhanced clock-like behavior. Our results demonstrate that combining different data sources increases topological accuracy and helps resolve conflicts between molecular and morphological data. Notably, we present a new hypothesis for the origin of sand dollars, and restructure the relationships between stem and crown echinoids in a way that implies a long stretch of undiscovered evolutionary history of the crown group in the late Paleozoic. Our efforts help bridge the gap between phylogenomics and phylogenetic paleontology, providing a model example of the benefits of combining the two. [Echinoidea; fossils; paleontology; phylogenomics; time calibration; total evidence.].

The first complete mitochondrial genome of the sand dollar Sinaechinocyamus mai (Echinoidea: Clypeasteroida)

Morphologic and molecular data often lead to different hypotheses of phylogenetic relationships. Such incongruence has been found in the echinoderm class Echinoidea. In particular, the phylogenetic status of the order Clypeasteroida is not well resolved. Complete mitochondrial genomes are currently available for 29 echinoid species, but no clypeasteroid had been sequenced to date. DNA extracted from a single live individual of Sinaechinocyamus mai was sequenced with 10× Genomics technology. This first complete mitochondrial genome (mitogenome) for the order Clypeasteroida is 15,756 base pairs in length. Phylogenomic analysis based on 34 ingroup taxa belonging to nine orders of the class Echinoidea show congruence between our new genetic inference and published trees based on morphologic characters, but also includes some intriguing differences that imply the need for additional investigation.

Trans-Species Polymorphism in Mitochondrial Genome of Camarodont Sea Urchins

Mitochondrial (mt) genomes of the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus demonstrate the identical patterns of intraspecific length variability of the ND6 gene, consisting of 489 bp (S variant) and 498 bp (L variant), respectively. For both species, the ND6 length difference is due to the 488A>G substitution, which changes the stop codon TAG in S variant for a tryptophan codon TGG in L variant and elongates the corresponding ND6 protein by three additional amino acids, Trp-Leu-Trp. The phylogenetic analysis based on mt genomes of sea urchins and related echinoderm groups from GenBank has shown the S and L ND6 variants as shared among the camarodont sea urchins; the rest of the echinoderms demonstrate the S variant only. The data suggest that the ND6 488A>G substitution can be the first example of the trans-species polymorphism in sea urchins, persisting at least since the time of the Odontophora diversification at the Eocene/Oligocene boundary (approximately 34 million years ago), which was characterized by an abrupt climate change and significant global ocean cooling. Alternative hypotheses, including the convergent RNA editing and/or codon reassignment, are not supported by direct comparisons of the ND6 gene sequences with the corresponding transcripts using the basic local alignment search tool (BLAST) of full sea urchin transcriptomes.