HpEts, an ets-related transcription factor implicated in primary mesenchyme cell differentiation in the sea urchin embryo

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.

CRISPR-Cas9-Mediated Gene Knockout in a Non-Model Sea Urchin, Heliocidaris crassispina

Sea urchins have been used as model organisms in developmental biology research and the genomes of several sea urchin species have been sequenced. Recently, genome editing technologies have become available for sea urchins, and methods for gene knockout using the CRISPRCas9 system have been established. Heliocidaris crassispina is an important marine fishery resource with edible gonads. Although H. crassispina has been used as a biological research material, its genome has not yet been published, and it is a non-model sea urchin for molecular biology research. However, as recent advances in genome editing technology have facilitated genome modification in non-model organisms, we applied genome editing using the CRISPR-Cas9 system to H. crassispina. In this study, we targeted genes encoding ETS transcription factor (HcEts) and pigmentation-related polyketide synthase (HcPks1). Gene fragments were isolated using primers designed by inter-specific sequence comparisons within Echinoidea. When Ets gene was targeted using two sgRNAs, one successfully introduced mutations and impaired skeletogenesis. In the Pks1 gene knockout, when two sgRNAs targeting the close vicinity of the site corresponding to the target site that showed 100% mutagenesis efficiency of the Pks1 gene in Hemicentrotus pulcherrimus, mutagenesis was not observed. However, two other sgRNAs targeting distant sites efficiently introduced mutations. In addition, Pks1 knockout H. crassispina exhibited an albino phenotype in the pluteus larvae and adult sea urchins after metamorphosis. This indicates that the CRISPRCas9 system can be used to modify the genome of the non-model sea urchin H. crassispina.

microRNA-1 regulates sea urchin skeletogenesis by directly targeting skeletogenic genes and modulating components of signaling pathways

microRNAs are evolutionarily conserved non-coding RNAs that direct post-transcriptional regulation of target transcripts. In vertebrates, microRNA-1 (miR-1) is expressed in muscle and has been found to play critical regulatory roles in vertebrate angiogenesis, a process that has been proposed to be analogous to sea urchin skeletogenesis. Results indicate that both miR-1 inhibitor and miR-1 mimic-injected larvae have significantly less F-actin enriched circumpharyngeal muscle fibers and fewer gut contractions. In addition, miR-1 regulates the positioning of skeletogenic primary mesenchyme cells (PMCs) and skeletogenesis of the sea urchin embryo. Interestingly, the gain-of-function of miR-1 leads to more severe PMC patterning and skeletal branching defects than its loss-of-function. The results suggest that miR-1 directly suppresses Ets1/2, Tbr, and VegfR7 of the skeletogenic gene regulatory network, and Nodal, and Wnt1 signaling components. This study identifies potential targets of miR-1 that impacts skeletogenesis and muscle formation and contributes to a deeper understanding of miR-1's function during development.

RAS-independent ERK activation by constitutively active KSR3 in non-chordate metazoa

During early development of the sea urchin embryo, activation of ERK signalling in mesodermal precursors is not triggered by extracellular RTK ligands but by a cell-autonomous, RAS-independent mechanism that was not understood. We discovered that in these cells, ERK signalling is activated through the transcriptional activation of a gene encoding a protein related to Kinase Suppressor of Ras, that we named KSR3. KSR3 belongs to a family of catalytically inactive allosteric activators of RAF. Phylogenetic analysis revealed that genes encoding kinase defective KSR3 proteins are present in most non-chordate metazoa but have been lost in flies and nematodes. We show that the structure of KSR3 factors resembles that of several oncogenic human RAF mutants and that KSR3 from echinoderms, cnidarians and hemichordates activate ERK signalling independently of RAS when overexpressed in cultured cells. Finally, we used the sequence of KSR3 factors to identify activating mutations of human B-RAF. These findings reveal key functions for this family of factors as activators of RAF in RAS-independent ERK signalling in invertebrates. They have implications on the evolution of the ERK signalling pathway and suggest a mechanism for its co-option in the course of evolution.

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.

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.

Global patterns of enhancer activity during sea urchin embryogenesis assessed by eRNA profiling

We used capped analysis of gene expression with sequencing (CAGE-seq) to profile eRNA expression and enhancer activity during embryogenesis of a model echinoderm: the sea urchin, Strongylocentrotus purpuratus. We identified more than 18,000 enhancers that were active in mature oocytes and developing embryos and documented a burst of enhancer activation during cleavage and early blastula stages. We found that a large fraction (73.8%) of all enhancers active during the first 48 h of embryogenesis were hyperaccessible no later than the 128-cell stage and possibly even earlier. Most enhancers were located near gene bodies, and temporal patterns of eRNA expression tended to parallel those of nearby genes. Furthermore, enhancers near lineage-specific genes contained signatures of inputs from developmental gene regulatory networks deployed in those lineages. A large fraction (60%) of sea urchin enhancers previously shown to be active in transgenic reporter assays was associated with eRNA expression. Moreover, a large fraction (50%) of a representative subset of enhancers identified by eRNA profiling drove tissue-specific gene expression in isolation when tested by reporter assays. Our findings provide an atlas of developmental enhancers in a model sea urchin and support the utility of eRNA profiling as a tool for enhancer discovery and regulatory biology. The data generated in this study are available at Echinobase, the public database of information related to echinoderm genomics.

The gene regulatory control of sea urchin gastrulation

