The micro1 gene is necessary and sufficient for micromere differentiation and mid/hindgut-inducing activity in the sea urchin embryo

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.

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.

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.

Diversification of spatiotemporal expression and copy number variation of the echinoid hbox12/pmar1/micro1 multigene family

Changes occurring during evolution in the cis-regulatory landscapes of individual members of multigene families might impart diversification in their spatiotemporal expression and function. The archetypal member of the echinoid hbox12/pmar1/micro1 family is hbox12-a, a homeobox-containing gene expressed exclusively by dorsal blastomeres, where it governs the dorsal/ventral gene regulatory network during embryogenesis of the sea urchin Paracentrotus lividus. Here we describe the inventory of the hbox12/pmar1/micro1 genes in P. lividus, highlighting that gene copy number variation occurs across individual sea urchins of the same species. We show that the various hbox12/pmar1/micro1 genes group into three subfamilies according to their spatiotemporal expression, which ranges from broad transcription throughout development to transient expression in either the animal hemisphere or micromeres of the early embryo. Interestingly, the promoter regions of those genes showing comparable expression patterns are highly similar, while differing from those of the other subfamilies. Strikingly, phylogenetic analysis suggests that the hbox12/pmar1/micro1 genes are species-specific, exhibiting extensive divergence in their noncoding, but not in their coding, sequences across three distinct sea urchin species. In spite of this, two micromere-specific genes of P. lividus possess a TCF/LEF-binding motif in a similar position, and their transcription relies on Wnt/β-catenin signaling, similar to the pmar1 and micro1 genes, which in other sea urchin species are involved in micromere specification. Altogether, our findings suggest that the hbox12/pmar1/micro1 gene family evolved rather rapidly, generating paralogs whose cis-regulatory sequences diverged following multiple rounds of duplication from a common ancestor.

Roles of hesC and gcm in echinoid larval mesenchyme cell development

To understand the roles of hesC and gcm during larval mesenchyme specification and differentiation in echinoids, we performed perturbation experiments for these genes in two distantly related euechinoids, Hemicentrotus pulcherrimus and Scaphechinus mirabilis. The number of larval mesenchyme cells increased when the translation of hesC was inhibited, thereby suggesting that hesC has a general role in larval mesenchyme development. We confirmed previous results by demonstrating that gcm is involved in pigment cell differentiation. Simultaneous inhibition of the translation of hesC and gcm induced a significant increase in the number of skeletogenic cells, which suggests that gcm functions in skeletogenic fate repression. Based on these observations, we suggest that: (i) hesC participates in some general aspects of mesenchymal cell development; and (ii) gcm is involved in the mechanism responsible for the binary specification of skeletogenic and pigment cell fates.

Large-scale gene expression study in the ophiuroid Amphiura filiformis provides insights into evolution of gene regulatory networks

Background

The evolutionary mechanisms involved in shaping complex gene regulatory networks (GRN) that encode for morphologically similar structures in distantly related animals remain elusive. In this context, echinoderm larval skeletons found in brittle stars and sea urchins provide an ideal system. Here, we characterize for the first time the development of the larval skeleton in the ophiuroid Amphiura filiformis and compare it systematically with its counterpart in sea urchin.

Results

We show that ophiuroids and euechinoids, that split at least 480 Million years ago (Mya), have remarkable similarities in tempo and mode of skeletal development. Despite morphological and ontological similarities, our high-resolution study of the dynamics of genetic regulatory states in A. filiformis highlights numerous differences in the architecture of their underlying GRNs. Importantly, the A.filiformispplx, the closest gene to the sea urchin double negative gate (DNG) repressor pmar1, fails to drive the skeletogenic program in sea urchin, showing important evolutionary differences in protein function. hesC, the second repressor of the DNG, is co-expressed with most of the genes that are repressed in sea urchin, indicating the absence of direct repression of tbr, ets1/2, and delta in A. filiformis. Furthermore, the absence of expression in later stages of brittle star skeleton development of key regulatory genes, such as foxb and dri, shows significantly different regulatory states.

Conclusion

Our data fill up an important gap in the picture of larval mesoderm in echinoderms and allows us to explore the evolutionary implications relative to the recently established phylogeny of echinoderm classes. In light of recent studies on other echinoderms, our data highlight a high evolutionary plasticity of the same nodes throughout evolution of echinoderm skeletogenesis. Finally, gene duplication, protein function diversification, and cis-regulatory element evolution all contributed to shape the regulatory program for larval skeletogenesis in different branches of echinoderms.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0039-x) contains supplementary material, which is available to authorized users.

The Maternal Maverick/GDF15-like TGF-β Ligand Panda Directs Dorsal-Ventral Axis Formation by Restricting Nodal Expression in the Sea Urchin Embryo

