Transfer of a large gene regulatory apparatus to a new developmental address in echinoid evolution

Antero-posterior patterning in the brittle star Amphipholis squamata and the evolution of echinoderm body plans

Although the adult pentaradial body plan of echinoderms evolved from a bilateral ancestor, identifying axial homologies between the morphologically divergent echinoderms and their bilaterian relatives has been an enduring problem in zoology. The expression of conserved bilaterian patterning genes in echinoderms provides a molecular framework for resolving this puzzle. Recent studies in juvenile asteroids suggest that the bilaterian antero-posterior axis maps onto the medio-lateral axis of the arms, perpendicular to the proximo-distal axis of each of the five rays of the pentaradial body plan. Here, we test this hypothesis in another echinoderm class, the ophiuroids, using the cosmopolitan brittle star Amphipholis squamata. Our results show that the general principles of axial patterning are similar to those described in asteroids, and comparisons with existing molecular data from other echinoderm taxa support the idea that medio-lateral deployment of the bilaterian AP patterning program across the rays predates the evolution of the asterozoans, and likely the echinoderm crown-group. Our data also reveal expression differences between A. squamata and asteroids, which we attribute to secondary modifications specific to ophiuroids. Together, this work provides important comparative data to reconstruct the evolution of axial properties in echinoderm body plans.

Co-option of the trichome-forming network initiated the evolution of a morphological novelty in Drosophila eugracilis

Identifying the molecular origins by which new morphological structures evolve is one of the long-standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through the partial co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR-Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, misexpression of Shavenbaby in the sheath of D. melanogaster, a naive species that lacks these projections, is sufficient to induce small trichomes. These induced projections rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its partial co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionary novelty, our work shows that the trichome-forming network is flexible enough that it can be partially co-opted in a new context and subsequently refined to produce unique apical projections that are barely recognizable compared with their simpler ancestral beginnings.

Contrasting the development of larval and adult body plans during the evolution of biphasic lifecycles in sea urchins

Biphasic lifecycles are widespread among animals, but little is known about how the developmental transition between larvae and adults is regulated. Sea urchins are a unique system for studying this phenomenon because of the stark differences between their bilateral larval and pentaradial adult body plans. Here, we use single-cell RNA sequencing to analyze the development of Heliocidaris erythrogramma (He), a sea urchin species with an accelerated, non-feeding mode of larval development. The sequencing time course extends from embryogenesis to roughly a day before the onset of metamorphosis in He larvae, which is a period that has not been covered by previous datasets. We find that the non-feeding developmental strategy of He is associated with several changes in the specification of larval cell types compared to sea urchins with feeding larvae, such as the loss of a larva-specific skeletal cell population. Furthermore, the development of the larval and adult body plans in sea urchins may utilize largely different sets of regulatory genes. These findings lay the groundwork for extending existing developmental gene regulatory networks to cover additional stages of biphasic lifecycles.

Localization and origins of juvenile skeletogenic cells in the sea urchin Lytechinuspictus

The development of the sea urchin larval body plan is well understood from extensive studies of embryonic patterning. However, fewer studies have investigated the late larval stages during which the unique pentaradial adult body plan develops. Previous work on late larval development highlights major tissue changes leading up to metamorphosis, but the location of specific cell types during juvenile development is less understood. Here, we improve on technical limitations by applying highly sensitive hybridization chain reaction fluorescent in situ hybridization (HCR-FISH) to the fast-developing and transparent sea urchin Lytechinus pictus, with a focus on skeletogenic cells. First, we show that HCR-FISH can be used in L. pictus to precisely localize skeletogenic cells in the rudiment. In doing so, we provide a detailed staging scheme for the appearance of skeletogenic cells around the rudiment prior to and during biomineralization and show that many skeletogenic cells unassociated with larval rods localize outside of the rudiment prior to localizing inside. Second, we show that downstream biomineralization genes have similar expression patterns during larval and juvenile skeletogenesis, suggesting some conservation of skeletogenic mechanisms during development between stages. Third, we find co-expression of blastocoelar and skeletogenic cell markers around juvenile skeleton located outside of the rudiment, which is consistent with data showing that cells from the non-skeletogenic mesoderm embryonic lineage contribute to the juvenile skeletogenic cell lineage. This work sets the foundation for subsequent studies of other cell types in the late larva of L. pictus to better understand juvenile body plan development, patterning, and evolution.

Shared regulatory function of non-genomic thyroid hormone signaling in echinoderm skeletogenesis

Thyroid hormones are crucial regulators of metamorphosis and development in bilaterians, particularly in chordate deuterostomes. Recent evidence suggests a role for thyroid hormone signaling, principally via 3,5,3',5'-Tetraiodo-L-thyronine (T4), in the regulation of metamorphosis, programmed cell death and skeletogenesis in echinoids (sea urchins and sand dollars) and sea stars. Here, we test whether TH signaling in skeletogenesis is a shared trait of Echinozoa (Echinoida and Holothouroida) and Asterozoa (Ophiourida and Asteroida). We demonstrate dramatic acceleration of skeletogenesis after TH treatment in three classes of echinoderms: sea urchins, sea stars, and brittle stars (echinoids, asteroids, and ophiuroids). Fluorescently labeled thyroid hormone analogues reveal thyroid hormone binding to cells proximal to regions of skeletogenesis in the gut and juvenile rudiment. We also identify, for the first time, a potential source of thyroxine during gastrulation in sea urchin embryos. Thyroxine-positive cells are present in tip of the archenteron. In addition, we detect thyroid hormone binding to the cell membrane and nucleus during metamorphic development in echinoderms. Immunohistochemistry of phosphorylated MAPK in the presence and absence of TH-binding inhibitors suggests that THs may act via phosphorylation of MAPK (ERK1/2) to accelerate initiation of skeletogenesis in the three echinoderm groups. Together, these results indicate that TH regulation of mesenchyme cell activity via integrin-mediated MAPK signaling may be a conserved mechanism for the regulation of skeletogenesis in echinoderm development. In addition, TH action via a nuclear thyroid hormone receptor may regulate metamorphic development. Our findings shed light on potentially ancient pathways of thyroid hormone activity in echinoids, ophiuroids, and asteroids, or on a signaling system that has been repeatedly co-opted to coordinate metamorphic development in bilaterians.

