Reorganization of sea urchin gene regulatory networks at least 268 million years ago as revealed by oldest fossil cidaroid echinoid

Single-Cell Transcriptomics Reveals Evolutionary Reconfiguration of Embryonic Cell Fate Specification in the Sea Urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of single-cell transcriptome analysis (scRNA-seq) developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states, respectively) reveal numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events and the primary signaling center are co-localized in the ancestral developmental gene regulatory network; remarkably, in H. erythrogramma, they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

Single-cell transcriptomics reveals evolutionary reconfiguration of embryonic cell fate specification in the sea urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of scRNA-seq developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states respectively) reveals numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events, and the primary signaling center are co-localized in the ancestral dGRN but remarkably, in H. erythrogramma they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

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.

Extreme phenotypic divergence and the evolution of development

As analyses of developmental mechanisms extend to ever more species, it becomes important to understand not just what is conserved or altered during evolution, but why. Closely related species that exhibit extreme phenotypic divergence can be uniquely informative in this regard. A case in point is the sea urchin genus Heliocidaris, which contains species that recently evolved a life history involving nonfeeding larvae following nearly half a billion years of prior evolution with feeding larvae. The resulting shift in selective regimes produced rapid and surprisingly extensive changes in developmental mechanisms that are otherwise highly conserved among echinoderm species. The magnitude and extent of these changes challenges the notion that conservation of early development in echinoderms is largely due to internal constraints that prohibit modification and instead suggests that natural selection actively maintains stability of inherently malleable trait developmental mechanisms over immense time periods. Knowing how and why natural selection changed during the evolution of nonfeeding larvae can also reveal why developmental mechanisms do and do not change in particular ways.

Perspectives on divergence of early developmental regulatory pathways: Insight from the evolution of echinoderm double negative gate

Evolution of gene regulatory networks (GRN) that orchestrate the highly coordinated course of development, is made possible by the network's robust nature for incorporating change without detrimental developmental outcome. It can be considered that the upstream network regulating early development, has immense influence over succeeding pathways thus may be less subjected to evolutionary modification. However, recent studies show incorporation of novel genes in such early developmental pathways such as the echinoderm pmar1 as evidence for drastic change occurring high in the GRN hierarchy. Here we discuss the mechanisms that underlie divergence of early developmental pathways utilizing promising insights from the evolution of echinoderm early mesoderm specification pathway of Pmar1-HesC double negative gate found solely in the euechinoid sea urchin lineage, as well as examples from other groups such as Spiralia and Drosophila.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates’ vascular tubulogenesis. On the other hand, vertebrates’ bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

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.

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

Inferred genetic architecture underlying evolution in a fossil stickleback lineage

Inferring the genetic architecture of evolution in the fossil record is difficult because genetic crosses are impossible, the acquisition of DNA is usually impossible and phenotype-genotype maps are rarely obvious. However, such inference is valuable because it reveals the genetic basis of microevolutionary change across many more generations than is possible in studies of extant taxa, thereby integrating microevolutionary process and macroevolutionary pattern. Here, we infer the genetic basis of pelvic skeleton reduction in Gasterosteus doryssus, a Miocene stickleback fish from a finely resolved stratigraphic sequence that spans nearly 17,000 years. Reduction in pelvic score, a categorical measure of pelvic structure, resulted primarily from reciprocal frequency changes of two discrete phenotypic classes. Pelvic vestiges also showed left-side larger asymmetry. These patterns implicate Pitx1, a large-effect gene whose deletion generates left-side larger asymmetry of pelvic vestiges in extant, closely related Gasterosteus aculeatus. In contrast, reductions in the length of the pelvic girdle and pelvic spines resulted from directional shifts of unimodal, continuous trait distributions, suggesting an additional suite of genes with minor, additive pelvic effects, again like G. aculeatus. Similar genetic architectures explain shared but phyletically independent patterns across 10 million years of stickleback evolution.

Evolution and Development at the Origin of a Phylum

Quantifying morphological evolution is key to determining the patterns and processes underlying the origin of phyla. We constructed a hierarchical morphological character matrix to characterize the radiation and establishment of echinoderm body plans during the early Paleozoic. This showed that subphylum-level clades diverged gradually through the Cambrian, and the distinctiveness of the resulting body plans was amplified by the extinction of transitional forms and obscured by convergent evolution during the Ordovician. Higher-order characters that define these body plans were not fixed at the origin of the phylum, countering hypotheses regarding developmental processes governing the early evolution of animals. Instead, these burdened characters were flexible, enabling continued evolutionary innovation throughout the clades' history.

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.