The cell behaviors associated with gastrulation in sea urchins have been well described. More recently, considerable progress has been made in elucidating gene regulatory networks (GRNs) that underlie the specification of early embryonic territories in this experimental model. This review integrates information from these two avenues of work. I discuss the principal cell movements that take place during sea urchin gastrulation, with an emphasis on molecular effectors of the movements, and summarize our current understanding of the gene regulatory circuitry upstream of those effectors. A case is made that GRN biology can provide a causal explanation of gastrulation, although additional analysis is needed at several levels of biological organization in order to provide a deeper understanding of this complex morphogenetic process.

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.

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.

Expression of exogenous mRNAs to study gene function in echinoderm embryos

With the completion of the genome sequencing projects, a new challenge for developmental biologists is to assign a function to the thousands of genes identified. Expression of exogenous mRNAs is a powerful, versatile and rapid technique that can be used to study gene function during development of the sea urchin. This chapter describes how this technique can be used to analyze gene function in echinoderm embryos, how it can be combined with cell transplantation to perform mosaic analysis and how it can be applied to identify downstream targets genes of transcription factors and signaling pathways. We describe specific examples of the use of overexpression of mRNA to analyze gene function, mention the benefits and current limitations of the technique and emphasize the importance of using different controls to assess the specificity of the effects observed. Finally, this chapter details the different steps, vectors and protocols for in vitro production of mRNA and phenotypic analysis.

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.

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.

Recent advances in functional perturbation and genome editing techniques in studying sea urchin development

Studies on the gene regulatory networks (GRNs) of sea urchin embryos have provided a basic understanding of the molecular mechanisms controlling animal development. The causal links in GRNs have been verified experimentally through perturbation of gene functions. Microinjection of antisense morpholino oligonucleotides (MOs) into the egg is the most widely used approach for gene knockdown in sea urchin embryos. The modification of MOs into a membrane-permeable form (vivo-MOs) has allowed gene knockdown at later developmental stages. Recent advances in genome editing tools, such as zinc-finger nucleases, transcription activator-like effector-based nucleases and the clustered regularly interspaced short palindromic repeat/clustered regularly interspaced short palindromic repeat-associated protein 9 (CRISPR/Cas9) system, have provided methods for gene knockout in sea urchins. Here, we review the use of vivo-MOs and genome editing tools in sea urchin studies since the publication of its genome in 2006. Various applications of the CRISPR/Cas9 system and their potential in studying sea urchin development are also discussed. These new tools will provide more sophisticated experimental methods for studying sea urchin development.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

Developmental effects of the protein kinase inhibitor kenpaullone on the sea urchin embryo

The selection and validation of bioactive compounds require multiple approaches, including in-depth analyses of their biological activity in a whole-animal context. We exploited the sea urchin embryo in a rapid, medium-scale range screening to test the effects of the small synthetic kinase inhibitor kenpaullone. We show that sea urchin embryos specifically respond to this molecule depending on both dose and timing of administration. Phenotypic effects of kenpaullone are not immediately visible, since this molecule affects neither the fertilization nor the spatial arrangement of blastomeres at early developmental stages. Nevertheless, kenpaullone exposure from the beginning of embryogenesis profoundly perturbs specification, detachment from the epithelium, and migration of the primary mesenchyme cells, thus affecting the whole embryonic epithelial mesenchymal transition process. Our results reaffirm the sea urchin embryo as an excellent and sensitive in vivo system, which provides straightforward and rapid response to external stimuli.

Comparative studies on the skeletogenic mesenchyme of echinoids

Skeletogenic mesenchyme cells in echinoids are suitable for studying developmental mechanisms, and have been used extensively. Most of these studies have been performed on species in the order Camarodonta, which are modern echinoids (subclass Euechinoidea) and are considered "model" echinoid species. In contrast, species belonging to other orders are studied less frequently, especially investigations of their molecular developmental biology such as gene regulatory networks. Recent studies on mesenchyme development in non-camarodont species suggest that these species are potential sources of comparative information to elucidate the mechanisms underlying skeletogenic mesenchyme development. In this review, the importance of using comparative data to understand development and evolution is discussed.

Mechanisms of the epithelial-to-mesenchymal transition in sea urchin embryos

Sea urchin mesenchyme is composed of the large micromere-derived spiculogenetic primary mesenchyme cells (PMC), veg2-tier macromere-derived non-spiculogenetic mesenchyme cells, the small micromere-derived germ cells, and the macro- and mesomere-derived neuronal mesenchyme cells. They are formed through the epithelial-to-mesenchymal transition (EMT) and possess multipotency, except PMCs that solely differentiate larval spicules. The process of EMT is associated with modification of epithelial cell surface property that includes loss of affinity to the apical and basal extracellular matrices, inter-epithelial cell adherens junctions and epithelial cell surface-specific proteins. These cell surface structures and molecules are endocytosed during EMT and utilized as initiators of cytoplasmic signaling pathways that often initiate protein phosphorylation to activate the gene regulatory networks. Acquisition of cell motility after EMT in these mesenchyme cells is associated with the expression of proteins such as Lefty, Snail and Seawi. Structural simplicity and genomic database of this model will further promote detailed EMT research.

Monohydroxylated polycyclic aromatic hydrocarbons influence spicule formation in the early development of sea urchins (Hemicentrotus pulcherrimus)