Gene regulatory network co-option is sufficient to induce a morphological novelty in Drosophila

Identifying the molecular origins by which new morphological structures evolve is one of the long standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR/Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, mis-expression of Shavenbaby in the sheath of D. melanogaster , a naïve species that lacks these extensions, is sufficient to induce small trichomes. These induced extensions rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionarily novelty, our work shows that the trichome-forming network is flexible enough that it can be co-opted in a new context, and subsequently refined to produce unique apical projections that are barely recognizable compared to their simpler ancestral beginnings.

Detection of TALEN-mediated genome cleavage during the early embryonic stage of the starfish Patiria pectinifera

BACKGROUND

Echinoderms have long been utilized as experimental materials to study the genetic control of developmental processes and their evolution. Among echinoderms, the molecular study of starfish embryos has received considerable attention across research topics such as gene regulatory network evolution and larval regeneration. Recently, experimental techniques to manipulate gene functions have been gradually established in starfish as the feasibility of genome editing methods was reported. However, it is still unclear when these techniques cause genome cleavage during the development of starfish, which is critical to understand the timeframe and applicability of the experiment during early development of starfish.

RESULTS

We herein reported that gene functions can be analyzed by the genome editing method TALEN in early embryos, such as the blastula of the starfish Patiria pectinifera. We injected the mRNA of TALEN targeting rar, which was previously constructed, into eggs of P. pectinifera and examined the efficiency of genome cleavage through developmental stages from 6 to 48 hours post fertilization.

CONCLUSION

The results will be key knowledge not only when designing TALEN-based experiments but also when assessing the results.

Tracing the history of cell types

A study of sea urchin and sea star larvae paves the way for understanding how cell types evolve and give rise to novel morphologies.

Thyroid hormone membrane receptor binding and transcriptional regulation in the sea urchin Strongylocentrotus purpuratus

Thyroid hormones (THs) are small amino acid derived signaling molecules with broad physiological and developmental functions in animals. Specifically, their function in metamorphic development, ion regulation, angiogenesis and many others have been studied in detail in mammals and some other vertebrates. Despite extensive reports showing pharmacological responses of invertebrate species to THs, little is known about TH signaling mechanisms outside of vertebrates. Previous work in sea urchins suggests that non-genomic mechanisms are activated by TH ligands. Here we show that several THs bind to sea urchin (Strongylocentrotus purpuratus) cell membrane extracts and are displaced by ligands of RGD-binding integrins. A transcriptional analysis across sea urchin developmental stages shows activation of genomic and non-genomic pathways in response to TH exposure, suggesting that both pathways are activated by THs in sea urchin embryos and larvae. We also provide evidence associating TH regulation of gene expression with TH response elements in the genome. In ontogeny, we found more differentially expressed genes in older larvae compared to gastrula stages. In contrast to gastrula stages, the acceleration of skeletogenesis by thyroxine in older larvae is not fully inhibited by competitive ligands or inhibitors of the integrin membrane receptor pathway, suggesting that THs likely activate multiple pathways. Our data confirms a signaling function of THs in sea urchin development and suggests that both genomic and non-genomic mechanisms play a role, with genomic signaling being more prominent during later stages of larval development.

Early expression onset of tissue-specific effector genes during the specification process in sea urchin embryos

BACKGROUND

In the course of animal developmental processes, various tissues are differentiated through complex interactions within the gene regulatory network. As a general concept, differentiation has been considered to be the endpoint of specification processes. Previous works followed this view and provided a genetic control scheme of differentiation in sea urchin embryos: early specification genes generate distinct regulatory territories in an embryo to express a small set of differentiation driver genes; these genes eventually stimulate the expression of tissue-specific effector genes, which provide biological identity to differentiated cells, in each region. However, some tissue-specific effector genes begin to be expressed in parallel with the expression onset of early specification genes, raising questions about the simplistic regulatory scheme of tissue-specific effector gene expression and the current concept of differentiation itself.

RESULTS

Here, we examined the dynamics of effector gene expression patterns during sea urchin embryogenesis. Our transcriptome-based analysis indicated that many tissue-specific effector genes begin to be expressed and accumulated along with the advancing specification GRN in the distinct cell lineages of embryos. Moreover, we found that the expression of some of the tissue-specific effector genes commences before cell lineage segregation occurs.

CONCLUSIONS

Based on this finding, we propose that the expression onset of tissue-specific effector genes is controlled more dynamically than suggested in the previously proposed simplistic regulation scheme. Thus, we suggest that differentiation should be conceptualized as a seamless process of accumulation of effector expression along with the advancing specification GRN. This pattern of effector gene expression may have interesting implications for the evolution of novel cell types.

