The evolutionary mechanisms involved in shaping complex gene regulatory networks (GRN) that encode for morphologically similar structures in distantly related animals remain elusive. In this context, echinoderm larval skeletons found in brittle stars and sea urchins provide an ideal system. Here, we characterize for the first time the development of the larval skeleton in the ophiuroid Amphiura filiformis and compare it systematically with its counterpart in sea urchin. We show that ophiuroids and euechinoids, that split at least 480 Million years ago (Mya), have remarkable similarities in tempo and mode of skeletal development. Despite morphological and ontological similarities, our high-resolution study of the dynamics of genetic regulatory states in A. filiformis highlights numerous differences in the architecture of their underlying GRNs. Importantly, the A.filiformis pplx, the closest gene to the sea urchin double negative gate (DNG) repressor pmar1, fails to drive the skeletogenic program in sea urchin, showing important evolutionary differences in protein function. hesC, the second repressor of the DNG, is co-expressed with most of the genes that are repressed in sea urchin, indicating the absence of direct repression of tbr, ets1/2, and delta in A. filiformis. Furthermore, the absence of expression in later stages of brittle star skeleton development of key regulatory genes, such as foxb and dri, shows significantly different regulatory states. Our data fill up an important gap in the picture of larval mesoderm in echinoderms and allows us to explore the evolutionary implications relative to the recently established phylogeny of echinoderm classes. In light of recent studies on other echinoderms, our data highlight a high evolutionary plasticity of the same nodes throughout evolution of echinoderm skeletogenesis. Finally, gene duplication, protein function diversification, and cis-regulatory element evolution all contributed to shape the regulatory program for larval skeletogenesis in different branches of echinoderms.

Cardiovascular lineages develop together with kidney, smooth muscle, and limb connective tissue progenitors from the lateral plate mesoderm (LPM). How the LPM initially emerges and how its downstream fates are molecularly interconnected remain unknown. Here, we isolate a pan-LPM enhancer in the zebrafish-specific draculin (drl) gene that provides specific LPM reporter activity from early gastrulation. In toto live imaging and lineage tracing of drl-based reporters captures the dynamic LPM emergence as lineage-restricted mesendoderm field. The drl pan-LPM enhancer responds to the transcription factors EomesoderminA, FoxH1, and MixL1 that combined with Smad activity drive LPM emergence. We uncover specific activity of zebrafish-derived drl reporters in LPM-corresponding territories of several chordates including chicken, axolotl, lamprey, Ciona, and amphioxus, revealing a universal upstream LPM program. Altogether, our work provides a mechanistic framework for LPM emergence as defined progenitor field, possibly representing an ancient mesodermal cell state that predates the primordial vertebrate embryo.

The brittle star Amphiura filiformis, which regenerates its arms post autotomy, is emerging as a useful model for studying the molecular underpinnings of regeneration, aided by the recent availability of some molecular resources. During regeneration a blastema initially is formed distally to the amputation site, and then a rapid rebuild is obtained by adding metameric units, which will eventually differentiate and become fully functional. In this work we first characterize the developmental process of the regenerating arms using two differentiation markers for muscle and skeletal structures - Afi-trop-1 and Afi-αcoll. Both genes are not expressed in the blastema and newly added undifferentiated metameric units. Their expression at different regenerating stages shows an early segregation of muscle and skeletal cells during the regenerating process, long before the metameric units become functional. We then studied the expression of a set of genes orthologous of the sea urchin transcription factors involved in the development of skeletal and non-skeletal mesoderm: Afi-ets1/2, Afi-alx1, Afi-tbr, Afi-foxB and Afi-gataC. We found that Afi-ets1/2, Afi-alx1, Afi-foxB and Afi-gataC are all expressed at the blastemal stage. As regeneration progresses those genes are expressed in a similar small undifferentiated domain beneath the distal growth cap, while in more advanced metameric units they become restricted to different skeletal domains. Afi-foxB becomes expressed in non-skeletal structures. This suggests that they might play a combinatorial role only in the early cell specification process and that subsequently they function independently in the differentiation of different structures. Afi-tbr is not present in the adult arm tissue at any stage of regeneration. In situ hybridization results have been confirmed with a new strategy for quantitative PCR (QPCR), using a subdivision of the three stages of regeneration into proximal (differentiated) and distal (undifferentiated) arm segments.