pmar1/phb homeobox genes and the evolution of the double-negative gate for endomesoderm specification in echinoderms

In several model animals, the earliest phases of embryogenesis are regulated by lineage-specific genes, such as Drosophila bicoid Sea urchin (echinoid) embryogenesis is initiated by zygotic expression of pmar1, a paired-class homeobox gene that has been considered to be present only in the lineage of modern urchins (euechinoids). In euechinoids, Pmar1 promotes endomesoderm specification by repressing the hairy and enhancer of split C (hesC) gene. Here, we have identified the basal echinoid (cidaroid) pmar1 gene, which also promotes endomesoderm specification but not by repressing hesC A further search for related genes demonstrated that other echinoderms have pmar1-related genes named phb Functional analyses of starfish Phb proteins indicated that, similar to cidaroid Pmar1, they promote activation of endomesoderm regulatory gene orthologs via an unknown repressor that is not HesC. Based on these results, we propose that Pmar1 may have recapitulated the regulatory function of Phb during the early diversification of echinoids and that the additional repressor HesC was placed under the control of Pmar1 in the euechinoid lineage. This case provides an exceptional model for understanding how early developmental processes diverge.

Echinoids from the Tesero Member (Werfen Formation) of the Dolomites (Italy): implications for extinction and survival of echinoids in the aftermath of the end-Permian mass extinction

The end-Permian mass extinction (∼252 Ma) was responsible for high rates of extinction and evolutionary bottlenecks in a number of animal groups. Echinoids, or sea urchins, were no exception, and the Permian to Triassic represents one of the most significant intervals of time in their macroevolutionary history. The extinction event was responsible for significant turnover, with the Permian–Triassic representing the transition from stem group echinoid-dominated faunas in the Palaeozoic to Mesozoic faunas dominated by crown group echinoids. This turnover is well-known, however, the environmental and taxonomic distribution of echinoids during the latest Permian and Early Triassic is not. Here we report on an echinoid fauna from the Tesero Member, Werfen Formation (latest Permian to Early Triassic) of the Dolomites (northern Italy). The fauna is largely known from disarticulated ossicles, but consists of both stem group taxa, and a new species of crown group echinoid, Eotiaris teseroensis n. sp. That these stem group echinoids were present in the Tesero Member indicates that stem group echinoids did not go extinct in the Dolomites coincident with the onset of extinction, further supporting other recent work indicating that stem group echinoids survived the end-Permian extinction. Furthermore, the presence of Eotiaris across a number of differing palaeoenvironments in the Early Triassic may have had implications for the survival of cidaroid echinoids during the extinction event.

Facies selectivity of benthic invertebrates in a Permian/Triassic boundary microbialite succession: Implications for the "microbialite refuge" hypothesis

Thrombolite and stromatolite habitats are becoming increasingly recognized as important refuges for invertebrates during Phanerozoic Oceanic Anoxic Events (OAEs); it is posited that oxygenic photosynthesis by cyanobacteria in these microbialites provided a refuge from anoxic conditions (i.e., the "microbialite refuge" hypothesis). Here, we test this hypothesis by investigating the distribution of ~34, 500 benthic invertebrate fossils found in ~100 samples from a microbialite succession that developed following the latest Permian mass extinction event on the Great Bank of Guizhou (South China), representing microbial (stromatolites and thrombolites) and non-microbial facies. The stromatolites were the least taxonomically diverse facies, and the thrombolites also recorded significantly lower diversities when compared to the non-microbial facies. Based on the distribution and ornamentation of the bioclasts within the thrombolites and stromatolites, the bioclasts are inferred to have been transported and concentrated in the non-microbial fabrics, that is, cavities around the microbial framework. Therefore, many of the identified metazoans from the post-extinction microbialites are not observed to have been living within a microbial mat. Furthermore, the lifestyle of many of the taxa identified from the microbialites was not suited for, or even amenable to, life within a benthic microbial mat. The high diversity of oxygen-dependent metazoans in the non-microbial facies on the Great Bank of Guizhou, and inferences from geochemical records, suggests that the microbialites and benthic communities developed in oxygenated environments, which disproves that the microbes were the source of the oxygenation. Instead, we posit that microbialite successions represent a taphonomic window for exceptional preservation of the biota, similar to a Konzentrat-Lagerstätte, which has allowed for diverse fossil assemblages to be preserved during intervals of poor preservation.

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.