An optimized Tet-On system for conditional control of gene expression in sea urchins

Sea urchins and other echinoderms are important experimental models for studying developmental processes. The lack of approaches for conditional gene perturbation, however, has made it challenging to investigate the late developmental functions of genes that have essential roles during early embryogenesis and genes that have diverse functions in multiple tissues. The doxycycline-controlled Tet-On system is a widely used molecular tool for temporally and spatially regulated transgene expression. Here, we optimized the Tet-On system to conditionally induce gene expression in sea urchin embryos. Using this approach, we explored the roles the MAPK signaling plays in skeletogenesis by expressing genes that perturb the pathway specifically in primary mesenchyme cells during later stages of development. We demonstrated the wide utility of the Tet-On system by applying it to a second sea urchin species and in cell types other than the primary mesenchyme cells. Our work provides a robust and flexible platform for the spatiotemporal regulation of gene expression in sea urchins, which will considerably enhance the utility of this prominent model system.

Thyroid hormone-induced cell death in sea urchin metamorphic development

Thyroid hormones (THs) are important regulators of development, metabolism and homeostasis in metazoans. Specifically, they have been shown to regulate the metamorphic transitions of vertebrates and invertebrates alike. Indirectly developing sea urchin larvae accelerate the formation of juvenile structures in response to thyroxine (T4) treatment, while reducing their larval arm length. The mechanisms underlying larval arm reduction are unknown and we hypothesized that programmed cell death (PCD) is linked to this process. To test this hypothesis, we measured larval arm retraction in response to different THs (T4, T3, rT3, Tetrac) and assessed cell death in larvae using three different methods (TUNEL, YO-PRO-1 and caspase-3 activity) in the sea urchin Strongylocentrotus purpuratus. We also compared the extent of PCD in response to TH treatment before and after the invagination of the larval ectoderm, which marks the initiation of juvenile development in larval sea urchin species. We found that T4 treatment results in the strongest reduction of larval arms but detected a significant increase of PCD in response to T4, T3 and Tetrac in post-ingression but not pre-ingression larvae. As post-ingression larvae have initiated metamorphic development and therefore allocate resources to both larval and the juvenile structures, these results provide evidence that THs regulate larval development differentially via PCD. PCD in combination with cell proliferation likely has a key function in sea urchin development.

Decoding the IGF1 signaling gene regulatory network behind alveologenesis from a mouse model of bronchopulmonary dysplasia

Lung development is precisely controlled by underlying gene regulatory networks (GRN). Disruption of genes in the network can interrupt normal development and cause diseases such as bronchopulmonary dysplasia (BPD) - a chronic lung disease in preterm infants with morbid and sometimes lethal consequences characterized by lung immaturity and reduced alveolarization. Here, we generated a transgenic mouse exhibiting a moderate severity BPD phenotype by blocking IGF1 signaling in secondary crest myofibroblasts (SCMF) at the onset of alveologenesis. Using approaches mirroring the construction of the model GRN in sea urchin's development, we constructed the IGF1 signaling network underlying alveologenesis using this mouse model that phenocopies BPD. The constructed GRN, consisting of 43 genes, provides a bird's eye view of how the genes downstream of IGF1 are regulatorily connected. The GRN also reveals a mechanistic interpretation of how the effects of IGF1 signaling are transduced within SCMF from its specification genes to its effector genes and then from SCMF to its neighboring alveolar epithelial cells with WNT5A and FGF10 signaling as the bridge. Consistently, blocking WNT5A signaling in mice phenocopies BPD as inferred by the network. A comparative study on human samples suggests that a GRN of similar components and wiring underlies human BPD. Our network view of alveologenesis is transforming our perspective to understand and treat BPD. This new perspective calls for the construction of the full signaling GRN underlying alveologenesis, upon which targeted therapies for this neonatal chronic lung disease can be viably developed.

Distinct regulatory states control the elongation of individual skeletal rods in the sea urchin embryo

Background

Understanding how gene regulatory networks (GRNs) control developmental progression is a key to the mechanistic understanding of morphogenesis. The sea urchin larval skeletogenesis provides an excellent platform to tackle this question. In the early stages of sea urchin skeletogenesis, skeletogenic genes are uniformly expressed in the skeletogenic lineage. Yet, during skeletal elongation, skeletogenic genes are expressed in distinct spatial sub‐domains. The regulation of differential gene expression during late skeletogenesis is not well understood.

Results

Here we reveal the dynamic expression of the skeletogenic regulatory genes that define a specific regulatory state for each pair of skeletal rods, in the sea urchin Paracentrotus lividus. The vascular endothelial growth factor (VEGF) signaling, essential for skeleton formation, specifically controls the migration of cells that form the postoral and distal anterolateral skeletogenic rods. VEGF signaling also controls the expression of regulatory genes in cells at the tips of the postoral rods, including the transcription factors Pitx1 and MyoD1. Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

Conclusions

Our study illuminates the fine‐tuning of the regulatory system during the transition from early to late skeletogenesis that gives rise to rod‐specific regulatory states.

The skeletogenic transcription factors form specific regulatory states in various skeletogenic sub‐populations.

Late VEGF signaling controls the regulatory states at the tips of the post‐oral and anterolateral skeletal rods.

VEGF signaling controls the expression of the transcription factors, MyoD1 and Pitx1.

Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