We previously demonstrated that monohydroxylated polycyclic aromatic hydrocarbons (OHPAHs), which are metabolites of polycyclic aromatic hydrocarbons (PAHs), act on calcified tissue and suppress osteoblastic and osteoclastic activity in the scales of teleost fish. The compounds may possibly influence other calcified tissues. Thus, the present study noted the calcified spicules in sea urchins and examined the effect of both PAHs and OHPAHs on spicule formation during the embryogenesis of sea urchins. After fertilization, benz[a]anthracene (BaA) and 4-hydroxybenz[a]anthracene (4-OHBaA) were added to seawater at concentrations of 10(-8) and 10(-7) M and kept at 18 °C. The influence of the compound was given at the time of the pluteus larva. At this stage, the length of the spicule was significantly suppressed by 4-OHBaA (10(-8) and 10(-7) M). BaA (10(-7) M) decreased the length of the spicule significantly, while the length did not change with BaA (10(-8) M). The expression of mRNAs (spicule matrix protein and transcription factors) in the 4-OHBaA (10(-7) M)-treated embryos was more strongly inhibited than were those in the BaA (10(-7) M)-treated embryos. This is the first study to demonstrate that OHPAHs suppress spicule formation in sea urchins.

microRNAs regulate β-catenin of the Wnt signaling pathway in early sea urchin development

Development of complex multicellular organisms requires careful regulation at both transcriptional and post-transcriptional levels. Post-transcriptional gene regulation is in part mediated by a class of non-coding RNAs of 21-25 nucleotides in length known as microRNAs (miRNAs). β-catenin, regulated by the canonical Wnt signaling pathway, has a highly evolutionarily conserved function in patterning early metazoan embryos, in forming the Anterior-Posterior axis, and in establishing the endomesoderm. Using reporter constructs and site-directed mutagenesis, we identified at least three miRNA binding sites within the 3' untranslated region (3'UTR) of the sea urchin β-catenin. Further, blocking these three miRNA binding sites within the β-catenin 3'UTR to prevent regulation of endogenous β-catenin by miRNAs resulted in a minor increase in β-catenin protein accumulation that is sufficient to induce aberrant gut morphology and circumesophageal musculature. These phenotypes are likely the result of increased transcript levels of Wnt responsive endomesodermal regulatory genes. This study demonstrates the importance of miRNA regulation of β-catenin in early development.

Regulatory logic and pattern formation in the early sea urchin embryo

We model the endomesoderm tissue specification process in the vegetal half of the early sea urchin embryo using Boolean models with continuous-time updating to represent the regulatory network that controls gene expression. Our models assume that the network interaction rules remain constant over time and the dynamics plays out on a predetermined program of cell divisions. An exhaustive search of two-node models, in which each node may represent a module of several genes in the real regulatory network, yields a unique network architecture that can accomplish the pattern formation task at hand--the formation of three latitudinal tissue bands from an initial state with only two distinct cell types. Analysis of an eight-gene model constructed from available experimental data reveals that it has a modular structure equivalent to the successful two-node case. Our results support the hypothesis that the gene regulatory network provides sufficient instructions for producing the correct pattern of tissue specification at this stage of development (between the fourth and tenth cleavages in the urchin embryo).

Larval mesenchyme cell specification in the primitive echinoid occurs independently of the double-negative gate

Echinoids (sea urchins) are divided into two major groups - cidaroids (a 'primitive' group) and euechinoids (a 'derived' group). The cidaroids are a promising model species for understanding the ancestral developmental mechanisms in echinoids, but little is known about the molecular mechanisms of cidaroid development. In euechinoids, skeletogenic mesenchyme cell specification is regulated by the double-negative gate (DNG), in which hesC represses the transcription of the downstream mesenchyme specification genes (alx1, tbr and ets1), thereby defining the prospective mesenchyme region. To estimate the ancestral mechanism of larval mesenchyme cell specification in echinoids, the expression patterns and roles of mesenchyme specification genes in the cidaroid Prionocidaris baculosa were examined. The present study reveals that the expression pattern and function of hesC in P. baculosa were inconsistent with the DNG model, suggesting that the euechinoid-type DNG is not utilized during cidaroid mesenchyme specification. In contrast with hesC, the expression patterns and functions of alx1, tbr and ets1 were similar between P. baculosa and euechinoids. Based on these results, we propose that the roles of alx1, tbr and ets1 in mesenchyme specification were established in the common ancestor of echinoids, and that the DNG system was acquired in the euechinoid lineage after divergence from the cidaroid ancestor. The evolutionary timing of the establishment of the DNG suggests that the DNG was originally related to micromere and/or primary mesenchyme cell formation but not to skeletogenic cell differentiation.

Genome-wide analysis of the skeletogenic gene regulatory network of sea urchins