Cell type phylogenetics informs the evolutionary origin of echinoderm larval skeletogenic cell identity

The multiplicity of cell types comprising multicellular organisms begs the question as to how cell type identities evolve over time. Cell type phylogenetics informs this question by comparing gene expression of homologous cell types in distantly related taxa. We employ this approach to inform the identity of larval skeletogenic cells of echinoderms, a clade for which there are phylogenetically diverse datasets of spatial gene expression patterns. We determined ancestral spatial expression patterns of alx1, ets1, tbr, erg, and vegfr, key components of the skeletogenic gene regulatory network driving identity of the larval skeletogenic cell. Here we show ancestral state reconstructions of spatial gene expression of extant eleutherozoan echinoderms support homology and common ancestry of echinoderm larval skeletogenic cells. We propose larval skeletogenic cells arose in the stem lineage of eleutherozoans during a cell type duplication event that heterochronically activated adult skeletogenic cells in a topographically distinct tissue in early development.

Conserved regulatory state expression controlled by divergent developmental gene regulatory networks in echinoids

Evolution of the animal body plan is driven by changes in developmental gene regulatory networks (GRNs), but how networks change to control novel developmental phenotypes remains, in most cases, unresolved. Here, we address GRN evolution by comparing the endomesoderm GRN in two echinoid sea urchins, Strongylocentrotus purpuratus and Eucidaris tribuloides, with at least 268 million years of independent evolution. We first analyzed the expression of twelve transcription factors and signaling molecules of the S. purpuratus GRN in E. tribuloides embryos, showing that orthologous regulatory genes are expressed in corresponding endomesodermal cell fates in the two species. However, perturbation of regulatory genes revealed that important regulatory circuits of the S. purpuratus GRN are significantly different in E. tribuloides For example, mesodermal Delta/Notch signaling controls exclusion of alternative cell fates in E. tribuloides but controls mesoderm induction and activation of a positive feedback circuit in S. purpuratus These results indicate that the architecture of the sea urchin endomesoderm GRN evolved by extensive gain and loss of regulatory interactions between a conserved set of regulatory factors that control endomesodermal cell fate specification.

A phylogenomic resolution of the sea urchin tree of life

BACKGROUND

Echinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved.

RESULTS

We performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We show that phylogenetic signal for novel resolutions of these lineages is strong and distributed throughout the genome, and fail to recover systematic biases as drivers of our results.

CONCLUSIONS

Our investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.

Paleogenomics of echinoids reveals an ancient origin for the double-negative specification of micromeres in sea urchins

Establishing a timeline for the evolution of novelties is a common, unifying goal at the intersection of evolutionary and developmental biology. Analyses of gene regulatory networks (GRNs) provide the ability to understand the underlying genetic and developmental mechanisms responsible for the origin of morphological structures both in the development of an individual and across entire evolutionary lineages. Accurately dating GRN novelties, thereby establishing a timeline for GRN evolution, is necessary to answer questions about the rate at which GRNs and their subcircuits evolve, and to tie their evolution to paleoenvironmental and paleoecological changes. Paleogenomics unites the fossil record and all aspects of deep time, with modern genomics and developmental biology to understand the evolution of genomes in evolutionary time. Recent work on the regulatory genomic basis of development in cidaroid echinoids, sand dollars, heart urchins, and other nonmodel echinoderms provides an ideal dataset with which to explore GRN evolution in a comparative framework. Using divergence time estimation and ancestral state reconstructions, we have determined the age of the double-negative gate (DNG), the subcircuit which specifies micromeres and skeletogenic cells in Strongylocentrotus purpuratus We have determined that the DNG has likely been used for euechinoid echinoid micromere specification since at least the Late Triassic. The innovation of the DNG thus predates the burst of post-Paleozoic echinoid morphological diversification that began in the Early Jurassic. Paleogenomics has wide applicability for the integration of deep time and molecular developmental data, and has wide utility in rigorously establishing timelines for GRN evolution.