Specification of the dorsal-ventral axis in the highly regulative sea urchin embryo critically relies on the zygotic expression of nodal, but whether maternal factors provide the initial spatial cue to orient this axis is not known. Although redox gradients have been proposed to entrain the dorsal-ventral axis by acting upstream of nodal, manipulating the activity of redox gradients only has modest consequences, suggesting that other factors are responsible for orienting nodal expression and defining the dorsal-ventral axis. Here we uncover the function of Panda, a maternally provided transforming growth factor beta (TGF-β) ligand that requires the activin receptor-like kinases (Alk) Alk3/6 and Alk1/2 receptors to break the radial symmetry of the embryo and orient the dorsal-ventral axis by restricting nodal expression. We found that the double inhibition of the bone morphogenetic protein (BMP) type I receptors Alk3/6 and Alk1/2 causes a phenotype dramatically more severe than the BMP2/4 loss-of-function phenotype, leading to extreme ventralization of the embryo through massive ectopic expression of nodal, suggesting that an unidentified signal acting through BMP type I receptors cooperates with BMP2/4 to restrict nodal expression. We identified this ligand as the product of maternal Panda mRNA. Double inactivation of panda and bmp2/4 led to extreme ventralization, mimicking the phenotype caused by inactivation of the two BMP receptors. Inhibition of maternal panda mRNA translation disrupted the early spatial restriction of nodal, leading to persistent massive ectopic expression of nodal on the dorsal side despite the presence of Lefty. Phylogenetic analysis indicates that Panda is not a prototypical BMP ligand but a member of a subfamily of TGF-β distantly related to Inhibins, Lefty, and TGF-β that includes Maverick from Drosophila and GDF15 from vertebrates. Indeed, overexpression of Panda does not appear to directly or strongly activate phosphoSmad1/5/8 signaling, suggesting that although this TGF-β may require Alk1/2 and/or Alk3/6 to antagonize nodal expression, it may do so by sequestering a factor essential for Nodal signaling, by activating a non-Smad pathway downstream of the type I receptors, or by activating extremely low levels of pSmad1/5/8. We provide evidence that, although panda mRNA is broadly distributed in the early embryo, local expression of panda mRNA efficiently orients the dorsal-ventral axis and that Panda activity is required locally in the early embryo to specify this axis. Taken together, these findings demonstrate that maternal panda mRNA is both necessary and sufficient to orient the dorsal-ventral axis. These results therefore provide evidence that in the highly regulative sea urchin embryo, the activity of spatially restricted maternal factors regulates patterning along the dorsal-ventral axis.

Early asymmetric cues triggering the dorsal/ventral gene regulatory network of the sea urchin embryo

Dorsal/ventral (DV) patterning of the sea urchin embryo relies on a ventrally-localized organizer expressing Nodal, a pivotal regulator of the DV gene regulatory network. However, the inceptive mechanisms imposing the symmetry-breaking are incompletely understood. In Paracentrotus lividus, the Hbox12 homeodomain-containing repressor is expressed by prospective dorsal cells, spatially facing and preceding the onset of nodal transcription. We report that Hbox12 misexpression provokes DV abnormalities, attenuating nodal and nodal-dependent transcription. Reciprocally, impairing hbox12 function disrupts DV polarity by allowing ectopic expression of nodal. Clonal loss-of-function, inflicted by blastomere transplantation or gene-transfer assays, highlights that DV polarization requires Hbox12 action in dorsal cells. Remarkably, the localized knock-down of nodal restores DV polarity of embryos lacking hbox12 function. Finally, we show that hbox12 is a dorsal-specific negative modulator of the p38-MAPK activity, which is required for nodal expression. Altogether, our results suggest that Hbox12 function is essential for proper positioning of the DV organizer.

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.

Evolutionary modification of T-brain (tbr) expression patterns in sand dollar

The sand dollars are a group of irregular echinoids that diverged from other regular sea urchins approximately 200 million years ago. We isolated two orthologs of T-brain (tbr), Smtbr and Pjtbr, from the indirect developing sand dollar Scaphechinus mirabilis and the direct developing sand dollar Peronella japonica, respectively. The expression patterns of Smtbr and Pjtbr during early development were examined by whole mount in situ hybridization. The expression of Smtbr was first detected in micromere descendants in early blastula stage, similar to tbr expression in regular sea urchins. However, unlike in regular sea urchin, Smtbr expression in middle blastula stage was detected in micromere-descendent cells and a subset of macromere-descendant cells. At gastrula stage, expression of Smtbr was detected in part of the archenteron as well as primary mesenchyme cells. A similar pattern of tbr expression was observed in early Peronella embryos. A comparison of tbr expression patterns between sand dollars and other echinoderm species suggested that broader expression in the endomesoderm is an ancestral character of echinoderms. In addition to the endomesoderm, Pjtbr expression was detected in the apical organ, the animal-most part of the ectoderm.

Evolutionary origins of blastoporal expression and organizer activity of the vertebrate gastrula organizer gene lhx1 and its ancient metazoan paralog lhx3