A central challenge of developmental and evolutionary biology is to understand the transformation of genetic information into morphology. Elucidating the connections between genes and anatomy will require model morphogenetic processes that are amenable to detailed analysis of cell/tissue behaviors and to systems-level approaches to gene regulation. The formation of the calcified endoskeleton of the sea urchin embryo is a valuable experimental system for developing such an integrated view of the genomic regulatory control of morphogenesis. A transcriptional gene regulatory network (GRN) that underlies the specification of skeletogenic cells (primary mesenchyme cells, or PMCs) has recently been elucidated. In this study, we carried out a genome-wide analysis of mRNAs encoded by effector genes in the network and uncovered transcriptional inputs into many of these genes. We used RNA-seq to identify >400 transcripts differentially expressed by PMCs during gastrulation, when these cells undergo a striking sequence of behaviors that drives skeletal morphogenesis. Our analysis expanded by almost an order of magnitude the number of known (and candidate) downstream effectors that directly mediate skeletal morphogenesis. We carried out genome-wide analysis of (1) functional targets of Ets1 and Alx1, two pivotal, early transcription factors in the PMC GRN, and (2) functional targets of MAPK signaling, a pathway that plays an essential role in PMC specification. These studies identified transcriptional inputs into >200 PMC effector genes. Our work establishes a framework for understanding the genomic regulatory control of a major morphogenetic process and has important implications for reconstructing the evolution of biomineralization in metazoans.

Developmental gene regulatory network evolution: insights from comparative studies in echinoderms

One of the central concerns of Evolutionary Developmental biology is to understand how the specification of cell types can change during evolution. In the last decade, developmental biology has progressed toward a systems level understanding of cell specification processes. In particular, the focus has been on determining the regulatory interactions of the repertoire of genes that make up gene regulatory networks (GRNs). Echinoderms provide an extraordinary model system for determining how GRNs evolve. This review highlights the comparative GRN analyses arising from the echinoderm system. This work shows that certain types of GRN subcircuits or motifs, i.e., those involving positive feedback, tend to be conserved and may provide a constraint on development. This conservation may be due to a required arrangement of transcription factor binding sites in cis regulatory modules. The review will also discuss ways in which novelty may arise, in particular through the co-option of regulatory genes and subcircuits. The development of the sea urchin larval skeleton, a novel feature that arose in echinoderms, has provided a model for study of co-option mechanisms. Finally, the types of GRNs that can permit the great diversity in the patterns of ciliary bands and their associated neurons found among these taxa are discussed. The availability of genomic resources is rapidly expanding for echinoderms, including genome sequences not only for multiple species of sea urchins but also a species of sea star, sea cucumber, and brittle star. This will enable echinoderms to become a particularly powerful system for understanding how developmental GRNs evolve.

Characterization and Endocytic Internalization of Epith-2 Cell Surface Glycoprotein during the Epithelial-to-Mesenchymal Transition in Sea Urchin Embryos

The epithelial cells of the sea urchin Hemicentrotus pulcherrimus embryo express an Epith-2, uncharacterized glycoprotein, on the lateral surface. Here, we describe internalization of Epith-2 during mesenchyme formation through the epithelial-to-mesenchymal transition (EMT). Epith-2 was first expressed on the entire egg surface soon after fertilization and on the blastomeres until the 4-cell stage, but was localized to the lateral surface of epithelial cells at and after the 16-cell stage throughout the later developmental period. However, primary mesenchyme cells (PMC) and secondary mesenchyme cells (SMC) that ingress by EMT lost Epith-2 from their cell surface by endocytosis during dissociation from the epithelium, which was associated with the appearance of cytoplasmic Epith-2 dots. The cytoplasmic Epith-2 retained a similar relative molecular mass to that of the cell surface immediately after ingression through the early period of the spreading to single cells. Then, Epith-2 was completely lost from the cytoplasm. Tyrosine residues of Epith-2 were phosphorylated. The endocytic retraction of Epith-2 was inhibited by herbimycin A (HA), a protein tyrosine kinase (PTK) inhibitor, and suramin, a growth factor receptor (GFR) inhibitor, suggesting the involvement of the GFR/PTK (GP) signaling pathway. These two GP inhibitors also inhibited PMC and SMC spreading to individual cells after ingression, but the dissociation of PMC and SMC from the epithelium was not inhibited. In suramin-treated embryos, dissociated mesenchyme cells migrated partially by retaining their epithelial morphology. In HA-treated embryos, no mesenchyme cells migrated. Thus, the EMT occurs in relation to internalization of Epith-2 from presumptive PMC and SMC.

Development of an embryonic skeletogenic mesenchyme lineage in a sea cucumber reveals the trajectory of change for the evolution of novel structures in echinoderms

BACKGROUND

The mechanisms by which the conserved genetic "toolkit" for development generates phenotypic disparity across metazoans is poorly understood. Echinoderm larvae provide a great resource for understanding how developmental novelty arises. The sea urchin pluteus larva is dramatically different from basal echinoderm larval types, which include the auricularia-type larva of its sister taxon, the sea cucumbers, and the sea star bipinnaria larva. In particular, the pluteus has a mesodermally-derived larval skeleton that is not present in sea star larvae or any outgroup taxa. To understand the evolutionary origin of this structure, we examined the molecular development of mesoderm in the sea cucumber, Parastichopus parvimensis.

RESULTS

By comparing gene expression in sea urchins, sea cucumbers and sea stars, we partially reconstructed the mesodermal regulatory state of the echinoderm ancestor. Surprisingly, we also identified expression of the transcription factor alx1 in a cryptic skeletogenic mesenchyme lineage in P. parvimensis. Orthologs of alx1 are expressed exclusively within the sea urchin skeletogenic mesenchyme, but are not expressed in the mesenchyme of the sea star, which suggests that alx1+ mesenchyme is a synapomorphy of at least sea urchins and sea cucumbers. Perturbation of Alx1 demonstrates that this protein is necessary for the formation of the sea cucumber spicule. Overexpression of the sea star alx1 ortholog in sea urchins is sufficient to induce additional skeleton, indicating that the Alx1 protein has not evolved a new function during the evolution of the larval skeleton.