Cambrian edrioasteroid reveals new mechanism for secondary reduction of the skeleton in echinoderms

Echinoderms are characterized by a distinctive high-magnesium calcite endoskeleton as adults, but elements of this have been drastically reduced in some groups. Herein, we describe a new pentaradial echinoderm, Yorkicystis haefneri n. gen. n. sp., which provides, to our knowledge, the oldest evidence of secondary non-mineralization of the echinoderm skeleton. This material was collected from the Cambrian Kinzers Formation in York (Pennsylvania, USA) and is dated as ca 510 Ma. Detailed morphological observations demonstrate that the ambulacra (i.e. axial region) are composed of flooring and cover plates, but the rest of the body (i.e. extraxial region) is preserved as a dark film and lacks any evidence of skeletal plating. Moreover, X-ray fluorescence analysis reveals that the axial region is elevated in iron. Based on our morphological and chemical data and on taphonomic comparisons with other fossils from the Kinzers Formation, we infer that the axial region was originally calcified, while the extraxial region was non-mineralized. Phylogenetic analyses recover Yorkicystis as an edrioasteroid, indicating that this partial absence of skeleton resulted from a secondary reduction. We hypothesize that skeletal reduction resulted from lack of expression of the skeletogenic gene regulatory network in the extraxial body wall during development. Secondary reduction of the skeleton in Yorkicystis might have allowed for greater flexibility of the body wall.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor–CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell–cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Lessons from a transcription factor: Alx1 provides insights into gene regulatory networks, cellular reprogramming, and cell type evolution

The skeleton-forming cells of sea urchins and other echinoderms have been studied by developmental biologists as models of cell specification and morphogenesis for many decades. The gene regulatory network (GRN) deployed in the embryonic skeletogenic cells of euechinoid sea urchins is one of the best understood in any developing animal. Recent comparative studies have leveraged the information contained in this GRN, bringing renewed attention to the diverse patterns of skeletogenesis within the phylum and the evolutionary basis for this diversity. The homeodomain-containing transcription factor, Alx1, was originally shown to be a core component of the skeletogenic GRN of the sea urchin embryo. Alx1 has since been found to be key regulator of skeletal cell identity throughout the phylum. As such, Alx1 is currently serving as a lens through which multiple developmental processes are being investigated. These include not only GRN organization and evolution, but also cell reprogramming, cell type evolution, and the gene regulatory control of morphogenesis. This review summarizes our current state of knowledge concerning Alx1 and highlights the insights it is yielding into these important developmental and evolutionary processes.

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.

Gene regulation of adult skeletogenesis in starfish and modifications during gene network co-option

The larval skeleton of the echinoderm is believed to have been acquired through co-option of a pre-existing gene regulatory network (GRN); that is, the mechanism for adult skeleton formation in the echinoderm was deployed in early embryogenesis during echinoderm diversification. To explore the evolutionary changes that occurred during co-option, we examined the mechanism for adult skeletogenesis using the starfish Patiria pectinifera. Expression patterns of skeletogenesis-related genes (vegf, vegfr, ets1/2, erg, alx1, ca1, and clect) suggest that adult skeletogenic cells develop from the posterior coelom after the start of feeding. Treatment with inhibitors and gene knockout using transcription activator-like effector nucleases (TALENs) suggest that the feeding-nutrient sensing pathway activates Vegf signaling via target of rapamycin (TOR) activity, leading to the activation of skeletogenic regulatory genes in starfish. In the larval skeletogenesis of sea urchins, the homeobox gene pmar1 activates skeletogenic regulatory genes, but in starfish, localized expression of the pmar1-related genes phbA and phbB was not detected during the adult skeleton formation stage. Based on these data, we provide a model for the adult skeletogenic GRN in the echinoderm and propose that the upstream regulatory system changed from the feeding-TOR-Vegf pathway to a homeobox gene-system during co-option of the skeletogenic GRN.

Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms

Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.

Homology of process: developmental dynamics in comparative biology

Comparative biology builds up systematic knowledge of the diversity of life, across evolutionary lineages and levels of organization, starting with evidence from a sparse sample of model organisms. In developmental biology, a key obstacle to the growth of comparative approaches is that the concept of homology is not very well defined for levels of organization that are intermediate between individual genes and morphological characters. In this paper, we investigate what it means for ontogenetic processes to be homologous, focusing specifically on the examples of insect segmentation and vertebrate somitogenesis. These processes can be homologous without homology of the underlying genes or gene networks, since the latter can diverge over evolutionary time, while the dynamics of the process remain the same. Ontogenetic processes like these therefore constitute a dissociable level and distinctive unit of comparison requiring their own specific criteria of homology. In addition, such processes are typically complex and nonlinear, such that their rigorous description and comparison requires not only observation and experimentation, but also dynamical modelling. We propose six criteria of process homology, combining recognized indicators (sameness of parts, morphological outcome and topological position) with novel ones derived from dynamical systems modelling (sameness of dynamical properties, dynamical complexity and evidence for transitional forms). We show how these criteria apply to animal segmentation and other ontogenetic processes. We conclude by situating our proposed dynamical framework for homology of process in relation to similar research programmes, such as process structuralism and developmental approaches to morphological homology.

FGF signalling plays similar roles in development and regeneration of the skeleton in the brittle star Amphiura filiformis