Brittle stars regenerate their whole arms post-amputation. Amphiura filiformis can now be used for molecular characterization of arm regeneration due to the availability of transcriptomic data. Previous work showed that specific developmental transcription factors known to take part in echinoderm skeletogenesis are expressed during adult arm regeneration in A. filiformis; however, the process of skeleton formation remained poorly understood. Here, we present the results of an in-depth microscopic analysis of skeletal morphogenesis during regeneration, using calcein staining, EdU labeling and in situ hybridization. To better compare different samples, we propose a staging system for the early A. filiformis arm regeneration stages based on morphological landmarks identifiable in living animals and supported by histological analysis. We show that the calcified spicules forming the endoskeleton first appear very early during regeneration in the dermal layer of regenerates. These spicules then mature into complex skeletal elements of the differentiated arm during late regeneration. The mesenchymal cells in the dermal area express the skeletal marker genes Afi-c-lectin, Afi-p58b and Afi-p19; however, EdU labeling shows that these dermal cells do not proliferate. A. filiformis arms regenerate through a consistent set of developmental stages using a distalization-intercalation mode, despite variability in regeneration rate. Skeletal elements form in a mesenchymal cell layer that does not proliferate and thus must be supplied from a different source. Our work provides the basis for future cellular and molecular studies of skeleton regeneration in brittle stars.

Limb development in salamanders differs from other tetrapods in that the first digits to form are the two most anterior (preaxial dominance). This has been proposed as a salamander novelty and its mechanistic basis is unknown. Salamanders are the only adult tetrapods able to regenerate the limb, and the contribution of preaxial dominance to limb regeneration is unclear. Here we show that during early outgrowth of the limb bud, a small cohort of cells express the orphan geneProd1together withBmp2, a critical player in digit condensation in amniotes. Disruption ofProd1with a gene-editing nuclease abrogates these cells, and blocks formation of the radius and ulna, and outgrowth of the anterior digits. Preaxial dominance is a notable feature of limb regeneration in the larval newt, but this changes abruptly after metamorphosis so that the formation of anterior and posterior digits occurs together within the autopodium resembling an amniote-like pattern.

Amongst the echinoderms the class Ophiuroidea is of particular interest for its phylogenetic position, ecological importance and developmental and regenerative biology. However, compared to other echinoderms, notably echinoids (sea urchins), relatively little is known about developmental changes in gene expression in ophiuroids. To address this issue, we have generated and assembled a large RNAseq data set of four key stages of development in the brittle star Amphiura filiformis and a de novo reference transcriptome of comparable quality to that of a model echinoderm—the sea urchin Strongylocentrotus purpuratus. Furthermore, we provide access to the new data via a web interface: http://www.echinonet.eu/shiny/Amphiura_filiformis/ . We have identified highly conserved genes associated with the development of a biomineralised skeleton. We also identify important class-specific characters, including the independent duplication of the msp130 class of genes in different echinoderm classes and the unique occurrence of spicule matrix (sm) genes in echinoids. Using a new quantification pipeline for our de novo transcriptome, validated with other methodologies, we find major differences between brittle stars and sea urchins in the temporal expression of many transcription factor genes. This divergence in developmental regulatory states is more evident in early stages of development when cell specification begins, rather than when cells initiate differentiation. Our findings indicate that there has been a high degree of gene regulatory network rewiring and clade-specific gene duplication, supporting the hypothesis of a convergent evolution of larval skeleton development in echinoderms.