CONCLUSIONS

The proposed echinoderm ancestral mesoderm state is highly conserved between the morphologically similar, but evolutionarily distant, auricularia and bipinnaria larvae. However, the auricularia, but not bipinnaria, also develops a simple skelotogenic cell lineage. Our data indicate that the first step in acquiring these novel cell fates was to re-specify the ancestral mesoderm into molecularly distinct territories. These new territories likely consisted of only a few cells with few regulatory differences from the ancestral state, thereby leaving the remaining mesoderm to retain its original function. The new territories were then free to take on a new fate. Partitioning of existing gene networks was a necessary pre-requisite to establish novelty in this system.

Zinc-finger nuclease-mediated targeted insertion of reporter genes for quantitative imaging of gene expression in sea urchin embryos

To understand complex biological systems, such as the development of multicellular organisms, it is important to characterize the gene expression dynamics. However, there is currently no universal technique for targeted insertion of reporter genes and quantitative imaging in multicellular model systems. Recently, genome editing using zinc-finger nucleases (ZFNs) has been reported in several models. ZFNs consist of a zinc-finger DNA-binding array with the nuclease domain of the restriction enzyme FokI and facilitate targeted transgene insertion. In this study, we successfully inserted a GFP reporter cassette into the HpEts1 gene locus of the sea urchin, Hemicentrotus pulcherrimus. We achieved this insertion by injecting eggs with a pair of ZFNs for HpEts1 with a targeting donor construct that contained ∼1-kb homology arms and a 2A-histone H2B-GFP cassette. We increased the efficiency of the ZFN-mediated targeted transgene insertion by in situ linearization of the targeting donor construct and cointroduction of an mRNA for a dominant-negative form of HpLig4, which encodes the H. pulcherrimus homolog of DNA ligase IV required for error-prone nonhomologous end joining. We measured the fluorescence intensity of GFP at the single-cell level in living embryos during development and found that there was variation in HpEts1 expression among the primary mesenchyme cells. These findings demonstrate the feasibility of ZFN-mediated targeted transgene insertion to enable quantification of the expression levels of endogenous genes during development in living sea urchin embryos.

The genomic regulatory control of skeletal morphogenesis in the sea urchin

A central challenge of developmental and evolutionary biology is to understand how anatomy is encoded in the genome. Elucidating the genetic mechanisms that control the development of specific anatomical features will require the analysis of model morphogenetic processes and an integration of biological information at genomic, cellular and tissue levels. The formation of the endoskeleton of the sea urchin embryo is a powerful experimental system for developing such an integrated view of the genomic regulatory control of morphogenesis. The dynamic cellular behaviors that underlie skeletogenesis are well understood and a complex transcriptional gene regulatory network (GRN) that underlies the specification of embryonic skeletogenic cells (primary mesenchyme cells, PMCs) has recently been elucidated. Here, we link the PMC specification GRN to genes that directly control skeletal morphogenesis. We identify new gene products that play a proximate role in skeletal morphogenesis and uncover transcriptional regulatory inputs into many of these genes. Our work extends the importance of the PMC GRN as a model developmental GRN and establishes a unique picture of the genomic regulatory control of a major morphogenetic process. Furthermore, because echinoderms exhibit diverse programs of skeletal development, the newly expanded sea urchin skeletogenic GRN will provide a foundation for comparative studies that explore the relationship between GRN evolution and morphological evolution.

Specific expression of a TRIM-containing factor in ectoderm cells affects the skeletal morphogenetic program of the sea urchin embryo

In the indirect developing sea urchin embryo, the primary mesenchyme cells (PMCs) acquire most of the positional and temporal information from the overlying ectoderm for skeletal initiation and growth. In this study, we characterize the function of the novel gene strim1, which encodes a tripartite motif-containing (TRIM) protein, that adds to the list of genes constituting the epithelial-mesenchymal signaling network. We report that strim1 is expressed in ectoderm regions adjacent to the bilateral clusters of PMCs and that its misexpression leads to severe skeletal abnormalities. Reciprocally, knock down of strim1 function abrogates PMC positioning and blocks skeletogenesis. Blastomere transplantation experiments establish that the defects in PMC patterning, number and skeletal growth depend upon strim1 misexpression in ectoderm cells. Furthermore, clonal expression of strim1 into knocked down embryos locally restores skeletogenesis. We also provide evidence that the Otp and Pax2/5/8 regulators, as well as FGFA, but not VEGF, ligand act downstream to strim1 in ectoderm cells, and that strim1 triggers the expression of the PMC marker sm30, an ectoderm-signaling dependent gene. We conclude that the strim1 function elicits specific gene expression both in ectoderm cells and PMCs to guide the skeletal biomineralization during morphogenesis.

Differential gene expression at coral settlement and metamorphosis--a subtractive hybridization study

Background

A successful metamorphosis from a planktonic larva to a settled polyp, which under favorable conditions will establish a future colony, is critical for the survival of corals. However, in contrast to the situation in other animals, e.g., frogs and insects, little is known about the molecular basis of coral metamorphosis. We have begun to redress this situation with previous microarray studies, but there is still a great deal to learn. In the present paper we have utilized a different technology, subtractive hybridization, to characterize genes differentially expressed across this developmental transition and to compare the success of this method to microarray.

Methodology/Principal Findings