Regeneration as an adult developmental process is in many aspects similar to embryonic development. Although many studies point out similarities and differences, no large-scale, direct and functional comparative analyses between development and regeneration of a specific cell type or structure in one animal exist. Here, we use the brittle star Amphiura filiformis to characterise the role of the FGF signalling pathway during skeletal development in embryos and arm regeneration. In both processes, we find ligands expressed in ectodermal cells that flank underlying skeletal mesenchymal cells, which express the receptors. Perturbation of FGF signalling showed inhibited skeleton formation in both embryogenesis and regeneration, without affecting other key developmental processes. Differential transcriptome analysis finds mostly differentiation genes rather than transcription factors to be downregulated in both contexts. Moreover, comparative gene analysis allowed us to discover brittle star-specific differentiation genes. In conclusion, our results show that the FGF pathway is crucial for skeletogenesis in the brittle star, as in other deuterostomes, and provide evidence for the re-deployment of a developmental gene regulatory module during regeneration.

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01787-9.

Transcription Factors of the Alx Family: Evolutionarily Conserved Regulators of Deuterostome Skeletogenesis

Members of the alx gene family encode transcription factors that contain a highly conserved Paired-class, DNA-binding homeodomain, and a C-terminal OAR/Aristaless domain. Phylogenetic and comparative genomic studies have revealed complex patterns of alx gene duplications during deuterostome evolution. Remarkably, alx genes have been implicated in skeletogenesis in both echinoderms and vertebrates. In this review, we provide an overview of current knowledge concerning alx genes in deuterostomes. We highlight their evolutionarily conserved role in skeletogenesis and draw parallels and distinctions between the skeletogenic gene regulatory circuitries of diverse groups within the superphylum.

The function of architecture and logic in developmental gene regulatory networks

An important contribution of systems biology is the insight that biological systems depend on the function of molecular interactions and not just on individual molecules. System level mechanisms are particularly important in the development of animals and plants which depends not just on transcription factors and signaling molecules, but also on regulatory circuits and gene regulatory networks (GRNs). However, since GRNs consist of transcription factors, it can be challenging to assess the function of regulatory circuits independently of the function of regulatory factors. The comparison of different GRNs offers a way to do so and leads to several observations. First, similar regulatory circuits operate in various developmental contexts and in different species, and frequently, these circuits are associated with similar developmental functions. Second, given regulatory circuits are often used at particular positions within the GRN hierarchy. Third, in some GRNs, regulatory circuits are organized in a particular order in respect to each other. And fourth, the evolution of GRNs occurs not just by co-option of regulatory genes but also by rewiring of regulatory linkages between conserved regulatory genes, indicating that the organization of interactions is important. Thus, even though in most instances the function of regulatory circuits remains to be discovered, it becomes evident that the architecture and logic of GRNs are functionally important for the control of genome activity and for the specification of the body plan.

Transcriptome analysis reveals wingless regulates neural development and signaling genes in the region of wing pigmentation of a polka-dotted fruit fly

How evolutionary novelties have arisen is one of the central questions in evolutionary biology. Preexisting gene regulatory networks or signaling pathways have been shown to be co-opted for building novel traits in several organisms. However, the structure of entire gene regulatory networks and evolutionary events of gene co-option for emergence of a novel trait are poorly understood. In this study, to explore the genetic and molecular bases of the novel wing pigmentation pattern of a polka-dotted fruit fly (Drosophila guttifera), we performed de novo genome sequencing and transcriptome analyses. As a result, we comprehensively identified the genes associated with the pigmentation pattern. Furthermore, we revealed that 151 of these associated genes were positively or negatively regulated by wingless, a master regulator of wing pigmentation. Genes for neural development, Wnt signaling, Dpp signaling, and effectors (such as enzymes) for melanin pigmentation were included among these 151 genes. None of the known regulatory genes that regulate pigmentation pattern formation in other fruit fly species were included. Our results suggest that the novel pigmentation pattern of a polka-dotted fruit fly might have emerged through multistep co-options of multiple gene regulatory networks, signaling pathways, and effector genes, rather than recruitment of one large gene circuit.

On the specificity of gene regulatory networks: How does network co-option affect subsequent evolution?

The process of multicellular organismal development hinges upon the specificity of developmental programs: for different parts of the organism to form unique features, processes must exist to specify each part. This specificity is thought to be hardwired into gene regulatory networks, which activate cohorts of genes in particular tissues at particular times during development. However, the evolution of gene regulatory networks sometimes occurs by mechanisms that sacrifice specificity. One such mechanism is network co-option, in which existing gene networks are redeployed in new developmental contexts. While network co-option may offer an efficient mechanism for generating novel phenotypes, losses of tissue specificity at redeployed network genes could restrict the ability of the affected traits to evolve independently. At present, there has not been a detailed discussion regarding how tissue specificity of network genes might be altered due to gene network co-option at its initiation, as well as how trait independence can be retained or restored after network co-option. A lack of clarity about network co-option makes it more difficult to speculate on the long-term evolutionary implications of this mechanism. In this review, we will discuss the possible initial outcomes of network co-option, outline the mechanisms by which networks may retain or subsequently regain specificity after network co-option, and comment on some of the possible evolutionary consequences of network co-option. We place special emphasis on the need to consider selectively-neutral outcomes of network co-option to improve our understanding of the role of this mechanism in trait evolution.

Evolution of wing pigmentation in Drosophila: Diversity, physiological regulation, and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo‐devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero‐distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis‐regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

The evolution of a new cell type was associated with competition for a signaling ligand