Notch-mediated lateral inhibition is an evolutionarily conserved mechanism patterning the ectoderm in echinoids

Notch signaling is a crucial cog in early development of euechinoid sea urchins, specifying both non-skeletogenic mesodermal lineages and serotonergic neurons in the apical neuroectoderm. Here, the spatial distributions and function of delta, gcm, and hesc, three genes critical to these processes in euechinoids, are examined in the distantly related cidaroid sea urchin Eucidaris tribuloides. Spatial distribution and experimental perturbation of delta and hesc suggest that the function of Notch signaling in ectodermal patterning in early development of E. tr ibuloides is consistent with canonical lateral inhibition. Delta transcripts were observed in t he archenteron, apical ectoderm, and lateral ectoderm in gastrulating e mbryos of E. tribuloides. Perturbation of Notch signaling by either delta morpholino or treatment of DAPT downregulated hesc and upregulated delta and gcm, resulting in ectopic expression of delta and gcm. Similarly, hesc perturbation mirrored the effects of delta perturbation. Interestingly, perturbation of delta or hesc resulted in more cells expressing gcm and supernumerary pigment cells, suggesting that pigment cell proliferation is regulated by Notch in E. tribuloides. These results are consistent with an evolutionary scenario whereby, in the echinoid ancestor, Notch signaling was deployed in the ectoderm to specify neurogenic progenitors and controlled pigment cell proliferation in the dorsal ectoderm.

A new stem group echinoid from the Triassic of China leads to a revised macroevolutionary history of echinoids during the end-Permian mass extinction

The Permian-Triassic bottleneck has long been thought to have drastically altered the course of echinoid evolution, with the extinction of the entire echinoid stem group having taken place during the end-Permian mass extinction. The Early Triassic fossil record of echinoids is, however, sparse, and new fossils are paving the way for a revised interpretation of the evolutionary history of echinoids during the Permian-Triassic crisis and Early Mesozoic. A new species of echinoid, Yunnanechinus luopingensis n. sp. recovered from the Middle Triassic (Anisian) Luoping Biota fossil Lagerstätte of South China, displays morphologies that are not characteristic of the echinoid crown group. We have used phylogenetic analyses to further demonstrate that Yunnanechinus is not a member of the echinoid crown group. Thus a clade of stem group echinoids survived into the Middle Triassic, enduring the global crisis that characterized the end-Permian and Early Triassic. Therefore, stem group echinoids did not go extinct during the Palaeozoic, as previously thought, and appear to have coexisted with the echinoid crown group for at least 23 million years. Stem group echinoids thus exhibited the Lazarus effect during the latest Permian and Early Triassic, while crown group echinoids did not.

Eric Davidson and deep time

Eric Davidson had a deep and abiding interest in the role developmental mechanisms played in generating evolutionary patterns documented in deep time, from the origin of the euechinoids to the processes responsible for the morphological architectures of major animal clades. Although not an evolutionary biologist, Davidson's interests long preceded the current excitement over comparative evolutionary developmental biology. Here I discuss three aspects at the intersection between his research and evolutionary patterns in deep time: First, understanding the mechanisms of body plan formation, particularly those associated with the early diversification of major metazoan clades. Second, a critique of early claims about ancestral metazoans based on the discoveries of highly conserved genes across bilaterian animals. Third, Davidson's own involvement in paleontology through a collaborative study of the fossil embryos from the Ediacaran Doushantuo Formation in south China.

A Conserved Role for VEGF Signaling in Specification of Homologous Mesenchymal Cell Types Positioned at Spatially Distinct Developmental Addresses in Early Development of Sea Urchins