Suppressive subtractive hybridization (SSH) was used to identify two pools of transcripts from the coral, Acropora millepora. One is enriched for transcripts expressed at higher levels at the pre-settlement stage, and the other for transcripts expressed at higher levels at the post-settlement stage. Virtual northern blots were used to demonstrate the efficacy of the subtractive hybridization technique. Both pools contain transcripts coding for proteins in various functional classes but transcriptional regulatory proteins were represented more frequently in the post-settlement pool. Approximately 18% of the transcripts showed no significant similarity to any other sequence on the public databases. Transcripts of particular interest were further characterized by in situ hybridization, which showed that many are regulated spatially as well as temporally. Notably, many transcripts exhibit axially restricted expression patterns that correlate with the pool from which they were isolated. Several transcripts are expressed in patterns consistent with a role in calcification.

Conclusions

We have characterized over 200 transcripts that are differentially expressed between the planula larva and post-settlement polyp of the coral, Acropora millepora. Sequence, putative function, and in some cases temporal and spatial expression are reported.

Regulative deployment of the skeletogenic gene regulatory network during sea urchin development

The well-known regulative properties of the sea urchin embryo, coupled with the recent elucidation of gene regulatory networks (GRNs) that underlie cell specification, make this a valuable experimental model for analyzing developmental plasticity. In the sea urchin, the primary mesenchyme cell (PMC) GRN controls the development of the embryonic skeleton. Remarkably, experimental manipulations reveal that this GRN can be activated in almost any cell of the embryo. Here, we focus on the activation of the PMC GRN during gastrulation by non-skeletogenic mesoderm (NSM) cells and by endoderm cells. We show that most transfating NSM cells are prospective blastocoelar cells, not prospective pigment cells, as was previously believed. Earlier work showed that the regulative deployment of the GRN, unlike its deployment in the micromere-PMC lineage, is independent of the transcriptional repressor Pmar1. In this work, we identify several additional differences in the upstream regulation of the GRN during normal and regulative development. We provide evidence that, despite these changes in the upstream regulation of the network, downstream regulatory genes and key morphoregulatory genes are deployed in transfating NSM cells in a fashion that recapitulates the normal deployment of the GRN, and which can account for the striking changes in migratory behavior that accompany NSM transfating. Finally, we report that mitotic cell division is not required for genomic reprogramming in this system, either within a germ layer (NSM transfating) or across a germ layer boundary (endoderm transfating).

The control of foxN2/3 expression in sea urchin embryos and its function in the skeletogenic gene regulatory network

Early development requires well-organized temporal and spatial regulation of transcription factors that are assembled into gene regulatory networks (GRNs). In the sea urchin, an endomesoderm GRN model explains much of the specification in the endoderm and mesoderm prior to gastrulation, yet some GRN connections remain incomplete. Here, we characterize FoxN2/3 in the primary mesenchyme cell (PMC) GRN state. Expression of foxN2/3 mRNA begins in micromeres at the hatched blastula stage and then is lost from micromeres at the mesenchyme blastula stage. foxN2/3 expression then shifts to the non-skeletogenic mesoderm and, later, to the endoderm. Here, we show that Pmar1, Ets1 and Tbr are necessary for activation of foxN2/3 in micromeres. The later endomesoderm expression of foxN2/3 is independent of the earlier expression of foxN2/3 in micromeres and is independent of signals from PMCs. FoxN2/3 is necessary for several steps in the formation of the larval skeleton. Early expression of genes for the skeletal matrix is dependent on FoxN2/3, but only until the mesenchyme blastula stage as foxN2/3 mRNA disappears from PMCs at that time and we assume that the protein is not abnormally long-lived. Knockdown of FoxN2/3 inhibits normal PMC ingression and foxN2/3 morphant PMCs do not organize in the blastocoel and fail to join the PMC syncytium. In addition, without FoxN2/3, the PMCs fail to repress the transfating of other mesodermal cells into the skeletogenic lineage. Thus, FoxN2/3 is necessary for normal ingression, for expression of several skeletal matrix genes, for preventing transfating and for fusion of the PMC syncytium.

P58-A and P58-B: novel proteins that mediate skeletogenesis in the sea urchin embryo

During sea urchin embryogenesis, the skeleton is produced by primary mesenchyme cells (PMCs). PMCs undergo a sequence of morphogenetic behaviors that includes ingression, directed migration, and cell-cell fusion. Ultimately, PMCs deposit the calcite-containing biomineral that forms the endoskeleton of the late embryo and early larva. The endoskeleton has a stereotypical structure and is the major determinant of the distinctive, angular shape of the larva. Although many candidate biomineralization proteins have been identified, functional data concerning these proteins are scant. Here, we identify and characterize two new biomineralization genes, p58-a and p58-b. We show that these two genes are highly conserved in Strongylocentrotus purpuratus and Lytechinus variegatus, two sea urchin species whose ancestors diverged approximately 100 mya. The p58-a and p58-b genes lie in tandem on the chromosome, suggesting that one of the two genes arose via a gene duplication event. The two genes encode closely related, type I transmembrane proteins. We have established by whole mount in situ hybridization that p58-a and p58-b are expressed specifically in the PMCs in both species. Knockdown of either gene by morpholino antisense oligonucleotides leads to profound defects in skeletogenesis, although skeletal elements are not completely eliminated. The P58-A and P58-B proteins do not appear to play a role in the specification, directed migration or differentiation of the PMCs, but most likely are directly involved in biomineralization during sea urchin embryonic development.