Expression of the LIM homeobox gene lhx1 (lim1) is specific to the vertebrate gastrula organizer. Lhx1 functions as a transcriptional regulatory core protein to exert `organizer' activity in Xenopus embryos. Its ancient paralog, lhx3 (lim3),is expressed around the blastopore in amphioxus and ascidian, but not vertebrate, gastrulae. These two genes are thus implicated in organizer evolution, and we addressed the evolutionary origins of their blastoporal expression and organizer activity. Gene expression analysis of organisms ranging from cnidarians to chordates suggests that blastoporal expression has its evolutionary root in or before the ancestral eumetazoan for lhx1,but possibly in the ancestral chordate for lhx3, and that in the ascidian lineage, blastoporal expression of lhx1 ceased, whereas endodermal expression of lhx3 has persisted. Analysis of organizer activity using Xenopus embryos suggests that a co-factor of LIM homeodomain proteins, Ldb, has a conserved function in eumetazoans to activate Lhx1, but that Lhx1 acquired organizer activity in the bilaterian lineage,Lhx3 acquired organizer activity in the deuterostome lineage and ascidian Lhx3 acquired a specific transactivation domain to confer organizer activity on this molecule. Knockdown analysis using cnidarian embryos suggests that Lhx1 is required for chordin expression in the blastoporal region. These data suggest that Lhx1 has been playing fundamental roles in the blastoporal region since the ancestral eumetazoan arose, that it contributed as an`original organizer gene' to the evolution of the vertebrate gastrula organizer, and that Lhx3 could be involved in the establishment of organizer gene networks.

Lessons from a gene regulatory network: echinoderm skeletogenesis provides insights into evolution, plasticity and morphogenesis

Significant new insights have emerged from the analysis of a gene regulatory network (GRN) that underlies the development of the endoskeleton of the sea urchin embryo. Comparative studies have revealed ways in which this GRN has been modified (and conserved) during echinoderm evolution, and point to mechanisms associated with the evolution of a new cell lineage. The skeletogenic GRN has also recently been used to study the long-standing problem of developmental plasticity. Other recent findings have linked this transcriptional GRN to morphoregulatory proteins that control skeletal anatomy. These new studies highlight powerful new ways in which GRNs can be used to dissect development and the evolution of morphogenesis.

Krüppel-like is required for nonskeletogenic mesoderm specification in the sea urchin embryo

The canonical Wnt pathway plays a central role in specifying vegetal cell fate in sea urchin embryos. SpKrl has been cloned as a direct target of nuclear beta-catenin. Using Hemicentrotus pulcherrimus embryos, here we show that HpKrl controls the specification of secondary mesenchyme cells (SMCs) through both cell-autonomous and non-autonomous means. Like SpKrl, HpKrl was activated in both micromere and macromere progenies. To examine the functions of HpKrl in each blastomere, we constructed chimeric embryos composed of blastomeres from control and morpholino-mediated HpKrl-knockdown embryos and analyzed the phenotypes of the chimeras. Micromere-swapping experiments showed that HpKrl is not involved in micromere specification, while micromere-deprivation assays indicated that macromeres require HpKrl for cell-autonomous specification. Transplantation of normal micromeres into a micromere-less host with morpholino revealed that macromeres are able to receive at least some micromere signals regardless of HpKrl function. From these observations, we propose that two distinct pathways of endomesoderm formation exist in macromeres, a Krl-dependent pathway and a Krl-independent pathway. The Krl-independent pathway may correspond to the Delta/Notch signaling pathway via GataE and Gcm. We suggest that Krl may be a downstream component of nuclear beta-catenin required by macromeres for formation of more vegetal tissues, not as a member of the Delta/Notch pathway, but as a parallel effector of the signaling (Krl-dependent pathway).

Gene regulatory networks and developmental plasticity in the early sea urchin embryo: alternative deployment of the skeletogenic gene regulatory network

Cell fates in the sea urchin embryo are remarkably labile, despite the fact that maternal polarity and zygotic programs of differential gene expression pattern the embryo from the earliest stages. Recent work has focused on transcriptional gene regulatory networks (GRNs) deployed in specific embryonic territories during early development. The micromere-primary mesenchyme cell(PMC) GRN drives the development of the embryonic skeleton. Although normally deployed only by presumptive PMCs, every lineage of the early embryo has the potential to activate this pathway. Here, we focus on one striking example of regulative activation of the skeletogenic GRN; the transfating of non-skeletogenic mesoderm (NSM) cells to a PMC fate during gastrulation. We show that transfating is accompanied by the de novo expression of terminal,biomineralization-related genes in the PMC GRN, as well as genes encoding two upstream transcription factors, Lvalx1 and Lvtbr. We report that Lvalx1, a key component of the skeletogenic GRN in the PMC lineage, plays an essential role in the regulative pathway both in NSM cells and in animal blastomeres. MAPK signaling is required for the expression of Lvalx1 and downstream skeletogenic genes in NSM cells, mirroring its role in the PMC lineage. We also demonstrate that Lvalx1 regulates the signal from PMCs that normally suppresses NSM transfating. Significantly,misexpression of Lvalx1 in macromeres (the progenitors of NSM cells)is sufficient to activate the skeletogenic GRN. We suggest that NSM cells normally deploy a basal mesodermal pathway and require only an Lvalx1-mediated sub-program to express a PMC fate. Finally, we provide evidence that, in contrast to the normal pathway, activation of the skeletogenic GRN in NSM cells is independent of Lvpmar1. Our studies reveal that, although most features of the micromere-PMC GRN are recapitulated in transfating NSM cells, different inputs activate this GRN during normal and regulative development.