There is presently a very limited understanding of the mechanisms that underlie the evolution of new cell types. The skeleton-forming primary mesenchyme cells (PMCs) of euechinoid sea urchins, derived from the micromeres of the 16-cell embryo, are an example of a recently evolved cell type. All adult echinoderms have a calcite-based endoskeleton, a synapomorphy of the Ambulacraria. Only euechinoids have a micromere-PMC lineage, however, which evolved through the co-option of the adult skeletogenic program into the embryo. During normal development, PMCs alone secrete the embryonic skeleton. Other mesoderm cells, known as blastocoelar cells (BCs), have the potential to produce a skeleton, but a PMC-derived signal ordinarily prevents these cells from expressing a skeletogenic fate and directs them into an alternative developmental pathway. Recently, it was shown that vascular endothelial growth factor (VEGF) signaling plays an important role in PMC differentiation and is part of a conserved program of skeletogenesis among echinoderms. Here, we report that VEGF signaling, acting through ectoderm-derived VEGF3 and its cognate receptor, VEGF receptor (VEGFR)-10-Ig, is also essential for the deployment of the skeletogenic program in BCs. This VEGF-dependent program includes the activation of aristaless-like homeobox 1 (alx1), a conserved transcriptional regulator of skeletogenic specification across echinoderms and an example of a “terminal selector” gene that controls cell identity. We show that PMCs control BC fate by sequestering VEGF3, thereby preventing activation of alx1 and the downstream skeletogenic network in BCs. Our findings provide an example of the regulation of early embryonic cell fates by direct competition for a secreted signaling ligand, a developmental mechanism that has not been widely recognized. Moreover, they reveal that a novel cell type evolved by outcompeting other embryonic cell lineages for an essential signaling ligand that regulates the expression of a gene controlling cell identity.

How do new cell types evolve? This study shows that mesoderm cells in sea urchin embryos diversified, at least in part, through a heterochronic shift in the expression of a key transcription factor, which led to competition for a signaling ligand and subsequent gene regulatory independence of the two cell types.

Introducing Trait Networks to Elucidate the Fluidity of Organismal Evolution Using Palaeontological Data

Explaining the evolution of animals requires ecological, developmental, paleontological, and phylogenetic considerations because organismal traits are affected by complex evolutionary processes. Modeling a plurality of processes, operating at distinct time-scales on potentially interdependent traits, can benefit from approaches that are complementary treatments to phylogenetics. Here, we developed an inclusive network approach, implemented in the command line software ComponentGrapher, and analyzed trait co-occurrence of rhinocerotoid mammals. We identified stable, unstable, and pivotal traits, as well as traits contributing to complexes, that may follow to a common developmental regulation, that point to an early implementation of the postcranial Bauplan among rhinocerotoids. Strikingly, most identified traits are highly dissociable, used repeatedly in distinct combinations and in different taxa, which usually do not form clades. Therefore, the genes encoding these traits are likely recruited into novel gene regulation networks during the course of evolution. Our evo-systemic framework, generalizable to other evolved organizations, supports a pluralistic modeling of organismal evolution, including trees and networks.

Genome-wide identification of binding sites and gene targets of Alx1, a pivotal regulator of echinoderm skeletogenesis

Alx1 is a conserved regulator of skeletogenesis in echinoderms and evolutionary changes in Alx1 sequence and expression have played a pivotal role in modifying programs of skeletogenesis within the phylum. Alx1 regulates a large suite of effector genes that control the morphogenetic behaviors and biomineral-forming activities of skeletogenic cells. To better understand the gene regulatory control of skeletogenesis by Alx1, we used genome-wide ChIP-seq to identify Alx1-binding sites and direct gene targets. Our analysis revealed that many terminal differentiation genes receive direct transcriptional inputs from Alx1. In addition, we found that intermediate transcription factors previously shown to be downstream of Alx1 all receive direct inputs from Alx1. Thus, Alx1 appears to regulate effector genes by indirect, as well as direct, mechanisms. We tested 23 high-confidence ChIP-seq peaks using GFP reporters and identified 18 active cis-regulatory modules (CRMs); this represents a high success rate for CRM discovery. Detailed analysis of a representative CRM confirmed that a conserved, palindromic Alx1-binding site was essential for expression. Our work significantly advances our understanding of the gene regulatory circuitry that controls skeletogenesis in sea urchins and provides a framework for evolutionary studies.

Pigeon foot feathering reveals conserved limb identity networks

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

Modularity, criticality, and evolvability of a developmental gene regulatory network

The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.

Possible cooption of a VEGF-driven tubulogenesis program for biomineralization in echinoderms

Biomineralization is the process by which living organisms use minerals to form hard structures that protect and support them. Biomineralization is believed to have evolved rapidly and independently in different phyla utilizing preexisting components. The mechanistic understanding of the regulatory networks that drive biomineralization and their evolution is far from clear. Sea urchin skeletogenesis is an excellent model system for studying both gene regulation and mineral uptake and deposition. The sea urchin calcite spicules are formed within a tubular cavity generated by the skeletogenic cells controlled by vascular endothelial growth factor (VEGF) signaling. The VEGF pathway is essential for biomineralization in echinoderms, while in many other phyla, across metazoans, it controls tubulogenesis and vascularization. Despite the critical role of VEGF signaling in sea urchin spiculogenesis, the downstream program it activates was largely unknown. Here we study the cellular and molecular machinery activated by the VEGF pathway during sea urchin spiculogenesis and reveal multiple parallels to the regulation of vertebrate vascularization. Human VEGF rescues sea urchin VEGF knockdown, vesicle deposition into an internal cavity plays a significant role in both systems, and sea urchin VEGF signaling activates hundreds of genes, including biomineralization and interestingly, vascularization genes. Moreover, five upstream transcription factors and three signaling genes that drive spiculogenesis are homologous to vertebrate factors that control vascularization. Overall, our findings suggest that sea urchin spiculogenesis and vertebrate vascularization diverged from a common ancestral tubulogenesis program, broadly adapted for vascularization and specifically coopted for biomineralization in the echinoderm phylum.

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.