Implication of HpEts in gene regulatory networks responsible for specification of sea urchin skeletogenic primary mesenchyme cells

The large micromeres of the 32-cell stage of sea urchin embryos are autonomously specified and differentiate into primary mesenchyme cells (PMCs), giving rise to the skeletogenic cells. We previously demonstrated that HpEts, an ets-related transcription factor, plays an essential role in the specification of PMCs in sea urchin embryos. In order to clarify the function of HpEts in the gene regulatory network involved in PMC specification, we analyzed the zygotic expression pattern and the cis-regulatory region of HpEts, and examined the activity of the HpEts protein as a transcription factor. Intron-based PCR reveals that zygotic expression of HpEts starts at the cleavage stage, and that the rate of transcription reaches maximum at the unhatched blastula stage. A series of progressive deletions of the fragments from -4.2 kbp to +1206 bp of the HpEts, which directs PMC-specific expression, caused a gradual decrease in the specificity, implying that coordination of several cis-regulatory elements regulates the expression in PMCs. A minimum cis-element required for the temporal expression is located within a 10 bp from -243 bp to -234 bp. The HpEts protein remains in the cytoplasm of entire embryonic cells in the cleavage stage. At the unhatched blastula stage, the HpEts protein translocates into the nucleus in presumptive PMCs. Transactivation assays demonstrate that the HpEts protein activates a promoter of Spicule Matrix Protein 50 (SM50), which is a target of HpEts, which binds to the regulatory region of SM50.

Functional evolution of Ets in echinoderms with focus on the evolution of echinoderm larval skeletons

Convergent evolution of echinoderm pluteus larva was examined from the standpoint of functional evolution of a transcription factor Ets1/2. In sea urchins, Ets1/2 plays a central role in the differentiation of larval skeletogenic mesenchyme cells. In addition, Ets1/2 is suggested to be involved in adult skeletogenesis. Conversely, in starfish, although no skeletogenic cells differentiate during larval development, Ets1/2 is also expressed in the larval mesoderm. Here, we confirmed that the starfish Ets1/2 is indispensable for the differentiation of the larval mesoderm. This result led us to assume that, in the common ancestors of echinoderms, Ets1/2 activates the transcription of distinct gene sets, one for the differentiation of the larval mesoderm and the other for the development of the adult skeleton. Thus, the acquisition of the larval skeleton involved target switching of Ets1/2. Specifically, in the sea urchin lineage, Ets1/2 activated a downstream target gene set for skeletogenesis during larval development in addition to a mesoderm target set. We examined whether this heterochronic activation of the skeletogenic target set was achieved by the molecular evolution of the Ets1/2 transcription factor itself. We tested whether starfish Ets1/2 induced skeletogenesis when injected into sea urchin eggs. We found that, in addition to ectopic induction of mesenchyme cells, starfish Ets1/2 can activate some parts of the skeletogenic pathway in these mesenchyme cells. Thus, we suggest that the nature of the transcription factor Ets1/2 did not change, but rather that some unidentified co-factor(s) for Ets1/2 may distinguish between targets for the larval mesoderm and for skeletogenesis. Identification of the co-factor(s) will be key to understanding the molecular evolution underlying the evolution of the pluteus larvae.

Activation of the skeletogenic gene regulatory network in the early sea urchin embryo

The gene regulatory network (GRN) that underlies the development of the embryonic skeleton in sea urchins is an important model for understanding the architecture and evolution of developmental GRNs. The initial deployment of the network is thought to be regulated by a derepression mechanism, which is mediated by the products of the pmar1 and hesC genes. Here, we show that the activation of the skeletogenic network occurs by a mechanism that is distinct from the transcriptional repression of hesC. By means of quantitative, fluorescent whole-mount in situ hybridization, we find that two pivotal early genes in the network, alx1 and delta, are activated in prospective skeletogenic cells prior to the downregulation of hesC expression. An analysis of the upstream regulation of alx1 shows that this gene is regulated by MAPK signaling and by the transcription factor Ets1; however, these inputs influence only the maintenance of alx1 expression and not its activation, which occurs by a distinct mechanism. By altering normal cleavage patterns, we show that the zygotic activation of alx1 and delta, but not that of pmar1, is dependent upon the unequal division of vegetal blastomeres. Based on these findings, we conclude that the widely accepted double-repression model is insufficient to account for the localized activation of the skeletogenic GRN. We postulate the existence of additional, unidentified repressors that are controlled by pmar1, and propose that the ability of pmar1 to derepress alx1 and delta is regulated by the unequal division of vegetal blastomeres.

Modularity and design principles in the sea urchin embryo gene regulatory network

The gene regulatory network (GRN) established experimentally for the pre-gastrular sea urchin embryo provides causal explanations of the biological functions required for spatial specification of embryonic regulatory states. Here we focus on the structure of the GRN which controls the progressive increase in complexity of territorial regulatory states during embryogenesis; and on the types of modular subcircuits of which the GRN is composed. Each of these subcircuit topologies executes a particular operation of spatial information processing. The GRN architecture reflects the particular mode of embryogenesis represented by sea urchin development. Network structure not only specifies the linkages constituting the genomic regulatory code for development, but also indicates the various regulatory requirements of regional developmental processes.

The cis-regulatory system of the tbrain gene: Alternative use of multiple modules to promote skeletogenic expression in the sea urchin embryo