Comparative studies of early development in echinoderms are revealing the tempo and mode of alterations to developmental gene regulatory networks and to the cell types they specify. In euechinoid sea urchins, skeletogenic mesenchyme (SM) ingresses prior to gastrulation at the vegetal pole and aligns into a ring-like array with two bilateral pockets of cells, the sites where spiculogenesis will later occur. In cidaroid sea urchins, the anciently diverged sister clade to euechinoid sea urchins, a homologous SM cell type ingresses later in development, after gastrulation has commenced, and consequently at a distinct developmental address. Thus, a heterochronic shift of ingression of the SM cell type occurred in one of the echinoid lineages. In euechinoids, specification and migration of SM are facilitated by vascular endothelial growth factor (VEGF) signaling. We describe spatiotemporal expression of vegf and vegfr and experimental manipulations targeting VEGF signaling in the cidaroid Eucidaris tribuloides. Spatially, vegf and vegfr mRNA localizes similarly as in euechinoids, suggesting conserved deployment in echinoids despite their spatially distinct development addresses of ingression. Inhibition of VEGF signaling in E. tribuloides suggests its role in SM specification is conserved in echinoids. Temporal discrepancies between the onset of vegf expression and SM ingression likely result in previous observations of SM "random wandering" behavior. Our results indicate that, although the SM cell type in echinoids ingresses into distinct developmental landscapes, it retains a signaling mechanism that restricts their spatial localization to a conserved developmental address where spiculogenesis later occurs.

IL17 factors are early regulators in the gut epithelium during inflammatory response to Vibrio in the sea urchin larva

IL17 cytokines are central mediators of mammalian immunity. In vertebrates, these factors derive from diverse cellular sources. Sea urchins share a molecular heritage with chordates that includes the IL17 system. Here, we characterize the role of epithelial expression of IL17 in the larval gut-associated immune response. The purple sea urchin genome encodes 10 IL17 subfamilies (35 genes) and 2 IL17 receptors. Most of these subfamilies are conserved throughout echinoderms. Two IL17 subfamilies are sequentially strongly upregulated and attenuated in the gut epithelium in response to bacterial disturbance. IL17R1 signal perturbation results in reduced expression of several response genes including an IL17 subtype, indicating a potential feedback. A third IL17 subfamily is activated in adult immune cells indicating that expression in immune cells and epithelia is divided among families. The larva provides a tractable model to investigate the regulation and consequences of gut epithelial IL17 expression across the organism.

DOI:http://dx.doi.org/10.7554/eLife.23481.001

To protect themselves from the constant invasion of harmful microbes, animals have evolved complex immune systems. The gut is one of the most active sites of the immune system and plays a key role in regulating immune responses. In mammals, cells lining the gut wall can sense the presence of harmful bacteria and communicate this information to tissues across the body by producing specialized proteins called Interleukin-17 (IL-17). IL-17 proteins are important for regulating inflammation and are thought to activate specific immune cells in an infected area.

Some aspects of immune systems are similar between different animal species, which can provide clues of how immunity evolved and how it is regulated. For example, sea urchins, which evolved 400-600 million years ago, begin life as simple larvae consisting of a few thousand cells. As oceans harbor a multitude of bacteria and viruses, sea urchin larvae need an efficient immune system to defend themselves. These larvae can respond to specific types of bacteria within a few hours after the microbes have entered their gut by modifying gene expression in distant cells. As these changes occur in cells that are removed from the bacteria, it is thought that the gut cells that initially sense the bacteria, somehow communicate this information.

Now, Buckley et al. exposed sea urchin larvae to a marine bacterium and measured the responses of the cells and their gene expression. The infection affected several types of cells, and in the first 24 hours, a subset of immune cells changed shape and started migrating to the gut wall. In addition, IL-17 gene expression changed significantly in gut cells in the early phases of the larval immune response. Buckley et al. identified three types of IL-17 proteins involved in sea urchin immunity: two that are important for the immune response in the gut during the larval stage, and a third that is only present in adults. These findings suggest that IL-17 signaling is an ancient and central element of gut-associated immune response, which even exists in animals that evolved long before humans.

These findings demonstrate that the sea urchin larva represents a unique and ideal experimental model to study immune responses in a living organism that is more closely related to mammals than some other models, like fruit flies or worms. By understanding the fundamental mechanisms that mediate gut health, this work may highlight new drug targets to treat conditions like Crohn’s disease and colon cancer.

DOI:http://dx.doi.org/10.7554/eLife.23481.002

Divergence of ectodermal and mesodermal gene regulatory network linkages in early development of sea urchins