A new ophiocistioid with soft-tissue preservation from the Silurian Herefordshire Lagerstätte, and the evolution of the holothurian body plan

Reconstructing the evolutionary assembly of animal body plans is challenging when there are large morphological gaps between extant sister taxa, as in the case of echinozoans (echinoids and holothurians). However, the inclusion of extinct taxa can help bridge these gaps. Here we describe a new species of echinozoan, Sollasina cthulhu, from the Silurian Herefordshire Lagerstätte, UK. Sollasina cthulhu belongs to the ophiocistioids, an extinct group that shares characters with both echinoids and holothurians. Using physical-optical tomography and computer reconstruction, we visualize the internal anatomy of S. cthulhu in three dimensions, revealing inner soft tissues that we interpret as the ring canal, a key part of the water vascular system that was previously unknown in fossil echinozoans. Phylogenetic analyses strongly suggest that Sollasina and other ophiocistioids represent a paraphyletic group of stem holothurians, as previously hypothesized. This allows us to reconstruct the stepwise reduction of the skeleton during the assembly of the holothurian body plan, which may have been controlled by changes in the expression of biomineralization genes.

Culture of and experiments with sea urchin embryo primary mesenchyme cells

Skeletogenesis in the sea urchin embryo gives rise to a pair of intricate endoskeletal spicules. Deposition of these skeletal elements in the early larva is the outcome of a morphogenetic program that begins with maternal inputs in the early zygote and results in the specification of the large micromere-primary mesenchyme cell (PMC) lineage. PMCs are of considerable interest as a model system, not only to dissect the mechanism of specific developmental processes, but also to investigate their evolution and the unrivaled level of control over the formation of a graded, mechanically robust, yet single crystalline biomineral. The ability to study gene regulatory circuits, cellular behavior, signaling pathways, and molecular players involved in biomineralization is significantly boosted by the high level of autonomy of PMCs. In fact, in the presence of horse serum, micromeres differentiate into PMCs and produce spicules in vitro, separated from the embryonic milieu. PMC culture eliminates indirect effects that can complicate the interpretation of experiments in vivo, offers superior spatiotemporal control, enables PMC-specific readouts, and is compatible with most imaging and characterization techniques. In this chapter, we provide an updated protocol, based on the pioneering work by Okazaki and Wilt, for the isolation of micromeres and subsequent culture of PMCs, as well as protocols for fixation and staining for fluorescent microscopy, preparation of cell cultures for electron microscopy, and the isolation of RNA.

Developmental effector gene regulation: Multiplexed strategies for functional analysis

The staggering complexity of the genome controls for developmental processes is revealed through massively parallel cis-regulatory analysis using new methods of perturbation and readout. The choice of combinations of these new methods is tailored to the system, question and resources at hand. Our focus is on issues that include the necessity or sufficiency of given cis-regulatory modules, cis-regulatory function in the normal spatial genomic context, and easily accessible high throughput and multiplexed analysis methods. In the sea urchin embryonic model, recombineered BACs offer new opportunities for consecutive modes of cis-regulatory analyses that answer these requirements, as we here demonstrate on a diverse suite of previously unstudied sea urchin effector genes expressed in skeletogenic cells. Positively active cis-regulatory modules were located in single Nanostring experiments per BAC containing the gene of interest, by application of our previously reported "barcode" tag vectors of which> 100 can be analyzed at one time. Computational analysis of DNA sequences that drive expression, based on the known skeletogenic regulatory state, then permitted effective identification of functional target site clusters. Deletion of these sub-regions from the parent BACs revealed module necessity, as simultaneous tests of the same regions in short constructs revealed sufficiency. Predicted functional inputs were then confirmed by site mutations, all generated and tested in multiplex formats. There emerged the simple conclusion that each effector gene utilizes a small subset of inputs from the skeletogenic GRN. These inputs may function to only adjust expression levels or in some cases necessary for expression. Since we know the GRN architecture upstream of the effector genes, we could then conceptually isolate and compare the wiring of the effector gene driver sub-circuits and identify the inputs whose removal abolish expression.

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.

From genome to anatomy: The architecture and evolution of the skeletogenic gene regulatory network of sea urchins and other echinoderms

The skeletogenic gene regulatory network (GRN) of sea urchins and other echinoderms is one of the most intensively studied transcriptional networks in any developing organism. As such, it serves as a preeminent model of GRN architecture and evolution. This review summarizes our current understanding of this developmental network. We describe in detail the most comprehensive model of the skeletogenic GRN, one developed for the euechinoid sea urchin Strongylocentrotus purpuratus, including its initial deployment by maternal inputs, its elaboration and stabilization through regulatory gene interactions, and its control of downstream effector genes that directly drive skeletal morphogenesis. We highlight recent comparative studies that have leveraged the euechinoid GRN model to examine the evolution of skeletogenic programs in diverse echinoderms, studies that have revealed both conserved and divergent features of skeletogenesis within the phylum. Last, we summarize the major insights that have emerged from analysis of the structure and evolution of the echinoderm skeletogenic GRN and identify key, unresolved questions as a guide for future work.