The genomic cis-regulatory systems controlling regulatory gene expression usually include multiple modules. The regulatory output of such systems at any given time depends on which module is directing the function of the basal transcription apparatus, and ultimately on the transcription factor inputs into that module. Here we examine regulation of the Strongylocentrotus purpuratus tbrain gene, a required activator of the skeletogenic specification state in the lineage descendant from the embryo micromeres. Alternate cis-regulatory modules were found to convey skeletogenic expression in reporter constructs. To determine their relative developmental functions in context, we made use of recombineered BAC constructs containing a GFP reporter and of derivatives from which specific modules had been deleted. The outputs of the various constructs were observed spatially by GFP fluorescence and quantitatively over time by QPCR. In the context of the complete genomic locus, early skeletogenic expression is controlled by an intron enhancer plus a proximal region containing a HesC site as predicted from network analysis. From ingression onward, however, a dedicated distal module utilizing positive Ets1/2 inputs contributes to definitive expression in the skeletogenic mesenchyme. This module also mediates a newly discovered negative Erg input which excludes non-skeletogenic mesodermal expression.

Global regulatory logic for specification of an embryonic cell lineage

Explanation of a process of development must ultimately be couched in the terms of the genomic regulatory code. Specification of an embryonic cell lineage is driven by a network of interactions among genes encoding transcription factors. Here, we present the gene regulatory network (GRN) that directs the specification of the skeletogenic micromere lineage of the sea urchin embryo. The GRN now includes all regulatory genes expressed in this lineage up to late blastula stage, as identified in a genomewide survey. The architecture of the GRN was established by a large-scale perturbation analysis in which the expression of each gene in the GRN was cut off by use of morpholinos, and the effects on all other genes were measured quantitatively. Several cis-regulatory analyses provided additional evidence. The explanatory power of the GRN suffices to provide a causal explanation for all observable developmental functions of the micromere lineage during the specification period. These functions are: (i) initial acquisition of identity through transcriptional interpretation of localized maternal cues; (ii) activation of specific regulatory genes by use of a double negative gate; (iii) dynamic stabilization of the regulatory state by activation of a feedback subcircuit; (iv) exclusion of alternative regulatory states; (v) presentation of a signal required by the micromeres themselves and of two different signals required for development of adjacent endomesodermal lineages; and (vi) lineage-specific activation of batteries of skeletogenic genes. The GRN precisely predicts gene expression responses and provides a coherent explanation of the biology of specification.

Twist is an essential regulator of the skeletogenic gene regulatory network in the sea urchin embryo

Recent work on the sea urchin endomesoderm gene regulatory network (GRN) offers many opportunities to study the specification and differentiation of each cell type during early development at a mechanistic level. The mesoderm lineages consist of two cell populations, primary and secondary mesenchyme cells (PMCs and SMCs). The micromere-PMC GRN governs the development of the larval skeleton, which is the exclusive fate of PMCs, and SMCs diverge into four lineages, each with its own GRN state. Here we identify a sea urchin ortholog of the Twist transcription factor, and show that it plays an essential role in the PMC GRN and later is involved in SMC formation. Perturbations of Twist either by morpholino knockdown or by overexpression result in defects in progressive phases of PMC development, including specification, ingression/EMT, differentiation and skeletogenesis. Evidence is presented that Twist expression is required for the maintenance of the PMC specification state, and a reciprocal regulation between Alx1 and Twist offers stability for the subsequent processes, such as PMC differentiation and skeletogenesis. These data illustrate the significance of regulatory state maintenance and continuous progression during cell specification, and the dynamics of the sequential events that depend on those earlier regulatory states.

FGF signals guide migration of mesenchymal cells, control skeletal morphogenesis and regulate gastrulation during sea urchin development

The sea urchin embryo is emerging as an attractive model to study morphogenetic processes such as directed migration of mesenchyme cells and cell sheet invagination, but surprisingly, few of the genes regulating these processes have yet been characterized. We present evidence that FGFA, the first FGF family member characterized in the sea urchin, regulates directed migration of mesenchyme cells, morphogenesis of the skeleton and gastrulation during early development. We found that at blastula stages, FGFA and a novel putative FGF receptor are expressed in a pattern that prefigures morphogenesis of the skeletogenic mesoderm and that suggests that FGFA is one of the elusive signals that guide migration of primary mesenchyme cells (PMCs). We first show that fgfA expression is correlated with abnormal migration and patterning of the PMCs following treatments that perturb specification of the ectoderm along the oral-aboral and animal-vegetal axes. Specification of the ectoderm initiated by Nodal is required to restrict fgfA to the lateral ectoderm, and in the absence of Nodal, fgfA is expressed ectopically throughout most of the ectoderm. Inhibition of either FGFA, FGFR1 or FGFR2 function severely affects morphogenesis of the skeleton. Furthermore,inhibition of FGFA and FGFR1 signaling dramatically delays invagination of the archenteron, prevents regionalization of the gut and abrogates formation of the stomodeum. We identified several genes acting downstream of fgfAin these processes, including the transcription factors pea3 and pax2/5/8 and the signaling molecule sprouty in the lateral ectoderm and SM30 and SM50 in the primary mesenchyme cells. This study identifies the FGF signaling pathway as an essential regulator of gastrulation and directed cell migration in the sea urchin embryo and as a key player in the gene regulatory network directing morphogenesis of the skeleton.