Developmental gene regulatory networks (GRNs) are assemblages of gene regulatory interactions that direct ontogeny of animal body plans. Studies of GRNs operating in the early development of euechinoid sea urchins have revealed that little appreciable change has occurred since their divergence ∼90 million years ago (mya). These observations suggest that strong conservation of GRN architecture was maintained in early development of the sea urchin lineage. Testing whether this holds for all sea urchins necessitates comparative analyses of echinoid taxa that diverged deeper in geological time. Recent studies highlighted extensive divergence of skeletogenic mesoderm specification in the sister clade of euechinoids, the cidaroids, suggesting that comparative analyses of cidaroid GRN architecture may confer a greater understanding of the evolutionary dynamics of developmental GRNs. Here I report spatiotemporal patterning of 55 regulatory genes and perturbation analyses of key regulatory genes involved in euechinoid oral-aboral patterning of nonskeletogenic mesodermal and ectodermal domains in early development of the cidaroid Eucidaris tribuloides These results indicate that developmental GRNs directing mesodermal and ectodermal specification have undergone marked alterations since the divergence of cidaroids and euechinoids. Notably, statistical and clustering analyses of echinoid temporal gene expression datasets indicate that regulation of mesodermal genes has diverged more markedly than regulation of ectodermal genes. Although research on indirect-developing euechinoid sea urchins suggests strong conservation of GRN circuitry during early embryogenesis, this study indicates that since the divergence of cidaroids and euechinoids, developmental GRNs have undergone significant, cell type-biased alterations.

Ancestral state reconstruction by comparative analysis of a GRN kernel operating in echinoderms

Diverse sampling of organisms across the five major classes in the phylum Echinodermata is beginning to reveal much about the structure and function of gene regulatory networks (GRNs) in development and evolution. Sea urchins are the most studied clade within this phylum, and recent work suggests there has been dramatic rewiring at the top of the skeletogenic GRN along the lineage leading to extant members of the euechinoid sea urchins. Such rewiring likely accounts for some of the observed developmental differences between the two major subclasses of sea urchins-cidaroids and euechinoids. To address effects of topmost rewiring on downstream GRN events, we cloned four downstream regulatory genes within the skeletogenic GRN and surveyed their spatiotemporal expression patterns in the cidaroid Eucidaris tribuloides. We performed phylogenetic analyses with homologs from other non-vertebrate deuterostomes and characterized their spatiotemporal expression by quantitative polymerase chain reaction (qPCR) and whole-mount in situ hybridization (WMISH). Our data suggest the erg-hex-tgif subcircuit, a putative GRN kernel, exhibits a mesoderm-specific expression pattern early in Eucidaris development that is directly downstream of the initial mesodermal GRN circuitry. Comparative analysis of the expression of this subcircuit in four echinoderm taxa allowed robust ancestral state reconstruction, supporting hypotheses that its ancestral function was to stabilize the mesodermal regulatory state and that it has been co-opted and deployed as a unit in mesodermal subdomains in distantly diverged echinoderms. Importantly, our study supports the notion that GRN kernels exhibit structural and functional modularity, locking down and stabilizing clade-specific, embryonic regulatory states.

Evolutionary rewiring of gene regulatory network linkages at divergence of the echinoid subclasses

Evolution of animal body plans occurs with changes in the encoded genomic programs that direct development, by alterations in the structure of encoded developmental gene-regulatory networks (GRNs). However, study of this most fundamental of evolutionary processes requires experimentally tractable, phylogenetically divergent organisms that differ morphologically while belonging to the same monophyletic clade, plus knowledge of the relevant GRNs operating in at least one of the species. These conditions are met in the divergent embryogenesis of the two extant, morphologically distinct, echinoid (sea urchin) subclasses, Euechinoidea and Cidaroidea, which diverged from a common late Paleozoic ancestor. Here we focus on striking differences in the mode of embryonic skeletogenesis in a euechinoid, the well-known model Strongylocentrotus purpuratus (Sp), vs. the cidaroid Eucidaris tribuloides (Et). At the level of descriptive embryology, skeletogenesis in Sp and Et has long been known to occur by distinct means. The complete GRN controlling this process is known for Sp. We carried out targeted functional analyses on Et skeletogenesis to identify the presence, or demonstrate the absence, of specific regulatory linkages and subcircuits key to the operation of the Sp skeletogenic GRN. Remarkably, most of the canonical design features of the Sp skeletogenic GRN that we examined are either missing or operate differently in Et. This work directly implies a dramatic reorganization of genomic regulatory circuitry concomitant with the divergence of the euechinoids, which began before the end-Permian extinction.