Developmental transcriptomics of the brittle star Amphiura filiformis reveals gene regulatory network rewiring in echinoderm larval skeleton evolution

Background

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

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1402-8) contains supplementary material, which is available to authorized users.

Thyroid Hormones Accelerate Initiation of Skeletogenesis via MAPK (ERK1/2) in Larval Sea Urchins (Strongylocentrotus purpuratus)

Thyroid hormones are important regulators of development and metabolism in animals. Their function via genomic and non-genomic actions is well-established in vertebrate species but remains largely elusive among invertebrates. Previous work suggests that thyroid hormones, principally 3,5,3',5'-Tetraiodo-L-thyronine (T4), regulate development to metamorphosis in sea urchins. Here we show that thyroid hormones, including T4, 3,5,3'-triiodo-l-thyronine (T3), and 3,5-Diiodothyronine (T2) accelerate initiation of skeletogenesis in sea urchin gastrulae and pluteus larvae of the sea urchin Strongylocentrotus purpuratus, as measured by skeletal spicule formation. Fluorescently conjugated hormones show T4 binding to primary mesenchyme cells in sea urchin gastrulae. Furthermore, our investigation of TH mediated skeletogenesis shows that Ets1, a transcription factor controlling initiation of skeletogenesis, is a target of activated (phosphorylated) mitogen-activated protein kinase [MAPK; extracellular signal-regulated kinase 1/2 (ERK1/2)]. As well, we show that PD98059, an inhibitor of ERK1/2 MAPK signaling, prevents the T4 mediated acceleration of skeletogenesis and upregulation of Ets1. In contrast, SB203580, an inhibitor of p38 MAPK signaling, did not inhibit the effect of T4. Immunohistochemistry revealed that T4 causes phosphorylation of ERK1/2 in presumptive primary mesenchyme cells and the basal membrane of epithelial cells in the gastrula. Pre-incubation of sea urchin gastrulae with RGD peptide, a competitive inhibitor of TH binding to integrins, inhibited the effect of T4 on skeletogenesis. Together, these experiments provide evidence that T4 acts via a MAPK- (ERK1/2) mediated integrin membrane receptor to accelerate skeletogenesis in sea urchin mesenchyme cells. These findings shed light, for the first time, on a putative non-genomic pathway of TH action in a non-chordate deuterostome and help elucidate the evolutionary history of TH signaling in animals.

Functional divergence of paralogous transcription factors supported the evolution of biomineralization in echinoderms

Alx1 is a pivotal transcription factor in a gene regulatory network that controls skeletogenesis throughout the echinoderm phylum. We performed a structure-function analysis of sea urchin Alx1 using a rescue assay and identified a novel, conserved motif (Domain 2) essential for skeletogenic function. The paralogue of Alx1, Alx4, was not functionally interchangeable with Alx1, but insertion of Domain 2 conferred robust skeletogenic function on Alx4. We used cross-species expression experiments to show that Alx1 proteins from distantly related echinoderms are not interchangeable, although the sequence and function of Domain 2 are highly conserved. We also found that Domain 2 is subject to alternative splicing and provide evidence that this domain was originally gained through exonization. Our findings show that a gene duplication event permitted the functional specialization of a transcription factor through changes in exon-intron organization and thereby supported the evolution of a major morphological novelty.

Wound healing, calcium signaling, and other novel pathways are associated with the formation of butterfly eyespots

Background

One hypothesis surrounding the origin of novel traits is that they originate from the co-option of pre-existing genes or larger gene regulatory networks into novel developmental contexts. Insights into a trait’s evolutionary origins can, thus, be gained via identification of the genes underlying trait development, and exploring whether those genes also function in other developmental contexts. Here we investigate the set of genes associated with the development of eyespot color patterns, a trait that originated once within the Nymphalid family of butterflies. Although several genes associated with eyespot development have been identified, the eyespot gene regulatory network remains largely unknown.

Results

In this study, next-generation sequencing and transcriptome analyses were used to identify a large set of genes associated with eyespot development of Bicyclus anynana butterflies, at 3-6 h after pupation, prior to the differentiation of the color rings. Eyespot-associated genes were identified by comparing the transcriptomes of homologous micro-dissected wing tissues that either develop or do not develop eyespots in wild-type and a mutant line of butterflies, Spotty, with extra eyespots. Overall, 186 genes were significantly up and down-regulated in wing tissues that develop eyespots compared to wing tissues that do not. Many of the differentially expressed genes have yet to be annotated. New signaling pathways, including the Toll, Fibroblast Growth Factor (FGF), extracellular signal–regulated kinase (ERK) and/or Jun N-terminal kinase (JNK) signaling pathways are associated for the first time with eyespot development. In addition, several genes involved in wound healing and calcium signaling were also found to be associated with eyespots.

Conclusions

Overall, this study provides the identity of many new genes and signaling pathways associated with eyespots, and suggests that the ancient wound healing gene regulatory network may have been co-opted to cells at the center of the pattern to aid in eyespot origins. New transcription factors that may be providing different identities to distinct wing sectors, and genes with sexually dimorphic expression in the eyespots were also identified.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4175-7) contains supplementary material, which is available to authorized users.

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.